An immersed boundary-lattice Boltzmann method for single- and multi-component fluid flows

International audienceThe paper presents a numerical method to simulate single-and multi-component fluid flows around moving/deformable solid boundaries, based on the coupling of Immersed Boundary (IB) and Lattice Boltzmann (LB) methods. The fluid domain is simulated with LB method using the single relaxation time BGK model, in which an interparticle potential model is applied for multi-component fluid flows. The IB-related force is directly calculated with the interpolated definition of the fluid macroscopic velocity on the Lagrangian points that define the immersed solid boundary. The present IB-LB method can better ensure the no-slip solid boundary condition, thanks to an improved spreading operator. The proposed method is validated through several 2D/3D single-and multi-component fluid test cases with a particular emphasis on wetting conditions on solid wall. Finally, a 3D two-fluid application case is given to show the feasibility of modeling the fluid transport via a cluster of beating cilia

Kinematic and dynamic forcing strategies for predicting the transport of inertial capsules via a combined lattice Boltzmann-Immersed Boundary method

Modeling the transport of deformable capsules under different flow regimens is crucial in a variety of fields, including oil rheology, blood flow and the dispersion of pollutants. The aim of this study is twofold. Firstly, a combined Lattice Boltzmann-Immersed Boundary (LBM-IB) approach is developed for predicting the transport of inertial deformable capsules. A Moving Least Squares (MLS) scheme has been implemented to correlate the pressure, velocity and force fields of the fluid domain with the capsule dynamics. This computational strategy has been named LBM Dynamic IB. Secondly, this strategy is directly compared with a more conventional approach, named LBM Kinematic IB, where capsules move with the same velocity of the surrounding fluid. Multiple test cases have been considered for assessing the accuracy and efficiency of the Dynamic over Kinematic IB scheme, including the stretching of circular capsules in shear flow, the transport in a plane Poiseuille flow of circular and biconcave capsules, with and without inertia. By monitoring the capsule geometry over time, the two schemes have been documented to be in excellent agreement, especially for low Capillary numbers (Ca $\leq$ 0.01), in the case of non-inertial capsules. Despite a moderate increase in computational burden, the presented LBM Dynamic IB scheme is the sole capable of predicting the dynamics of both non-inertial and inertial deformable capsules. The proposed approach can be efficiently employed for studying the transport of blood cells, cancer cells and nano/micro capsules within a capillary flow

Simulation of incompressible viscous flows around moving objects by a variant of immersed boundary-Lattice Boltzmann method

A variant of immersed boundary-lattice Boltzmann method (IB-LBM) is presented in this paper to simulate incompressible viscous flows around moving objects. As compared with the conventional IB-LBM where the force density is computed explicitly by Hook's law or the direct forcing method and the non-slip condition is only approximately satisfied, in the present work, the force density term is considered as the velocity correction which is determined by enforcing the non-slip condition at the boundary. The lift and drag forces on the moving object can be easily calculated via the velocity correction on the boundary points. The capability of the present method for moving objects is well demonstrated through its application to simulate flows around a moving circular cylinder, a rotationally oscillating cylinder, and an elliptic flapping wing. Furthermore, the simulation of flows around a flapping flexible airfoil is carried out to exhibit the ability of the present method for implementing the elastic boundary condition. It was found that under certain conditions, the flapping flexible airfoil can generate larger propulsive force than the flapping rigid airfoil

Recent advances in the simulation of particle-laden flows

A substantial number of algorithms exists for the simulation of moving particles suspended in fluids. However, finding the best method to address a particular physical problem is often highly non-trivial and depends on the properties of the particles and the involved fluid(s) together. In this report we provide a short overview on a number of existing simulation methods and provide two state of the art examples in more detail. In both cases, the particles are described using a Discrete Element Method (DEM). The DEM solver is usually coupled to a fluid-solver, which can be classified as grid-based or mesh-free (one example for each is given). Fluid solvers feature different resolutions relative to the particle size and separation. First, a multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine resolution) is presented to study the behavior of particle stabilized fluid interfaces and second, a Smoothed Particle Hydrodynamics implementation (mesh-free, meso-scale resolution, similar to the particle size) is introduced to highlight a new player in the field, which is expected to be particularly suited for flows including free surfaces.Comment: 16 pages, 4 figure

Bridging the computational gap between mesoscopic and continuum modeling of red blood cells for fully resolved blood flow

We present a computational framework for the simulation of blood flow with fully resolved red blood cells (RBCs) using a modular approach that consists of a lattice Boltzmann solver for the blood plasma, a novel finite element based solver for the deformable bodies and an immersed boundary method for the fluid-solid interaction. For the RBCs, we propose a nodal projective FEM (npFEM) solver which has theoretical advantages over the more commonly used mass-spring systems (mesoscopic modeling), such as an unconditional stability, versatile material expressivity, and one set of parameters to fully describe the behavior of the body at any mesh resolution. At the same time, the method is substantially faster than other FEM solvers proposed in this field, and has an efficiency that is comparable to the one of mesoscopic models. At its core, the solver uses specially defined potential energies, and builds upon them a fast iterative procedure based on quasi-Newton techniques. For a known material, our solver has only one free parameter that demands tuning, related to the body viscoelasticity. In contrast, state-of-the-art solvers for deformable bodies have more free parameters, and the calibration of the models demands special assumptions regarding the mesh topology, which restrict their generality and mesh independence. We propose as well a modification to the potential energy proposed by Skalak et al. 1973 for the red blood cell membrane, which enhances the strain hardening behavior at higher deformations. Our viscoelastic model for the red blood cell, while simple enough and applicable to any kind of solver as a post-convergence step, can capture accurately the characteristic recovery time and tank-treading frequencies. The framework is validated using experimental data, and it proves to be scalable for multiple deformable bodies

Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [ Phys. Rev. E 47 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed