[论文解读] A Unified Gas-kinetic Particle Method for Multiscale Photon Transport

Weiming Li, Chang Liu+3 more|arXiv (Cornell University)|Oct 14, 2018

Gas Dynamics and Kinetic Theory参考文献 30被引用 5

一句话总结

该论文提出了一种统一气体动力学粒子（UGKP）方法，用于多尺度光子输运，能够在单一基于粒子的框架中无缝捕捉粒子输运、碰撞和扩散区域。通过将粒子输运与碰撞通过源自动力学方程的时间演化通量函数相耦合，该方法实现了渐近保持性，准确恢复了光学厚极限下的扩散方程，并在自由输运区域实现了精确的粒子追踪，且在所有区域均通过基准问题验证了其有效性。

ABSTRACT

In this work, we present a unified gas-kinetic particle (UGKP) method for the simulation of multiscale photon transport. The multiscale nature of the particle method mainly comes from the recovery of the time evolution flux function in the unified gas-kinetic scheme (UGKS) through a coupled dynamic process of particle transport and collision. This practice improves the original operator splitting approach in the Monte Carlo method, such as the separated treatment of particle transport and collision. As a result, with the variation of the ratio between numerical time step and local photon's collision time, different transport physics can be fully captured in a single computation. In the diffusive limit, the UGKP method could recover the solution of the diffusion equation with the cell size and time step being much larger than the photon's mean free path and the mean collision time. In the free transport limit, it presents an exact particle tracking process as the original Monte Carlo method. In the transition regime, the weights of particle free transport and collision are determined by the ratio of local numerical time step to the photon's collision time. Several one-dimensional numerical examples covering all transport regimes from the optically thin to optically thick are computed to validate the accuracy and efficiency of the current scheme. In comparison with the $S_N$ discrete ordinate method, the UGKP method is based on particles and avoids the discretization of particle velocity space, which does not suffer from the ray effect.

研究动机与目标

• 开发一种统一的数值方法，能够在不需域分解或独立求解器的情况下，准确模拟从自由输运到扩散的所有区域的光子输运。
• 克服传统蒙特卡罗方法在扩散区域中的局限性，避免因过度追踪碰撞而导致计算成本过高。
• 将统一气体动力学格式（UGKS）扩展为基于粒子的公式，保留蒙特卡罗方法的优点（无射线效应），同时实现渐近保持性。
• 在耦合辐射-材料系统上验证该方法，确保在所有输运区域中能量交换和平衡行为的准确性。

提出的方法

• UGKP方法采用基于粒子的角向方向离散化，避免了确定性方法所需的速度空间离散化。
• 通过在时间上对动力学方程进行积分，构建多尺度通量函数，动态耦合粒子输运与碰撞过程。
• 通量的时间演化取决于数值时间步长与局部光子平均自由程的比值，从而实现区域间的自动过渡。
• 该方法使用反映自由输运与碰撞相对贡献的粒子权重，其由局部碰撞时间尺度决定。
• 采用统一的粒子演化格式替代算子分裂，确保在所有区域中对输运与碰撞的一致处理。
• 通过引入辐射与物质之间的能量交换项，将方法扩展至辐射-材料耦合，材料能量通过热弛豫模型演化。

实验结果

研究问题

• RQ1 基于粒子的方法能否在多尺度光子输运中实现渐近保持性，在光学厚极限下恢复扩散方程，且不受时间步长限制？
• RQ2 如何在单一算法中耦合粒子输运与碰撞，以在所有输运区域保持精度？
• RQ3 UGKP方法是否在自由输运区域保持了蒙特卡罗方法的精确粒子追踪行为，同时在扩散区域保持稳定？
• RQ4 该方法能否在所有区域（包括扩散极限下的平衡态）准确模拟耦合辐射-材料系统？
• RQ5 在多尺度输运问题中，UGKP方法与离散 ordinates 方法（DOM）和隐式蒙特卡罗方法相比，在精度和效率方面表现如何？

主要发现

• 在光学厚区域，UGKP方法即使在时间步长和单元尺寸远大于光子平均自由程和碰撞时间的情况下，仍能准确恢复扩散方程。
• 在自由输运区域（ε = 1），UGKP解与通过离散 ordinates 方法获得的参考解高度一致，证实了精确的粒子追踪能力。
• 在中间区域（ε = 10⁻²），UGKP方法与参考解表现出极佳的一致性，展示了在过渡尺度上的鲁棒性。
• 在扩散区域（ε = 10⁻⁴），UGKP解收敛至平衡扩散解，验证了渐近保持性。
• 在辐射-材料耦合中，UGKP方法在扩散区域实现了辐射与物质之间的能量平衡，t = 0.03时E与cu_r紧密匹配。
• 该方法避免了确定性离散 ordinates 方法固有的射线效应，并消除了在不同区域间使用域分解或混合求解器的需求。