An SPH pressure correction algorithm for multiphase flows with large density ratio

In this paper, we propose an interfacial pressure correction algorithm for smoothed particle hydrodynamics (SPH) simulation of multiphase flows with large density ratios. This correction term is based on the assumption of small deformation of the interface, and derived from perturbation expansion analysis. It is also proven to be applicable in cases with complex interfaces. This correction algorithm helps to overcome the discontinuities of the pressure gradient over the interfaces, which may cause unphysical gap between different phases. This proposed correction algorithm is implemented on a recent multiphase SPH model, which is based on the assumption of pressure continuity over the interfaces. The coupled dynamic solid boundary treatment is used to simulate solid walls; and a cut‐off pressure is applied to avoid negative particle pressure, which may cause computational instabilities in SPH. Three numerical examples of air–water flows, including sloshing, dam breaking, and water entry, are presented and compared with experimental data, indicating the robustness of our pressure correction algorithm in multiphase simulations with large density ratios. Copyright © 2015 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.     

一种基于黎曼解处理大密度比多相流SPH的改进算法

多相流界面存在密度、黏性等物理场间断,直接采用传统光滑粒子水动力学(smoothedparticle hydrodynamics,SPH)方法进行数值模拟,界面附近的压力和速度存在震荡.一套基于黎曼解能够处理大密度比的多相流SPH计算模型被提出,该模型利用黎曼解在处理接触间断问题方面的优势,将黎曼解引入到SPH多相流计算模型中,为了能够准确求解多相流体物理黏性、减小黎曼耗散,对黎曼形式的SPH动量方程进行了改进,又将Adami固壁边界与黎曼单侧问题相结合来施加多相流SPH固壁边界,同时模型中考虑了表面张力对小尺度异相界面的影响,该模型没有添加任何人工黏性、人工耗散和非物理人工处理技术,能够反应多相流真实物理黏性和物理演变状态.采用该模型首先对三种不同粒子间距离散下方形液滴震荡问题进行了数值模拟,验证了该模型在处理异相界面的正确性和模型本身的收敛性;后又通过对Rayleigh--Taylor不稳定、单气泡上浮、双气泡上浮问题进行了模拟计算,结果与文献对比吻合度高,异相界面捕捉清晰,结果表明,本文改进的多相流SPH模型能够稳定、有效的模拟大密度比和黏性比的多相流问题.      

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.     

SPH固壁边界处理方法的改进

提出一种适用于光滑质点水动力学（SPH）方法的改进的边界处理方法。在这种方法中,边界粒子的压力可通过其周围的流体粒子的压力插值得到,从而改进了耦合边界法在边界上压力不准的问题。运用这种改进的边界处理方法模拟了二维方形水箱中的非线性晃荡问题以及二维楔形体自由入水问题。模拟结果与实验结果吻合较好,证明了此改进的边界处理方法是有效的。     

多相流动的光滑粒子流体动力学方法研究综述

光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)具有粒子方法的无网格和全拉格朗日特征, 适用于具有界面大变形、不连续性和多物理场的多相流的高精度模拟. SPH方法模拟多相流已有大量报道, 具体的实现方式也大不相同. 本文首先阐述了采用SPH方法模拟流体的基本控制方程, 以及求解过程中需要考虑的流体压力求解、表面张力、固体边界等问题. 整理和总结了基于SPH方法进行多相流模拟的主要实现方式: (1)双流体模型的拉格朗日求解器: 两相离散为两组独立SPH粒子, 并用显式相间作用耦合两相; (2)多相SPH方法: SPH方法对多相流模拟的自然延伸, 相间作用由SPH参数隐式描述; (3) SPH与其他离散方法的耦合: 差异较大的两相各自采用不同离散方法, 发挥不同拉格朗日方法的优点; (4) SPH和基于网格方法的耦合: 网格方法处理简单的单相流动主体, 获得精度和效率间的平衡. 另外, 还在模拟参数物理化等方面论述了与SPH方法模拟多相流相关的一些改进和修正方法, 并在最后讨论和建议了提高多相流SPH模拟效率和精度的措施.      

Comparative study of two corrective gradient models in the simulation of multiphase flows using moving particle semi-implicit method

