An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

A kernel gradient free (KGF) SPH method

The finite particle method (FPM) is a modified SPH method with high order accuracy while retaining the advantages of SPH in modeling problems with free surfaces, moving interfaces, and large deformations. In both SPH and FPM, kernel gradient is necessary in kernel and particle approximation of a field function and its derivatives. In this paper, a new FPM is presented, which only involves kernel function itself in kernel and particle approximation. The kernel gradient is not necessary in the whole computation, and this approach is thus referred to as a kernel gradient free (KGF) SPH method. This is helpful when a kernel function is not differentiable or the resultant kernel gradient is not sufficiently smooth, and thus it is more general in selecting a kernel function. Moreover, different from the original FPM with an asymmetric corrective matrix, in the new FPM, the resultant corrective matrix is symmetric, and this is advantageous in particle approximations. A series of numerical examples have been conducted to show the efficiencies of KGF‐SPH including one‐dimensional mathematical tests of polynomial functions with equal or variable smoothing length and two‐dimensional incompressible fluid flow of shear cavity. It is found that KGF‐SPH is comparable with FPM in accuracy and is flexible as SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Development of highly accurate interpolation method for mesh‐free flow simulations II. Application of CIVA method to incompressible flow simulations

This is the second report on the development of a highly accurate interpolation method, which is called cubic interpolation with volume/area (CIVA) co‐ordinates, for mesh‐free flow simulations. In this paper, the method of determining the c‐parameter of CIVA using a constant curvature condition is first considered for the two‐ and three‐dimensional cases. A computation of a three‐dimensional passive scalar advection problem is performed for accuracy verification and for comparison with widely used methods. Then, an application algorithm of the CIVA method respecting incompressible fluid simulation is presented. As the incompressible condition based on Lagrangian approaches causes problems, in this paper we consider the condition based on the conventional Eulerian approach. The CIVA‐based incompressible flow simulation algorithm enables a highly accurate simulation of many kinds of problems that have complicated geometries and involve complicated phenomena. To confirm the facts, numerical analyzes are executed for some benchmark problems, namely flow in a square cavity, free surface sloshing and moving boundary problems in complex geometries. The results show that the method achieves high accuracy and has high flexibility, even for the flows involving high Reynolds number, complicated geometries, moving boundaries and free surfaces. Copyright © 2000 John Wiley & Sons, Ltd.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Enhancement of the accuracy of the finite volume particle method for the simulation of incompressible flows

A finite volume particle (FVP) method for simulation of incompressible flows that provides enhanced accuracy is proposed. In this enhanced FVP method, a dummy neighbor particle is introduced for each particle in the calculation and used for the discretization of the gradient model and Laplacian model. The error‐compensating term produced by introducing the dummy neighbor particle enables higher order terms to be calculated. The proposed gradient model and Laplacian model are applied in both pressure and pressure gradient calculations. This enhanced FVP scheme provides more accurate simulations of incompressible flows. Several 2‐dimensional numerical simulations are given to confirm its enhanced performance.     

Consistent inlet and outlet boundary conditions for particle methods

In this paper, simple and consistent open boundary conditions are presented for the numerical simulation of viscous incompressible laminar flows. The present approach is based on an arbitrary Lagrangian-Eulerian particle method using upwind interpolation. Three kinds of inlet/outlet boundary conditions are proposed for particle methods, a pressure specified inlet/outlet condition, a velocity profile specified inlet/outlet condition, and a fully developed flow outlet condition. These inlet/outlet conditions are realized by using boundary particles and modification to the physical value such as velocity. Poiseuille flows, flows over a backward-facing step, and flows in a T-shape branch are calculated. The results are compared with those of mesh-based methods such as the finite volume method. The method presented herein exhibits accuracy and numerical stability.     

A finite element method for compressible viscous fluid flows

This work presents a mixed three‐dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a ‘stable’ numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf–sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady‐state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright © 2010 John Wiley & Sons, Ltd.     

A stable and accurate projection method on a locally refined staggered mesh

