粒子法中压力振荡的机理研究

移动粒子半隐式法(moving particle semi-implicit method, MPS)是一种适用于不可压缩流体的无网格方法, MPS方法常应用于自由表面大变形问题.MPS 方法提出至今一直存在着严重的压力振荡问题. 本研究针对MPS 方法中存在的压力振荡现象, 首先将实际的物理问题简化为一维模型, 并从粒子之间相互位置关系的角度说明了MPS 方法中压力波动产生的原因.在采用MPS方法进行模拟时, 加入了粒子碰撞模型, 通过对碰撞系数的选择从而控制粒子之间的相互位置关系.并且对经典的溃坝问题进行了模拟, 结果表明随着碰撞系数的增加, 粒子数密度偏差的波动幅度都会减小, 从而压力振荡的幅度得到了有效的抑制.并且对比了两种不同核函数对压力振荡的影响, 结果表明: 采用高斯核函数时, 压力振荡的幅度更小, 这是因为采用高斯核函数时, 相同的粒子位置波动幅度将会得到较小的粒子数密度偏差的波动.由于在模拟过程中粒子运动的随机性, 这将导致粒子数密度偏差产生随机的波动, 从而产生压力振荡, 因此粒子法中的压力振荡很难彻底消除.      

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.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

Improvement of stability in moving particle semi‐implicit method

Moving particle semi‐implicit (MPS) method is one of the particle methods, which can be used to analyze incompressible free surface flow without surface tracking by a mesh or a scalar quantity. However, MPS causes unphysical numerical oscillation of pressure with high frequencies. We proposed a new formulation for the source term of Poisson equation of pressure. The proposed source term consists of three parts, one main part and two error‐compensating parts. With proper selection of the coefficients for the error‐compensating parts, we can suppress the unphysical pressure oscillation. Smoother pressure distributions are obtained in hydrostatic pressure and dam break problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Improved MPS method with variable‐size particles

Using variable‐size particles in the moving particle semi‐implicit method (MPS) could lead to inaccurate predictions and/or numerical instability. In this paper, a variable‐size particle moving particle semi‐implicit method (VSP‐MPS) scheme is proposed for the MPS method to achieve more reliable simulations with variable‐size particles. To improve stability and accuracy, a new gradient model is developed based on a previously developed MPS scheme that requires no surface detection MPS. The dynamic particle coalescing and splitting algorithm is revised to achieve dynamic multi‐resolution. A cubic spline function with additional function is employed as the kernel function. The effectiveness of the VSP‐MPS method is demonstrated by three verification examples, that is, a hydrostatic pressure problem, a complicated free surface flow problem with large deformation, and a dynamic impact problem. The new VSP‐MPS scheme with variable‐size particles is found to have balanced efficiency and accuracy that is suitable for simulating large systems with complex flow patterns. Copyright © 2015 John Wiley & Sons, Ltd.     

A moving particle semi‐implicit method for free surface flow: Improvement in inter‐particle force stabilization and consistency restoring

The moving particle semi‐implicit (MPS) method has been widely applied in free surface flows. However, the implementation of MPS remains limited because of compressive instability occurred when the particles are under compressive stress states. This study proposed an inter‐particle force stabilization and consistency restoring MPS (IFS‐CR‐MPS) method to overcome this numerical instability. For inter‐particle force stabilization, a hyperbolic‐shaped quintic kernel function is developed with a non‐negative and smooth second order derivative to satisfy the stability criterion under compressive stress state. Then, a contrastive study is conducted on the contradiction between the common understanding of the conventional MPS hyperbolic‐shaped kernel function and its performance. The result shows that the conventional MPS hyperbolic‐shaped kernel function can easily cause violent repulsive inter‐particle force and then lead to the compressive instability. Therefore, the first order derivative of the modified hyperbolic‐shaped quintic kernel function is recommended as the form of the contribution of the neighbor particles to achieve a more stable inter‐particle repulsive force. For consistency restoring, the Taylor series expansion and the hyperbolic‐shaped quintic kernel are combined to improve the accuracy of the viscosity and pressure calculation. The IFS‐CR‐MPS algorithm is subsequently verified by the inviscid hydrostatic pressure, jet impacting, and viscous droplet impacting problems. These results can be used for choosing kernel function and the contribution of neighbor particles in particle methods. Copyright © 2016 John Wiley & Sons, Ltd.     

粒子法中压力振荡问题的研究

