A finite element immersed boundary method for fluid flow around rigid objects

This paper proposes a new immersed boundary (IB) method for solving fluid flow problems in the presence of rigid objects which are not represented by the mesh. Solving the flow around objects with complex shapes may involve extensive meshing work that has to be repeated each time a change in the geometry is needed. Important benefit would be reached if we are able to solve the flow without the need of generating a mesh that fits the shape of the immersed objects. This work presents a finite element IB method using a discretization covering the entire domain of interest, including the volume occupied by immersed objects, and which produces solutions of the flow satisfying accurately the boundary conditions at the surface of immersed bodies. In other words the finite element solution represents accurately the presence of immersed bodies while the mesh does not. This is done by including additional degrees of freedom on interface cut elements which are then eliminated at element level. The boundary of immersed objects is defined using a level set function. Solutions are shown for various flow problems and the accuracy of the present approach is measured with respect to solutions obtained on body‐fitted meshes. Copyright © 2010 Crown in the right of Canada.     

An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

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.     

An immersed smoothed point interpolation method (IS‐PIM) for fluid‐structure interaction problems

An immersed smoothed point interpolation method using 3‐node triangular background cells is proposed to solve 2D fluid‐structure interaction problems for solids with large deformation/displacement placed in incompressible viscous fluid. In the framework of immersed‐type method, the governing equations can be decomposed into 3 parts on the basis of the fictitious fluid assumption. The incompressible Navier‐Stokes equations are solved using the semi‐implicit characteristic‐based split scheme, and solids are simulated using the newly developed edge‐based smoothed point interpolation method. The fictitious fluid domain can be used to calculate the coupling force. The numerical results show that immersed smoothed point interpolation method can avoid remeshing for moving solid based on immersed operation and simulate the contact phenomenon without an additional treatment between the solid and the fluid boundary. The influence from information transfer between solid domain and fluid domain on fluid‐structure interaction problems has been investigated. The numerical results show that the proposed interpolation schemes will generally improve the accuracy for simulating both fluid flows and solid structures.     

A ghost-cell immersed boundary method for large eddy simulation of flows in complex geometries

An efficient ghost-cell immersed boundary (IB) method is proposed for large eddy simulations of three-dimensional incompressible flow in complex geometries. In the framework of finite volume method, the Navier–Stokes equations are integrated using an explicit time advancement scheme on a collocated mesh. Since the IB method is known to generate an unphysical velocity field inside the IB that violates the mass conservation of the cells near the IB, a new IB treatment is devised to eliminate the unphysical velocity generated near the IB and to improve the pressure distribution on the body surface. To validate the proposed method, both laminar and turbulent flow cases are presented. In particular, large eddy simulations were performed to simulate the turbulent flows over a circular cylinder and a sphere at subcritical Reynolds numbers. The computed results show good agreements with the published numerical and experimental data.     

Numerical modeling of flow past a volumeless and thin rigid body using direct forcing immersed boundary method

The new capability has been added as the numerical method for modeling volumeless and thin rigid bodies to the direct forcing immersed boundary (DFIB) method. The DFIB approach is based on adding a virtual force to the Navier–Stokes equations of incompressible flow to account for the interaction between the fluid and structures. The volume of a solid function (VOS) identifies the stationary or moving solid structures in a given fluid domain. A new VOS-based algorithm was developed to identify thin, rigid structure boundary points in fluid flow and ensure that the fluid cannot cross through the boundary of a thin rigid structure while moving or stationary. The DFIB method was first validated in a three-dimensional (3D) turbulent flow over a circular cylinder. The large-eddy simulation simulated the turbulent flow scales. The proposed algorithm was tested using a 3D turbulent flow past a stationary and rotating Savonius wind turbine that functions as a thin, rigid body. The validation results showed that the selected DFIB approach, combined with the novel algorithm, could simulate a thin, volumeless, rigid structure that is stationary and rotating in incompressible turbulent flows. The current method is also applicable for two-way fluid-structure interaction problems.     

Moving immersed boundary method