In this study, two corrective gradient models (CGMs) are compared in the simulation of multiphase flows. Linear consistency of the gradient model of moving particle semi-implicit (MPS) method has been recovered by introducing corrective matrix. However, it is found that particles tend to disperse along the streamline while using the CGM proposed in a previous study. Particle shifting (PS) schemes are necessary to reduce the irregularity of particle distribution to stabilize the calculation. To enhance the accuracy and stability, another CGM with dummy particle (CGMD) was proposed in our previous study. This enhanced CGM is featured by linear consistency and purely repulsive pressure gradient force. In this study, this enhanced CGM is modified and applied to multiphase flow simulation. Comparative study suggests that the modified CGM with PS scheme is capable of calculating various multiphase flows and predicting the interface evolution both clearly and accurately.     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

An enhanced treatment of boundary conditions in implicit smoothed particle hydrodynamics简

In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

An SPH multi‐fluid model based on quasi buoyancy for interface stabilization up to high density ratios and realistic wave speed ratios

We introduce a smoothed particle hydrodynamics (SPH) concept for the stabilization of the interface between 2 fluids. It is demonstrated that the change in the pressure gradient across the interface leads to a force imbalance. This force imbalance is attributed to the particle approximation implicit to SPH. To stabilize the interface, a pressure gradient correction is proposed. In this approach, the multi‐fluid pressure gradients are related to the (gravitational and fluid) accelerations. This leads to a quasi‐buoyancy correction for hydrostatic (stratified) flows, which is extended to nonhydrostatic flows. The result is a simple density correction that involves no parameters or coefficients. This correction is included as an extra term in the SPH momentum equation. The new concept for the stabilization of the interface is explored in 5 case studies and compared with other multi‐fluid models. The first case is the stagnant flow in a tank: The interface remains stable up to density ratios of 1:1000 (typical for water and air), in combination with artificial wave speed ratios up to 1:4. The second and third cases are the Rayleigh‐Taylor instability and the rising bubble, where a reasonable agreement between SPH and level‐set models is achieved. The fourth case is an air flow across a water surface up to density ratios of 1:100, artificial wave speed ratios of 1:4, and high air velocities. The fifth case is about the propagation of internal gravity waves up to density ratios of 1:100 and artificial wave speed ratios of 1:4. It is demonstrated that the quasi‐buoyancy model may be used to stabilize the interface between 2 fluids up to high density ratios, with real (low) viscosities and more realistic wave speed ratios than achieved by other weakly compressible SPH multi‐fluid models. Real wave speed ratios can be achieved as long as the fluid velocities are not very high. Although the wave speeds may be artificial in many cases, correct and realistic wave speed ratios are essential in the modelling of heat transfer between 2 fluids (eg, in engineering applications such as gas turbines).     

Dynamic simulation of liquid-filled flexible multibody systems via absolute nodal coordinate formulation and SPH method

A new computation method is proposed to study the coupled dynamics of a partially liquid-filled flexible multibody system, where the liquid is modeled by using the Smoothed Particle Hydrodynamics (SPH) method and the flexible bodies are described by using the Absolute Nodal Coordinate Formulation (ANCF). Extra virtual particles are introduced and embedded in the liquid neighboring the rigid or flexible boundaries in order to prevent field particles from penetrating the boundary and force them to follow the deformation of flexible boundary. The interaction forces between the liquid and the flexible multibody system are transmitted by the virtual particles. The domain decomposition is used to improve the efficiency of interaction detection in SPH computation. A predictor-corrector scheme is used to solve the governing equations of liquid discretized by SPH particles. The generalized-alpha method based on sparse matrix storage skill is used to solve a huge set of dynamic equations of the multibody system. The OpenMP+OpenACC based parallel computation skills are embedded in the iteration processes to speed up the computation efficiency. Finally, three numerical examples are given to validate the proposed computation method.     

Simulating surface tension with smoothed particle hydrodynamics

A method for simulating two‐phase flows including surface tension is presented. The approach is based upon smoothed particle hydrodynamics (SPH). The fully Lagrangian nature of SPH maintains sharp fluid–fluid interfaces without employing high‐order advection schemes or explicit interface reconstruction. Several possible implementations of surface tension force are suggested and compared. The numerical stability of the method is investigated and optimal choices for numerical parameters are identified. Comparisons with a grid‐based volume of fluid method for two‐dimensional flows are excellent. The methods presented here apply to problems involving interfaces of arbitrary shape undergoing fragmentation and coalescence within a two‐phase system and readily extend to three‐dimensional problems. Boundary conditions at a solid surface, high viscosity and density ratios, and the simulation of free‐surface flows are not addressed. Copyright © 2000 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 comparative study of direct‐forcing immersed boundary‐lattice Boltzmann methods for stationary complex boundaries

