[论文解读] A fast algorithm for radiative transport in isotropic media

Kui Ren, Rongting Zhang+1 more|arXiv (Cornell University)|Oct 4, 2016

Numerical methods in inverse problems参考文献 63被引用 4

一句话总结

该论文提出了一种快速算法，用于在各向同性介质中求解稳态辐射转移方程（ERT），通过将该方程重新表述为角向平均解的弗雷德霍姆积分方程，随后利用与GMRES结合的核无关快速多极方法（FMM）高效求解。该方法实现了线性计算复杂度，并且无需修改即可处理均匀和非均匀介质，通过避免在主要计算阶段进行角度离散化，相较于传统方法显著提升了计算速度。

ABSTRACT

We propose in this work a fast numerical algorithm for solving the equation of radiative transfer (ERT) in isotropic media. The algorithm has two steps. In the first step, we derive an integral equation for the angularly averaged ERT solution by taking advantage of the isotropy of the scattering kernel, and solve the integral equation with a fast multipole method (FMM). In the second step, we solve a scattering-free transport equation to recover the original ERT solution. Numerical simulations are presented to demonstrate the performance of the algorithm for both homogeneous and inhomogeneous media.

研究动机与目标

• 开发一种高效的数值方法，用于求解各向同性介质中的辐射转移方程（ERT），该方程由于高维性而仍具计算挑战性。
• 克服依赖角度离散化的传统确定性求解器所面临的高计算成本和内存需求。
• 通过统一框架实现在均匀和非均匀介质中快速求解ERT。
• 利用积分方程公式化和快速多极方法加速求解过程，同时保持精度。

提出的方法

• 将ERT重述为角向平均解 ∫Φ(x,v)d𝑣 的第二类弗雷德霍姆积分方程。
• 应用GMRES等迭代求解器求解所得的体积分方程。
• 使用基于切比雪夫插值的核无关快速多极方法（FMM）加速GMRES中的矩阵-向量乘积。
• 采用低秩近似技术压缩FMM中的相互作用项，降低计算成本。
• 在获得角向平均解后，通过求解无散射输运方程来恢复解的完整角依赖性。
• 通过缓存总吸收系数的线积分计算结果，减少在非均匀介质中重复计算的开销。

实验结果

研究问题

• RQ1 是否可以将各向同性介质中的ERT重述为一种避免在主要计算阶段进行角度离散化的体积分方程？
• RQ2 如何有效将快速多极方法应用于ERT的积分公式化，以实现与问题规模成线性关系的缩放？
• RQ3 该算法是否能在不改变结构的前提下，保持在均匀和非均匀介质中的精度与效率？
• RQ4 使用少量切比雪夫插值节点对FMM加速解的精度和收敛性有何影响？
• RQ5 该方法是否可作为复杂或各向异性介质中更通用输运求解器的预条件子？

主要发现

• 所提出的算法在空间问题规模上实现了线性计算复杂度，表现出在大规模问题上的可扩展性。
• 使用少量切比雪夫插值节点的FMM可对真实解提供高度精确的近似，表明其具有作为预条件子的潜力。
• 该方法在均匀和非均匀介质中均表现稳健，切换介质类型时无需修改算法结构。
• 积分方程中的核在均匀和非均匀介质中保持结构相同，简化了实现过程。
• 数值模拟证实了该方法的高效性和准确性，其性能在包括复杂散射和吸收分布在内的多种场景下得到验证。
• 该方法在计算最密集的阶段避免了角度离散化，显著降低了内存使用量，使得更大规模问题的求解成为可能。