Semi-Implicit Lagrangian Voronoi Approximation for Compressible Viscous Fluid Flows

TL;DR

This paper introduces SILVA, a semi-implicit Lagrangian Voronoi scheme for simulating complex compressible and multi-phase fluid flows, combining accuracy with computational stability across different Mach regimes.

Contribution

It presents a novel semi-implicit Lagrangian Voronoi method that improves stability and accuracy in simulating compressible viscous flows without grid deterioration.

Findings

Effective in low-Mach number flows

Handles shocks and multi-phase flows accurately

Reduces time step restrictions compared to explicit methods

Abstract

This paper contributes to the recent investigations of Lagrangian methods based on Voronoi meshes. The aim is to design a new conservative numerical scheme that can simulate complex flows and multi-phase problems with more accuracy than SPH (Smoothed Particle Hydrodynamics) methods but, unlike diffuse interface models on fixed grid topology, does not suffer from the deteriorating quality of the computational grid. The numerical solution is stored at particles, which move with the fluid velocity and also play the role of the generators of the computational mesh, that is efficiently re-constructed at each time step. The main novelty stems from combining a Lagrangian Voronoi scheme with a semi-implicit integrator for compressible flows. This allows to model low-Mach number flows without the extremely stringent stability constraint on the time step and with the correct scaling of numerical viscosity. The implicit linear system for the unknown pressure is obtained by splitting the reversible from the irreversible (viscous) part of the dynamics, and then using entropy conservation of the reversible sub-system to derive an auxiliary elliptic equation. The method, called SILVA (Semi-Implicit Lagrangian Voronoi Approximation), is validated in a variety of test cases that feature diverse Mach numbers, shocks and multi-phase flows.