传统移动粒子半隐式法MPS(Moving Particle Semi-implicit Method)中一直存在压力振荡问题,针对此问题对MPS方法进行改进。改进的MPS方法,采用一种新型抑制压力振荡的压力泊松方程离散格式;在核函数的选择方面,采用能够增加计算稳定性的二次样条核函数;并且针对MPS方法中粒子插值不完整问题,对粒子插值不完整性进行了修正。应用改进的MPS方法对溃坝问题进行数值模拟验证。结果表明,应用改进的MPS方法能够得到更为光滑的压力场空间分布。对模拟过程中的检测点压力进行采集,并且与实验值进行对比分析,发现改进的MPS方法能够有效地抑制模拟过程中的压力振荡,而且与实验值接近。同时应用改进的MPS方法对静水问题进行验证模拟,发现改进的MPS方法能够有效地抑制模拟过程中的压力振荡,而且监测点的压力与理论解接近。改进的MPS方法对今后应用MPS方法模拟实际工程问题,并且获得准确稳定的压力值有着重要的意义。     

A pressure correction method for fluid–particle interaction flow: Direct‐forcing method and sedimentation flow

A direct‐forcing pressure correction method is developed to simulate fluid–particle interaction problems. In this paper, the sedimentation flow is investigated. This method uses a pressure correction method to solve incompressible flow fields. A direct‐forcing method is introduced to capture the particle motions. It is found that the direct‐forcing method can also be served as a wall‐boundary condition. By applying Gauss's divergence theorem, the formulas for computing the hydrodynamic force and torque acting on the particle from flows are derived from the volume integral of the particle instead of the particle surface. The order of accuracy of the present method is demonstrated by the errors of velocity, pressure, and wall stress. To demonstrate the efficiency and capability of the present method, sedimentations of many spherical particles in an enclosure are simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of plunging wave breaking by the weakly compressible smoothed particle hydrodynamic method

In this paper, a numerical meshless method called the weakly compressible smoothed particle hydrodynamic (WCSPH) method, which is a two-dimensionalmodel of a weakly compressible fluid, is applied to simulate the plunging wave breaking process. This model solves the viscous fluid equations to obtain the velocity and density fields and also solves the equation of state to obtain the pressure field. The WCSPH method is demonstrated to have a higher computational efficiency than the basic SPH model. To simulate the turbulent behavior of the fluid flow in the wave breaking procedure, a sub particle scale (SPS) model is used, which is obtained from the Large eddy simulation (LES) theory. To consider the accuracy of the standard WCSPH model (WCSPH model without considering the turbulent effect), a dam break test is performed, and model results are compared with available experimental data.     

Some improvements on free surface simulation by the particle finite element method

The Lagrangian method has become increasingly popular in numerical simulation of free surface problems. In this paper, after a brief review of a recent Lagrangian method, namely the particle finite element method, some issues are discussed and some improvements are made. The least‐square finite element method is adopted to simplify the solving of the Navier–Stokes equations. An adaptive time method is derived to obtain suitable time steps. A mass correction procedure is imported to improve the mass conservation in long time calculations and time discretization scheme is adopted to decrease the pressure oscillations during the calculations. Finally, the method is used to simulate a series of examples and the results are compared with the commercial FLOW3D code. Copyright © 2008 John Wiley & Sons, Ltd.     

人工黏性增强三维移动粒子半隐算法稳定性

通过综合比较移动粒子半隐方法(moving-particle semi-implicitmethod, MPS)各种稳定性计算方法, 考虑了多种可能导致计算不稳定的因素, 并首次提出加入人工黏性来抑制非物理压力振荡, 得到了较为稳定的三维MPS算法.采用各种稳定性方案对比模拟了三维立方流体旋转状态下的变形及三维静水压问题, 并进一步与商业软件的流体体积函数方法(volumeof fluid, VOF)模型计算结果对比, 验证了其正确性.结果表明:对三维情形, 单纯应用已有的稳定性算法仍难以满足MPS计算稳定性的要求, 而进一步辅以该文提出的人工黏性方法则可以在准确性及稳定性方面均取得较好的效果, 且计算结果显示, 通过该文提出的开关模式增加人工黏性并不影响压强泊松方程的求解精度, 因而不会对流场的求解产生负面作用.     

A particle–gridless hybrid method for incompressible flows

