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

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

An improved Rhie–Chow interpolation scheme for the smoothed‐interface immersed boundary method

Rhie–Chow interpolation is a commonly used method in CFD calculations on a co‐located mesh in order to suppress non‐physical pressure oscillations arising from chequerboard effects. A fully parallelized smoothed‐interface immersed boundary method on a co‐located grid is described in this paper. We discuss the necessity of modifications to the original Rhie–Chow interpolation in order to deal with a locally refined mesh. Numerical simulation with the modified scheme of Choi shows that numerical dissipation due to Rhie–Chow interpolation introduces significant errors at the immersed boundary. To address this issue, we develop an improved Rhie–Chow interpolation scheme that is shown to increase the accuracy in resolving the flow near the immersed boundary. We compare our improved scheme with the modified scheme of Choi by parallel simulations of benchmark flows: (i) flow past a stationary cylinder; (ii) flow past an oscillating cylinder; and (iii) flow past a stationary elliptical cylinder, where Reynolds numbers are tested in the range 10–200. Our improved scheme is significantly more accurate and compares favourably with a staggered grid algorithm. We also develop a scheme to compute the boundary force for the direct‐forcing immersed boundary method efficiently. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Computations of flow over a flexible plate using the hybrid Cartesian/immersed boundary method

A hybrid Cartesian/immersed boundary code is developed and applied to interactions between a flexible plate and a surrounding fluid. The velocities at the immersed boundary (IB) nodes are reconstructed by interpolations along local normal lines to an interface. A new criterion is suggested to distribute the IB nodes near an interface. The suggested criterion guarantees a closed fluid domain by a set of the IB nodes and it is applicable to a zero‐thickness body. To eliminate the pressure interpolation at the IB nodes, the hybrid staggered/non‐staggered grid method is adapted. The developed code is validated by comparisons with other experimental and computational results of flow around an in‐line oscillating cylinder. Good agreements are achieved for velocity profiles and vorticity and pressure contours. As applications to the fluid–structure interaction, oscillations of flexible plate in a resting fluid and flow over a flexible plate are simulated. The elastic deformations of the flexible plate are modelled based on the equations of motion for plates considering the fluid pressure as the external load on the plate. Two non‐dimensional parameters are identified and their effects on the damping of the plate motion are examined. Grid convergence tests are carried out for both cases. Copyright © 2007 John Wiley & Sons, Ltd.     

基于浸入边界法的鱼体自主游动的数值模拟

对游动或飞行生物自主运动特性的深入研究,可促进仿生学的进一步发展。本文以"C"型游动鱼作为研究对象,建立了自主游动的柔性鱼模型。此模型较为真实地反映了鱼自主游动时鱼体内力(由鱼体肌肉收缩提供)、鱼体运动和外界流体之间的耦合作用。基于传统的反馈力方法和混合有限元浸入边界法对鱼的自主游动进行了数值模拟。分析了鱼自主游动启动阶段和巡游阶段流场特性及鱼体运动特征。模拟结果表明,受到鱼体自身组织结构和外界流场作用,鱼游动时通过呈"C"型和类"S"型的不断转换,以获取能量,实现鱼体自主游动。     

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.     

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.     

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.     

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.     

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.     

An improved ghost‐cell immersed boundary method for compressible flow simulations

This study presents an improved ghost‐cell immersed boundary approach to represent a solid body in compressible flow simulations. In contrast to the commonly used approaches, in the present work, ghost cells are mirrored through the boundary described using a level‐set method to farther image points, incorporating a higher‐order extra/interpolation scheme for the ghost‐cell values. A sensor is introduced to deal with image points near the discontinuities in the flow field. Adaptive mesh refinement is used to improve the representation of the geometry efficiently in the Cartesian grid system. The improved ghost‐cell method is validated against four test cases: (a) double Mach reflections on a ramp, (b) smooth Prandtl–Meyer expansion flows, (c) supersonic flows in a wind tunnel with a forward‐facing step, and (d) supersonic flows over a circular cylinder. It is demonstrated that the improved ghost‐cell method can reach the accuracy of second order in L¹ norm and higher than first order in L^∞ norm. Direct comparisons against the cut‐cell method demonstrate that the improved ghost‐cell method is almost equally accurate with better efficiency for boundary representation in high‐fidelity compressible flow simulations. Copyright © 2016 John Wiley & Sons, Ltd.     

An immersed boundary method for the incompressible Navier–Stokes equations in complex geometry

