Three‐dimensional immersed boundary conditions for moving solids in the lattice‐Boltzmann method

This paper establishes the range of validity for a previously published three‐dimensional moving solid boundary condition for the lattice‐Boltzmann method. This method was reasonably formulated from a mass and momentum balance perspective, but was only verified for a small range of (primarily two‐dimensional) problems. One of the advantages of this boundary condition is that it offers resolution at the sub‐grid scale, allowing for accurate and stable calculation of the force and torque for solids which are moving through a lattice, even for small solid sizes relative to the computational grid size. We verify the boundary condition for creeping flows by comparison to analytical solutions that include both the force and the torque on fixed and moving spheres, and then follow this with comparisons to experimental and empirical results for both fixed as well moving spheres in inertial flows. Finally, we compare simulation results to numerical results of other investigators for the settling of an offset sphere and the drafting–kissing–tumbling of two sedimenting spheres. We found that an accurate calculation of the collision‐operator weighting used to obtain sub‐grid‐scale resolution was necessary in order to prevent spikes in the velocities, forces, and moments when solid objects cross‐computational cells. The wide range of comparisons collected and presented in this paper can be used to establish the validity of other numerical models, in addition to the one examined here. Published in 2007 by John Wiley & Sons, Ltd.     

On the immersed boundary‐lattice Boltzmann simulations of incompressible flows with freely moving objects

For simulating freely moving problems, conventional immersed boundary‐lattice Boltzmann methods encounter two major difficulties of an extremely large flow domain and the incompressible limit. To remove these two difficulties, this work proposes an immersed boundary‐lattice Boltzmann flux solver (IB‐LBFS) in the arbitrary Lagragian–Eulerian (ALE) coordinates and establishes a dynamic similarity theory. In the ALE‐based IB‐LBFS, the flow filed is obtained by using the LBFS on a moving Cartesian mesh, and the no‐slip boundary condition is implemented by using the boundary condition‐enforced immersed boundary method. The velocity of the Cartesian mesh is set the same as the translational velocity of the freely moving object so that there is no relative motion between the plate center and the mesh. This enables the ALE‐based IB‐LBFS to study flows with a freely moving object in a large open flow domain. By normalizing the governing equations for the flow domain and the motion of rigid body, six non‐dimensional parameters are derived and maintained to be the same in both physical systems and the lattice Boltzmann framework. This similarity algorithm enables the lattice Boltzmann equation‐based solver to study a general freely moving problem within the incompressible limit. The proposed solver and dynamic similarity theory have been successfully validated by simulating the flow around an in‐line oscillating cylinder, single particle sedimentation, and flows with a freely falling plate. The obtained results agree well with both numerical and experimental data. Copyright © 2016 John Wiley & Sons, Ltd.     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

A hybrid immersed‐boundary and multi‐block lattice Boltzmann method for simulating fluid and moving‐boundaries interactions

A numerical method is developed for modelling the interactions between incompressible viscous fluid and moving boundaries. The principle of this method is introducing the immersed‐boundary concept in the framework of the lattice Boltzmann method, and improving the accuracy and efficiency of the simulation by refining the mesh near moving boundaries. Besides elastic boundary with a constitutive law, the method can also efficiently simulate solid moving‐boundary interacting with fluid by employing the direct forcing technique. The method is validated by the simulations of flow past a circular cylinder, two cylinders moving with respect to each other and flow around a hovering wing. The versatility of the method is demonstrated by the numerical studies including elastic filament flapping in the wake of a cylinder and fish‐like bodies swimming in quiescent fluid. Copyright © 2006 John Wiley & Sons, Ltd.     

An improved hybrid Cartesian/immersed boundary method for fluid–solid flows

An improved hybrid Cartesian/immersed boundary method is proposed based on ghost point treatment. A second‐order Taylor series expansion is used to evaluate the values at the ghost points, and an inverse distance weighting method to interpolate the values due to its properties of preserving local extrema and smooth reconstruction. The present method effectively eliminates numerical instabilities caused by matrix inversion and flexibly adopts the interpolation in the vicinity of the boundary. Some typical fluid–solid flows, including viscous flow past a circular cylinder, a sphere, two cylinders in a side‐by‐side arrangement, and an array of 18 staggered cylinders, are examined. These benchmark simulations reasonably indicate the reliability and capability of the present method. Copyright © 2007 John Wiley & Sons, Ltd.     

A divergence‐free interpolation scheme for the immersed boundary method

The immersed boundary approach for the modeling of complex geometries in incompressible flows is examined critically from the perspective of satisfying boundary conditions and mass conservation. It is shown that the system of discretized equations for mass and momentum can be inconsistent, if the velocity is used in defining the force density to satisfy the boundary conditions. As a result, the velocity is generally not divergence free and the pressure at locations in the vicinity of the immersed boundary is not physical. However, the use of the pseudo‐velocities in defining the force density, as frequently done when the governing equations are solved using a fractional step or projection method, combined with the use of the specified velocity on the immersed boundary, is shown to result in a consistent set of equations which allows a divergence‐free velocity but, depending on the time step, is shown to have the undesirable effects of inaccurately satisfying the boundary conditions and allowing a significant permeability of the immersed boundary. If the time step is reduced sufficiently, the boundary conditions on the immersed boundary can be satisfied. However, this entails an unacceptable increase in computational expense. Two new methods that satisfy the boundary conditions and allow a divergence‐free velocity while avoiding the increased computational expense are presented and shown to be second‐order accurate in space. The first new method is based on local time step reduction. This method is suitable for problems where the immersed boundary does not move. For these problems, the first new method is shown to be closely related to the second new method. The second new method uses an optimization scheme to minimize the deviation from the interpolation stencil used to represent the immersed boundary while ensuring a divergence‐free velocity. This method performs well for all problems, including those where the immersed boundary moves relative to the grid. Additional results include showing that the force density that is added to satisfy the boundary conditions at the immersed boundary is unbounded as the time step is reduced and that the pressure in the vicinity of the immersed boundary is unphysical, being strongly a function of the time step. A method of computing the total force on an immersed boundary which takes into account the specifics of the numerical solver used in the iterative process and correctly computes the total force irrespective of the residual level is also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

A comparative study of direct‐forcing immersed boundary‐lattice Boltzmann methods for stationary complex boundaries

In this study, we assess several interface schemes for stationary complex boundary flows under the direct‐forcing immersed boundary‐lattice Boltzmann methods (IB‐LBM) based on a split‐forcing lattice Boltzmann equation (LBE). Our strategy is to couple various interface schemes, which were adopted in the previous direct‐forcing immersed boundary methods (IBM), with the split‐forcing LBE, which enables us to directly use the direct‐forcing concept in the lattice Boltzmann calculation algorithm with a second‐order accuracy without involving the Navier–Stokes equation. In this study, we investigate not only common diffuse interface schemes but also a sharp interface scheme. For the diffuse interface scheme, we consider explicit and implicit interface schemes. In the calculation of velocity interpolation and force distribution, we use the 2‐ and 4‐point discrete delta functions, which give the second‐order approximation. For the sharp interface scheme, we deal with the exterior sharp interface scheme, where we impose the force density on exterior (solid) nodes nearest to the boundary. All tested schemes show a second‐order overall accuracy when the simulation results of the Taylor–Green decaying vortex are compared with the analytical solutions. It is also confirmed that for stationary complex boundary flows, the sharper the interface scheme, the more accurate the results are. In the simulation of flows past a circular cylinder, the results from each interface scheme are comparable to those from other corresponding numerical schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparative study of lattice Boltzmann methods using bounce‐back schemes and immersed boundary ones for flow acoustic problems

