First-order positivity-preserving entropy stable spectral collocation scheme for the 3-D compressible Navier-Stokes equations

In this paper, we extend the positivity–preserving, entropy stable first–order finite volume-type scheme developed for the one-dimensional compressible Navier-Stokes equations in [1] to three spatial dimensions. The new firstorder scheme is provably entropy stable, design–order accurate for smooth solutions, and guarantees the pointwise positivity of thermodynamic variables for 3-D compressible viscous flows. Similar to the 1-D counterpart, the proposed scheme for the 3-D Navier-Stokes equations is discretized on Legendre-Gauss-Lobatto grids used for high-order spectral collocation methods. The positivity of density is achieved by adding an artificial dissipation in the form of the first-order Brenner-Navier-Stokes diffusion operator. Another distinctive feature of the proposed scheme is that the Navier–Stokes viscous terms are discretized by high–order spectral collocation summationby-parts operators. To eliminate time step stiffness caused by the high-order approximation of the viscous terms, the velocity and temperature limiters developed for the 1-D compressible Navier-Stokes equations in [1] are generalized to three spatial dimensions. These limiters bound the magnitude of velocity and temperature gradients and preserve the entropy stability and positivity properties of the baseline scheme. Numerical results are presented to demonstrate design-order accuracy and positivity-preserving properties of the new first-order scheme for 2-D and 3-D inviscid and viscous flows with strong shocks and contact discontinuities. Corresponding author. Department of Mathematics and Statistics, Tel.: +1 757 683 3423. E-mail address: nyamalee@odu.edu Preprint submitted to Elsevier November 8, 2021 ar X iv :2 11 1. 03 23 9v 1 [ m at h. N A ] 5 N ov 2 02 1