A new algorithm for two-dimensional transport for astrophysical simulations. I. General formulation and tests for the one-dimensional spherical case

Ivan - Hubeny, Adam Burrows

Research output: Contribution to journalArticle

35 Citations (Scopus)

Abstract

We derive new equations using the mixed-frame approach tor one- and two-dimensional (axisymmetnc) time-dependent radiation transport and the associated couplings with matter. Our formulation is multigroup and multiangle and includes anisotropic scattering, frequency (energy)-dependent scattering and absorption, complete velocity dependence to order v/c, rotation, and energy redistribution due to inelastic scattering. Hence, the "2D" realization is actually "6 1/2"-dimensional. The effects of radiation viscosity are automatically incorporated. Moreover, we develop accelerated lambda iteration, Krylov subspace (GMRES), discontinuous finite element, and Feautrier numerical methods for solving the equations and present the results of one-dimensional numerical tests of the new formalism. The virtues of the mixed-frame approach include simple velocity dependence with no velocity derivatives, straight characteristics, simple physical interpretation, and clear generalization to higher dimensions. Our treatment can be used for both photon and neutrino transport, but we focus on neutrino transport and applications to core-collapse supernova theory in the discussions and examples.

Original languageEnglish (US)
Pages (from-to)1458-1487
Number of pages30
JournalAstrophysical Journal
Volume659
Issue number2 I
DOIs
StatePublished - Apr 20 2007

Fingerprint

astrophysics
scattering
formulations
neutrinos
simulation
radiation transport
numerical method
iteration
supernovae
energy
inelastic scattering
viscosity
formalism
photons
radiation
test
effect

Keywords

  • Hydrodynamics
  • Methods: numerical
  • Neutrinos
  • Supernovae: general

ASJC Scopus subject areas

  • Space and Planetary Science

Cite this

@article{291606eb2544435188d635f4ee30e2a5,
title = "A new algorithm for two-dimensional transport for astrophysical simulations. I. General formulation and tests for the one-dimensional spherical case",
abstract = "We derive new equations using the mixed-frame approach tor one- and two-dimensional (axisymmetnc) time-dependent radiation transport and the associated couplings with matter. Our formulation is multigroup and multiangle and includes anisotropic scattering, frequency (energy)-dependent scattering and absorption, complete velocity dependence to order v/c, rotation, and energy redistribution due to inelastic scattering. Hence, the {"}2D{"} realization is actually {"}6 1/2{"}-dimensional. The effects of radiation viscosity are automatically incorporated. Moreover, we develop accelerated lambda iteration, Krylov subspace (GMRES), discontinuous finite element, and Feautrier numerical methods for solving the equations and present the results of one-dimensional numerical tests of the new formalism. The virtues of the mixed-frame approach include simple velocity dependence with no velocity derivatives, straight characteristics, simple physical interpretation, and clear generalization to higher dimensions. Our treatment can be used for both photon and neutrino transport, but we focus on neutrino transport and applications to core-collapse supernova theory in the discussions and examples.",
keywords = "Hydrodynamics, Methods: numerical, Neutrinos, Supernovae: general",
author = "Hubeny, {Ivan -} and Adam Burrows",
year = "2007",
month = "4",
day = "20",
doi = "10.1086/512179",
language = "English (US)",
volume = "659",
pages = "1458--1487",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "2 I",

}

TY - JOUR

T1 - A new algorithm for two-dimensional transport for astrophysical simulations. I. General formulation and tests for the one-dimensional spherical case

AU - Hubeny, Ivan -

AU - Burrows, Adam

PY - 2007/4/20

Y1 - 2007/4/20

N2 - We derive new equations using the mixed-frame approach tor one- and two-dimensional (axisymmetnc) time-dependent radiation transport and the associated couplings with matter. Our formulation is multigroup and multiangle and includes anisotropic scattering, frequency (energy)-dependent scattering and absorption, complete velocity dependence to order v/c, rotation, and energy redistribution due to inelastic scattering. Hence, the "2D" realization is actually "6 1/2"-dimensional. The effects of radiation viscosity are automatically incorporated. Moreover, we develop accelerated lambda iteration, Krylov subspace (GMRES), discontinuous finite element, and Feautrier numerical methods for solving the equations and present the results of one-dimensional numerical tests of the new formalism. The virtues of the mixed-frame approach include simple velocity dependence with no velocity derivatives, straight characteristics, simple physical interpretation, and clear generalization to higher dimensions. Our treatment can be used for both photon and neutrino transport, but we focus on neutrino transport and applications to core-collapse supernova theory in the discussions and examples.

AB - We derive new equations using the mixed-frame approach tor one- and two-dimensional (axisymmetnc) time-dependent radiation transport and the associated couplings with matter. Our formulation is multigroup and multiangle and includes anisotropic scattering, frequency (energy)-dependent scattering and absorption, complete velocity dependence to order v/c, rotation, and energy redistribution due to inelastic scattering. Hence, the "2D" realization is actually "6 1/2"-dimensional. The effects of radiation viscosity are automatically incorporated. Moreover, we develop accelerated lambda iteration, Krylov subspace (GMRES), discontinuous finite element, and Feautrier numerical methods for solving the equations and present the results of one-dimensional numerical tests of the new formalism. The virtues of the mixed-frame approach include simple velocity dependence with no velocity derivatives, straight characteristics, simple physical interpretation, and clear generalization to higher dimensions. Our treatment can be used for both photon and neutrino transport, but we focus on neutrino transport and applications to core-collapse supernova theory in the discussions and examples.

KW - Hydrodynamics

KW - Methods: numerical

KW - Neutrinos

KW - Supernovae: general

UR - http://www.scopus.com/inward/record.url?scp=34248328362&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34248328362&partnerID=8YFLogxK

U2 - 10.1086/512179

DO - 10.1086/512179

M3 - Article

VL - 659

SP - 1458

EP - 1487

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 2 I

ER -