A conservative Godunov-type mesh-free hydrodynamics scheme for weakly compressible multiphase flows

Abstract

We present a Godunov-type mesh-free hydrodynamics scheme for the numerical simulation of weakly compressible single- and multiphase flows. The method employs a consistent second-order spatial reconstruction of the primitive variables at the inter-particle interface, applies slope limiting, and solves the moving Riemann problem using a subsonic HLLC approximate Riemann solver. Particles are advected with a quasi-Lagrangian material velocity incorporating an inherent self-relaxation mechanism that suppresses the formation of anisotropic particle distributions. Near material interfaces, one-sided renormalized gradients are used to handle density discontinuities and slope limiting is applied for each phase independently. The Lagrangian character of material interfaces is preserved by constraining the relaxation mechanism to act tangentially to the interface. For same-phase particle interactions, the Riemann problem is centered at the arithmetic mean of the two material velocities. For cross-phase interactions, it is centered at the contact-wave velocity. This choice naturally eliminates spurious mass fluxes across phase boundaries. The methodology is validated on a set of challenging benchmarks, including the inviscid Taylor-Green vortex, material interface advection, surface gravity waves near hydrostatic equilibrium, the Kelvin-Helmholtz instability, and the Rayleigh-Taylor instability under gravitational stratification.