An Anti-Clustering Model for Stability Enhancement of a 3D Moving Particle Semi-Implicit Method and Two-Phase Coupling between MPS and Euler Grids

As a Lagrangian gridless particle method, the MPS (Moving Particle Semi-implicit) method has a wide engineering application. However, for complex 3D flows, unphysical pressure oscillations often occur and result in the failure of simulations. This paper compares the stability enhancement methods proposed by different researchers to develop a 3D, stable MPS method. The results indicate that the proposed methods are incapable of eliminating the particle clustering that leads to instability as the main source in coarser particle spacing cases. An anti-clustering model, referring to the SPH (Smoothed Particle Hydrodynamics) artificial viscosity model, is proposed to further reduce instability. Combining various proposed methods and models, several typical examples are simulated comparatively. The results are compared with those of the VOF (Volume of Fluid) model using commercial software to validate the accuracy and stability of the combination of the proposed methods and models. It is concluded that (1) 3D cases that adopt a high-order Laplacian model and high-order source terms in PPE are more accurate than those adopting the low-order operators; (2) the proposed anti-clustering model can produce a tuned interparticle force to prevent particle clustering and introduce no additional viscosity effects in the flow of the normal state, which plays a very positive role for further stability enhancement of MPS; (3) particle resolution significantly maintains simulation accuracy given the stable algorithms by the combination of stability enhancement methods. The 3D MPS method is coupled with the Euler grid (FLUENT V17 software, ANSYS, Pittsburgh, PA, USA) in two phases. In particular, the 3D MPS algorithm is used to calculate the liquid-phase change from the continuous to the dispersed, and the finite volume method based on the Euler grid is adopted to measure the corresponding gas-phase motion. The atomization of the liquid jet under static air flow is calculated and compared with the results of the VOF method, which can capture the continuous interface.