2024/02/23 by Chen Liu, Liu, Chen, Buzzard, Gregery T. +1 · 1 citation
Engineering · #65M12 #65M60 #65N30 #90C25 #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #Elasticity and Material Modeling #FOS: Mathematics #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2402.15645
openalex publication_date 2024/02/23 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We consider an optimization based limiter for enforcing positivity of internal energy in a semi-implicit scheme for solving gas dynamics equations. With Strang splitting, the compressible Navier-Stokes system is splitted into the compressible Euler equations, solved by the positivity-preserving Runge-Kutta discontinuous Galerkin (DG) method, and the parabolic subproblem, solved by Crank-Nicolson method with interior penalty DG method. Such a scheme is at most second order accurate in time, high order accurate in space, conservative, and preserves positivity of density. To further enforce the positivity of internal energy, we impose an optimization based limiter for the total energy variable to post process DG polynomial cell averages. The optimization based limiter can be efficiently implemented by the popular first order convex optimization algorithms such as the Douglas-Rachford splitting method if using the optimal algorithm parameters. Numerical tests suggest that the DG method with ℚk basis and the optimization-based limiter is robust for demanding low pressure problems such as high speed flows.