In this paper, a robust projection method on a locally refined mesh is proposed for two‐ and three‐dimensional viscous incompressible flows. The proposed method is robust not only when the interface between two meshes is located in a smooth flow region but also when the interface is located in a flow region with large gradients and/or strong unsteadiness. In numerical simulations, a locally refined mesh saves many grid points in regions of relatively small gradients compared with a uniform mesh. For efficiency and ease of implementation, we consider a two‐level blocked structure, for which both of the coarse and fine meshes are uniform Cartesian ones individually. Unfortunately, the introduction of the two‐level blocked mesh results in an important but difficult issue: coupling of the coarse and fine meshes. In this paper, by properly addressing the issue of the coupling, we propose a stable and accurate projection method on a locally refined staggered mesh for both two‐ and three‐dimensional viscous incompressible flows. The proposed projection method is based on two principles: the linear interpolation technique and the consistent discretization of both sides of the pressure Poisson equation. The proposed algorithm is straightforward owing to the linear interpolation technique, is stable and accurate, is easy to extend from two‐ to three‐dimensional flows, and is valid even when flows with large gradients cross the interface between the two meshes. The resulting pressure Poisson equation is non‐symmetric on a locally refined mesh. The numerical results for a series of exact solutions for 2D and 3D viscous incompressible flows verify the stability and accuracy of the proposed projection method. The method is also applied to some challenging problems, including turbulent flows around particles, flows induced by impulsively started/stopped particles, and flows induced by particles near solid walls, to test the stability and accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

Incompressible smoothed particle hydrodynamics simulation of multifluid flows

This paper presents an incompressible SPH (ISPH) technique to simulate multifluid flows. The SPH method is a mesh‐free particle modeling approach that can treat free surfaces and multi‐interfaces in a simple and efficient manner. The ISPH method employs an incompressible hydrodynamic formulation to solve the fluid pressure that ensures a stable pressure field. Two multifluid ISPH models are proposed following different interface treatments: the coupled ISPH model does not distinguish the different fluid phases and applies the standard ISPH technique across the interface, whereas the decoupled ISPH model first treats each fluid phase separately and then couples the different phases by applying pressure and shear stress continuities across the interface. The two proposed models were used to investigate a gravity underflow with a low density ratio in a Generalized Reservoir Hydrodynamics (GRH) flume and a horizontal lock exchange flow with a high density ratio. Comparisons with data and relevant numerical error analysis indicated that the decoupled model performed well in cases of both low and high density ratios because of the accurate treatment of interface boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

A volume-of-fluid method for incompressible free surface flows

This paper proposes a hybrid volume-of-fluid (VOF) level-set method for simulating incompressible two-phase flows. Motion of the free surface is represented by a VOF algorithm that uses high resolution differencing schemes to algebraically preserve both the sharpness of interface and the boundedness of volume fraction. The VOF method is specifically based on a simple order high resolution scheme lower than that of a comparable method, but still leading to a nearly equivalent order of accuracy. Retaining the mass conservation property, the hybrid algorithm couples the proposed VOF method with a level-set distancing algorithm in an implicit manner when the normal and the curvature of the interface need to be accurate for consideration of surface tension. For practical purposes, it is developed to be efficiently and easily extensible to three-dimensional applications with a minor implementation complexity. The accuracy and convergence properties of the method are verified through a wide range of tests: advection of rigid interfaces of different shapes, a three-dimensional air bubble's rising in viscous liquids, a two-dimensional dam-break, and a three-dimensional dam-break over an obstacle mounted on the bottom of a tank. The standard advection tests show that the volume advection algorithm is comparable in accuracy with geometric interface reconstruction algorithms of higher accuracy than other interface capturing-based methods found in the literature. The numerical results for the remainder of tests show a good agreement with other numerical solutions or available experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

A pressure-based Mach-uniform method for viscous fluid flows

A pressure-based, Mach-uniform method is developed by combining the SLAU2 numerical scheme and the higher temporal order pressure-based algorithm. This hybrid combination compensates the limitation of the SLAU2 numerical scheme in the low-Mach number regime and deficiencies of the pressure-based method in the high-Mach number regime. A momentum interpolation method is proposed to replace the Rhie-Chow interpolation for accurate shock-capturing and to alleviate the carbuncle phenomena. The momentum interpolation method is consistent in addition to preserving pressure–velocity coupling in the incompressible limit . The postulated pressure equation allows the algorithm to compute the subsonic flows without empirical scaling of numerical dissipation at low-Mach number computation. Several test cases involving a broad range of Mach number regimes are presented. The numerical results demonstrate that the present algorithm is remarkable for the calculation of viscous fluid flows at arbitrary Mach number including shock wave/laminar boundary layer interaction and aerodynamics heating problem.     

