Computational Eulerian hydrodynamics and Galilean invariance

Brant E Robertson, Andrey V. Kravtsov, Nickolay Y. Gnedin, Tom Abel, Douglas H. Rudd

Research output: Contribution to journalArticle

51 Citations (Scopus)

Abstract

Eulerian hydrodynamical simulations are a powerful and popular tool for modelling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

Original languageEnglish (US)
Pages (from-to)2463-2476
Number of pages14
JournalMonthly Notices of the Royal Astronomical Society
Volume401
Issue number4
DOIs
StatePublished - Feb 2010
Externally publishedYes

Fingerprint

invariance
hydrodynamics
Kelvin-Helmholtz instability
simulation
astrophysics
spatial resolution
method
fluids
high resolution
fluid
modeling
effect

Keywords

  • Hydrodynamics
  • Instabilities
  • Methods: numerical

ASJC Scopus subject areas

  • Space and Planetary Science
  • Astronomy and Astrophysics

Cite this

Computational Eulerian hydrodynamics and Galilean invariance. / Robertson, Brant E; Kravtsov, Andrey V.; Gnedin, Nickolay Y.; Abel, Tom; Rudd, Douglas H.

In: Monthly Notices of the Royal Astronomical Society, Vol. 401, No. 4, 02.2010, p. 2463-2476.

Research output: Contribution to journalArticle

Robertson, Brant E ; Kravtsov, Andrey V. ; Gnedin, Nickolay Y. ; Abel, Tom ; Rudd, Douglas H. / Computational Eulerian hydrodynamics and Galilean invariance. In: Monthly Notices of the Royal Astronomical Society. 2010 ; Vol. 401, No. 4. pp. 2463-2476.
@article{cff100a87f27402b944cb042212498cb,
title = "Computational Eulerian hydrodynamics and Galilean invariance",
abstract = "Eulerian hydrodynamical simulations are a powerful and popular tool for modelling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.",
keywords = "Hydrodynamics, Instabilities, Methods: numerical",
author = "Robertson, {Brant E} and Kravtsov, {Andrey V.} and Gnedin, {Nickolay Y.} and Tom Abel and Rudd, {Douglas H.}",
year = "2010",
month = "2",
doi = "10.1111/j.1365-2966.2009.15823.x",
language = "English (US)",
volume = "401",
pages = "2463--2476",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "4",

}

TY - JOUR

T1 - Computational Eulerian hydrodynamics and Galilean invariance

AU - Robertson, Brant E

AU - Kravtsov, Andrey V.

AU - Gnedin, Nickolay Y.

AU - Abel, Tom

AU - Rudd, Douglas H.

PY - 2010/2

Y1 - 2010/2

N2 - Eulerian hydrodynamical simulations are a powerful and popular tool for modelling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

AB - Eulerian hydrodynamical simulations are a powerful and popular tool for modelling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

KW - Hydrodynamics

KW - Instabilities

KW - Methods: numerical

UR - http://www.scopus.com/inward/record.url?scp=74549179771&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=74549179771&partnerID=8YFLogxK

U2 - 10.1111/j.1365-2966.2009.15823.x

DO - 10.1111/j.1365-2966.2009.15823.x

M3 - Article

VL - 401

SP - 2463

EP - 2476

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 4

ER -