VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

Abstract

We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation.

arXiv comment23 pages, 15 figures. Accepted for publication in MNRAS

References in corpus (64)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• Mass Loss: Its Effect on the Evolution and Fate of High-Mass Stars
• Theory of Core-Collapse Supernovae
• Theoretical Challenges in Galaxy Formation
• Star Clusters Across Cosmic Time
• Rosseland and Planck Mean Opacities for Protoplanetary Discs
• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• The RAGE radiation-hydrodynamic code
• RAMSES-RT: Radiation hydrodynamics in the cosmological context
• Multi--dimensional Cosmological Radiative Transfer with a Variable Eddington Tensor Formalism
• H II regions: Witnesses to massive star formation
• A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
• Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks
• A Module for Radiation Hydrodynamic Calculations With ZEUS-2D Using Flux-Limited Diffusion
• Stellar irradiated discs and implications on migration of embedded planets I: equilibrium discs
• Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods
• A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
• Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets
• Magnetohydronamic Evolution of HII Regions in Molecular Clouds: Simulation Methodology, Tests, and Uniform Media
• An Unstable Truth: How Massive Stars get their Mass
• A Radiation Transfer Solver for Athena using Short Characteristics
• AREPO-RT: Radiation hydrodynamics on a moving mesh
• Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation
• Beyond Flux-Limited Diffusion: Parallel Algorithms for Multidimensional Radiation Hydrodynamics
• A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
• Crash: A Block-Adaptive-Mesh Code for Radiative Shock Hydrodynamics - Implementation and Verification
• Hybrid Characteristics: 3D radiative transfer for parallel adaptive mesh refinement hydrodynamics
• A New Numerical Scheme for Cosmic Ray Transport
• Numerical Simulations of Turbulent Molecular Clouds Regulated by Reprocessed Radiation Feedback from Nascent Super Star Clusters
• Atmospheric Regimes and Trends on Exoplanets and Brown Dwarfs
• Transport and Accretion in Planet-Forming Disks
• CASTRO: A New Compressible Astrophysical Solver. II. Gray Radiation Hydrodynamics
• Radiation Magnetohydrodynamics In Global Simulations Of Protoplanetary Disks
• Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?
• A numerical model for multigroup radiation hydrodynamics
• On the stability of radiation-pressure-dominated cavities
• Formation and Evolution of Disks around Young Stellar Objects
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• A fast, robust, and simple implicit method for adaptive time-stepping on adaptive mesh-refinement grids
• Zooming in on Individual Star Formation: Low- and High-mass Stars
• Radiation-hydrodynamical simulations of massive star formation using Monte Carlo radiative transfer: I. Algorithms and numerical methods
• Radiation hydrodynamics integrated in the code PLUTO
• Radiation pressure driving of a dusty atmosphere
• Hybrid Adaptive Ray-Moment Method (HARM${}^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
• The reliability of approximate radiation transport methods for irradiated disk studies
• Saturation of the MRI in Strongly Radiation Dominated Accretion Disks
• Radiation pressure in super star cluster formation
• A Hybrid Godunov Method for Radiation Hydrodynamics
• A new hybrid radiative transfer method for massive star formation
• A general hybrid radiation transport scheme for star formation simulations on an adaptive grid
• Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics
• Radiation-Hydrodynamics with MPI-AMRVAC: Flux-Limited Diffusion
• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations
• Arepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
• Accretion Disk Coronae in High Luminosity Systems
• Radiation hydrodynamics including irradiation and adaptive mesh refinement with AZEuS. I. Methods
• Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
• A high-performance and portable asymptotic preserving radiation hydrodynamics code with the M1 model
• Gray Radiation Hydrodynamics with the FLASH Code for Astrophysical Applications

Cited by in corpus (15)

• Nitrogen-enriched, highly-pressurized nebular clouds surrounding a super star cluster at cosmic noon
• Outflows Driven by Direct and Reprocessed Radiation Pressure in Massive Star Clusters
• Infrared Radiation Feedback Does Not Regulate Star Cluster Formation
• How magnetic field and stellar radiative feedback influences the collapse and the stellar mass spectrum of a massive star forming clump
• The role of the turbulence driving mode for the Initial Mass Function
• Three Dimensional Natures of Massive Star Envelopes
• An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics
• Tree-based solvers for adaptive mesh refinement code FLASH -- III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources
• Population III star formation in the presence of turbulence, magnetic fields and ionizing radiation feedback
• Magnetic fields limit the mass of Population III stars even before the onset of protostellar radiation feedback
• Starlight-driven flared-staircase geometry in radiation hydrodynamic models of protoplanetary disks
• Protostellar disc structure and dynamics during star formation from cloud-scale initial conditions
• Radiation magnetohydrodynamics simulations of Population III star formation during the Epoch of Reionization
• SPRAY: A smoothed particle radiation hydrodynamics code for modeling high intensity laser-plasma interactions
• A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure