流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

A combined level set/ghost cell immersed boundary representation for floating body simulations

A six degrees of freedom (6DOF) algorithm is implemented in the open‐source CFD code REEF3D. The model solves the incompressible Navier–Stokes equations. Complex free surface dynamics are modeled with the level set method based on a two‐phase flow approach. The convection terms of the velocities and the level set method are treated with a high‐order weighted essentially non‐oscillatory discretization scheme. Together with the level set method for the free surface capturing, this algorithm can model the movement of rigid floating bodies and their interaction with the fluid. The 6DOF algorithm is implemented on a fixed grid. The solid‐fluid interface is represented with a combination of the level set method and ghost cell immersed boundary method. As a result, re‐meshing or overset grids are not necessary. The capability, accuracy, and numerical stability of the new algorithm is shown through benchmark applications for the fluid‐body interaction problem. Copyright © 2016 John Wiley & Sons, Ltd.     

A mass‐conservative staggered immersed boundary model for solving the shallow water equations on complex geometries

In this work, an approach is proposed for solving the 3D shallow water equations with embedded boundaries that are not aligned with the underlying horizontal Cartesian grid. A hybrid cut‐cell/ghost‐cell method is used together with a direction‐splitting implicit solver: Ghost cells are used for the momentum equations in order to prescribe the correct boundary condition at the immersed boundary, while cut cells are used in the continuity equation in order to conserve mass. The resulting scheme is robust, does not suffer any time step limitation for small cut cells, and conserves fluid mass up to machine precision. Moreover, the solver displays a second‐order spatial accuracy, both globally and locally. Comparisons with analytical solutions and reference numerical solutions on curvilinear grids confirm the quality of the method. Copyright © 2015 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.     

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.     

Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

Efficient coupling of direct forcing immersed boundary-lattice Boltzmann method and finite element method to simulate fluid-structure interactions

Fluid-structure interactions (FSI) of rigid and flexible bodies are simulated in this article. For the fluid flow, multidirect forcing immersed boundary method (IBM) is adopted to capture the moving boundary, and lattice Boltzmann method (LBM) is used to evolve the flow field. Compared with our previous no-penetration IBM, less iterations are required in this work. In addition, larger velocity in lattice units can be used and the nonphysical force oscillations are suppressed due to the C³ 6-point kernel. Multi-relaxation-time collision operator and local grid refinement are also adopted in LBM to enhance the numerical stability and efficiency. To improve the efficiency of the FSI coupling algorithm, the mesh of the deformable structure can be coarser than the Lagrangian mesh using Newton-Cotes formulas to integrate the traction on the structure surface. A variety of benchmarks, including flow around a circular cylinder with Reynold number ranging from 20 to 200, forced oscillation of a circular cylinder, vortex-induced vibration (VIV) of an elastically mounted circular cylinder in two degrees of freedom, and VIV of an elastic cantilever beam attached to a circular cylinder, are carried out to evaluate the accuracy and stability of the present coupling algorithm. For the benchmark FSI problem considered in this article, a reduction of 54% of the calculation time is achieved using coarser structure mesh. As an application of the FSI coupling algorithm, the mechanism of an elastic beam in the wake of a circular cylinder is discussed.     

Semi‐coupled air/water immersed boundary approach for curvilinear dynamic overset grids with application to ship hydrodynamics

