BETHE-Hydro: An Arbitrary Lagrangian-Eulerian Multi-dimensional Hydrodynamics Code for Astrophysical Simulations
arXiv:0805.3356 · doi:10.1086/591272
Abstract
In this paper, we describe a new hydrodynamics code for 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 2-order convergence.
51 pages in emulateapj, including 25 figures, replace with version accepted to ApJS, corrected typos and included minor referee's comments
References in corpus (23)
- Relativistic Equation of State of Nuclear Matter for Supernova and Neutron Star
- On the Nature of Core Collapse Supernova Explosions
- Inside the Supernova: A Powerful Convective Engine
- Fundamental differences between SPH and grid methods
- Explosions of O-Ne-Mg Cores, the Crab Supernova, and Subluminous Type II-P Supernovae
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport
- Simulations of Magnetically-Driven Supernova and Hypernova Explosions in the Context of Rapid Rotation
- Two-dimensional hydrodynamic core-collapse supernova simulations with spectral neutrino transport. I. Numerical method and results for a 15 M_sun star
- Radiation hydrodynamics with neutrinos: Variable Eddington factor method for core-collapse supernova simulations
- A New Mechanism for Core-Collapse Supernova Explosions
- Shock Breakout in Core-Collapse Supernovae and its Neutrino Signature
- Probing the gravitational well: No supernova explosion in spherical symmetry with general relativistic Boltzmann neutrino transport
- Supernova Simulations with Boltzmann Neutrino Transport: A Comparison of Methods
- Improved Models of Stellar Core Collapse and Still no Explosions: What is Missing?
- Modeling Core-Collapse Supernovae in 3-Dimensions
- Features of the Acoustic Mechanism of Core-Collapse Supernova Explosions
- Two-dimensional, Time-dependent, Multi-group, Multi-angle Radiation Hydrodynamics Test Simulation in the Core-Collapse Supernova Context
- Conservative General Relativistic Radiation Hydrodynamics in Spherical Symmetry and Comoving Coordinates
- Simulating Self-gravitating Hydrodynamic Flows
- The First Second of a Type-II Supernova: Convection, Accretion, and Shock Propagation
- 3D Numerical Experimentation on the Core Helium Flash of low-mass Red Giants
- Block Structured Adaptive Mesh and Time Refinement for Hybrid, Hyperbolic + N-body Systems
- A New Algorithm for 2-D Transport for Astrophysical Simulations: I. General Formulation and Tests for the 1-D Spherical Case
Cited by in corpus (12)
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- A Model for Gravitational Wave Emission from Neutrino-Driven Core-Collapse Supernovae
- Criteria for Core-Collapse Supernova Explosions by the Neutrino Mechanism
- A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry
- Gravitational Wave Signals from 3D Neutrino Hydrodynamics Simulations of Core-Collapse Supernovae
- Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
- A Global Turbulence Model for Neutrino-Driven Convection in Core-Collapse Supernovae
- Using the PPML approach for constructing a low-dissipation, operator-splitting scheme for numerical simulations of hydrodynamic flows
- A Moving Mesh Hydrodynamic Solver for ChaNGa
- The origin of tail-like structures around protoplanetary disks
- Measuring the Effects of Artificial Viscosity in SPH Simulations of Rotating Fluid Flows