PRECONDITIONED BI-CONJUGATE GRADIENT METHOD FOR RADIATIVE TRANSFER IN SPHERICAL MEDIA

被引:18
作者
Anusha, L. S. [1 ]
Nagendra, K. N. [1 ]
Paletou, F. [2 ]
Leger, L. [2 ]
机构
[1] Indian Inst Astrophys, Bangalore 560034, Karnataka, India
[2] Univ Toulouse, Lab Astrophys Toulouse Tarbes, CNRS, F-31400 Toulouse, France
关键词
line: formation; methods: numerical; radiative transfer; scattering; APPROXIMATE LAMBDA-OPERATORS; LINE FORMATION; SYMMETRIC SYSTEMS; ALI METHOD; FORMULATION; SCHEME;
D O I
10.1088/0004-637X/704/1/661
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
A robust numerical method called the Preconditioned Bi-Conjugate Gradient (Pre-BiCG) method is proposed for the solution of the radiative transfer equation in spherical geometry. A variant of this method called Stabilized Preconditioned Bi-Conjugate Gradient (Pre-BiCG-STAB) is also presented. These are iterative methods based on the construction of a set of bi-orthogonal vectors. The application of the Pre-BiCG method in some benchmark tests shows that the method is quite versatile, and can handle difficult problems that may arise in astrophysical radiative transfer theory.
引用
收藏
页码:661 / 671
页数:11
相关论文
共 30 条
[1]  
[Anonymous], 2000, Iterative Methods for Sparse Linear Systems
[2]   Very efficient methods for multilevel radiative transfer in atomic and molecular lines [J].
Asensio Ramos, A. ;
Trujillo-Bueno, J. .
RADIATIVE TRANSFER AND APPLICATIONS TO VERY LARGE TELESCOPES, 2006, 18 :25-+
[3]  
AUER L, 1991, NATO ADV SCI I C-MAT, V341, P9
[4]  
AUER L, 1994, ASTRON ASTROPHYS, V292, P599
[5]  
Auer L. H., 1984, Methods in radiative transfer, P237
[6]   ANGULAR QUADRATURE PERTURBATIONS IN RADIATIVE TRANSFER THEORY [J].
CANNON, CJ .
JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER, 1973, 13 (07) :627-633
[7]  
CHANDRASEKHAR S, 1934, MNRAS, V94, P522
[8]   On the accuracy of the ALI method for solving the radiative transfer equation [J].
Chevallier, L ;
Paletou, F ;
Rutily, B .
ASTRONOMY & ASTROPHYSICS, 2003, 411 (02) :221-227
[9]   Solving radiative transfer with line overlaps using Gauss-Seidel algorithms [J].
Daniel, F. ;
Cernicharo, J. .
ASTRONOMY & ASTROPHYSICS, 2008, 488 (03) :1237-1247
[10]   An implicit integral method to solve selected radiative transfer problems .4. The case of spherical geometry [J].
Gros, M ;
Crivellari, L ;
Simonneau, E .
ASTROPHYSICAL JOURNAL, 1997, 489 (01) :331-345