Introduction of turbulent model in a mixed finite volume/finite element method

The aim of this work is to present a new numerical method to compute turbulent flows in complex configurations. With this in view, a k-? model with wall functions has been introduced in a mixed finite volume/finite element method. The numerical method has been developed to deal with compressible flows but is also able to compute nearly incompressible flows. The physical model and the numerical method are first described, then validation results for an incompressible flow over a backward-facing step and for a supersonic flow over a compression ramp are presented. Comparisons are performed with experimental data and with other numerical results. These simulations show the ability of the present method to predict turbulent flows, and this method will be applied to simulate complex industrial flows (flow inside the combustion chamber of gas turbine engines). The main goal of this paper is not to test turbulence models, but to show that this numerical method is a solid base to introduce more sophisticated turbulence model.     

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.     

Development of a hybrid particle‐mesh method for two‐phase flow simulations

A hybrid particle‐mesh method was developed for efficient and accurate simulations of two‐phase flows. In this method, the main component of the flow is solved using the constrained interpolated profile/multi‐moment finite volumemethod; the two‐phase interface is rendered using the finite volume particle (FVP) method. The effect of surface tension is evaluated using the continuum surface force model. Numerical particles in the FVP method are distributed only on the surface of the liquid in simulating the interface between liquid and gas; these particles are used to determine the density of each mesh grid. An artificial term was also introduced to mitigate particle clustering in the direction of maximum compression and sparse discretization errors in the stretched direction. This enables accurate interface tracking without diminishing numerical efficiency. Two benchmark simulations are used to demonstrate the validity of the method developed and its numerical stability. Copyright © 2016 John Wiley & Sons, Ltd.     

一种考虑界面不连续的改进的有限粒子法

有限粒子法（finite particle method，FPM）作为SPH（smoothed particle hydrodynamics）方法的重要改进，有效提高了边界区域粒子的近似精度，但是当FPM处理多物理场时，在不连续界面附近的计算精度会大大降低，并且FPM必须满足的矩阵非奇异性也提高了对界面处理的要求。本文中基于DSPH（discontinuous SPH）方法，提出了一种考虑界面不连续的改进FPM—DSFPM（discontinuous special FPM）法，旨在改善FPM在界面不连续处的计算精度，从而进一步提高其计算效率和稳定性。首先，分析了DSFPM的核近似精度。其次，根据不同的工程问题，给出DSFPM处理小变形和大变形问题的算法流程。利用DSFPM、DSPH和FPM等3种方法对弹性铝块小变形碰撞冲击算例进行了模拟，通过对比分析铝块的速度和应力以及计算时间验证了DSFPM算法在非连续界面处计算精度和计算效率的优势。最后，通过结合DSFPM和DFPM（discontinuous FPM）实现了对于大变形问题的模拟。     

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 stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

A conservative flux-corrected continuous finite element method for fluid interface dynamics

This paper presents a continuous finite element solution for fluid flows with interfaces. The method is founded on the sign-preserving flux correction transport methodology and extends nonoscillatory finite element algorithm capabilities to predict interface motion efficiently. The procedure is composed of three main stages, along the lines of the conservative level set method: transport of phase function, reconstruction of phase function, and solution of equations of motion of two incompressible fluids. The flux correction technique takes action on the three steps. Limiting process incorporates a straightforward refinement to remove global mass residuals present in the earliest version of the algorithm. This is of particular importance in the transport step. Moreover, new method retains the efficacy of the original. To reconstruct the phase function after transport, a novel nonlinear (and conservative) streamlined diffusion equation is proposed, with an anisotropic diffusivity comprising artificial compression and diffusive fluxes along interface displacements direction. A substantial reduction of unphysical overshoots along the interface is reached by an improved bound estimation that includes interface information. Complete operation of the correction algorithm for two incompressible fluids flows requires two pressure solutions. We explore a reduced form to circumvent this extra burden. Numerical experiments verify the formulation by reproducing stringent benchmarks both for transport/reinitialization and for two-fluid interface propagation.