A simple SPH algorithm for multi‐fluid flow with high density ratios

In this paper, we describe an SPH algorithm for multi‐fluid flow, which is efficient, simple and robust. We derive the inviscid equations of motion from a Lagrangian together with the constraint provided by the continuity equation. The viscous flow equations then follow by adding a viscous term. Rigid boundaries are simulated using boundary force particles in a manner similar to the immersed boundary method. Each fluid is approximated as weakly compressible with a speed of sound sufficiently large to guarantee that the relative density variations are typically 1%. When the SPH force interaction is between two particles of different fluids, we increase the pressure terms. This simple procedure stabilizes the interface between the fluids. The equations of motion are integrated using a time stepping rule based on a second‐order symplectic integrator. When linear and angular momentum should be conserved exactly, they are conserved to within round‐off errors. We test the algorithm by simulating a variety of problems involving fluids with a density ratio in the range 1–1000. The first of these is a free surface problem with no rigid boundaries. It involves the flow of an elliptical distribution with one fluid inside the other. We show that the simulations converge as the particle spacing decreases, and the results are in good agreement with the exact inviscid, incompressible theory. The second test is similar to the first but involves the nonlinear oscillation of the fluids. As in the first test, the agreement with theory is very good, and the method converges. The third test is the simulation of waves at the interface between two fluids. The method is shown to converge, and the agreement with theory is satisfactory. The fourth test is the Rayleigh–Taylor instability for a configuration considered by other authors. Key parameters are shown to converge, and the agreement with other authors is good. The fifth and final test is how well the SPH method simulates gravity currents with density ratios in the range 2–30. The results of these simulations are in very good agreement with those of other authors and in satisfactory agreement with experimental results.Copyright © 2012 John Wiley & Sons, Ltd.     

The capturing of free surfaces in incompressible multi‐fluid flows

By treating it as a contact discontinuity in the density field, a free surface between two immiscible fluids can be automatically ‘captured’ by the enforcement of conservation laws. A surface‐capturing method of this kind requires no special tracking or fitting treatment for the free surface, thereby offering the advantage of algorithm simplicity over the surface‐tracking or the surface‐fitting method. A surface‐capturing method based on a new multi‐fluid incompressible Navier–Stokes formulation is developed. It is applied to a variety of free‐surface flows, including the Rayleigh–Taylor instability problem, the ship waves around a Wigley hull and a model bubble‐rising problem to demonstrate the validity and versatility of the present method. Copyright © 2000 John Wiley & Sons, Ltd.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

Buoyancy modelling with incompressible SPH for laminar and turbulent flows

This work aims to model buoyant, laminar or turbulent flows, using a two‐dimensional incompressible smoothed particle hydrodynamics model with accurate wall boundary conditions. The buoyancy effects are modelled through the Boussinesq approximation coupled to a heat equation, which makes it possible to apply an incompressible algorithm to compute the pressure field from a Poisson equation. Based on our previous work [1], we extend the unified semi‐analytical wall boundary conditions to the present model. The latter is also combined to a Reynolds‐averaged Navier–Stokes approach to treat turbulent flows. The k ? ? turbulence model is used, where buoyancy is modelled through an additional term in the k ? ? equations like in mesh‐based methods. We propose a unified framework to prescribe isothermal (Dirichlet) or to impose heat flux (Neumann) wall boundary conditions in incompressible smoothed particle hydrodynamics. To illustrate this, a theoretical case is presented (laminar heated Poiseuille flow), where excellent agreement with the theoretical solution is obtained. Several benchmark cases are then proposed: a lock‐exchange flow, two laminar and one turbulent flow in differentially heated cavities, and finally a turbulent heated Poiseuille flow. Comparisons are provided with a finite volume approach using an open‐source industrial code. Copyright © 2015 John Wiley & Sons, Ltd.     

A remedy for numerical oscillations in weakly compressible smoothed particle hydrodynamics

Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) can lead to non‐physical oscillations in the pressure and density fields when simulating incompressible flow problems. This in turn may result in tensile instability and sometimes divergence. In this paper, it is shown that this difficulty originates from the specific form of spatial discretization used for the pressure term when solving the mass conservation equation. After describing the pressure–velocity decoupling problem associated with the so‐called colocated grid methods, a modified approach is presented that overcomes this problem using a different discretization scheme for the second derivative of pressure. The modified scheme is employed for solving a number of benchmark problems including both single‐phase and two‐phase test cases. Copyright © 2010 John Wiley & Sons, Ltd.     

