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 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.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Supplement Series
Medium: X Size: Article No. 4
["Article No. 4"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the implementation of a two-moment-based general-relativistic multigroup radiation transport module in theGeneral-relativisticmultigridnumerical (Gmunu) code. On top of solving the general-relativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth- and first-order moments of the radiations in the Eulerian-frame. An analytic closure relation is used to obtain the higher order moments and close the system. The finite-volume discretization has been adopted for the radiation moments. The advection in spatial space and frequency-space are handled explicitly. In addition, the radiation–matter interaction terms, which are very stiff in the optically thick region, are solved implicitly. The implicit–explicit Runge–Kutta schemes are adopted for time integration. We test the implementation with a number of numerical benchmarks from frequency-integrated to frequency-dependent cases. Furthermore, we also illustrate the astrophysical applications in hot neutron star and core-collapse supernovae modelings, and compare with other neutrino transport codes.

    more » « less
  2. This paper develops a tree-topological local mesh refinement (TLMR) method on Cartesian grids for the simulation of bio-inspired flow with multiple moving objects. The TLMR nests refinement mesh blocks of structured grids to the target regions and arrange the blocks in a tree topology. The method solves the time-dependent incompressible flow using a fractional-step method and discretizes the Navier-Stokes equation using a finite-difference formulation with an immersed boundary method to resolve the complex boundaries. When iteratively solving the discretized equations across the coarse and fine TLMR blocks, for better accuracy and faster convergence, the momentum equation is solved on all blocks simultaneously, while the Poisson equation is solved recursively from the coarsest block to the finest ones. When the refined blocks of the same block are connected, the parallel Schwarz method is used to iteratively solve both the momentum and Poisson equations. Convergence studies show that the algorithm is second-order accurate in space for both velocity and pressure, and the developed mesh refinement technique is benchmarked and demonstrated by several canonical flow problems. The TLMR enables a fast solution to an incompressible flow problem with complex boundaries or multiple moving objects. Various bio-inspired flows of multiple moving objects show that the solver can save over 80% computational time, proportional to the grid reduction when refinement is applied. 
    more » « less
  3. 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, 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. 
    more » « less
  4. This paper focuses on the laminar boundary layer startup process (momentum and thermal) in incompressible flows. The unsteady boundary layer equations can be solved via similarity analysis by normalizing the stream-wise (x), wall-normal (y) and time (t) coordinates by a variable η and τ, respectively. The resulting ODEs are solved by a finite difference explicit algorithm. This can be done for two cases: flat plate flow where the change in pressure are zero (Blasius solution) and wedge or Falkner-Skan flow where the changes in pressure can be favorable (FPG) or adverse (APG). In addition, transient passive scalar transport is examined by setting several Prandtl numbers in the governing equation at two different wall thermal conditions: isothermal and isoflux. Numerical solutions for the transient evolution of the momentum and thermal boundary layer profiles are compared with analytical approximations for both small times (unsteady flow) and large (steady-state flow) times. 
    more » « less
  5. Context. High-frequency very-long-baseline interferometry (VLBI) observations can now resolve the event-horizon-scale emission from sources in the immediate vicinity of nearby supermassive black holes. Future space-VLBI observations will access highly lensed features of black hole images – photon rings – that will provide particularly sharp probes of strong-field gravity. Aims. Focusing on the particular case of the supermassive black hole M 87*, our goal is to explore a wide variety of accretion flows onto a Kerr black hole and to understand their corresponding images and visibilities. We are particularly interested in the visibility on baselines to space, which encodes the photon ring shape and whose measurement could provide a stringent test of the Kerr hypothesis. Methods. We developed a fully analytical model of stationary, axisymmetric accretion flows with a variable disk thickness and a matter four-velocity that can smoothly interpolate between purely azimuthal rotation and purely radial infall. To determine the observational appearance of such flows, we numerically integrated the general-relativistic radiative transfer equation in the Kerr spacetime, taking care to include the effects of thermal synchrotron emission and absorption. We then Fourier transformed the resulting images and analyzed their visibility amplitudes along the directions parallel and orthogonal to the black hole spin projected on the observer sky. Results. Our images generically display a wedding cake structure composed of discrete, narrow photon rings ( n  = 1, 2, …) stacked on top of broader primary emission that surrounds a central brightness depression of model-dependent size. At 230 GHz, the n  = 1 ring is always visible, but the n  = 2 ring is sometimes suppressed due to absorption. At 345 GHz, the medium is optically thinner and the n  = 2 ring displays clear signatures in both the image and visibility domains. We also examine the thermal synchrotron emissivity in the equatorial plane and show that it exhibits an exponential dependence on the radius for the preferred M 87* parameters. Conclusions. The black hole shadow is a model-dependent phenomenon – even for diffuse, optically thin sources – and should not be regarded as a generic prediction of general relativity. Observations at 345 GHz are promising for future space-VLBI measurements of the photon ring shape, since at this frequency the signal of the n  = 2 ring persists despite the disk thickness and nonzero absorption featured in our models. Future work is needed to investigate whether this conclusion holds in a larger variety of reasonable models. 
    more » « less