A particle–gridless hybrid method for the analysis of incompressible flows is presented. The numerical scheme consists of Lagrangian and Eulerian phases as in an arbitrary Lagrangian–Eulerian (ALE) method, where a new‐time physical property at an arbitrary position is determined by introducing an artificial velocity. For the Lagrangian calculation, the moving‐particle semi‐implicit (MPS) method is used. Diffusion and pressure gradient terms of the Navier–Stokes equation are calculated using the particle interaction models of the MPS method. As an incompressible condition, divergence of velocity is used while the particle number density is kept constant in the MPS method. For the Eulerian calculation, an accurate and stable convection scheme is developed. This convection scheme is based on a flow directional local grid so that it can be applied to multi‐dimensional convection problems easily. A two‐dimensional pure convection problem is calculated and a more accurate and stable solution is obtained compared with other schemes. The particle–gridless hybrid method is applied to the analysis of sloshing problems. The amplitude and period of sloshing are predicted accurately by the present method. The range of the occurrence of self‐induced sloshing predicted by the present method shows good agreement with the experimental data. Calculations have succeeded even for the higher injection velocity range, where the grid method fails to simulate. Copyright © 1999 John Wiley & Sons, Ltd.     

An improved SPH model for multiphase flows with large density ratios

This paper presents a new smoothed particle hydrodynamics (SPH) model for simulating multiphase fluid flows with large density ratios. The new SPH model consists of an improved discretization scheme, an enhanced multiphase interface treatment algorithm, and a coupled dynamic boundary treatment technique. The presented SPH discretization scheme is developed from Taylor series analysis with kernel normalization and kernel gradient correction and is then used to discretize the Navier‐Stokes equation to obtain improved SPH equations of motion for multiphase fluid flows. The multiphase interface treatment algorithm involves treating neighboring particles from different phases as virtual particles with specially updated density to maintain pressure consistency and a repulsive interface force between neighboring interface particles into the pressure gradient to keep sharp interface. The coupled dynamic boundary treatment technique includes a soft repulsive force between approaching fluid and solid particles while the information of virtual particles are approximated using the improved SPH discretization scheme. The presented SPH model is applied to 3 typical multiphase flow problems including dam breaking, Rayleigh‐Taylor instability, and air bubble rising in water. It is demonstrated that inherent multiphase flow physics can be well captured while the dynamic evolution of the complex multiphase interfaces is sharp with consistent pressure across the interfaces.     

MPS与GPU结合数值模拟LNG液舱晃荡

液舱晃荡是一种在外部激励作用下部分装载的液舱内液体的波动现象,它会对液舱结构强度和运输船舶稳性产生危害.移动粒子半隐式法(moving particle semi-implicit,MPS)是一种典型的无网格粒子类方法,可以有效地模拟剧烈的液舱晃荡问题.但MPS方法存在计算效率低的缺点,难以模拟大规模三维问题,而GPU并行加速技术已广泛应用于科学计算领域.因此,本文将MPS方法与GPU并行加速技术相结合,采用CUDA程序语言编写,自主开发了MPSGPU-SJTU求解器,对三维液化天然气(liquefiednatural gas, LNG)型液舱晃荡进行了数值模拟.通过三种不同粒子间距的数值模拟,验证了求解器的收敛性,其中最大计算粒子数达到了200多万.与其他研究结果相比,MPSGPU-SJTU求解器能够准确地预测壁面砰击压力,并且捕捉晃荡过程中自由面的大幅度变形和强非线性破碎现象.相比CPU求解器的计算时间,GPU并行加速技术可以大幅度地减小计算时长,提高MPS方法的计算效率.本文将LNG型液舱与方型液舱的晃荡进行对比,结果表明在高充液率下LNG型液舱可以有效地减小晃荡幅值和壁面砰击压力.但在中低充液率下,LNG型液舱则会加剧晃荡,自由面呈现明显的三维特征.本文还进一步研究了水和LNG两种不同介质的液舱晃荡现象,数值模拟结果表明二者的流场基本相似,砰击压力则正比于液体密度.      

Turbulence model and numerical scheme assessment for buffet computations

The prediction of shock‐induced oscillations over transonic rigid airfoils is important for a better understanding of the buffeting phenomenon. The unsteady resolution of the Navier–Stokes equations is performed with various transport‐equation turbulence models in which corrections are added for non‐equilibrium flows. The lack of numerical efficiency due to the CFL stability condition is circumvented by the use of a wall law approach and a dual time stepping method. Moreover, various numerical schemes are used to try and be independent of the numerical discretization. Comparisons are made with the experimental results obtained for the supercritical RA16SC1 airfoil. They show the interest in using the SST correction or realizability conditions to get correct predictions of the frequency, amplitude and pressure fluctuations over the airfoil. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

A simple incompressible flux splitting for sharp free surface capturing

