skip to main content

Title: Multigroup Radiation Magnetohydrodynamics Based on Discrete Ordinates including Compton Scattering

We present a formulation and numerical algorithm to extend the scheme for gray radiation magnetohydrodynamics (MHD) developed by Jiang to include the frequency dependence via the multigroup approach. The entire frequency space can be divided into an arbitrary number of groups in the lab frame, and we follow the time-dependent evolution of frequency-integrated specific intensities along discrete rays inside each group. Spatial transport of photons is done in the lab frame while all the coupling terms are solved in the fluid rest frame. Lorentz transformation is used to connect different frames. The radiation transport equation is solved fully implicitly in time while the MHD equations are evolved explicitly so that time step is not limited by the speed of light. A finite volume approach is used for transport in both spatial and frequency spaces to conserve the radiation energy density and momentum. The algorithm includes photon absorption, electron scattering, as well as Compton scattering, which is calculated by solving the Kompaneets equation. The algorithm is accurate for a wide range of optical depth conditions and can handle both radiation-pressure- and gas-pressure-dominated flows. It works for both Cartesian and curvilinear coordinate systems with adaptive mesh refinement. We provide a variety more » of test problems including a radiating sphere, shadow test, absorption of a moving gas, Bondi-type flows, as well as a collection of test problems for thermal and bulk Compton scattering. We also discuss examples where frequency dependence can make a big difference compared with the gray approach.

« less
Publication Date:
Journal Name:
The Astrophysical Journal Supplement Series
Page Range or eLocation-ID:
Article No. 4
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract We present STARFORGE (STAR FORmation in Gaseous Environments): a new numerical framework for 3D radiation MHD simulations of star formation that simultaneously follow the formation, accretion, evolution, and dynamics of individual stars in massive giant molecular clouds (GMCs) while accounting for stellar feedback, including jets, radiative heating and momentum, stellar winds, and supernovae. We use the GIZMO code with the MFM mesh-free Lagrangian MHD method, augmented with new algorithms for gravity, timestepping, sink particle formation and accretion, stellar dynamics, and feedback coupling. We survey a wide range of numerical parameters/prescriptions for sink formation and accretion and find very small variations in star formation history and the IMF (except for intentionally-unphysical variations). Modules for mass-injecting feedback (winds, SNe, and jets) inject new gas elements on-the-fly, eliminating the lack of resolution in diffuse feedback cavities otherwise inherent in Lagrangian methods. The treatment of radiation uses GIZMO’s radiative transfer solver to track 5 frequency bands (IR, optical, NUV, FUV, ionizing), coupling direct stellar emission and dust emission with gas heating and radiation pressure terms. We demonstrate accurate solutions for SNe, winds, and radiation in problems with known similarity solutions, and show that our jet module is robust to resolution and numerical details,more »and agrees well with previous AMR simulations. STARFORGE can scale up to massive (>105M⊙) GMCs on current supercomputers while predicting the stellar (≳ 0.1M⊙) range of the IMF, permitting simulations of both high- and low-mass cluster formation in a wide range of conditions.« less
  2. ABSTRACT We derive a consistent set of moment equations for cosmic ray (CR)-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts, we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly isotropic DF and/or strong-scattering regimes, and the arbitrarily anisotropic DF or free-streaming regimes, as well as allowing for anisotropic scattering and transport/magnetic field structure. We present the appropriate two-moment closure and equations for various choices of evolved variables, including the CR phase space DF f, number density n, total energy e, kinetic energy ϵ, and their fluxes or higher moments, and the appropriate coupling terms to the gas. We show that this naturally includes and generalizes a variety of terms including convection/fluid motion, anisotropic CR pressure, streaming, diffusion, gyro-resonant/streaming losses, and re-acceleration. We discuss how this extends previous treatments of CR transport including diffusion and moment methods and popular forms of the Fokker–Planck equation, as well as how this differs from the analogous M1-RHD equations. We also present two different methods for incorporating a reduced speed of light (RSOL) to reduce time-step limitations: In both, we carefully address where the RSOL (versus true c) must appear for themore »correct behaviour to be recovered in all interesting limits, and show how current implementations of CRs with an RSOL neglect some additional terms.« less

    Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of ‘resonant drag instabilities’ (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the non-linear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper, we consider conditions resembling giant molecular clouds (GMCs), H ii regions, and distributed starbursts, where optical depths are modest (≲1), single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing non-linear on time-scales that are much shorter than the characteristic times of the outflow.more »The instabilities produce filamentary and plume-like or ‘horsehead’ nebular morphologies that are remarkably similar to observed dust structures in GMCs and H ii regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas ‘behind,’ an order-unity fraction of the gas is always efficiently entrained by dust.

    « less
  4. Abstract Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly-understood. We explore this using high-resolution, multi-frequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultra-faint dwarf through Milky Way mass scales, including H+He photo-ionization; photo-electric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically-thin) and moments-based M1 (exact when optically-thick). The most important RFB channels on galaxy scales are photo-ionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (∼1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photo-ionization heating from the UV background suppressing accretion. In MW-mass galaxies, meta-galactic backgrounds have negligible effects; but local photo-ionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other “rapid” FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyper-dense, bound star clusters. This makes star formation more violent and “bursty” when SNe explode in these hyper-clustered objects: thus, including RFB “smoothes” SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IRmore »multiple-scattering is rare (negligible in dwarfs, $\sim 10\%$ of RP in massive galaxies): absorption occurs primarily in “normal” GMCs with AV ∼ 1.« less
  5. Abstract The acceleration and transport of energetic electrons during solar flares is one of the outstanding topics in solar physics. Recent X-ray and radio imaging and spectroscopy observations have provided diagnostics of the distribution of nonthermal electrons and suggested that, in certain flare events, electrons are primarily accelerated in the loop top and likely experience trapping and/or scattering effects. By combining the focused particle transport equation with magnetohydrodynamic (MHD) simulations of solar flares, we present a macroscopic particle model that naturally incorporates electron acceleration and transport. Our simulation results indicate that physical processes such as turbulent pitch-angle scattering can have important impacts on both electron acceleration in the loop top and transport in the flare loop, and their influences are highly energy-dependent. A spatial-dependent turbulent scattering with enhancement in the loop top can enable both efficient electron acceleration to high energies and transport of abundant electrons to the footpoints. We further generate spatially resolved synthetic hard X-ray (HXR) emission images and spectra, revealing both the loop-top and footpoint HXR sources. Similar to the observations, we show that the footpoint HXR sources are brighter and harder than the loop-top HXR source. We suggest that the macroscopic particle model provides new insightsmore »into understanding the connection between the observed loop-top and footpoint nonthermal emission sources by combining the particle model with dynamically evolving MHD simulations of solar flares.« less