[论文解读] A 3D radiative transfer framework: XI. multi-level NLTE

P. H. Hauschildt, E. Baron|arXiv (Cornell University)|Apr 16, 2014

Atmospheric Ozone and Climate参考文献 22被引用 11

一句话总结

本文在PHOENIX代码中提出了一套完整的3D多级非局部热动平衡（NLTE）辐射转移框架，通过分层区域分解和MPI并行化，将成熟的1D方法扩展至3D球坐标系。该方法在各类测试案例中与1D结果高度一致，包括一个复杂的Fe I–III模型，证明了在现代超级计算机上实现具有详细原子物理过程的全3D NLTE计算现在在计算上是可行的。

ABSTRACT

Multi-level non-local thermodynamic equilibrium (NLTE) radiation transfer calculations have become standard throughout the stellar atmospheres community and are applied to all types of stars as well as dynamical systems such as novae and supernovae. Even today spherically symmetric 1D calculations with full physics are computationally intensive. We show that full NLTE calculations can be done with fully 3 dimensional (3D) radiative transfer. With modern computational techniques and current massive parallel computational resources, full detailed solution of the multi-level NLTE problem coupled to the solution of the radiative transfer scattering problem can be solved without sacrificing the micro physics description. We extend the use of a rate operator developed to solve the coupled NLTE problem in spherically symmetric 1D systems. In order to spread memory among processors we have implemented the NLTE/3D module with a hierarchical domain decomposition method that distributes the NLTE levels, radiative rates, and rate operator data over a group of processes so that each process only holds the data for a fraction of the voxels. Each process in a group holds all the relevant data to participate in the solution of the 3DRT problem so that the 3DRT solution is parallelized within a domain decomposition group. We solve a spherically symmetric system in 3D spherical coordinates in order to directly compare our well-tested 1D code to the 3D case. We compare three levels of tests: a) a simple H+He test calculation, b) H+He+CNO+Mg, c) H+He+Fe. The last test is computationally large and shows that realistic astrophysical problems are solvable now, but they do require significant computational resources. With presently available computational resources it is possible to solve the full 3D multi-level problem with the same detailed micro-physics as included in 1D modeling.

研究动机与目标

• 将成熟的1D多级NLTE辐射转移方法扩展至完整的3D球坐标系，并保持精确的微观物理细节。
• 实现对恒星大气、新星、超新星及其他动态系统中高保真度3D NLTE建模，其中3D结构与非LTE效应至关重要。
• 通过采用分层区域分解与MPI的可扩展、内存高效的并行化方法，克服3D NLTE计算的计算挑战。
• 在从H+He到Fe I–III的多种原子复杂度下，将3D NLTE框架与经过充分验证的1D结果进行对比验证。
• 证明利用当前超级计算资源求解具有详细模型原子的现实、大规模3D NLTE问题的可行性。

提出的方法

• 将1D速率算子形式化方法（Hauschildt 1993）推广至3D球坐标系，对每个体素求解耦合的NLTE速率方程。
• 实施分层区域分解策略，将NLTE能级、辐射速率及速率算子数据分布至MPI进程，最大限度降低每个进程的内存占用。
• 采用非局部算子分裂法求解包含线和连续谱散射的3D辐射转移方程，确保数值稳定性和准确性。
• 采用Menzel定义计算LTE粒子数密度，并求解包含电荷与电离平衡约束的完整速率方程组。
• 在每个区域分解组内并行化3DRT求解过程，使每个进程能够持有本地求解更新所需的所有数据。
• 利用高性能计算资源（Cray XE30、JUQUEEN、NERSC），最大测试案例使用49,248个MPI进程和2.6TB存储。

实验结果

研究问题

• RQ1 能否在与1D模型相同的微观物理细节下，实现完整的3D多级NLTE辐射转移框架？
• RQ2 在相同的物理设置和原子模型下，3D NLTE解与经过充分验证的1D解相比如何？
• RQ3 对于Fe I–III等复杂离子，3D NLTE计算的计算需求与扩展行为如何？
• RQ4 大气结构中的3D效应相较于1D近似，对光谱合成的影响有多大？
• RQ5 在当前超级计算基础设施下，对超新星或活跃恒星等天体物理系统进行真实3D NLTE建模是否可行？

主要发现

• 3D NLTE框架在所有测试案例中与1D结果高度一致，包括H+He、H+He+CNO+Mg以及H+He+Fe I–III，差异低于数值分辨率极限。
• 在H+He+Fe I–III测试案例中，模型包含57,232条谱线和24,100个束缚-自由跃迁，每轮迭代需约2.6TB存储，耗时2680秒（使用49,248个MPI进程）。
• 1D与3D通量谱之间的相对差异在大多数波长范围内低于2%，仅因有限角分辨率产生微小偏差。
• 包含50万个波长点的3D模型需约2.6TB存储以完整保存光谱数据，证明该框架可扩展至大规模原子模型。
• 该方法成功处理了Fe I–III等具有数千个能级和谱线的复杂离子，证明了对天体物理相关系统进行3D NLTE建模在计算上已可行。
• 分层区域分解与基于MPI的并行化实现了高效的内存分布，使现代超级计算机能够求解大规模3D NLTE问题。