A consistent incompressible SPH method for internal flows with fixed and moving boundaries

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first‐order consistent discretization schemes both for the first‐order and second‐order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two‐dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data. It is shown that using smaller smoothing lengths, the proposed method can provide desirable accuracies with relatively less computational cost for two‐dimensional problems. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled WC-TL SPH method for simulation of hydroelastic problems

A coupled weakly compressible (WC) and total Lagrangian (TL) smoothed particle hydrodynamics (SPH) method is developed for simulating hydroelastic problems. The fluid phase is simulated using WCSPH method, while the structural dynamics are solved using TLSPH method. Fluid and solid components of the method are validated separately. A sloshing water tank problem is solved to test the WCSPH method while oscillation of a thin plate and large deformation of a cantilever beam are simulated to test the TLSPH method. After validating each component, the coupled WC-TL SPH scheme is used to simulate two benchmark hydroelastic problems. The first test case shows the evolution of water column with an elastic boundary gate, and the second one investigates the breaking water column impact on elastic structures. The agreement between WC-TL SPH results and literature data shows the ability of the proposed method in simulating hydroelastic phenomena.     

A consistent and fast weakly compressible smoothed particle hydrodynamics with a new wall boundary condition

A modified weakly compressible smoothed particle hydrodynamics (WCSPH) is presented, which utilizes consistent discretization schemes for spatial derivatives in the flow equations. Here, each SPH particle is considered as a computational point that represents a specific part of the fluid. To overcome non‐physical oscillations that usually arise in standard WCSPH, we modified the mass conservation equation by using a numerical filter. This modification is based on the difference between two discretization schemes used for the term . Furthermore, a new implementation of wall boundary condition in SPH is introduced. This condition is imposed on the pressure of wall boundary particles to ensure that the acceleration of each boundary particle in normal direction to the wall is zero. Thus, no penetration through walls will occur. To examine the performance of the modified method, we solved a series of two‐dimensional incompressible internal flow benchmark problems. By comparing the result with analytical solutions and the results of the standard WCSPH, we show that the use of consistent schemes in conjunction with the proposed numerical filter improves both accuracy and speed of the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of a vortex method to free surface flows

Vortex methods have found wide applications in various practical problems. The use of vortex methods in free surface flow problems, however, is still very limited. This paper demonstrates a vortex method for practical computation of non-linear free surface flows produced by moving bodies. The method is a potential flow formulation which uses the exact non-linear free surface boundary condition at the exact location of the instantaneous free surface. The position of the free surface, on which vortices are distributed, is updated using a Lagrangian scheme following the fluid particles on the free surface. The vortex densities are updated by the non-linear dynamic boundary condition, derived from the Euler equations, with an iterative Lagrangian numerical scheme. The formulation is tested numerically for a submerged circular cylinder in unsteady translation. The iteration is shown to converge for all cases. The results of the unsteady simulations agree well with classical linearized solutions. The stability of the method is also discussed.     

A new finite element formulation for both compressible and nearly incompressible fluid dynamics

A new finite element formulation designed for both compressible and nearly incompressible viscous flows is presented. The formulation combines conservative and non‐conservative dependent variables, namely, the mass–velocity (density * velocity), internal energy and pressure. The central feature of the method is the derivation of a discretized equation for pressure, where pressure contributions arising from the mass, momentum and energy balances are taken implicitly in the time discretization. The method is applied to the analysis of laminar flows governed by the Navier–Stokes equations in both compressible and nearly incompressible regimes. Numerical examples, covering a wide range of Mach number, demonstrate the robustness and versatility of the new method. Copyright © 2000 John Wiley & Sons, Ltd.     

A numerical investigation of 2d,steady free surface flows

Systematic tests have been performed to study the behaviour of a numerical method developed to calculate 2D, steady free surface flows. The Reynolds equations are solved in the physical space by employing a non–orthogonal staggered grid, while the k-ε model is adopted to approximate the Reynolds stresses. The free surface is calculated following an iterative procedure and various parameters that affect convergence and accuracy of the numerical solution have been examined. Calculated results are compared with measured data for two cases, i.e. the wave generation above a bottom topography at various Froude numbers and the free surface formation above a submerged hydrofoil. © 1997 John Wiley & Sons, Ltd.     

The robustness issue on multigrid schemes applied to the Navier–Stokes equations for laminar and turbulent,incompressible and compressible flows