In order to find applicable treatments of moving boundary conditions based on the lattice Boltzmann method in flow acoustic problems, three bounce‐back (BB) methods and four kinds of immersed boundary (IB) methods are compared. We focused on fluid–solid boundary conditions for flow acoustic problems especially the simulations of sound waves from moving boundaries. BB methods include link bounce‐back, interpolation bounce‐back and unified interpolation bounce‐back methods. Five IB methods are explicit and implicit direct‐forcing (Explicit‐IB and Implicit‐IB), two kinds of partially saturated computational methods and ghost fluid method. In order to reduce the spurious pressure generated by the fresh grid node changing from solid domain to fluid domain for BB methods and sharp IB methods, we proposed two new kinds of treatments and compared them with two existing ones. Simulations of the benchmark problems prove that the local evolutionary iteration (LI) is the best one in treatments of the fresh nodes. In addition, for standing boundary problems, although BB methods have a little higher accuracy, all the methods have similar accuracy. However, for moving boundary problems, IB methods are more appropriate than BB methods, because IB methods' smooth interpolation of pressure eld produces less disturbing spurious pressure waves. With improved treatments of fresh nodes, BB methods are also acceptable for moving boundary acoustic problems. In comparative tests in respective type, unified interpolation bounce‐back with LI, Implicit‐IB, and ghost fluid with LI are the best choices. Copyright © 2013 John Wiley & Sons, Ltd.     

A multi‐energy‐level lattice Boltzmann model for the compressible Navier–Stokes equations

In this paper, we propose a new lattice Boltzmann model for the compressible Navier–Stokes equations. The new model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. As the 25‐bit model, we obtained the equilibrium distribution functions and the compressible Navier–Stokes equations with the second accuracy of the truncation errors. The numerical examples show that the model can be used to simulate the shock waves, contact discontinuities and supersonic flows around circular cylinder. The numerical results are compared with those obtained by traditional method. Copyright © 2007 John Wiley & Sons, Ltd.     

A hybrid immersed boundary and material point method for simulating 3D fluid–structure interaction problems

A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd.     

Analytical and numerical studies of the boundary slip in the immersed boundary‐thermal lattice Boltzmann method

We analytically and numerically investigate the boundary slip, including the velocity slip and the temperature jump, in immersed boundary‐thermal lattice Boltzmann methods (IB‐TLBMs) with the two‐relaxation‐time collision operator. We derive the theoretical equation for the relaxation parameters considering the effect of the advection velocity on the temperature jump of the IB‐TLBMs. The analytical and numerical solutions demonstrate that the proposed iterative correction methods without the computational cost of the sparse matrix solver reduce the boundary slip and boundary‐value deviation as effectively as the implicit correction method for any relaxation time. Because the commonly used multi‐direct forcing method does not consider the contributions of the body force to the momentum flux, it cannot completely eliminate the boundary slip because of the numerical instability for a long relaxation time. Both types of proposed iterative correction methods are more numerically stable than the implicit correction method. In simulations of flow past a circular cylinder and of natural convection, the present iterative correction methods yield adequate results without the errors of the velocity slip, the temperature jump, and the boundary‐value deviation for any relaxation time parameters and for any number of Lagrangian points per length. The combination of the present methods and the two‐relaxation‐time collision operator is suitable for simulating fluid flow with thermal convection in the multiblock method in which the relaxation time increases in inverse proportion to the grid size.     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems

Embedded boundary methods for CFD (computational fluid dynamics) simplify a number of issues. These range from meshing the fluid domain, to designing and implementing Eulerian‐based algorithms for fluid–structure applications featuring large structural motions and/or deformations. Unfortunately, embedded boundary methods also complicate other issues such as the treatment of the wall boundary conditions in general, and fluid–structure transmission conditions in particular. This paper focuses on this aspect of the problem in the context of compressible flows, the finite volume method for the fluid, and the finite element method for the structure. First, it presents a numerical method for treating simultaneously the fluid pressure and velocity conditions on static and dynamic embedded interfaces. This method is based on the exact solution of local, one‐dimensional, fluid–structure Riemann problems. Next, it describes two consistent and conservative approaches for computing the flow‐induced loads on rigid and flexible embedded structures. The first approach reconstructs the interfaces within the CFD solver. The second one represents them as zero level sets, and works instead with surrogate fluid/structure interfaces. For example, the surrogate interfaces obtained simply by joining contiguous segments of the boundary surfaces of the fluid control volumes that are the closest to the zero level sets are explored in this work. All numerical algorithms presented in this paper are applicable with any embedding CFD mesh, whether it is structured or unstructured. Their performance is illustrated by their application to the solution of three‐dimensional fluid–structure interaction problems associated with the fields of aeronautics and underwater implosion. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐dimensional compact finite difference immersed boundary method

