REMIX SPH -- improving mixing in smoothed particle hydrodynamics simulations using a generalised, material-independent approach
arXiv:2407.18587 · doi:10.1016/j.jcp.2025.113907
Abstract
We present REMIX, a smoothed particle hydrodynamics (SPH) scheme designed to alleviate effects that typically suppress mixing and instability growth at density discontinuities in SPH simulations. We approach this problem by directly targeting sources of kernel smoothing error and discretisation error, resulting in a generalised, material-independent formulation that improves the treatment both of discontinuities within a single material, for example in an ideal gas, and of interfaces between dissimilar materials. This approach also leads to improvements in capturing wider hydrodynamic behaviour unrelated to mixing. We demonstrate marked improvements in three-dimensional test scenarios, focusing on cases with particles of equal mass across the simulation. This choice is particularly relevant for use cases in astrophysics and engineering -- specifically those in which particles are free to evolve over a large range of density scales -- where bespoke choices of unequal particle masses in the initial conditions cannot easily be used to address emergent and evolving density discontinuities. We achieve these improvements while maintaining sharp discontinuities; without introducing additional equation of state dependence in, for example, particle volume elements; and without contrived or targeted corrections. Our methods build upon a fully compressible and thermodynamically consistent core-SPH construction, retaining Galilean invariance as well as conservation of mass, momentum, and energy. REMIX is integrated in the open-source, state-of-the-art SWIFT code and is designed with computational efficiency also in mind, meaning that its improved hydrodynamic treatment can be used for high-resolution simulations without prohibitive cost to run-speed.
55 pages, 36 figures, published in Journal of Computational Physics
References in corpus (43)
- The cosmological simulation code GADGET-2
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- Common Envelope Evolution: Where we stand and how we can move forward
- Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
- Smoothed Particle Hydrodynamics and Magnetohydrodynamics
- Fundamental differences between SPH and grid methods
- Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
- A General Class of Lagrangian Smoothed Particle Hydrodynamics Methods and Implications for Fluid Mixing Problems
- Smoothed Particle Hydrodynamics in Astrophysics
- Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
- Inviscid SPH
- Gasoline2: A Modern SPH Code
- Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver
- The EAGLE simulations of galaxy formation: the importance of the hydrodynamics scheme
- A Density Independent Formulation of Smoothed Particle Hydrodynamics
- Multiphase Smoothed-Particle Hydrodynamics
- Subsonic turbulence in smoothed particle hydrodynamics and moving-mesh simulations
- Resolving mixing in Smoothed Particle Hydrodynamics
- AQUA: A Collection of HO Equations of State for Planetary Models
- SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch
- A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
- SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
- Astrophysical Weighted Particle Magnetohydrodynamics
- Computational Eulerian Hydrodynamics and Galilean Invariance
- CRKSPH - A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme
- Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
- A new equation of state for dense hydrogen-helium mixtures II: taking into account hydrogen-helium interactions
- Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations
- Improving smoothed particle hydrodynamics with an integral approach to calculating gradients
- Bifurcation in the history of Uranus and Neptune: the role of giant impacts
- Numerical aspects of Giant Impact simulations
- SWIFT: Using task-based parallelism, fully asynchronous communication, and graph partition-based domain decomposition for strong scaling on more than 100,000 cores
- Planetary Giant Impacts: Convergence of High-Resolution Simulations using Efficient Spherical Initial Conditions and SWIFT
- Atmospheric Erosion by Giant Impacts onto Terrestrial Planets: A Scaling Law for any Speed, Angle, Mass, and Density
- The Lagrangian hydrodynamics code MAGMA2
- Immediate origin of the Moon as a post-impact satellite
- A new insight into the consistency of smoothed particle hydrodynamics
- Enhanced mixing in giant impact simulations with a new Lagrangian method
- An efficient and generalized consistency correction method for weakly-compressible SPH
- The effect of pre-impact spin on the Moon-forming collision
- Density Independent Smoothed Particle Hydrodynamics for Non-Ideal Equation of State
- Sensitivity of non-radiative cloud-wind interactions to the hydrodynamics solver
Cited by in corpus (5)
- On the origin of Jupiter's fuzzy core: constraints from N-body, impact and evolution simulations
- No dilute core produced in simulations of giant impacts on to Jupiter
- On the consistency of jet feedback modelling across different astrophysics hydrodynamical codes
- Optimized smoothing kernels for SPH
- Shock-capturing particle hydrodynamics with reproducing kernels