A numerical method for low Mach number compressible flows by simultaneous relaxation of dependent variables

TL;DR

This paper introduces a new numerical method for low Mach number compressible flows that simultaneously relaxes multiple variables and conserves energy, improving accuracy in various flow scenarios.

Contribution

The study proposes a conservative finite difference scheme that enhances energy conservation and accurately models low Mach number flows without approximations.

Findings

Excellent energy conservation in sound wave analysis

Discrete conservation of momentum, energy, and entropy in 3D flows

Accurate turbulence and convection modeling in complex flows

Abstract

Density varies spatiotemporally in low Mach number flows. Hence, incompressibility cannot be assumed, and the density must be accurately solved. Various methods have been proposed to analyze low Mach number flows, but their energy conservation properties have not been investigated in detail. This study proposes a new method for simultaneously relaxing velocity, pressure, density, and internal energy using a conservative finite difference scheme with excellent energy conservation properties to analyze low Mach number flows. In the analysis for sound wave propagation in an inviscid compressible flow, the amplitude amplification ratio and frequency of sound wave obtained by this numerical method agree well with the theoretical values. In the analysis for a three-dimensional periodic inviscid compressible flow, each total amount for the momentum, total energy, and entropy are discretely conserved. When no approximation, such as low Mach number approximations, is applied to the fundamental equations, the excellent conservation properties of momentum, total energy, and entropy are achieved. In decaying compressible isotropic turbulence, this computational method can capture turbulence fluctuations. Analyzing the Taylor-Green decaying vortex, we confirmed the validity of this computational scheme for compressible viscous flows. In the calculation for the natural convection in a cavity, the validity of this numerical method was presented even in incompressible flows considering density variation. In a three-dimensional Taylor decaying vortex problem, it was shown that this numerical method can accurately calculate incompressible flows. We clarified the accuracy and validity of the present numerical method by analyzing various flow models and demonstrated the possibility of applying this method to complex flow fields.