The paper's leitmotiv is condensed in one word: robustness. This is a real hindrance for the successful implementation of any multigrid scheme for solving the Navier–Stokes set of equations. In this paper, many hints are given to improve this issue. Instead of looking for the best possible speed‐up rate for a particular set of problems, at a given regime and in a given condition, the authors propose some ideas pursuing reasonable speed‐up rates in any situation. In a previous paper, the authors presented a multigrid method for solving the incompressible turbulent RANS equations, with particular care in the robustness and flexibility of the solution scheme. Here, these concepts are further developed and extended to compressible laminar and turbulent flows. This goal is achieved by introducing a non‐linear multigrid scheme for compressible laminar (NS equations) and turbulent flow (RANS equations), taking benefit of a convenient master–slave implementation strategy. Copyright © 2004 John Wiley & Sons, Ltd.     

A finite element model for free surface flows on fixed meshes

In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

SPH方法Delaunay三角刨分与自由液面重构

光滑粒子法(SPH)作为一种拉格朗日型无网格方法,兼具欧拉网格方法和拉格朗日网格方法的优势,已经成功应用于科学和工程的众多领域。SPH方法后处理一般基于无规则分布的粒子,不如网格类方法后处理简便、直接。另外,SPH方法模拟自由液面流动等问题时,通过粒子位置难以重构自由液面的准确位置。发展一种基于Delaunay三角刨分的SPH后处理方法,即先基于SPH粒子位置利用Delaunay三角刨分建立三角网格,然后将粒子信息转化成网格单元/节点信息,从而可以在三角网格上进行后处理,实现基于网格方法的后处理功能,并可以在三角网格上直接提取或重构自由液面。将本文的方法应用到液滴碰撞和溃坝流SPH模拟结果的后处理中,得到了非常好的结果,表明本文的方法有效可靠。     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.     

A corrected symmetric SPH method to simulate viscoelastic free surface flows based on the PTT model

In this work, a corrected symmetric and periodic density reinitialized SPH (CSPDR‐SPH) method is proposed and extended to simulate the viscoelastic free surface flows based on the Phan–Thien–Tanner model. The improvements mainly lie in deriving a corrected symmetric kernel gradient, and combining it with a periodic density reinitialization procedure. In addition, a simple artificial viscosity and a simple artificial stress form are adopted. Thus, the CSPDR‐SPH method has higher accuracy and better stability than the SPH method, and conserves both linear and angular momentums. The consistency and convergence of the CSPDR‐SPH method are justified by approximating a function in one and two dimensions. The merits of CSPDR‐SPH method are demonstrated by several benchmarks. The simple flow in a two‐dimensional channel is investigated to show the capability of the CSPDR‐SPH method to simulate the viscoelastic free surface flow. Then the CSPDR‐SPH method is extended to simulate the impacting drop problem. Numerical results show that the CSPDR‐SPH method can precisely capture the viscoelastic free surface. The Reynolds number, Weissenberg number and elongation parameter have remarkable effect on the flows. Copyright © 2012 John Wiley & Sons, Ltd.     

An unsteady single‐phase level set method for viscous free surface flows

The single‐phase level set method for unsteady viscous free surface flows is presented. In contrast to the standard level set method for incompressible flows, the single‐phase level set method is concerned with the solution of the flow field in the water (or the denser) phase only. Some of the advantages of such an approach are that the interface remains sharp, the computation is performed within a fluid with uniform properties and that only minor computations are needed in the air. The location of the interface is determined using a signed distance function, and appropriate interpolations at the fluid/fluid interface are used to enforce the jump conditions. A reinitialization procedure has been developed for non‐orthogonal grids with large aspect ratios. A convective extension is used to obtain the velocities at previous time steps for the grid points in air, which allows a good estimation of the total derivatives. The method was applied to three unsteady tests: a plane progressive wave, sloshing in a two‐dimensional tank, and the wave diffraction problem for a surface ship, and the results compared against analytical solutions or experimental data. The method can in principle be applied to any problem in which the standard level set method works, as long as the stress on the second phase can be specified (or neglected) and no bubbles appear in the flow during the computation. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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

A locally DIV‐free projection scheme for incompressible flows based on non‐conforming finite elements

We present a projection scheme whose end‐of‐step velocity is locally pointwise divergence free, using a continuous ?₁ approximation for the velocity in the momentum equation, a first‐order Crouzeix–Raviart approximation at the projection step, and a ?₀ approximation for the pressure in both steps. The analysis of the scheme is done only for grids that guarantee the existence of a divergence free conforming ?₁ interpolant for the velocity. Optimal estimates for the velocity error in L²‐ and H¹‐norms are deduced. The numerical results demonstrate that these estimates should also hold on grids on which the continuous ?₁ approximation for the velocity locks. Since the end‐of‐step velocity is locally solenoidal, the scheme is recommendable for problems requiring good mass conservation. Copyright © 2005 John Wiley & Sons, Ltd.