Skip to main content
×
Home

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

Abstract
Abstract

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.

Copyright
Corresponding author
*Corresponding author.Email:praveen@math.tifrbng.res.in
References
Hide All
[1]Barth T. and Ohlberger M.Finite volume methods: Foundation and analysis. In Stein E., Borst R. D., and Hughes T. J., editors, Encyclopedia of Computational Mechanics, volume 1, chapter 15. Wiley, 2004.
[2]Blaisdell G., Spyropoulos E., and Qin J.The effect of the formulation of nonlinear terms on aliasing errors in spectral methods. Applied Numerical Mathematics, 21(3):207219, 1996.
[3]Dumbser M. and Toro E. F.On universal Osher-type schemes for general nonlinear hyperbolic conservation laws. Commun. Comput. Phys., 10:635671, 2011.
[4]Fjordholm U., Mishra S., and Tadmor E.Arbitrarily high-order accurate entropy stable essentially nonoscillatory schemes for systems of conservation laws. SIAM Journal on Numerical Analysis, 50(2):544573, 2012.
[5]Gerritsen M. and Olsson P.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]Geuzaine C. and Remacle J.-F.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]Harten A.On the symmetric form of systems of conservation laws with entropy. Journal of Computational Physics, 49(1):151164, 1983.
[8]Honein A. E. and Moin P.Higher entropy conservation and numerical stability of compressible turbulence simulations. Journal of Computational Physics, 201(2):531545, 2004.
[9]Hughes T., Franca L., and Mallet M.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]Ismail F. and Roe P. L.Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks. J. Comput. Phys., 228(15):54105436, Aug. 2009.
[11]Ismail F., Roe P. L., and Nishikawa H.A proposed cure to the carbuncle phenomenon. In Computational Fluid Dynamics (2006): Proc. of Fourth International Conference on CFD, Ghent, Belgium, July 2006.
[12]Jameson A.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]Jameson A., Schmidt W., and Turkel E.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]Kennedy C. A. and Gruber A.Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid. Journal of Computational Physics, 227(3):16761700, 2008.
[15]Le P.Floch, Mercier J., and Rohde C.Fully discrete, entropy conservative schemes of arbitraryorder. SIAM Journal on Numerical Analysis, 40(5):19681992, 2002.
[16]Morinishi Y.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.
[17]Muller J.-D.On Triangles and flow. PhD thesis, The University of Michigan, 1996.
[18]Pirozzoli S.Numerical methods for high-speed flows. Annual Review of Fluid Mechanics, 43(1):163194, 2011.
[19]Quirk J. J.A contribution to the great Riemann solver debate. International Journal for Numerical Methods in Fluids, 18(6):555574, 1994.
[20]Roe P. L.Approximate Riemann solvers, parameter vectors, and differenceschemes. Journal of Computational Physics, 43(2):357372, 1981.
[21]Roe P. L.Affordable, entropy consistent flux functions. In Eleventh International Conference on Hyperbolic Problems: Theory, Numerics and Applications, Lyon, 2006.
[22]Scandaliato A. L. and Liou M.-S.AUSM-based high-order solution for Euler equations. Commun. Comput. Phys., 12:10961120, 2012.
[23]Shoeybi M., Svard M., Ham F. E., and Moin P.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]Shu C.-W. and Osher S.Efficient implementation of essentially non-oscillatory shock-capturing schemes. Journal of Computational Physics, 77(2):439471, 1988.
[25]Subbareddy P. K. and Candler G. V.A fully discrete, kinetic energy consistent finite-volume scheme for compressible flows. Journal of Computational Physics, 228(5):13471364, 2009.
[26]Tadmor E.The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Mathematics of Computation, 49(179):91103,1987.
[27]Tadmor E.Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems. Acta Numerica, pp. 451512, 2004.
[28] B.Thornber, Drikakis D., Williams R., and Youngs D.On entropy generation and dissipation of kinetic energy in high-resolution shock-capturing schemes. Journal of Computational Physics, 227(10):48534872, 2008.
[29]Leer B. van. Upwind and high-resolution methods for compressible flow: From donor cell to residual-distribution schemes. Commun. Comput. Phys., 1:192206, 2006.
[30]Woodward P. and Colella P.The numerical simulation of two-dimensional fluid flow with strong shocks. Journal of Computational Physics, 54(1):115173, 1984.
[31]Wu H., Shen L., and Shen Z.A hybrid numerical method to cure numerical shock instability. Commun. Comput. Phys., 8:12641271, 2010.
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? *
×

Keywords:

Metrics

Full text views

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

Abstract views

Total abstract views: 344 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 22nd November 2017. This data will be updated every 24 hours.