This paper first applies a flux vector‐type splitting method based on the numerical speed of sound for computing incompressible single and multifluid flows. Here, a preconditioning matrix based on Chorin's artificial compressibility concept is used to modify the incompressible multifluid Navier–Stokes equations to be hyperbolic and density or volume fraction‐independent. The current approach can reduce eigenvalues disparity induced from density or volume fraction ratios and enhance numerical stability. Also, a simple convection‐pressure flux‐splitting method with high‐order essentially nonoscillatory‐type primitive variable extrapolations coupled with monotone upstream‐centered schemes for conservation laws‐type volume fraction recompressed reconstruction is used to maintain the preservation of sharp interface evolutions in multifluid flow simulations. Benchmark tests including a solid rotation test of a notched two‐dimensional cylinder, the evolution of spiral and rotational shapes of deformable circles, a dam breaking problem, and the Rayleigh–Taylor instability were chosen to validate the current incompressible multifluid methodology. An incompressible driven cavity was also chosen to check the robustness of the proposed method on the computation of single fluid incompressible flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

A pressure correction‐volume of fluid method for simulation of two‐phase flows

A pressure correction method coupled with the volume of fluid (VOF) method is developed to simulate two‐phase flows. A volume fraction function is introduced in the VOF method and is governed by an advection equation. A modified monotone upwind scheme for a conservation law (modified MUSCL) is used to solve the solution of the advection equation. To keep the initial sharpness of an interface, a slope modification scheme is introduced. The continuum surface tension (CST) model is used to calculate the surface tension force. Three schemes, central‐upwind, Parker–Youngs, and mixed schemes, are introduced to compute the interface normal vector and the gradient of the volume fraction function. Moreover, a height function technique is applied to compute the local curvature of the interface. Several basic test problems are performed to check the order of accuracy of the present numerical schemes for computing the interface normal vector and the gradient of the volume fraction function. Three physical problems, two‐dimensional broken dam problem, static drop, and spurious currents, and three‐dimensional rising bubble, are performed to demonstrate the efficiency and accuracy of the pressure correction method. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparative study of truly incompressible and weakly compressible SPH methods for free surface incompressible flows

In this paper, the performance of the incompressible SPH (ISPH) method and an improved weakly compressible SPH (IWCSPH) method for free surface incompressible flows are compared and analyzed. In both methods, the Navier–Stokes equations are solved, and no artificial viscosity is used. The ISPH algorithm in this paper is based on the classical SPH projection method with common treatments on solid boundaries and free surfaces. The IWCSPH model includes some advanced corrective algorithms in density approximation and solid boundary treatment (SBT). In density approximation, the moving least squares (MLS) approach is applied to re‐initialize density every several steps to obtain smoother and more stable pressure fields. An improved coupled dynamic SBT algorithm is implemented to obtain stable pressure values near solid wall areas and, thus, to minimize possible numerical oscillations brought in by the solid boundaries. Three representative numerical examples, including a benchmark test for hydrostatic pressure, a dam breaking problem and a liquid sloshing problem, are comparatively analyzed with ISPH and IWCSPH. It is demonstrated that the present IWCSPH is more attractive than ISPH in modeling free surface incompressible flows as it is more accurate and more stable with comparable or even less computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

An extension of the SIMPLE based discontinuous Galerkin solver to unsteady incompressible flows

In this paper, we present a SIMPLE based algorithm in the context of the discontinuous Galerkin method for unsteady incompressible flows. Time discretization is done fully implicit using backward differentiation formulae (BDF) of varying order from 1 to 4. We show that the original equation for the pressure correction can be modified by using an equivalent operator stemming from the symmetric interior penalty (SIP) method leading to a reduced stencil size. To assess the accuracy as well as the stability and the performance of the scheme, three different test cases are carried out: the Taylor vortex flow, the Orr‐Sommerfeld stability problem for plane Poiseuille flow and the flow past a square cylinder. (1) Simulating the Taylor vortex flow, we verify the temporal accuracy for the different BDF schemes. Using the mixed‐order formulation, a spatial convergence study yields convergence rates of k + 1 and k in the L²‐norm for velocity and pressure, respectively. For the equal‐order formulation, we obtain approximately the same convergence rates, while the absolute error is smaller. (2) The stability of our method is examined by simulating the Orr–Sommerfeld stability problem. Using the mixed‐order formulation and adjusting the penalty parameter of the symmetric interior penalty method for the discretization of the viscous part, we can demonstrate the long‐term stability of the algorithm. Using pressure stabilization the equal‐order formulation is stable without changing the penalty parameter. (3) Finally, the results for the flow past a square cylinder show excellent agreement with numerical reference solutions as well as experiments. Copyright © 2015 John Wiley & Sons, Ltd.