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.

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

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

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

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

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

The Finite Volume Scheme Preserving Maximum Principle for Two-Dimensional Time-Fractional Fokker–Planck Equations on Distorted Meshes

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

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

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

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

A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure

A Positivity-Preserving and Robust Fast Solver for Time-Fractional Convection–Diffusion Problems

A positive meshless finite difference scheme for scalar conservation laws with adaptive artificial viscosity driven by fault detection

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

Positivity-preserving high-order schemes for conservation laws on arbitrarily distributed point clouds with a simple weno limiter

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

A high-order compact finite difference method on nonuniform time meshes for variable coefficient reaction–subdiffusion problems with a weak initial singularity

High Order Finite Volume Schemes for Solving the Non-Conservative Convection Equations on the Unstructured Grids.

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

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations