Skip to main content
    • Aa
    • Aa

Kinetic Energy Preserving and Entropy Stable Finite Volume Schemes for Compressible Euler and Navier-Stokes Equations


Centered numerical fluxes can be constructed for compressible Euler equations which preserve kinetic energy in the semi-discrete finite volume scheme. The essential feature is that the momentum flux should be of the form where and are any consistent approximations to the pressure and the mass flux. This scheme thus leaves most terms in the numerical flux unspecified and various authors have used simple averaging. Here we enforce approximate or exact entropy consistency which leads to a unique choice of all the terms in the numerical fluxes. As a consequence novel entropy conservative flux that also preserves kinetic energy for the semi-discrete finite volume scheme has been proposed. These fluxes are centered and some dissipation has to be added if shocks are present or if the mesh is coarse. We construct scalar artificial dissipation terms which are kinetic energy stable and satisfy approximate/exact entropy condition. Secondly, we use entropy-variable based matrix dissipation flux which leads to kinetic energy and entropy stable schemes. These schemes are shown to be free of entropy violating solutions unlike the original Roe scheme. For hypersonic flows a blended scheme is proposed which gives carbuncle free solutions for blunt body flows. Numerical results for Euler and Navier-Stokes equations are presented to demonstrate the performance of the different schemes.

Corresponding author
Linked references
Hide All

This list contains references from the content that can be linked to their source. For a full set of references and notes please see the PDF or HTML where available.

[1] T. Barth and M. Ohlberger Finite volume methods: Foundation and analysis. In E. Stein , R. D. Borst , and T. J. Hughes , editors, Encyclopedia of Computational Mechanics, volume 1, chapter 15. Wiley, 2004.

[2] G. Blaisdell , E. Spyropoulos , and J. Qin The effect of the formulation of nonlinear terms on aliasing errors in spectral methods. Applied Numerical Mathematics, 21(3):207219, 1996.

[4] U. Fjordholm , S. Mishra , and E. Tadmor Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM Journal on Numerical Analysis, 50(2):544573, 2012.

[5] M. Gerritsen and P. Olsson Designing an efficient solution strategy for fluid flows: 1. A stable high orderfinite differencescheme and sharp shock resolution for the Euler equations. Journal of Computational Physics, 129(2):245262, 1996.

[6] C. Geuzaine and J.-F. Remacle Gmsh: A 3-d finite element mesh generator with built-in preand post-processing facilities. International Journal for Numerical Methods in Engineering, 79(11):13091331, 2009.

[7] A. Harten On the symmetric form of systems of conservation laws with entropy. Journal of Computational Physics, 49(1):151164, 1983.

[8] A. E. Honein and P. Moin Higher entropy conservation and numerical stability of compressible turbulence simulations. Journal of Computational Physics, 201(2):531545, 2004.

[9] T. Hughes , L. Franca , and M. Mallet A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier-Stokes equations and the second law of thermodynamics. Computer Methods in Applied Mechanics and Engineering, 54(2):223234, 1986.

[10] F. Ismail and P. L. Roe Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks. J. Comput. Phys., 228(15):54105436, Aug. 2009.

[12] A. Jameson Formulation of kinetic energy preserving conservative schemes for gas dynamics and direct numerical simulation of one-dimensional viscous compressible flow in a shock tube using entropy and kinetic energy preserving schemes. J. Sci. Comput., 34(2):188208, Feb. 2008.

[13] A. Jameson , W. Schmidt , and E. Turkel Numerical solution of the Euler equations by finite volume methods using Runge-Kutta time-stepping schemes. In AIAA 14th Fluid and Plasma Dynamic Conference, Palo Alto, California, June 1981.

[14] C. A. Kennedy and A. Gruber Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid. Journal of Computational Physics, 227(3):16761700, 2008.

[16] Y. Morinishi Skew-symmetric form of convective terms and fully conservative finite difference schemes for variable density low-mach number flows. Journal of Computational Physics, 229(2):276300, 2010.

[18] S. Pirozzoli Numerical methods for high-speed flows. Annual Review of Fluid Mechanics, 43(1):163194, 2011.

[19] J. J. Quirk A contribution to the great Riemann solver debate. International Journal for Numerical Methods in Fluids, 18(6):555574, 1994.

[20] P. L. Roe Approximate Riemann solvers, parameter vectors, and differenceschemes. Journal of Computational Physics, 43(2):357372, 1981.

[23] M. Shoeybi , M. Svard , F. E. Ham , and P. Moin An adaptive implicit-explicit scheme for the DNS and LES of compressible flows on unstructured grids. Journal of Computational Physics, 229(17):59445965, 2010.

[24] C.-W. Shu and S. Osher Efficient implementation of essentially non-oscillatory shock-capturing schemes. Journal of Computational Physics, 77(2):439471, 1988.

[25] P. K. Subbareddy and G. V. Candler A fully discrete, kinetic energy consistent finite-volume scheme for compressible flows. Journal of Computational Physics, 228(5):13471364, 2009.

[26] E. Tadmor The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Mathematics of Computation, 49(179):91103,1987.

[30] P. Woodward and P. Colella The numerical simulation of two-dimensional fluid flow with strong shocks. Journal of Computational Physics, 54(1):115173, 1984.

Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 19 *
Loading metrics...

Abstract views

Total abstract views: 177 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 25th June 2017. This data will be updated every 24 hours.