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 positivity-preserving edge-centred finite volume scheme for heterogeneous and anisotropic diffusion problems on polygonal meshes

An Efficient Positivity-Preserving Finite Volume Scheme for the Nonequilibrium Three-Temperature Radiation Diffusion Equations 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 Positivity-Preserving and Free Energy Dissipative Hybrid Scheme for the Poisson-Nernst-Planck Equations on Polygonal and Polyhedral Meshes

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

An enhanced flux continuity three-dimensional finite element method for heterogeneous and anisotropic diffusion problems on general meshes

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

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

An asymptotic preserving scheme for the model on polygonal and conical meshes

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

A virtual element method for the steady-state Poisson-Nernst-Planck equations on polygonal meshes

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

Composite finite volume schemes for the diffusion equation on unstructured meshes

A locally conservative staggered least squares method on polygonal meshes

Edge multiscale finite element methods for semilinear parabolic problems with heterogeneous coefficients

A Graded Mesh Refinement for 2D Poisson's Equation on Non Convex Polygonal Domains

Simulating anisotropic diffusion processes with smoothed particle hydrodynamics

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

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids