A Hybrid Advection Scheme for Conserving Angular Momentum on a Refined Cartesian Mesh
arXiv:1404.5942 · doi:10.1088/0067-0049/212/2/23
Abstract
We test a new "hybrid" scheme for simulating dynamical fluid flows in which cylindrical components of the momentum are advected across a rotating Cartesian coordinate mesh. This hybrid scheme allows us to conserve angular momentum to machine precision while capitalizing on the advantages offered by a Cartesian mesh, such as a straightforward implementation of mesh refinement. Our test focuses on measuring the real and imaginary parts of the eigenfrequency of unstable axisymmetric modes that naturally arise in massless polytropic tori having a range of different aspect ratios, and quantifying the uncertainty in these measurements. Our measured eigenfrequencies show good agreement with the results obtained from the linear stability analysis of Kojima (1986) and from nonlinear hydrodynamic simulations performed on a cylindrical coordinate mesh by Woodward et al. (1994). When compared against results conducted with a traditional Cartesian advection scheme, the hybrid scheme achieves qualitative convergence at the same or, in some cases, much lower grid resolutions and conserves angular momentum to a much higher degree of precision. As a result, this hybrid scheme is much better suited for simulating astrophysical fluid flows, such as accretion disks and mass-transferring binary systems.
arXiv comment47 pages, 19 figures, 3 tables, accepted for publication in ApJS
References in corpus (20)
- PLUTO: a Numerical Code for Computational Astrophysics
- Compact Remnant Mass Function: Dependence on the Explosion Mechanism and Metallicity
- V1309 Scorpii: merger of a contact binary
- Population synthesis for double white dwarfs I.Close detached systems
- Simulating Radiating and Magnetized Flows in Multi-Dimensions with ZEUS-MP
- Remnant evolution after a carbon-oxygen white dwarf merger
- Collisions of white dwarfs as a new progenitor channel for type Ia supernovae
- Testing Convergence for Global Accretion Disks
- The hot subdwarf B + white dwarf binary KPD 1930+2752 - a Supernova Type Ia progenitor candidate
- Remnants of Binary White Dwarf Mergers
- SNSPH: A Parallel 3-D Smoothed Particle Radiation Hydrodynamics Code
- A conservative orbital advection scheme for simulations of magnetized shear flows with the PLUTO code
- Numerical Simulations of the Onset and Stability of Dynamical Mass Transfer in Binaries
- Spectra of Type Ia Supernovae from Double Degenerate Mergers
- Application of the Cubed-Sphere Grid to Tilted Black-Hole Accretion Disks
- The Stability of Double White Dwarf Binaries Undergoing Direct Impact Accretion
- Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?
- The Relative Stability Against Merger of Close, Compact Binaries
- Constructing Synchronously Rotating Double White Dwarf Binaries
- A Numerical Method for Studying Super-Eddington Mass Transfer in Double White Dwarf Binaries
Cited by in corpus (6)
- GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
- FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
- How To Model Supernovae in Simulations of Star and Galaxy Formation
- White Dwarf Mergers on Adaptive Meshes I. Methodology and Code Verification
- Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for the Einstein Toolkit
- Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations