Thumbnail
Access Restriction
Open

Author Dykema, P. G. ♦ Klein, R. I. ♦ Castor, J. I.
Source United States Department of Energy Office of Scientific and Technical Information
Content type Text
Language English
Subject Keyword PHYSICS ♦ STARS ♦ RADIANT HEAT TRANSFER ♦ TWO-DIMENSIONAL CALCULATIONS ♦ MANY-DIMENSIONAL CALCULATIONS ♦ TENSORS ♦ NUMERICAL SOLUTION ♦ TIN ♦ FINITE ELEMENT METHOD ♦ LTE ♦ COMPUTER CODES ♦ THREE-DIMENSIONAL CALCULATIONS ♦ ITERATIVE METHODS ♦ ATOMS ♦ LINE FORMATION ♦ TWO-LEVEL SYSTEMS
Abstract We describe a new code{emdash}ALTAIR{emdash}that has been developed to solve the time-dependent non-LTE problem for a multilevel atom in two-dimensional axisymmetric geometry on a Lagrangian mesh of arbitrary complexity. The method and results are presented here for a two-level atom. The extension to a multilevel atom is made by the equivalent two-level atom (ETLA) methods we have recently described in Paper II of this series. The source function of the line formed by the two-level atom is iterated to self-consistency with the radiation field, which is obtained from a system of coupled moment equations that employs an Eddington tensor closure using the double-splitting method described in Paper I of this series. The Eddington tensor itself is derived from a formal solution of the photon transport equation, which is based on a new discontinuous finite-element technique; the tensor moment system uses a continuous representation. The Eddington tensor is updated in an outer iteration loop, within which the double-splitting iteration is used to find the self-consistent source function for a given tensor. The resulting method shows very rapid convergence of the scattering iteration, even for optically thick regions with scattering albedo near unity. In addition, the spatial structure of the solution is rather insensitive to the shape of the spatial zones; the calculation of problems with zone aspect ratios of 1000 to 1 causes no difficulty. Application to several benchmark calculations, including radiative transfer in a nonorthogonal mesh, is discussed. We find that for two-dimensional problems of astrophysical interest, many angle rays are required to ensure an accurate solution. {copyright} {ital 1996 The American Astronomical Society.}
ISSN 0004637X
Educational Use Research
Learning Resource Type Article
Publisher Date 1996-02-01
Publisher Department Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publisher Place United States
Journal Astrophysical Journal
Volume Number 457
Issue Number 2
Organization Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)


Open content in new tab

   Open content in new tab