Skip to main content

A new way to conserve total energy for Eulerian hydrodynamic simulations with self-gravity

Author(s): Jiang, Yan-Fei; Belyaev, Mikhail; Goodman, Jeremy J.; Stone, James M.

To refer to this page use:
Abstract: We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow. (C) 2012 Elsevier B.V. All rights reserved.
Publication Date: Feb-2013
Electronic Publication Date: 9-Aug-2012
Citation: Jiang, Yan-Fei, Belyaev, Mikhail, Goodman, Jeremy, Stone, James M. (2013). A new way to conserve total energy for Eulerian hydrodynamic simulations with self-gravity. NEW ASTRONOMY, 19 (48 - 55. doi:10.1016/j.newast.2012.08.002
DOI: doi:10.1016/j.newast.2012.08.002
ISSN: 1384-1076
Pages: 48 - 55
Type of Material: Journal Article
Journal/Proceeding Title: NEW ASTRONOMY
Version: Author's manuscript

Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.