Numerical Methods for the Simulation of Dynamical Mass Transfer in Binaries

Abstract

We describe computational tools that have been developed to simulate dynamical mass transfer in semi-detached, polytropic binaries that are initially executing synchronous rotation upon circular orbits. Initial equilibrium models are generated with a self-consistent field algorithm; models are then evolved in time with a parallel, explicit, Eulerian hydrodynamics code with no assumptions made about the symmetry of the system. Poisson's equation is solved along with the equations of ideal fluid mechanics to allow us to treat the nonlinear tidal distortion of the components in a fully self-consistent manner. We present results from several standard numerical experiments that have been conducted to assess the general viability and validity of our tools, and from benchmark simulations that follow the evolution of two detached systems through five full orbits (up to approximately 90 stellar dynamical times). These benchmark runs allow us to gauge the level of quantitative accuracy with which simulations of semi-detached systems can be performed using presently available computing resources. We find that we should be able to resolve mass transfer at levels $\dot{M}/M>fewx10^{-}5$ per orbit through approximately 20 orbits with each orbit taking about 30 hours of computing time on parallel computing platforms.

34 pages, 20 eps figures, submitted to ApJS

References in corpus (2)

• Gamma-Ray Bursts: Afterglow Implications, Progenitor Clues and Prospects
• Testing Hydrodynamic Models of LMC X-4 with UV and X-ray Spectra

Cited by in corpus (29)

• Prelude to a double degenerate merger: the onset of mass transfer and its impact on gravitational waves and surface detonations
• Gravitational Waves from Gravitational Collapse
• One-armed Spiral Instability in a Low T/|W| Postbounce Supernova Core
• SNSPH: A Parallel 3-D Smoothed Particle Radiation Hydrodynamics Code
• Unexpected Dynamical Instabilities In Differentially Rotating Neutron Stars
• Numerical Simulations of the Onset and Stability of Dynamical Mass Transfer in Binaries
• Runaway Coalescence at the Onset of Common Envelope Episodes
• Evolution of Close White Dwarf Binaries
• The Stability of Double White Dwarf Binaries Undergoing Direct Impact Accretion
• Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?
• Disk-Planet Interaction Simulations: (I) Baroclinic Generation of Vortensity and Non-Axisymmetric Rossby-Wave-Instability
• White Dwarf Mergers on Adaptive Meshes I. Methodology and Code Verification
• Nonlinear Development of the Secular Bar-mode Instability in Rotating Neutron Stars
• Octo-Tiger: A New, 3D Hydrodynamic Code for Stellar Mergers that uses HPX Parallelisation
• Numerical evolutions of nonlinear r-modes in neutron stars
• Twin Binaries: Studies of Stability, Mass Transfer, and Coalescence
• From Piz Daint to the Stars: Simulation of Stellar Mergers using High-Level Abstractions
• A Numerical Method for Generating Rapidly Rotating Bipolytropic Structures in Equilibrium
• A Very Fast And Angular Momentum Conserving Tree Code
• Mass Transfer in Binary Stars using SPH. I. Numerical Method
• A Comparison of Grid-Based and SPH Binary Mass-Transfer and Merger Simulations
• Neutron stars formation and Core Collapse Supernovae
• Numerical Simulations of Mass Transfer in Binaries with Bipolytropic Components
• A Numerical Method for Studying Super-Eddington Mass Transfer in Double White Dwarf Binaries
• Further Evidence for an Elliptical Instability in Rotating Fluid Bars and Ellipsoidal Stars
• Hydrodynamic simulations of WD-WD mergers and the origin of RCB stars
• An approximate solver for Riemann and Riemann-like Ellipsoidal Configurations
• A multi-dimensional, adiabatic, hydrodynamics code for studying tidal excitation
• Halted-Pendulum Relaxation: Application to White Dwarf Binary Initial Data