REMIX SPH -- improving mixing in smoothed particle hydrodynamics simulations using a generalised, material-independent approach

• Published in: arXiv.org
• Main Authors: Sandnes, Thomas D, Eke, Vincent R, Kegerreis, Jacob A, Massey, Richard J, Ruiz-Bonilla, Sergio, Schaller, Matthieu, Teodoro, Luis F A
• Format: Paper
• Language: English
• Published: Ithaca Cornell University Library, arXiv.org 26.07.2024
• Subjects: Compressibility; Density; Discontinuity; Dissimilar materials; Equations of state; Fluid mechanics; Ideal gas; Initial conditions; Interface stability; Simulation; Smooth particle hydrodynamics; Source code
• ISSN: 2331-8422

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 hydrodynamic behaviour unrelated to mixing, such as in shocks. We demonstrate marked improvements in three-dimensional test scenarios, focusing on more challenging cases with particles of equal mass across the simulation. This validates our methods for use-cases relevant across applications spanning astrophysics and engineering, where particles are free to evolve over a large range of density scales, or where emergent and evolving density discontinuities cannot easily be corrected by choosing bespoke particle masses in the initial conditions. 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 in mind, which means that its improved hydrodynamic treatment can be used for high-resolution simulations without significant cost to run-speed.