We present a compact finite differences method for the calculation of two‐dimensional viscous flows in biological fluid dynamics applications. This is achieved by using body‐forces that allow for the imposition of boundary conditions in an immersed moving boundary that does not coincide with the computational grid. The unsteady, incompressible Navier–Stokes equations are solved in a Cartesian staggered grid with fourth‐order Runge–Kutta temporal discretization and fourth‐order compact schemes for spatial discretization, used to achieve highly accurate calculations. Special attention is given to the interpolation schemes on the boundary of the immersed body. The accuracy of the immersed boundary solver is verified through grid convergence studies. Validation of the method is done by comparison with reference experimental results. In order to demonstrate the application of the method, 2D small insect hovering flight is calculated and compared with available experimental and computational results. Copyright © 2009 John Wiley & Sons, Ltd.     

An immersed boundary method for fluid flows around rigid objects

In this paper, we present an immersed boundary method for solving fluid flow problems in the presence of static and moving rigid objects. A FEM is used starting from a base mesh that does not represent exactly rigid objects (non?body?conforming mesh). At each time step, the base mesh is locally modified to provide a new mesh fitting the boundary of the rigid objects. The mesh is also locally improved using edge swapping to enhance the quality of the elements. The Navier–Stokes equations are then solved on this new mesh. The velocity of moving objects is imposed through standard Dirichlet boundary conditions. We consider a number of test problems and compare the numerical solutions with those obtained on classical body?fitted meshes whenever possible. Copyright © 2014 John Wiley & Sons, Ltd.     

An improved immersed boundary‐lattice Boltzmann method based on force correction technique

In this paper, an improved immersed boundary‐lattice Boltzmann method based on the force correction technique is presented for fluid‐structure interaction problems including the moving boundary interfaces. By introducing a force correction coefficient, the non‐slip boundary conditions are much better enforced compared with the conventional immersed boundary‐lattice Boltzmann methods. In addition, the implicit and iterative calculations are avoided; thus, the computational cost is reduced dramatically. Several numerical experiments are carried out to test the efficiency of the method. It is found that the method has the second‐order accuracy, and the non‐slip boundary conditions are enforced indeed. The numerical results also show that the present method is a suitable tool for fluid‐structure interaction problems involving complex moving boundaries.     

Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods

A robust, accurate, and computationally efficient interface tracking algorithm is a key component of an embedded computational framework for the solution of fluid–structure interaction problems with complex and deformable geometries. To a large extent, the design of such an algorithm has focused on the case of a closed embedded interface and a Cartesian computational fluid dynamics grid. Here, two robust and efficient interface tracking computational algorithms capable of operating on structured as well as unstructured three‐dimensional computational fluid dynamics grids are presented. The first one is based on a projection approach, whereas the second one is based on a collision approach. The first algorithm is faster. However, it is restricted to closed interfaces and resolved enclosed volumes. The second algorithm is therefore slower. However, it can handle open shell surfaces and underresolved enclosed volumes. Both computational algorithms exploit the bounding box hierarchy technique and its parallel distributed implementation to efficiently store and retrieve the elements of the discretized embedded interface. They are illustrated, and their respective performances are assessed and contrasted, with the solution of three‐dimensional, nonlinear, dynamic fluid–structure interaction problems pertaining to aeroelastic and underwater implosion applications. Copyright © 2012 John Wiley & Sons, Ltd.