BETHE-hydro

An arbitrary Lagrangian-Eulerian multidimensional hydrodynamics code for astrophysical simulations

Jeremiah W. Murphy, Adam S. Burrows

Research output: Contribution to journalArticle

27 Citations (Scopus)

Abstract

In this paper, we describe a new hydrodynamics code for one- and two-dimensional (1D and 2D) astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian (ALE) approach, in which the gradient and divergence operators are made compatible using the support-operator method. We present 1D and 2D gravity solvers that are finite differenced using the support-operator technique, and the resulting system of linear equations are solved using the tridiagonal method for 1D simulations and an iterative multigrid-preconditioned conjugate-gradient method for 2D simulations. Rotational terms are included for 2D calculations using cylindrical coordinates. We document an incompatibility between a subcell pressure algorithm to suppress hourglass motions, and the subcell remapping algorithm and present a modified subcell pressure scheme that avoids this problem. Strengths of this code include a straightforward structure, enabling simple inclusion of additional physics packages, the ability to use a general equation of state, and most importantly, the ability to solve self-gravitating hydrodynamic flows on time-dependent, arbitrary grids. In what follows, we describe in detail the numerical techniques employed and,with a large suite of tests, demonstrate that BETHE-hydro finds accurate solutions with second-order convergence.

Original languageEnglish (US)
Pages (from-to)209-241
Number of pages33
JournalAstrophysical Journal, Supplement Series
Volume179
Issue number1
DOIs
StatePublished - Nov 1 2008

Fingerprint

astrophysics
hydrodynamics
operators
simulation
conjugate gradient method
incompatibility
cylindrical coordinates
linear equations
equation of state
divergence
equations of state
physics
grids
inclusions
gravity
gravitation
gradients
method
code
document

All Science Journal Classification (ASJC) codes

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

@article{be59e86339634f1a80dc9f8fd99529f2,
title = "BETHE-hydro: An arbitrary Lagrangian-Eulerian multidimensional hydrodynamics code for astrophysical simulations",
abstract = "In this paper, we describe a new hydrodynamics code for one- and two-dimensional (1D and 2D) astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian (ALE) approach, in which the gradient and divergence operators are made compatible using the support-operator method. We present 1D and 2D gravity solvers that are finite differenced using the support-operator technique, and the resulting system of linear equations are solved using the tridiagonal method for 1D simulations and an iterative multigrid-preconditioned conjugate-gradient method for 2D simulations. Rotational terms are included for 2D calculations using cylindrical coordinates. We document an incompatibility between a subcell pressure algorithm to suppress hourglass motions, and the subcell remapping algorithm and present a modified subcell pressure scheme that avoids this problem. Strengths of this code include a straightforward structure, enabling simple inclusion of additional physics packages, the ability to use a general equation of state, and most importantly, the ability to solve self-gravitating hydrodynamic flows on time-dependent, arbitrary grids. In what follows, we describe in detail the numerical techniques employed and,with a large suite of tests, demonstrate that BETHE-hydro finds accurate solutions with second-order convergence.",
author = "Murphy, {Jeremiah W.} and Burrows, {Adam S.}",
year = "2008",
month = "11",
day = "1",
doi = "https://doi.org/10.1086/591272",
language = "English (US)",
volume = "179",
pages = "209--241",
journal = "Astrophysical Journal, Supplement Series",
issn = "0067-0049",
publisher = "IOP Publishing Ltd.",
number = "1",

}

BETHE-hydro : An arbitrary Lagrangian-Eulerian multidimensional hydrodynamics code for astrophysical simulations. / Murphy, Jeremiah W.; Burrows, Adam S.

In: Astrophysical Journal, Supplement Series, Vol. 179, No. 1, 01.11.2008, p. 209-241.

Research output: Contribution to journalArticle

TY - JOUR

T1 - BETHE-hydro

T2 - An arbitrary Lagrangian-Eulerian multidimensional hydrodynamics code for astrophysical simulations

AU - Murphy, Jeremiah W.

AU - Burrows, Adam S.

PY - 2008/11/1

Y1 - 2008/11/1

N2 - In this paper, we describe a new hydrodynamics code for one- and two-dimensional (1D and 2D) astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian (ALE) approach, in which the gradient and divergence operators are made compatible using the support-operator method. We present 1D and 2D gravity solvers that are finite differenced using the support-operator technique, and the resulting system of linear equations are solved using the tridiagonal method for 1D simulations and an iterative multigrid-preconditioned conjugate-gradient method for 2D simulations. Rotational terms are included for 2D calculations using cylindrical coordinates. We document an incompatibility between a subcell pressure algorithm to suppress hourglass motions, and the subcell remapping algorithm and present a modified subcell pressure scheme that avoids this problem. Strengths of this code include a straightforward structure, enabling simple inclusion of additional physics packages, the ability to use a general equation of state, and most importantly, the ability to solve self-gravitating hydrodynamic flows on time-dependent, arbitrary grids. In what follows, we describe in detail the numerical techniques employed and,with a large suite of tests, demonstrate that BETHE-hydro finds accurate solutions with second-order convergence.

AB - In this paper, we describe a new hydrodynamics code for one- and two-dimensional (1D and 2D) astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian (ALE) approach, in which the gradient and divergence operators are made compatible using the support-operator method. We present 1D and 2D gravity solvers that are finite differenced using the support-operator technique, and the resulting system of linear equations are solved using the tridiagonal method for 1D simulations and an iterative multigrid-preconditioned conjugate-gradient method for 2D simulations. Rotational terms are included for 2D calculations using cylindrical coordinates. We document an incompatibility between a subcell pressure algorithm to suppress hourglass motions, and the subcell remapping algorithm and present a modified subcell pressure scheme that avoids this problem. Strengths of this code include a straightforward structure, enabling simple inclusion of additional physics packages, the ability to use a general equation of state, and most importantly, the ability to solve self-gravitating hydrodynamic flows on time-dependent, arbitrary grids. In what follows, we describe in detail the numerical techniques employed and,with a large suite of tests, demonstrate that BETHE-hydro finds accurate solutions with second-order convergence.

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

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

U2 - https://doi.org/10.1086/591272

DO - https://doi.org/10.1086/591272

M3 - Article

VL - 179

SP - 209

EP - 241

JO - Astrophysical Journal, Supplement Series

JF - Astrophysical Journal, Supplement Series

SN - 0067-0049

IS - 1

ER -