Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA parallel fluid solid coupling model using LAMMPS and Palabos based on the immersed boundary method
92
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
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).
We study experimentally the collision between a sphere falling through a viscous fluid, and a solid plate below. It is known that there is a well-defined threshold Stokes number above which the sphere rebounds from such a collision. Our experiment tests for direct contact between the colliding bodies, and contrary to prior theoretical predictions, shows that solid-on-solid contact occurs even for Stokes numbers just above the threshold for rebounding. The dissipation is fluid-dominated, though details of the contact mechanics depend on the surface and bulk properties of the solids. Our experiments and a model calculation indicate that mechanical contact between the two colliding objects is generic and will occur for any realistic surface roughness.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
