[论文解读] An efficient algorithm for two-dimensional radiative transfer in axisymmetric circumstellar envelopes and disks

C. P. Dullemond, R. Turolla|arXiv (Cornell University)|Mar 30, 2000

Atmospheric Ozone and Climate参考文献 1被引用 9

一句话总结

本文提出扩展短特征（ESC）方法，一种在球坐标系下用于轴对称原行星环境的高效二维辐射转移算法，可实现从恒星表面到秒差距量级的精确建模。该方法通过专用的角向离散化和光束捆绑实现通量守恒，并利用加速λ迭代（ALI）和Ng算法实现快速收敛，其有效性已在星周包夹星云中的尘埃散射和坍缩原恒星云中非局域热动平衡（non-LTE）HCO+谱线形成问题中得到验证。

ABSTRACT

We present an algorithm for two-dimensional radiative transfer in axisymmetric, circumstellar media. The formal integration of the transfer equation is performed by a generalization of the short characteristics (SC) method to spherical coordinates. Accelerated Lambda Iteration (ALI) and Ng's algorithm are used to converge towards a solution. By taking a logarithmically spaced radial coordinate grid, the method has the natural capability of treating problems that span several decades in radius, in the most extreme case from the stellar radius up to parsec scale. Flux conservation is guaranteed in spherical coordinates by a particular choice of discrete photon directions and a special treatment of nearly-radially outward propagating radiation. The algorithm works well from zero up to very high optical depth, and can be used for a wide variety of transfer problems, including non-LTE line formation, dust continuum transfer and high temperature processes such as compton scattering. In this paper we focus on multiple scattering off dust grains and on non-LTE transfer in molecular and atomic lines. Line transfer is treated according to an ALI scheme for multi-level atoms/molecules, and includes both random and systematic velocity fields. The algorithms are implemented in a multi-purpose user-friendly radiative transfer program named RADICAL. We present two example computations: one of dust scattering in the Egg Nebula, and one of non-LTE line formation in rotational transitions of HCO$^{+}$ in a flattened protostellar collapsing cloud.

研究动机与目标

• 开发一种鲁棒、精确且计算高效的二维辐射转移算法，用于覆盖多个空间尺度的轴对称星周结构。
• 通过将短特征方法适配至球坐标系并采用专门的角向分箱处理径向传播的光线，确保在半径跨越多个数量级时实现通量守恒。
• 实现高动态范围建模——从恒星表面到秒差距量级包层——而无需网格加密或精度损失。
• 支持多种物理过程，包括非局域热动平衡谱线传输、尘埃散射和连续谱发射，适用于复杂非球对称几何结构。
• 提供一个用户友好、多功能的工具（RADICAL），用于原恒星、渐近巨星分支后星及活动星系核的天体物理建模。

提出的方法

• 该算法将短特征方法推广至球坐标系，实现对二维轴对称几何中辐射转移方程的精确形式解。
• 采用对数间距的径向网格，自然适应从恒星半径到秒差距尺度的问题。
• 对几乎径向向外传播的光线实施特殊处理，将其归入单一角向分箱，以保持径向通量守恒。
• 采用加速λ迭代（ALI）和Ng算法，加速非局域热动平衡及多能级谱线传输问题中的收敛速度。
• 针对谱线传输，实现了适用于多能级原子/分子的ALI方案，涵盖随机与系统性速度场。
• 核心算法已实现于RADICAL代码中，该程序为多用途、用户友好的天体物理辐射转移程序。

实验结果

研究问题

• RQ1 在球坐标系下，二维辐射转移算法能否在跨越多个数量级半径的高光学厚度与低光学厚度区域均实现精确建模？
• RQ2 在建模具有大径向动态范围的延伸星周包层时，如何在球坐标系中保持通量守恒？
• RQ3 扩展短特征方法能否在旋转、坍缩的原恒星云中实现非局域热动平衡谱线传输的快速且稳定的收敛？
• RQ4 该算法在建模复杂非球对称结构（如具有强散射与自吸收特征的双极星云）时表现如何？
• RQ5 该方法在多大程度上可推广以包含多样化的辐射过程，如尘埃发射、康普顿散射和热轫致辐射？

主要发现

• ESC方法在几个百分点内重现了1D球对称解，证实其精度与现有1D代码的一致性。
• 通过将几乎径向的光线捆绑至单一角向分箱，该算法在半径跨越多个数量级时实现了通量守恒，消除了径向通量误差。
• 在星云Egg Nebula的模拟中，模型成功再现了观测到的不对称双峰谱线轮廓及自吸收特征，这归因于环状结构遮挡了中心辐射。
• 在HCO+谱线传输案例中，模型显示侧向视线的谱线轮廓比极向视线更宽，这是由于赤道区域密度和激发温度更高所致。
• 坍缩云模型中的非局域热动平衡能级布居与激发温度表现出强烈的径向与角向依赖性，其中离心力半径位于1.5×10^15 cm处，标志着盘状结构的外缘。
• RADICAL代码成功以高保真度模拟了尘埃散射与非局域热动平衡谱线传输，展示了该方法在多样化天体物理问题中的通用性。