For many problems in ship hydrodynamics, the effects of air flow on the water flow are negligible (the frequently called free surface conditions), but the air flow around the ship is still of interest. A method is presented where the water flow is decoupled from the air solution, but the air flow uses the unsteady water flow as a boundary condition. The authors call this a semi‐coupled air/water flow approach. The method can be divided into two steps. At each time step the free surface water flow is computed first with a single‐phase method assuming constant pressure and zero stress on the interface. The second step is to compute the air flow assuming the free surface as a moving immersed boundary (IB). The IB method developed for Cartesian grids (Annu. Rev. Fluid Mech. 2005; 37 :239–261) is extended to curvilinear grids, where no‐slip and continuity conditions are used to enforce velocity and pressure boundary conditions for the air flow. The forcing points close to the IB can be computed and corrected under a sharp interface condition, which makes the computation very stable. The overset implementation is similar to that of the single‐phase solver (Comput. Fluids 2007; 36 :1415–1433), with the difference that points in water are set as IB points even if they are fringe points. Pressure–velocity coupling through pressure implicit with splitting of operators or projection methods is used for water computations, and a projection method is used for the air. The method on each fluid is a single‐phase method, thus avoiding ill‐conditioned numerical systems caused by large differences of fluid properties between air and water. The computation is only slightly slower than the single‐phase version, with complete absence of spurious velocity oscillations near the free surface, frequently present in fully coupled approaches. Validations are performed for laminar Couette flow over a wavy boundary by comparing with the analytical solution, and for the surface combatant model David Taylor Model Basin (DTMB) 5512 by comparing with Experimental Fluid Dynamics (EFD) and the results of two‐phase level set computations. Complex flow computations are demonstrated for the ONR Tumblehome DTMB 5613 with superstructure subject to waves and wind, including 6DOF motions and broaching in SS7 irregular waves and wind. Copyright © 2008 John Wiley & Sons, Ltd.     

Adaptive Runge–Kutta discontinuous Galerkin method for complex geometry problems on Cartesian grid

A Cartesian grid method using immersed boundary technique to simulate the impact of body in fluid has become an important research topic in computational fluid dynamics because of its simplification, automation of grid generation, and accuracy of results. In the frame of Cartesian grid, one often uses finite volume method with second order accuracy or finite difference method. In this paper, an h‐adaptive Runge–Kutta discontinuous Galerkin (RKDG) method on Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is developed. A ghost cell immersed boundary treatment with the modification of normal velocity is presented. The method is validated versus well documented test problems involving both steady and unsteady compressible flows through complex bodies over a wide range of Mach numbers. The numerical results show that the present boundary treatment to some extent reduces the error of entropy and demonstrate the efficiency, robustness, and versatility of the proposed approach. Copyright © 2013 John Wiley & Sons, Ltd.     

浸入边界法及其应用

浸入边界法主要用于模拟存在复杂外形结构的流场的运动情况和处理各种动边界问题,目前已广泛应用于计算流体力学领域.浸入边界法既是数学建模方法又是数值离散方法,它将物体边界与流体的相互作用通过在流体运动方程中加体积力项来体现,并在数值计算中采用简单的笛卡尔网格,避免了按照物体边界形状生成贴体网格时所遇到的各种问题.浸入边界法分为连续力法和离散力法:连续力法主要用于处理弹性边界问题,它的力源项满足特定的力学关系式;离散力法主要用于处理固体界面问题,它的力源项由边界条件推导得到.着重阐述了浸入边界法的基本原理和数学构造,对目前已有的几种不同的浸入边界法做了简单地介绍,并给出了一些应用实例,最后提出了浸入边界法未来的发展方向.     

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.     

Comparison of body‐fitted,embedded and immersed solutions of low Reynolds‐number 3‐D incompressible flows

The solutions obtained for low Reynolds‐number incompressible flows using the same flow solver and solution technique on body‐fitted, embedded surface and immersed body grids of similar size are compared. The cases considered are a sphere at Re = 100 and an idealized stented aneurysm. It is found that the solutions using all these techniques converge to the same grid‐independent solution. On coarser grids, the effect of higher‐order boundary conditions is noticeable. Therefore, if the manual labor required to set up a body‐fitted domain is excessive (as is often the case for patient‐specific geometries with medical devices), and/or computing resources are plentiful, the embedded surface and immersed body approaches become very attractive options. Copyright © 2007 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.     

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.     

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.     

Numerical simulation of unsteady multidimensional free surface motions by level set method

