A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
arXiv:1306.0010 · doi:10.1088/0067-0049/206/2/21
Abstract
We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas and radiation variables including the non-stiff source terms, and a local implicit method to integrate the stiff source terms. We adopt the M1 closure relation and include all leading source terms. We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. Our code is dimensionally unsplit in one, two, and three space dimensions and is parallelized using MPI. The streaming and diffusion limits are well-described by the M1 closure model, and our implementation shows excellent behavior for a problem with a concentrated radiation source containing both regimes simultaneously. Our operator-split method is ideally suited for problems with a slowly varying radiation field and dynamical gas flows, in which the effect of the RSLA is minimal. We present an analysis of the dispersion relation of RHD linear waves highlighting the conditions of applicability for the RSLA. To demonstrate the accuracy of our method, we utilize a suite of radiation and RHD tests covering a broad range of regimes, including RHD waves, shocks, and equilibria, which show second-order convergence in most cases. As an application, we investigate radiation-driven ejection of a dusty, optically thick shell in the interstellar medium (ISM). Finally, we compare the timing of our method with other well-known iterative schemes for the RHD equations. Our code implementation, Hyperion, is suitable for a wide variety of astrophysical applications and will be made freely available on the Web.
30 pages, 29 figures, accepted for publication in ApJS
References in corpus (10)
- Athena: A New Code for Astrophysical MHD
- The RAGE radiation-hydrodynamic code
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
- Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks
- Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
- A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
- A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor
- A Radiation Transfer Solver for Athena using Short Characteristics
- Implementation of Sink Particles in the Athena Code
- A radiation-hydrodynamics scheme valid from the transport to the diffusion limit
Cited by in corpus (76)
- Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers
- The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
- Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
- A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
- Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo
- When Feedback Fails: The Scaling and Saturation of Star Formation Efficiency
- A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
- The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert
- A New Multi-Energy Neutrino Radiation-Hydrodynamics Code in Full General Relativity and Its Application to Gravitational Collapse of Massive Stars
- AREPO-RT: Radiation hydrodynamics on a moving mesh
- An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
- A New Numerical Scheme for Cosmic Ray Transport
- Numerical Simulations of Turbulent Molecular Clouds Regulated by Reprocessed Radiation Feedback from Nascent Super Star Clusters
- Numerical Simulations of Turbulent, Molecular Clouds Regulated by Radiation Feedback Forces I: Star Formation Rate and Efficiency
- Outflows Driven by Quasars in High-Redshift Galaxies with Radiation Hydrodynamics
- Disruption of Molecular Clouds by Expansion of Dusty H II Regions
- Global Hydromagnetic Simulations of Protoplanetary Disks with Stellar Irradiation and Simplified Thermochemistry
- Monte Carlo Radiation Hydrodynamics with Implicit Methods
- Neutrino transport in general relativistic neutron star merger simulations
- Should One Use the Ray-by-Ray Approximation in Core-Collapse Supernova Simulations?
- An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
- Radiation hydrodynamics simulations of massive star cluster formation in giant molecular clouds
- Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields
- Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Method and Tests
- Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations
- First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
- Predicting FIR lines from simulated galaxies
- Radiation-driven outflows from and radiative support in dusty tori of active galactic nuclei
- Outflows Driven by Direct and Reprocessed Radiation Pressure in Massive Star Clusters
- Geometrically thick obscuration by radiation-driven outflow from magnetized tori of active galactic nuclei
- Cosmic-Ray Transport in Simulations of Star-forming Galactic Disks
- A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
- RIGEL: Simulating dwarf galaxies at solar mass resolution with radiative transfer and feedback from individual massive stars
- Numerical Simulations of Turbulent Molecular Clouds Regulated by Radiation Feedback Forces II: Radiation-Gas Interactions and Outflows
- Dusty Cloud Acceleration by Radiation Pressure in Rapidly Star-Forming Galaxies
- A general hybrid radiation transport scheme for star formation simulations on an adaptive grid
- Resolution Requirements and Resolution Problems in Simulations of Radiative Feedback in Dusty Gas
- The Elephant in the Room: The Importance of Where and When Massive Stars Form in Molecular Clouds
- Global 3D Radiation Hydrodynamic Simulations of Proto-Jupiter's Convective Envelope
- Infrared Radiation Feedback Does Not Regulate Star Cluster Formation
- Radiation-Hydrodynamics with MPI-AMRVAC: Flux-Limited Diffusion
- Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs
- A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes
- A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
- Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay
- Three-temperature radiation hydrodynamics with PLUTO: Tests and applications to protoplanetary disks
- VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
- From Piz Daint to the Stars: Simulation of Stellar Mergers using High-Level Abstractions
- Rank-3 Moment Closures in General Relativistic Neutrino Transport
- Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
- General-relativistic radiation transport scheme in I: Implementation of two-moment based multi-frequency radiative transfer and code tests
- Arepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
- 3D simulations of AGB stellar winds -- I. Steady winds and dust formation
- The Launching of Cosmic Ray Driven Outflows
- No self-shadowing instability in 2D radiation-hydrodynamical models of irradiated protoplanetary disks
- Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating
- 3D Radiation-Hydrodynamic Simulations Resolving Interior of Rapidly Accreting Primordial Protostar
- Radiation transport methods in star formation simulations
- Time Dependent Radiation Hydrodynamics on a Moving Mesh
- Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
- Three-temperature radiation hydrodynamics with PLUTO: Thermal and kinematic signatures of accreting protoplanets
- Toward realistic simulations of magneto-thermal winds from weakly-ionized protoplanetary disks
- CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
- An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics
- Non-LTE radiation hydrodynamics in PLUTO
- Dusty Cloud Acceleration with Multiband Radiation
- On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
- Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
- The time step constraint in radiation hydrodynamics
- A Reduced-Speed-of-Light Formulation of the Magnetohydrodynamic-Particle-in-Cell Method
- A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure based on Data-Driven Shape Function
- Radiation-magnetohydrodynamics with MPI-AMRVAC using flux-limited diffusion
- Star Formation for Predictive Primordial Galaxy Formation
- Numerical Methods for Simulating Star Formation
- Octo-Tiger's New Hydro Module and Performance Using HPX+CUDA on ORNL's Summit
- The 21-cm signature of X-ray heated halos around galaxies during cosmic dawn