Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients

Abstract

The correct evaluation of gradients is at the cornerstone of the smoothed particle hydrodynamics (SPH) technique. Using an integral approach to estimate gradients has proven to enhance accuracy substantially. Such approach retains the Lagrangian structure of SPH equations and is fully conservative. But, in practice, it is difficult to make the Lagrangian formulation totally consistent to an exact partition of the unity. In this paper we study, among other things, the connection between the choice of the volume elements (VEs), which enters in the SPH summations, and the accuracy in the gradient estimation within the integral approach scheme (ISPH). A new variant of VEs are proposed which improve the partition of the unity and is fully compatible with the Lagrangian formulation of SPH, including the grad-h corrections. Using analytic considerations, simple static toy models in 1D, and a few full 3D test cases, we show that any improvement in the partition of the unity also leads to a better calculation of gradients when the integral approach is used jointly. Additionally, we propose an easy-to-implement modification of the ISPH scheme, which makes it more flexible and better suited to handle sharp density contrasts. The ISPH code built with the proposed scheme has been validated with a good number of standard tests, some of them involving contact discontinuities. The performance of the code was excellent in all of them, showing that an improvement in the partition of the unity is not detrimental of the good conservation of energy, momentum, and entropy typical of Lagrangian schemes.

arXiv comment20 pages, 20 figures, published in Astronomy and Astrophysics

References in corpus (23)

• Cosmological SPH simulations: The entropy equation
• 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
• A Density Independent Formulation of Smoothed Particle Hydrodynamics
• Multiphase Smoothed-Particle Hydrodynamics
• A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
• Numerical Convergence in Smoothed Particle Hydrodynamics
• CRKSPH - A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme
• Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests
• A one-parameter family of interpolating kernels for Smoothed Particle Hydrodynamics studies
• Improving smoothed particle hydrodynamics with an integral approach to calculating gradients
• The Lagrangian hydrodynamics code MAGMA2
• SPHYNX: an accurate density-based SPH method for astrophysical applications
• The Kelvin-Helmholtz instability and smoothed particle hydrodynamics
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Testing the concept of integral approach to derivatives within the smoothed particle hydrodynamics technique in astrophysical scenarios
• A simple, entropy-based dissipation trigger for SPH
• Smoothed particle magnetohydrodynamics with the geometric density average force expression

Cited by in corpus (5)

• Volume and Mass Conservation in Lagrangian Meshfree Methods
• Sub-stellar engulfment by a main sequence star: where is the lithium?
• Adaptive Compressible Smoothed Particle Hydrodynamics
• Axisymmetric smoothed particle magneto-hydrodynamics
• Mixing Sinc kernels to improve interpolations in smoothed particle hydrodynamics without pairing instability