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.

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

A Positivity-preserving High Order Finite Volume Compact-WENO Scheme for Compressible Euler Equations

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

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

A New High-Order Compact Scheme Of Unstructured Finite Volume Method

On High Order Positivity-Preserving Well-Balanced Finite Volume Methods for the Euler Equations with Gravitation

High-order accurate structure-preserving finite volume schemes on adaptive moving meshes for shallow water equations: Well-balancedness and positivity

A positivity-preserving implicit-explicit scheme with high order polynomial basis for compressible Navier-Stokes equations

Outflow Positivity Limiting for Hyperbolic Conservation Laws. Part I: Framework and Recipe

High Order Well-Balanced Positivity-Preserving Scale-invariant AWENO Scheme for Euler Systems with Gravitational Field

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 low dissipation finite difference nested multi-resolution WENO scheme for Euler/Navier-Stokes equations

A new discretely divergence-free positivity-preserving high-order finite volume method for ideal MHD equations

Positivity-preserving discontinuous Galerkin methods with Lax-Wendroff time discretizations

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

Uniformly High-Order Structure-Preserving Discontinuous Galerkin Methods for Euler Equations with Gravitation: Positivity and Well-Balancedness

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

Finite difference alternative unequal-sized weighted essentially non-oscillatory schemes for hyperbolic conservation laws

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