Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe divergence-conforming immersed boundary method: Application to vesicle and capsule dynamics
110
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic knot vectors to obtain discretizations of closed curves with C^2 inter-element continuity. In three-dimensional settings, we use analysis-suitable bi-cubic T-splines to obtain discretizations of closed surfaces with at least C^1 inter-element continuity. Large spurious changes of the fluid volume inside closed co-dimension one solids is a well-known issue for IB methods. The DCIB method results in volume changes orders of magnitude lower than conventional IB methods. This is a byproduct of discretizing the velocity-pressure pair with divergence-conforming B-splines, which lead to negligible incompressibility errors at the Eulerian level. The higher inter-element continuity of divergence-conforming B-splines is also crucial to avoid the quadrature/interpolation errors of IB methods becoming the dominant discretization error. Benchmark and application problems of vesicle and capsule dynamics are solved, including mesh-independence studies and comparisons with other numerical methods.
rate research
Read More
The study of viscous fluid flow coupled with rigid or deformable solids has many applications in biological and engineering problems, e.g., blood cell transport, drug delivery, and particulate flow. We developed a partitioned approach to solve this coupled Multiphysics problem. The fluid motion was solved by Palabos (Parallel Lattice Boltzmann Solver), while the solid displacement and deformation was simulated by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The coupling was achieved through the immersed boundary method (IBM). The code modeled both rigid and deformable solids exposed to flow. The code was validated with the Jeffery orbits of an ellipsoid particle in shear flow, red blood cell stretching test, and effective blood viscosity flowing in tubes. It demonstrated essentially linear scaling from 512 to 8192 cores for both strong and weak scaling cases. The computing time for the coupling increased with the solid fraction. An example of the fluid-solid coupling was given for flexible filaments (drug carriers) transport in a flowing blood cell suspensions, highlighting the advantages and capabilities of the developed code.
We present a novel moving immersed boundary method (IBM) and employ it in direct numerical simulations (DNS) of the closed-vessel swirling von Karman flow in laminar and turbulent regimes. The IBM extends direct-forcing approaches by leveraging a time integration scheme, that embeds the immersed boundary forcing step within a semi-implicit iterative Crank-Nicolson scheme. The overall method is robust, stable, and yields excellent results in canonical cases with static and moving boundaries. The moving IBM allows us to reproduce the geometry and parameters of the swirling von Karman flow experiments in (F. Ravelet, A. Chiffaudel, and F. Daviaud, JFM 601, 339 (2008)) on a Cartesian grid. In these DNS, the flow is driven by two-counter rotating impellers fitted with curved inertial stirrers. We analyze the transition from laminar to turbulent flow by increasing the rotation rate of the counter-rotating impellers to attain the four Reynolds numbers 90, 360, 2000, and 4000. In the laminar regime at Reynolds number 90 and 360, we observe flow features similar to those reported in the experiments and in particular, the appearance of a symmetry-breaking instability at Reynolds number 360. We observe transitional turbulence at Reynolds number 2000. Fully developed turbulence is achieved at Reynolds number 4000. Non-dimensional torque computed from simulations matches correlations from experimental data. The low Reynolds number symmetries, lost with increasing Reynolds number, are recovered in the mean flow in the fully developed turbulent regime, where we observe two tori symmetrical about the mid-height plane. We note that turbulent fluctuations in the central region of the device remain anisotropic even at the highest Reynolds number 4000, suggesting that isotropization requires significantly higher Reynolds numbers.
The dissolution of solids has created spectacular geomorphologies ranging from centimeter-scale cave scallops to the kilometer-scale stone forests of China and Madagascar. Mathematically, dissolution processes are modeled by a Stefan problem, which describes how the motion of a phase-separating interface depends on local concentration gradients, coupled to a fluid flow. Simulating these problems is challenging, requiring the evolution of a free interface whose motion depends on the normal derivatives of an external field in an ever-changing domain. Moreover, density differences created in the fluid domain induce self-generated convecting flows that further complicate the numerical study of dissolution processes. In this contribution, we present a numerical method for the simulation of the Stefan problem coupled to a fluid flow. The scheme uses the Immersed Boundary Smooth Extension method to solve the bulk advection-diffusion and fluid equations in the complex, evolving geometry, coupled to a {theta}-L scheme that provides stable evolution of the boundary. We demonstrate third-order temporal and pointwise spatial convergence of the scheme for the classical Stefan problem, and second-order temporal and pointwise spatial convergence when coupled to flow. Examples of dissolution of solids that result in high-Rayleigh number convection are numerically studied, and qualitatively reproduce the complex morphologies observed in recent experiments.
An important class of fluid-structure problems involve the dynamics of ordered arrays of immersed, flexible fibers. While specialized numerical methods have been developed to study fluid-fiber systems, they become infeasible when there are many, rather than a few, fibers present, nor do these methods lend themselves to analytical calculation. Here, we introduce a coarse-grained continuum model, based on local-slender body theory, for elastic fibers immersed in a viscous Newtonian fluid. It takes the form of an anisotropic Brinkman equation whose skeletal drag is coupled to elastic forces. This model has two significant benefits: (1) the density effects of the fibers in a suspension become analytically manifest, and (2) it allows for the rapid simulation of dense suspensions of fibers in regimes inaccessible to standard methods. As a first validation, without fitting parameters, we achieve very reasonable agreement with 3D Immersed Boundary simulations of a bed of anchored fibers bent by a shear flow. Secondly, we characterize the effect of density on the relaxation time of fiber beds under oscillatory shear, and find close agreement to results from full numerical simulations. We then study buckling instabilities in beds of fibers, using our model both numerically and analytically to understand the role of fiber density and the structure of buckling transitions. We next apply our model to study the flow-induced bending of inclined fibers in a channel, as has been recently studied as a flow rectifier, examining the nature of the internal flows within the bed, and the emergence of inhomogeneous permeability. Finally, we extend the method to study a simple model of metachronal waves on beds of actuated fibers, as a model for ciliary beds. Our simulations reproduce qualitatively the pumping action of coordinated waves of compression through the bed.
A conjugate heat transfer (CHT) immersed boundary (IB and CHTIB) method is developed for use with laminar and turbulent flows with low to moderate Reynolds numbers. The method is validated with the canonical flow of two co-annular rotating cylinders at $Re=50$ which shows second order accuracy of the $L_{2}$ and $L_{infty}$ error norms of the temperature field over a wide rage of solid to fluid thermal conductivities, $kappa_{s}/kappa_{f} = left(9-100right)$. To evaluate the CHTIBM with turbulent flow a fully developed, heated, turbulent channel $left(Re_{u_{tau}}=150text{ and } kappa_{s}/kappa_{f}=4 right)$ is used which shows near perfect correlation to previous direct numerical simulation (DNS) results. The CHTIB method is paired with a momentum IB method (IBM), both of which use a level set field to define the wetted boundaries of the fluid/solid interfaces and are applied to the flow solver implicitly with rescaling of the difference operators of the finite volume (FV) method (FVM).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
