Abstract
The gas dynamics equations, coupled with a static gravitational field, admit the hydrostatic balance where the flux produced by the pressure is exactly canceled by the gravitational source term. Many astrophysical problems involve the hydrodynamical evolution in a gravitational field, therefore it is essential to correctly capture the effect of gravitational force in the simulations. Improper treatment of the gravitational force can lead to a solution which either oscillates around the equilibrium, or deviates from the equilibrium after a long time run. In this paper we design high order well-balanced finite difference WENO schemes to this system, which can preserve the hydrostatic balance state exactly and at the same time can maintain genuine high order accuracy. Numerical tests are performed to verify high order accuracy, well-balanced property, and good resolution for smooth and discontinuous solutions. The main purpose of the well-balanced schemes designed in this paper is to well resolve small perturbations of the hydrostatic balance state on coarse meshes. The more difficult issue of convergence towards such hydrostatic balance state from an arbitrary initial condition is not addressed in this paper.
Similar content being viewed by others
References
Audusse, E., Bouchut, F., Bristeau, M.-O., Klein, R., Perthame, B.: A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows. SIAM J. Sci. Comput. 25, 2050–2065 (2004)
Bale, D.S., LeVeque, R.J., Mitran, S., Rossmanith, J.A.: A wave propagation method for conservation laws and balance laws with spatially varying flux functions. SIAM J. Sci. Comput. 24, 955–978 (2002)
Balsara, D.S., Shu, C.-W.: Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J. Comput. Phys. 160, 405–452 (2000)
Bermudez, A., Vazquez, M.E.: Upwind methods for hyperbolic conservation laws with source terms. Comput. Fluids 23, 1049–1071 (1994)
Botta, N., Klein, R., Langenberg, S., Lützenkirchen, S.: Well-balanced finite volume methods for nearly hydrostatic flows. J. Comput. Phys. 196, 539–565 (2004)
Jiang, G., Shu, C.-W.: Efficient implementation of weighted ENO schemes. J. Comput. Phys. 126, 202–228 (1996)
LeVeque, R.J.: Balancing source terms and flux gradients on high-resolution Godunov methods: the quasi-steady wave-propagation algorithm. J. Comput. Phys. 146, 346–365 (1998)
LeVeque, R.J., Bale, D.S.: Wave propagation methods for conservation laws with source terms. In: Proceedings of the 7th International Conference on Hyperbolic Problems, pp. 609–618 (1998)
Li, G., Lu, C., Qiu, J.: Hybrid well-balanced WENO schemes with different indicators for shallow water equations. J. Sci. Comput. In press (2011). doi:10.1007/s10915-011-9520-4
Liu, X.-D., Osher, S., Chan, T.: Weighted essentially nonoscillatory schemes. J. Comput. Phys. 115, 200–212 (1994)
Luo, J., Xu, K., Liu, N.: A well-balanced symplecticity-preserving gas-kinetic scheme for hydrodynamic equations under gravitational field. SIAM J. Sci. Comput. 33, 2356–2381 (2011)
Perthame, B., Simeoni, C.: A kinetic scheme for the Saint-Venant system with a source term. Calcolo 38, 201–231 (2001)
Roe, P.L.: Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys. 43, 357–372 (1981)
Shu, C.-W.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Quarteroni, A. (ed.) Advanced Numerical Approximation of Nonlinear Hyperbolic Equations. Lecture Notes in Mathematics, vol. 1697, pp. 325–432. Springer, Berlin (1998)
Shu, C.-W.: High order weighted essentially non-oscillatory schemes for convection dominated problems. SIAM Rev. 51, 82–126 (2009)
Shu, C.-W., Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys. 77, 439–471 (1988)
Slyz, A., Prendergast, K.H.: time independent gravitational fields in the BGK scheme for hydrodynamics. Astron. Astrophys. Suppl. Ser. 139, 199–217 (1999)
Tian, C.T., Xu, K., Chan, K.L., Deng, L.C.: A three-dimensional multidimensional gas-kinetic scheme for the Navier-Stokes equations under gravitational fields. J. Comput. Phys. 226, 2003–2027 (2007)
Xing, Y., Shu, C.-W.: High order finite difference WENO schemes with the exact conservation property for the shallow water equations. J. Comput. Phys. 208, 206–227 (2005)
Xing, Y., Shu, C.-W.: High order well-balanced finite difference WENO schemes for a class of hyperbolic systems with source terms. J. Sci. Comput. 27, 477–494 (2006)
Xu, K.: A well-balanced gas-kinetic scheme for the shallow-water equations with source terms. J. Comput. Phys. 178, 533–562 (2002)
Xu, K., Luo, J., Chen, S.: A well-balanced kinetic scheme for gas dynamic equations under gravitational field. Adv. Appl. Math. Mech. 2, 200–210 (2010)
Zingale, M., Dursi, L.J., ZuHone, J., Calder, A.C., Fryxell, B., Plewa, T., Truran, J.W., Caceres, A., Olson, K., Ricker, P.M., Riley, K., Rosner, R., Siegel, A., Timmes, F.X., Vladimirova, N.: Mapping initial hydrostatic models in Godunov codes. Astrophys. J. Suppl. Ser. 143, 539–565 (2002)
Author information
Authors and Affiliations
Corresponding author
Additional information
Dedicated to Professor Stanley Osher on the occasion of his 70th birthday.
Y. Xing research was sponsored by the Office of Advanced Scientific Computing Research; U.S. Department of Energy. The work was partially performed at the ORNL, which was managed by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725.
C.-W. Shu research was supported by ARO grant W911NF-11-1-0091 and NSF grant DMS-1112700.
Rights and permissions
About this article
Cite this article
Xing, Y., Shu, CW. High Order Well-Balanced WENO Scheme for the Gas Dynamics Equations Under Gravitational Fields. J Sci Comput 54, 645–662 (2013). https://doi.org/10.1007/s10915-012-9585-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-012-9585-8