Abstract:Various formulations of smooth-particle hydrodynamics (SPH) have been proposed, intended to resolve certain difficulties in the treatment of fluid mixing instabilities. Most have involved changes to the algorithm which either introduce artificial correction terms or violate arguably the greatest advantage of SPH over other methods: manifest conservation of energy, entropy, momentum, and angular momentum. Here, we show how a class of alternative SPH equations of motion (EOM) can be derived self-consistently from a discrete particle Lagrangian (guaranteeing manifest conservation) in a manner which tremendously improves treatment of instabilities and contact discontinuities. Saitoh & Makino recently noted that the volume element used to discretize the EOM does not need to explicitly invoke the mass density (as in the 'standard' approach); we show how this insight, and the resulting degree of freedom, can be incorporated into the rigorous Lagrangian formulation that retains ideal conservation properties and includes the 'Grad-h' terms that account for variable smoothing lengths. We derive a general EOM for any choice of volume element (particle 'weights') and method of determining smoothing lengths. We then specify this to a 'pressure-entropy formulation' which resolves problems in the traditional treatment of fluid interfaces. Implementing this in a new version of the GADGET code, we show it leads to good performance in mixing experiments (e.g. Kelvin-Helmholtz & blob tests). And conservation is maintained even in strong shock/blastwave tests, where formulations without manifest conservation produce large errors. This also improves the treatment of sub-sonic turbulence, and lessens the need for large kernel particle numbers. The code changes are trivial and entail no additional numerical expense. This provides a general framework for self-consistent derivation of different 'flavors' of SPH.

A general computational framework for Lagrangian hydrodynamic scheme. I: Unification of staggered‐grid and cell‐centered methods

Unification of Lagrangian staggered-grid hydrodynamics and cell-centered hydrodynamics in one dimension

A New Coupled Computational Method in Conjunction with Three-Dimensional Finite Volume Schemes for Nonlinear Coupled Constitutive Relations

A Computational Strategy for Nonlinear Coupled Constitutive Relations of Rarefied Nonequilibrium Flows on Unstructured Grids

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

Second-order decoupled energy-stable schemes for Cahn-Hilliard-Navier-Stokes equations

A new thermodynamically compatible finite volume scheme for Lagrangian gas dynamics

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system

A Linear Unconditionally Stable Scheme for the Incompressible Cahn–Hilliard–Navier–Stokes Phase-Field Model

Linear and Energy-Stable Method with Enhanced Consistency for the Incompressible Cahn–Hilliard–Navier–Stokes Two-Phase Flow Model

A high‐order pure streamfunction method in general curvilinear coordinates for unsteady incompressible viscous flow with complex geometry

A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems

The stabilization-free HDG method for fluid-structure interaction in a unified mixed formulation on Alfeld splits

A fully-decoupled second-order-in-time and unconditionally energy stable scheme for a phase-field model of two phase flow with variable density

A Finite Volume Scheme for Savage-Hutter Equations on Unstructured Grids.

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

A high order conservative semi-Lagrangian discontinuous Galerkin method for two-dimensional transport simulations

A learned conservative semi-Lagrangian finite volume scheme for transport simulations

Preliminary Extension of the Unified Coordinate System Approach to Computation of Viscous Flows

A Finite Volume Method for 2D In viscous Lagrangian Hydrodynamics Based on Characteristics Theory

A Unified Computational Fluid Dynamics Framework from Rarefied to Continuum Regimes