Abstract:We present a new hydrodynamic scheme named Godunov Density-Independent Smoothed Particle Hydrodynamics (GDISPH), that can accurately handle shock waves and contact discontinuities without any manually tuned parameters. This is in contrast to the standard formulation of smoothed particle hydrodynamics (SSPH), which requires the parameters for an artificial viscosity term to handle the shocks and struggles to accurately handle the contact discontinuities due to unphysical repulsive forces, resulting in surface tension that disrupts pressure equilibrium and suppresses fluid instabilities. While Godunov SPH (GSPH) can handle the shocks without the parameters by using solutions from a Riemann solver, it still cannot fully handle the contact discontinuities. Density-Independent Smoothed Particle Hydrodynamics (DISPH), one of several schemes proposed to handle contact discontinuities more effectively than SSPH, demonstrates superior performance in our tests involving strong shocks and contact discontinuities. However, DISPH still requires the artificial viscosity term. We integrate the Riemann solver into DISPH in several ways, yielding some patterns of GDISPH. The results of standard tests such as the one-dimensional Riemann problem, pressure equilibrium, Sedov-Taylor, and Kelvin-Helmholtz tests are favourable to GDISPH Case 1 and GDISPH Case 2, as well as DISPH. We conclude that GDISPH Case 1 has an advantage over GDISPH Case 2, effectively handling shocks and contact discontinuities without the need for specific parameters or introducing any additional numerical diffusion.
Instrumentation and Methods for Astrophysics,Cosmology and Nongalactic Astrophysics,Astrophysics of Galaxies,Computational Physics,Fluid Dynamics
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve two key problems encountered in fluid dynamics simulations: **the accurate treatment of shock waves and contact discontinuities**. Specifically, the paper proposes a new method named Godunov Density - Independent Smoothed Particle Hydrodynamics (GDISPH) to overcome the deficiencies in existing methods.
#### 1. Treatment of shock waves
In the traditional Smoothed Particle Hydrodynamics (SPH) method, in order to correctly handle shock waves, artificial viscosity terms usually need to be introduced, and these parameters need to be manually adjusted. However, this method has the following problems:
- **Parameter sensitivity**: Under different physical conditions, the optimal parameter values will be different, resulting in unstable results.
- **Additional numerical diffusion**: The artificial viscosity not only acts in the shock wave region but also wrongly affects the shear flow region, resulting in unnecessary viscous effects.
Godunov SPH (GSPH) has solved these problems by using the Riemann solver, and can handle shock waves without relying on manual parameter adjustment. But GSPH is still difficult to fully handle contact discontinuities.
#### 2. Treatment of contact discontinuities
Contact discontinuities refer to interfaces in the fluid where the density or pressure suddenly changes but the velocity is continuous. The traditional SPH method has the following problems when handling contact discontinuities:
- **Non - physical repulsive force**: Since SPH assumes that the local density distribution is differentiable, which is not valid at the contact discontinuity, resulting in non - physical surface tension effects, suppressing the development of fluid instabilities.
- **Pressure imbalance**: The pressure calculation at the contact discontinuity is inaccurate, resulting in pressure fluctuations.
To solve these problems, Saitoh and Makino (2013) proposed Density - Independent SPH (DISPH), which assumes that the pressure is continuous and differentiable, thus handling contact discontinuities more accurately. However, DISPH still needs to introduce artificial viscosity to handle shock waves.
#### Proposal of GDISPH
In order to solve the above two problems simultaneously, this paper proposes the GDISPH method, which combines the advantages of GSPH and DISPH:
- **Automatically handle shock waves**: By integrating the Riemann solver, GDISPH can handle shock waves without relying on manual parameter adjustment.
- **Accurately handle contact discontinuities**: It inherits the contact discontinuity handling ability of DISPH and avoids non - physical surface tension effects.
In this way, GDISPH can effectively handle shock waves and contact discontinuities without additional parameters, providing more accurate fluid dynamics simulation results.
### Summary
The main objective of this paper is to develop a new fluid dynamics scheme - GDISPH, to overcome the limitations of existing methods in handling shock waves and contact discontinuities, thereby providing a more accurate and stable fluid dynamics simulation tool.