A new scheme for multidimensional line transfer. III. A two-dimensional Lagrangian variable tensor method with discontinuous finite-element Sn transport
- University of California, Lawrence Livermore National Laboratory, Livermore, California 94550 (United States)
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.}
- Research Organization:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 279894
- Journal Information:
- Astrophysical Journal, Vol. 457, Issue 2; Other Information: PBD: Feb 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
A RADIATION TRANSFER SOLVER FOR ATHENA USING SHORT CHARACTERISTICS
Krylov Iterative Methods and the Degraded Effectiveness of Diffusion Synthetic Acceleration for Multidimensional S{sub N} Calculations in Problems with Material Discontinuities