In this study, we assess several interface schemes for stationary complex boundary flows under the direct‐forcing immersed boundary‐lattice Boltzmann methods (IB‐LBM) based on a split‐forcing lattice Boltzmann equation (LBE). Our strategy is to couple various interface schemes, which were adopted in the previous direct‐forcing immersed boundary methods (IBM), with the split‐forcing LBE, which enables us to directly use the direct‐forcing concept in the lattice Boltzmann calculation algorithm with a second‐order accuracy without involving the Navier–Stokes equation. In this study, we investigate not only common diffuse interface schemes but also a sharp interface scheme. For the diffuse interface scheme, we consider explicit and implicit interface schemes. In the calculation of velocity interpolation and force distribution, we use the 2‐ and 4‐point discrete delta functions, which give the second‐order approximation. For the sharp interface scheme, we deal with the exterior sharp interface scheme, where we impose the force density on exterior (solid) nodes nearest to the boundary. All tested schemes show a second‐order overall accuracy when the simulation results of the Taylor–Green decaying vortex are compared with the analytical solutions. It is also confirmed that for stationary complex boundary flows, the sharper the interface scheme, the more accurate the results are. In the simulation of flows past a circular cylinder, the results from each interface scheme are comparable to those from other corresponding numerical schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulations of wave interactions with vertical wave barriers using the SPH method

This paper focuses on the fluid boundary separation problem of the conventional dynamic solid boundary treatment (DSBT) and proposes a modified DSBT (MDSBT). Classic 2D free dam break flows and 3D dam break flows against a rectangular box are used to assess the performance of this MDSBT in free surface flow and violent fluid–structure interaction, respectively. Another test, water column oscillations in a U‐tube, is specially designed to reveal the applicability of dealing with two types of particular boundaries: the wet–dry solid boundary and the large‐curvature solid boundary. A comparison between the numerical results and the experimental data shows that the MDSBT is capable of eliminating the fluid boundary separation, improving the accuracy of the solid boundary pressure calculations and preventing the unphysical penetration of fluid particles. Using a 2D SPH numerical wave tank with MDSBT, the interactions between regular waves and a simplified vertical wave barrier are simulated. The numerical results reveal that the maximum horizontal force occurs at the endpoint of the vertical board, and with the enlargement of the relative submerged board length, the maximum moment grows linearly; furthermore, the relative average mass transportation under the breakwater initially increases to 11.14 per wave strike but is later reduced. The numerical simulation of a full‐scale 3D wave barrier with two vertical boards shows that the wave and structure interactions in the practical project are far more complicated than in the simplified 2D models. The SPH model using the MDSBT is capable of providing a reference for engineering designs. Copyright © 2014 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

改进的SPH边界处理方法与土体大变形模拟

沙土滑坡往往会造成重大的人身财产损失,研究这类土体大变形问题对防灾工程具有指导意义。光滑粒子流体动力学SPH（Smoothed Particle Hydrodynamics）方法是一种拉格朗日型无网格粒子法,十分适用于模拟大变形问题。在SPH方法中,合适的边界处理方法一直是个难点,传统的边界虚粒子法或排斥力法较难模拟复杂边界。本文引入了一种能处理任意形状边界的方法——统一半解析壁面边界条件处理方法USAW（unified semi-analytical wall boundary conditions）,通过在控制方程中引入修正因子并保留边界面积分项来弥补边界缺失。为了更准确模拟问题域边界,提出无质量边界粒子的新概念。利用该方法成功模拟了土体滑坡算例,验证了方法的可靠性,并避免了边界零粒子层问题。通过数值模拟,分析了内摩擦角和黏聚力等土体物性参数对滑坡过程的影响。最后,应用该方法研究了滑坡冲击楔形体时的压力响应。     

An integrated smoothed particle hydrodynamics model for complex interfacial flows with large density ratios

In this paper, an integrated smoothed particle hydrodynamics (SPH) model for complex interfacial flows with large density ratios is developed. The discrete continuity equation and acceleration equation are obtained by considering the time derivative of the volume of particle and Eckart's continuum Lagrangian equation. A continuum surface force model is used to meet the fact that surface force may not be distributed uniformly on each side of the interface. An improved boundary condition is imposed to model wall free-slip and no-slip condition for interfacial flows with large density ratios. Particle shifting algorithm (PSA) is added for interfacial flows by imposing the normal correction near the interface, called as Interface-PSA. Then four representative numerical examples, including droplet deformation, Rayleigh-Taylor instability, dam breaking, and bubble rising, are presented and compared well with reference data. It is demonstrated that inherent interfacial flow physics can be well captured, including surface tension and the dynamic evolution of the complex interfaces.