A nodally exact convection–diffusion–reaction scheme developed in Cartesian grids is applied to solve the flow equations in irregular domains within the framework of immersed boundary (IB) method. The artificial momentum forcing term applied at certain points in the flow and inside the body of any shape allows the imposition of no‐slip velocity condition to account for the body of complex boundary. Development of an interpolation scheme that can accurately lead to no‐slip velocity condition along the IB is essential since Cartesian grid lines generally do not coincide with the IB. The results simulated from the proposed IB method agree well with other numerical and experimental results for several chosen benchmark problems. The accuracy and fidelity of the IB flow solver to predict flows with irregular IBs are therefore demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.     

Adaptive mesh refinement for simulating fluid-structure interaction using a sharp interface immersed boundary method

In this article, adaptive mesh refinement (AMR) is performed to simulate flow around both stationary and moving boundaries. The finite-difference approach is applied along with a sharp interface immersed boundary (IB) method. The Lagrangian polynomial is employed to facilitate the interpolation from a coarse to a fine grid level, while a weighted-average formula is used to transfer variables inversely. To save memory, the finest grid is only generated in the local areas close to the wall boundary, and the mesh is dynamically reconstructed based on the location of the wall boundary. The Navier-Stokes equations are numerically solved through the second-order central difference scheme in space and the third-order Runge-Kutta time integration. Flow around a circular cylinder rotating in a square domain is firstly simulated to examine the accuracy and convergence rate. Then three cases are investigated to test the validity of the present method: flow past a stationary circular cylinder at low Reynolds numbers, flow past a forced oscillating circular cylinder in the transverse direction at various frequencies, and a free circular cylinder subjected to vortex-induced vibration in two degrees of freedom. Computational results agree well with these in the literature and the flow fields are smooth around the interface of different refinement levels. The effect of refinement level has also been evaluated. In addition, a study for the computational efficiency shows that the AMR approach is helpful to reduce the total node number and speed up the time integration, which could prompt the application of the IB method when a great near-wall spatial resolution is required.     

Improvement of mass source/sink for an immersed boundary method

An improved immersed boundary method using a mass source/sink as well as momentum forcing is developed for simulating flows over or inside complex geometries. The present method is based on the Navier–Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A more accurate formulation of the mass source/sink is derived by considering mass conservation of the virtual cells in the fluid crossed by the immersed boundary. Two flow problems (the decaying vortex problem and uniform flow past a circular cylinder) are used to validate the proposed formulation. The results indicate that the accuracy near the immersed boundary is improved by introducing the accurate mass source/sink. Copyright © 2006 John Wiley & Sons, Ltd.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Computation of incompressible flows with immersed bodies by a simple ghost cell method

The incompressible Navier–Stokes equations are solved by an implicit pressure correction method on Cartesian meshes with local refinement. A simple and stable ghost cell method is developed to treat the boundary condition for the immersed bodies in the flow field. Multigrid methods are developed for both velocity and pressure correction to enhance the stability and convergence of the solution process. It is shown that the spatial accuracy of the method is second order in L₂ norm for both velocity and pressure. Various steady and unsteady flows over a 2D circular cylinder and a 3D sphere are computed to validate the present method. The capability of the present method to treat a moving body is also demonstrated. Copyright © 2008 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.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

Simulation of two‐phase flow with moving immersed boundaries

A two‐dimensional multi‐phase model for immiscible binary fluid flow including moving immersed objects is presented. The fluid motion is described by the incompressible Navier–Stokes equation coupled with a phase‐field model based on van der Waals' free energy density and the Cahn–Hilliard equation. A new phase‐field boundary condition was implemented with minimization of the free energy in a direct way, to specifically improve the physical behavior of the contact line dynamics for moving immersed objects. Numerical stability and execution time were significantly improved by the use of the new boundary condition. Convergence toward the analytical solution was demonstrated for equilibrium contact angle, the Lucas–Washburn theory and Stefan's problem. The proposed model may be used for multi‐phase flow problems with moving boundaries of complex geometry, such as the penetration of fluid into a deformable, porous medium. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of a three-dimensional solver based on the local domain-free discretization and immersed boundary method and its application for incompressible flow problems

Recently, the author and two other coauthors have proposed a two-dimensional hybrid local domain-free discretization and immersed boundary method (LDFD-IBM), which can be used to solve the flow problem with complex geometries. In this paper, the LDFD-IBM is extended to solve a three-dimensional unsteady incompressible flow with the complex computational domain. The technical issues related to the implementation of the LDFD-IBM in three-dimensional problems are discussed in detail, particularly for the discretization of Navier-Stokes equations, mesh strategies for a three-dimensional flow, and the fast algorithm on the identification of the status of mesh nodes (ie, to identify if the mesh node is located in the solid domain, in the fluid domain, or near the immersed boundary). Numerical tests show that the LDFD-IBM can accurately solve three-dimensional incompressible problems with ease.