skip to main content

Title: AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-knownmore » shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.« less
Authors:
 [1] ;  [2] ;  [3]
  1. Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 (United States)
  2. Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 (United States)
  3. Canadian Institute for Theoretical Astrophysics, Toronto, ON M5S3H4 (Canada)
Publication Date:
OSTI Identifier:
22340206
Resource Type:
Journal Article
Resource Relation:
Journal Name: Astrophysical Journal, Supplement Series; Journal Volume: 213; Journal Issue: 1; Other Information: Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ABSORPTION; ALGORITHMS; ANALYTICAL SOLUTION; COMPUTERS; CONVERGENCE; DIFFUSION; INTERPOLATION; MAGNETOHYDRODYNAMICS; PARALLEL PROCESSING; QUADRATURES; RADIANT HEAT TRANSFER; SCATTERING; SOURCE TERMS; TIME DEPENDENCE; VELOCITY; VISIBLE RADIATION