A novel implicit immersed boundary method of high accuracy and efficiency is presented for the simulation of incompressible viscous flow over complex stationary or moving solid boundaries. A boundary force is often introduced in many immersed boundary methods to mimic the presence of solid boundary, such that the overall simulation can be performed on a simple Cartesian grid. The current method inherits this idea and considers the boundary force as a Lagrange multiplier to enforce the no‐slip constraint at the solid boundary, instead of applying constitutional relations for rigid bodies. Hence excessive constraint on the time step is circumvented, and the time step only depends on the discretization of fluid Navier‐Stokes equations, like the CFL condition in present work. To determine the boundary force, an additional moving force equation is derived. The dimension of this derived system is proportional to the number of Lagrangian points describing the solid boundaries, which makes the method very suitable for moving boundary problems since the time for matrix update and system solving is not significant. The force coefficient matrix is made symmetric and positive definite so that the conjugate gradient method can solve the system quickly. The proposed immersed boundary method is incorporated into the fluid solver with a second‐order accurate projection method as a plug‐in. The overall scheme is handled under an efficient fractional step framework, namely, prediction, forcing, and projection. Various simulations are performed to validate current method, and the results compare well with previous experimental and numerical studies.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

A new modification of the immersed‐boundary method for simulating flows with complex moving boundaries

In this paper, a new immersed‐boundary method for simulating flows over complex immersed, moving boundaries is presented. The flow is computed on a fixed Cartesian mesh and the solid boundaries are allowed to move freely through the mesh. The present method is based on a finite‐difference approach on a staggered mesh together with a fractional‐step method. It must be noted that the immersed boundary is generally not coincident with the position of the solution variables on the grid, therefore, an appropriate strategy is needed to construct a relationship between the curved boundary and the grid points nearby. Furthermore, a momentum forcing is added on the body boundaries and also inside the body to satisfy the no‐slip boundary condition. The immersed boundary is represented by a series of interfacial markers, and the markers are also used as Lagrangian forcing points. A linear interpolation is then used to scale the Lagrangian forcing from the interfacial markers to the corresponding grid points nearby. This treatment of the immersed‐boundary is used to simulate several problems, which have been validated with previous experimental results in the open literature, verifying the accuracy of the present method. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order stable interpolations for immersed boundary methods

The analysis and improvement of an immersed boundary method (IBM) for simulating turbulent flows over complex geometries are presented. Direct forcing is employed. It consists in interpolating boundary conditions from the solid body to the Cartesian mesh on which the computation is performed. Lagrange and least squares high‐order interpolations are considered. The direct forcing IBM is implemented in an incompressible finite volume Navier–Stokes solver for direct numerical simulations (DNS) and large eddy simulations (LES) on staggered grids. An algorithm to identify the body and construct the interpolation schemes for arbitrarily complex geometries consisting of triangular elements is presented. A matrix stability analysis of both interpolation schemes demonstrates the superiority of least squares interpolation over Lagrange interpolation in terms of stability. Preservation of time and space accuracy of the original solver is proven with the laminar two‐dimensional Taylor–Couette flow. Finally, practicability of the method for simulating complex flows is demonstrated with the computation of the fully turbulent three‐dimensional flow in an air‐conditioning exhaust pipe. Copyright © 2006 John Wiley & Sons, Ltd.     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

An immersed boundary vector potential-vorticity meshless solver of the incompressible Navier–Stokes equation

We present a strong form meshless solver for numerical solution of the nonstationary, incompressible, viscous Navier–Stokes equations in two (2D) and three dimensions (3D). We solve the flow equations in their stream function-vorticity (in 2D) and vector potential-vorticity (in 3D) formulation, by extending to 3D flows the boundary condition-enforced immersed boundary method, originally introduced in the literature for 2D problems. We use a Cartesian grid, uniform or locally refined, to discretize the spatial domain. We apply an explicit time integration scheme to update the transient vorticity equations, and we solve the Poisson type equation for the stream function or vector potential field using the meshless point collocation method. Spatial derivatives of the unknown field functions are computed using the discretization-corrected particle strength exchange method. We verify the accuracy of the proposed numerical scheme through commonly used benchmark and example problems. Excellent agreement with the data from the literature was achieved. The proposed method was shown to be very efficient, having relatively large critical time steps.     

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.     

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.     

Solution of flow around complex‐shaped surfaces by an immersed boundary‐body conformal enrichment method

Solving the flow around objects with complex shapes may involve extensive meshing work that has to be repeated each time a change in the geometry is needed. Time consuming meshing can be avoided when the solution algorithm can tackle grids that do not fit the shape of immersed objects. This work presents applications of a recently proposed immersed boundary—body conformal enrichment method to the solution of the flow around complex shaped surfaces such as those of a metallic foam matrix. The method produces solutions of the flow satisfying accurately Dirichlet boundary conditions imposed on the immersed fluid/solid interface. The boundary of immersed objects is defined using a level‐set function, and the finite element discretization of interface elements is enriched with additional degrees of freedom, which are eliminated at element level. The method is first validated in the case of flow problems for which reference solutions on body‐conformal grids can be obtained: flow around an array of spheres and flow around periodic arrays of cylinders. Then, solutions are shown for the more complex flow inside a metallic foam matrix. A multiscale approach combining the solution at the pore level by the immersed boundary method and the macro‐scale solution with simulated permeability is used to solve actual experimental configurations. The computed pressure drop as a function of the flow rate on the macro scale configuration replicating two experimental setups is compared with the experimental data for various foam thicknesses. Copyright © 2011 National Research Council Canada     

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.     

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.     

Application of local DFD method to simulate unsteady flows around an oscillating circular cylinder

In this paper, the recently proposed local domain‐free discretization (DFD) method is applied to simulate incompressible flows around an oscillating circular cylinder. It is found that it is very easy for the local DFD method to handle such moving boundary flow problems. This is because it does not need to move the mesh, which is indeed needed in traditional methods. Numerical experiments show that the present numerical results agree very well with the available data in the literature, and that the local DFD method is an effective tool for the computation of moving boundary flow problems. Copyright © 2008 John Wiley & Sons, Ltd.     

An immersed boundary method for unstructured meshes in depth averaged shallow water models

The representation of geometries as buildings, flood barriers or dikes in free surface flow models implies tedious and time‐consuming operations in order to define accurately the shape of these objects when using a body fitted numerical mesh. The immersed boundary method is an alternative way to define solid bodies inside the computational domain without the need of fitting the mesh boundaries to the shape of the object. In the direct forcing immersed boundary method, a solid body is represented by a grid of Lagrangian markers, which define its shape and which are independent from the fluid Eulerian mesh. This paper presents a new implementation of the immersed boundary method in an unstructured finite volume solver for the 2D shallow water equations. Moving least‐squares is used to transmit information between the grid of Lagrangian markers and the fluid Eulerian mesh. The performance of the proposed implementation is analysed in three test cases involving different flow conditions: the flow around a spur dike, a dam break flow with an isolated obstacle and the flow around an array of obstacles. A very good agreement between the classic body fitted approach and the immersed boundary method was found. The differences between the results obtained with both methods are less relevant than the errors because of the intrinsic shallow water assumptions. Copyright © 2015 John Wiley & Sons, Ltd.     

Explicit formulations for evaluation of velocity gradients using boundary‐domain integral equations in 2D and 3D viscous flows

In this paper, explicit boundary‐domain integral equations for evaluating velocity gradients are derived from the basic velocity integral equations. A free term is produced in the new strongly singular integral equation, which is not included in recent formulations using the complex variable differentiation method (CVDM) to compute velocity gradients (Int. J. Numer. Meth. Fluids 2004; 45 :463–484; Int. J. Numer. Meth. Fluids 2005; 47 :19–43). The strongly singular domain integrals involved in the new integral equations are accurately evaluated using the radial integration method (RIM). Considerable computational time for evaluating integrals of velocity gradients can be saved by using present formulation than using CVDM. The formulation derived in this paper together with those presented in reference (Int. J. Numer. Meth. Fluids 2004; 45 :463–484) for 2D and in (Int. J. Numer. Meth. Fluids 2005; 47 :19–43) for 3D problems constitutes a complete boundary‐domain integral equation system for solving full Navier–Stokes equations using primitive variables. Three numerical examples for steady incompressible viscous flow are given to validate the derived formulations. Copyright © 2007 John Wiley & Sons, Ltd.