Difference between smoothed particle hydrodynamics and moving particle semi-implicit operators

Smoothed particle hydrodynamics (SPH) and moving particle semi-implicit (MPS) methods are representative meshfree particle methods used to compute Lagrangian mechanics. The approximations of differential operators in the SPH and MPS methods have several similarities, but the theoretical discussion of the difference between them is limited. This study mathematically describes the difference via a comprehensive derivation of the first- and second-order derivative operators for each method. The comprehensive derivation indicates that the SPH and MPS operators are consistent with the pressure Poisson equation and moving least-squares approximation, respectively. The variation in consistency can explain the difference in the schemes of the incompressible flow problem. Additionally, the comprehensive derivation of the MPS operators can result in novel second-order and anisotropic operators. This study strengthens the theoretical understanding of the SPH and MPS methods and facilitates the selection of the appropriate method by users. Furthermore, the proposed MPS operators contribute to developing methods with adaptive or multiscale particle distributions

A comparison between weakly-compressible smoothed particle hydrodynamics (WCSPH) and moving particle semi-implicit (MPS) methods for 3d dam-break flows

Lagrangian particle-based methods have opened new perspectives for the investigation of complex problems with large free-surface deformation. Some well-known particle-based methods adopted to solve non-linear hydrodynamics problems are the smoothed particle hydrodynamics (SPH) and the moving particle semi-implicit (MPS). Both methods modeled the continuum by a system of Lagrangian particles (points) but adopting distinct approaches for the numerical operators, pressure calculation, and boundary conditions. Despite the ability of the particle-based methods in modeling highly nonlinear hydrodynamics, some shortcomings, such as unstable pressure computation and high computational cost remains. In order to assess the performance of these two methods, the weaklycompressible SPH (WCSPH) parallel solver, DualSPHysics, and an in-house incompressible MPS solver are adopted in this work. Two test cases consisting of threedimensional (3D) dam-break problems are simulated, and wave heights, pressures and forces are compared with available experimental data. The influence of the artificial viscosity on the accuracy of WCSPH is investigated. Computational times of both solvers are also compared. Finally, the relative benefits of the methods for solving free-surface problems are discussed, therefore providing directions of their applicability.Comment: 26 pages, 12 figure

Development of an incompressible smoothed particle hydrodynamics method for electrohydrodynamics of immiscible fluids and rigid particles

An incompressible smoothed particle hydrodynamics method for modeling immiscible and isothermal flow of two- and three-phase Newtonian fluids and solid particles subject to an external electric field has been developed. Continuum surface force method is used to calculate the surface tension forces on fluid-fluid interfaces. The materials are assumed to be either perfect or leaky dielectrics. Solid particles are modeled using viscous penalty method coupled with rigidity constraints. The equations are discretized using corrected derivatives and artificial particle displacement is used to ensure homogeneous particle distribution. The projection method is used to advance the governing equations of the flow and electric field in time. The components of the scheme are tested in three stages of two- and three-phase hydrodynamics, multiphase electrohydrodynamics and fluid-structure/solid interaction. The results of each stage is compared to experimental and numerical data available in literature and their validity is established. The combination of the individual elements of the numerical method is used to simulate the motion of rigid particles submerged in Newtonian fluids subject to an external electric field. The behavior of the particles are found to be in agreement with experimental and numerical observations found in the literature. This shows the applicability of the proposed incompressible smoothed particle hydrodynamics scheme in simulating such complex and relatively unexplored phenomena

Numerical simulation of multiphase jet fragmentation using Smoothed Particle Hydrodynamics

This thesis is devoted to the study of multiphase jet fragmentation using Smoothed Particle Hydrodynamics (SPH). The theoretical aspects of three hydrodynamic instabilities, namely the Kelvin-Helmholtz instability (KHI), Rayleigh-Taylor instability (RTI), and Rayleigh Plateau instability (RPI) are reviewed. The linear growth rate of the combined KHI and RTI are derived by means of linear perturbation in chapter 2. The linear growth rate of the multiphase RPI is presented in chapter 7. An overview of the Smoothed Particle Hydrodynamics is given in chapter 3. A pseudo-consistent SPH scheme is presented for the simulation of multiphase flow problems. Additionally, two interface stabilisation models are presented: quasi-buoyancy model and gas-repulsion model. When used in combination with the pseudo-consistent SPH scheme, these models are found to be superior than those presented in the weakly-compressible SPH literature and allows for the simulations for density ratio up to three-magnitudes. The development of an idealised KHI and a KHI subjected to constant gravitational acceleration (stratified shear instability) is examined in chapter 5. The extracted linear growth rate are compared with the theoretical growth rate presented both in the literature and in chapter 2 for the purpose of validation. The development of a single- and multi-mode RTI are studied by means of SPH in chapter 6. Chapter 7 presents the results for the three-dimensional RPI occurring between two fluids. Based on the knowledge acquired in chapter 5-7, the multiphase jet fragmentation driven by the previously mentioned hydrodynamic instabilities are presented in chapter 8. Finally, the major research findings and recommendations are summarised in chapter 9

Numerical simulation of multiphase jet fragmentation using Smoothed Particle Hydrodynamics

This thesis is devoted to the study of multiphase jet fragmentation using Smoothed Particle Hydrodynamics (SPH). The theoretical aspects of three hydrodynamic instabilities, namely the Kelvin-Helmholtz instability (KHI), Rayleigh-Taylor instability (RTI), and Rayleigh Plateau instability (RPI) are reviewed. The linear growth rate of the combined KHI and RTI are derived by means of linear perturbation in chapter 2. The linear growth rate of the multiphase RPI is presented in chapter 7. An overview of the Smoothed Particle Hydrodynamics is given in chapter 3. A pseudo-consistent SPH scheme is presented for the simulation of multiphase flow problems. Additionally, two interface stabilisation models are presented: quasi-buoyancy model and gas-repulsion model. When used in combination with the pseudo-consistent SPH scheme, these models are found to be superior than those presented in the weakly-compressible SPH literature and allows for the simulations for density ratio up to three-magnitudes. The development of an idealised KHI and a KHI subjected to constant gravitational acceleration (stratified shear instability) is examined in chapter 5. The extracted linear growth rate are compared with the theoretical growth rate presented both in the literature and in chapter 2 for the purpose of validation. The development of a single- and multi-mode RTI are studied by means of SPH in chapter 6. Chapter 7 presents the results for the three-dimensional RPI occurring between two fluids. Based on the knowledge acquired in chapter 5-7, the multiphase jet fragmentation driven by the previously mentioned hydrodynamic instabilities are presented in chapter 8. Finally, the major research findings and recommendations are summarised in chapter 9