The role of time integration in energy conservation in Smoothed Particle Hydrodynamics fluid dynamics simulations

TL;DR

This paper investigates how different time integration schemes affect energy conservation in SPH fluid simulations, providing explicit formulas to measure deviations and demonstrating the impact through impact and vortex simulations.

Contribution

It introduces a posteriori expressions to evaluate energy deviations caused by time integration schemes in SPH, including the first stable non-dissipative simulation.

Findings

Implicit schemes improve energy conservation in viscous flows.

First stable non-dissipative SPH simulation achieved.

Explicit formulas for energy deviation assessment provided.

Abstract

The choice of a time integration scheme is a crucial aspect of any transient fluid simulation, and Smoothed-Particle Hydrodynamics (SPH) is no exception. The influence of the time integration scheme on energy balance is here addressed. To do so, explicit expressions allowing to compute the deviations from the energy balance, induced by the time integration scheme, are provided. These expressions, computed \textit{a posteriori}, are valid for different integration methods. Such formulation is tested with the simulation of a two-dimensional non-viscous impact of two jets, with no artificial dissipation terms. To the best of our knowledge, this is the first stable simulation of a non-dissipative system with a weakly-compressible SPH method. A viscous case, the Taylor-Green vortex, has also been simulated. Results show that an implicit time integration scheme also behaves better in a viscous context.