This paper presents a numerical method that couples the incompressible Navier–Stokes equations with the level set method in a curvilinear co‐ordinate system for study of free surface flows. The finite volume method is used to discretize the governing equations on a non‐staggered grid with a four‐step fractional step method. The free surface flow problem is converted into a two‐phase flow system on a fixed grid in which the free surface is implicitly captured by the zero level set. We compare different numerical schemes for advection of the level set function in a generalized curvilinear format, including the third order quadratic upwind interpolation for convective kinematics (QUICK) scheme, and the second and third order essentially non‐oscillatory (ENO) schemes. The level set equations of evolution and reinitialization are validated with benchmark cases, e.g. a stationary circle, a rotating slotted disk and stretching of a circular fluid element. The coupled system is then applied to a travelling solitary wave, and two‐ and three‐dimensional dam breaking problems. Some interesting free surface phenomena are revealed by the computational results, such as, the large free surface vortices, air entrapment and splashing of the water surge front. The computational results are in excellent agreement with theoretical predictions and experimental data, where they are available. Copyright © 2003 John Wiley & Sons, Ltd.     

An embedded boundary framework for compressible turbulent flow and fluid–structure computations on structured and unstructured grids

The finite volume method with exact two‐phase Riemann problems (FIVER) is a two‐faceted computational method for compressible multi‐material (fluid–fluid, fluid–structure, and multi‐fluid–structure) problems characterized by large density jumps, and/or highly nonlinear structural motions and deformations. For compressible multi‐phase flow problems, FIVER is a Godunov‐type discretization scheme characterized by the construction and solution at the material interfaces of local, exact, two‐phase Riemann problems. For compressible fluid–structure interaction (FSI) problems, it is an embedded boundary method for computational fluid dynamics (CFD) capable of handling large structural deformations and topological changes. Originally developed for inviscid multi‐material computations on nonbody‐fitted structured and unstructured grids, FIVER is extended in this paper to laminar and turbulent viscous flow and FSI problems. To this effect, it is equipped with carefully designed extrapolation schemes for populating the ghost fluid values needed for the construction, in the vicinity of the fluid–structure interface, of second‐order spatial approximations of the viscous fluxes and source terms associated with Reynolds averaged Navier–Stokes (RANS)‐based turbulence models and large eddy simulation (LES). Two support algorithms, which pertain to the application of any embedded boundary method for CFD to the robust, accurate, and fast solution of FSI problems, are also presented in this paper. The first one focuses on the fast computation of the time‐dependent distance to the wall because it is required by many RANS‐based turbulence models. The second algorithm addresses the robust and accurate computation of the flow‐induced forces and moments on embedded discrete surfaces, and their finite element representations when these surfaces are flexible. Equipped with these two auxiliary algorithms, the extension of FIVER to viscous flow and FSI problems is first verified with the LES of a turbulent flow past an immobile prolate spheroid, and the computation of a series of unsteady laminar flows past two counter‐rotating cylinders. Then, its potential for the solution of complex, turbulent, and flexible FSI problems is also demonstrated with the simulation, using the Spalart–Allmaras turbulence model, of the vertical tail buffeting of an F/A‐18 aircraft configuration and the comparison of the obtained numerical results with flight test data. Copyright © 2014 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.     

Assessment of regularized delta functions and feedback forcing schemes for an immersed boundary method

We present an improved immersed boundary method for simulating incompressible viscous flow around an arbitrarily moving body on a fixed computational grid. To achieve a large Courant–Friedrichs–Lewy number and to transfer quantities between Eulerian and Lagrangian domains effectively, we combined the feedback forcing scheme of the virtual boundary method with Peskin's regularized delta function approach. Stability analysis of the proposed method was carried out for various types of regularized delta functions. The stability regime of the 4‐point regularized delta function was much wider than that of the 2‐point delta function. An optimum regime of the feedback forcing is suggested on the basis of the analysis of stability limits and feedback forcing gains. The proposed method was implemented in a finite‐difference and fractional‐step context. The proposed method was tested on several flow problems, including the flow past a stationary cylinder, inline oscillation of a cylinder in a quiescent fluid, and transverse oscillation of a circular cylinder in a free‐stream. The findings were in excellent agreement with previous numerical and experimental results. Copyright © 2008 John Wiley & Sons, Ltd.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.