A 3D radiative transfer framework I. Non-local operator splitting and continuum scattering problems

被引:53
作者
Hauschildt, PH [1 ]
Baron, E
机构
[1] Hamburger Sternwarte, Gojenbergsweg 112, D-21029 Hamburg, Germany
[2] Univ Oklahoma, Homer L Dodge Dept Phys & Astron, Norman, OK 73019 USA
[3] Lawrence Berkeley Natl Lab, CRD, Berkeley, CA 94720 USA
基金
美国国家科学基金会;
关键词
radiative transfer;
D O I
10.1051/0004-6361:20053846
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
We describe a highly flexible framework to solve 3D radiation transfer problems in scattering dominated environments based on a long characteristics piece-wise parabolic formal solution and an operator splitting method. We find that the linear systems are efficiently solved with iterative solvers such as Gauss-Seidel and Jordan techniques. We use a sphere-in-a-box test model to compare the 3D results to 1D solutions in order to assess the accuracy of the method. We have implemented the method for static media, however, it can be used to solve problems in the Eulerian-frame for media with low velocity fields.
引用
收藏
页码:273 / 284
页数:12
相关论文
共 31 条
[1]  
ADAM J, 1990, ASTRON ASTROPHYS, V240, P541
[2]  
Anderson E., 1992, LAPACK User's Guide
[3]  
Auer L, 2003, ASTR SOC P, V288, P3
[4]  
AUER L, 1994, ASTRON ASTROPHYS, V292, P599
[5]  
Bendicho PF, 1997, ASTRON ASTROPHYS, V324, P161
[6]   ANGULAR QUADRATURE PERTURBATIONS IN RADIATIVE TRANSFER THEORY [J].
CANNON, CJ .
JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER, 1973, 13 (07) :627-633
[7]  
DURMUHL R, 1986, MATRIZEN ANWENDUNGEN, V2
[8]   A new scheme for multidimensional line transfer .3. A two-dimensional Lagrangian variable tensor method with discontinuous finite-element Sn transport [J].
Dykema, PG ;
Klein, RI ;
Castor, JI .
ASTROPHYSICAL JOURNAL, 1996, 457 (02) :892-921
[9]  
Golub GH, 1989, MATRIX COMPUTATIONS
[10]  
HAMANN WR, 1987, NUMERICAL RAD TRANSF, P35