Abstract:This paper presents a general positivity-preserving algorithm for implicit high-order finite volume schemes that solve compressible Euler and Navier-Stokes equations to ensure the positivity of density and internal energy (or pressure). Previous positivity-preserving algorithms are mainly based on the slope limiting or flux limiting technique, which rely on the existence of low-order positivity-preserving schemes. This dependency poses serious restrictions on extending these algorithms to temporally implicit schemes since it is difficult to know if a low-order implicit scheme is positivity-preserving. In the present paper, a new positivity-preserving algorithm is proposed in terms of the flux correction technique. And the factors of the flux correction are determined by a residual correction procedure. For a finite volume scheme that is capable of achieving a converged solution, we show that the correction factors are in the order of unity with additional high-order terms corresponding to the spatial and temporal rates of convergence. Therefore, the proposed positivity-preserving algorithm is accuracy-reserving and asymptotically consistent. The notable advantage of this method is that it does not rely on the existence of low-order positivity-preserving baseline schemes. Therefore, it can be applied to the implicit schemes solving Euler and especially Navier-Stokes equations. In the present paper, the proposed technique is applied to an implicit dual time-stepping finite volume scheme with temporal second-order and spatial high-order accuracy. The present positivity-preserving algorithm is implemented in an iterative manner to ensure that the dual time-stepping iteration will converge to the positivity-preserving solution. Another similar correction technique is also proposed to ensure that the solution remains positivity-preserving at each sub-iteration. Numerical results demonstrate that the proposed algorithm preserves positive density and pressure in all test cases and significantly improves the robustness of the numerical schemes.

An Efficient Positivity-Preserving Finite Volume Scheme for the Nonequilibrium Three-Temperature Radiation Diffusion Equations on Polygonal Meshes

A Positivity-Preserving and Free Energy Dissipative Hybrid Scheme for the Poisson-Nernst-Planck Equations on Polygonal and Polyhedral Meshes

A positivity-preserving edge-centred finite volume scheme for heterogeneous and anisotropic diffusion problems on polygonal meshes

Simple Positivity Preserving Nonlinear Finite Volume Scheme for Subdiffusion Equations on General Non-conforming Distorted Meshes

A symmetric and coercive finite volume scheme preserving the discrete maximum principle for anisotropic diffusion equations on star-shaped polygonal meshes

On positivity preservation in nonlinear finite volume method for multi-term fractional subdiffusion equation on polygonal meshes

A New Positivity-Preserving Technique for High-Order Schemes to Solve Extreme Problems of Euler Equations on Structured Meshes

A pyramid scheme for three-dimensional diffusion equations on polyhedral meshes

A general positivity-preserving algorithm for implicit high-order finite volume schemes solving the Euler and Navier-Stokes equations

A robust, interpolation‐free and monotone finite volume scheme for diffusion equations on arbitrary quadrilateral meshes

A Monotone Finite Volume Scheme for Diffusion Equations on General Non-Conforming Meshes

Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential

A monotone diffusion scheme for 3D general meshes: Application to radiation hydrodynamics in the equilibrium diffusion limit

A positivity-preserving, linear, energy stable and convergent numerical scheme for the Poisson–Nernst–Planck (PNP) system

Nonorthogonal Hexahedral Mesh Finite Volume Difference Method for the 3-D Diffusion Equation

High-order Conservative Positivity-Preserving DG-interpolation for Deforming Meshes and Application to Moving Mesh DG Simulation of Radiative Transfer.

Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

A novel formulation for heat conduction using non-convex meshes based on smoothed finite element method

High order conservative Lagrangian schemes for two-dimensional radiation hydrodynamics equations in the equilibrium-diffusion limit

A combined finite volume–nonconforming/mixed-hybrid finite element scheme for degenerate parabolic problems