一种动态截断的积分传输矩阵低阶方法用于求解中性粒子输运问题

《Annals of Nuclear Energy》：A dynamically-truncated integral transport matrix low-order method for solving neutral particle transport problems

　　

一种计算效率高的低阶中性粒子输运方法。一种动态截断的积分传输矩阵方法(ITMM)。一种动态截断算法，用于构建低阶积分传输矩阵(LITM)。

**论文解读：动态截断积分传输矩阵低阶方法（LITMM）用于中性粒子输运问题**

**研究背景与问题**
中性粒子输运方程（NTE）的求解在核工程等领域中至关重要，但其高维相空间离散导致计算量巨大。传统迭代方法在低粒子耗散问题（如高散射比、大光学厚度）中收敛缓慢，尤其在特征值问题中，随着优势比接近1，幂迭代效率显著下降。扩散近似作为低阶方法虽能降低计算成本，但在边界附近、吸收主导介质及强各向异性散射配置中精度有限。不连续因子（DF）方法通过引入界面通量不连续性来改进精度，但其生成依赖于计算密集的蒙特卡罗模拟或多群确定论计算，且在不同问题状态（如泄漏、能谱、控制棒布局）下需重新生成，导致计算成本高昂。因此，研究人员开发了一种新的低阶方法——低阶积分传输矩阵方法（LITMM），旨在避免现有方法的缺陷，同时保持计算效率与精度。该研究发表在《Annals of Nuclear Energy》，由G?khan Pediz和Yousry Y. Azmy完成。

**主要技术方法**
为开展研究，作者采用了以下关键技术方法：
1. **动态截断算法**：基于积分传输矩阵（ITM）元素随光学距离增大而衰减的特性，在构建过程中动态截断低于用户指定阈值的元素，生成稀疏的低阶积分传输矩阵（LITM）。该算法通过位掩码列表（SBML）实现高效扫描，避免全ITM构建的二次复杂度。
2. **AHOT-N0空间离散**：采用任意高阶节点输运方法（零阶局部展开）在二维笛卡尔几何下离散NTE，导出ITM的表达式。
3. **BiCGSTAB求解器**：用于求解LITM线性系统，替代传统源迭代。
4. **Arnoldi迭代**：用于临界问题中的特征值求解，通过投影到Hessenberg矩阵获取最大特征值（即有效增殖因子k_eff）。
5. **基准问题**：使用BIBLIS-2D和IAEA-2D压水堆基准问题（转换为输运模型）进行验证，并采用均匀S₆角求积组和真空边界条件。

**研究结果**
**4.1 固定源问题**
在均质和非均质方形网格上（40×40至160×160）测试，改变散射比（c）和总截面（Σ_t）。结果表明：LITM密度随光学厚度增加而降低（从1%至10%），验证了截断的合理性；误差（MARE）在10^-8–10^-7阈值下通常小于1%，但呈现非单调性，归因于AHOT-N0产生的负标量通量导致ITM元素符号变化。条件数随散射比增加而增大，导致误差上升。计算时间对比显示，LITMM在光学厚度足够大且问题规模较大时优于矩阵自由Krylov（MFK）方法。

**4.2 符号变化ITM元素对误差趋势的影响**
通过对比Step法（无负通量）与AHOT-N0，证实了非单调误差行为源于符号变化的ITM元素。Step法下误差随截断阈值减小单调递减，而AHOT-N0则非单调。

**4.3 LITM构建成本缩放**
在固定光学厚度下，构建时间随网格数N呈近似线性缩放（O(N^1.05)），远优于全ITM的O(N²)，证明了动态截断算法的可扩展性。

**4.4 内存需求评估**
LITMM需存储稀疏矩阵（COO格式），内存随问题规模线性增长（每非零元素约24字节）。与MFK相比，内存占用更高，但通过稀疏矩阵-向量乘积实现了计算加速。

**4.5 屏蔽输运问题**
采用屏蔽问题（源区与屏蔽区），在48×48至192×192网格上测试。LITMM能准确估计源远区通量数量级（即使通量低至10^-14），相对误差约14%，且在源区和屏蔽区内精度高。误差随网格细化而增大，但可通过减小截断阈值缓解。

**4.6 临界问题**
在BIBLIS-2D和IAEA-2D基准上，采用不同网格和截断阈值。LITMM计算的k_eff与高阶参考解误差在10^-7阈值下低于10 pcm，群通量分布MARE在芯区小于1%。计算时间对比表明，在光学厚度较大的网格（34×34、68×68）上，LITMM速度提升显著（最高8.4倍），且随着网格细化，优势减弱，但在136×136网格上仍优于MFK（部分情况）。LITM密度在快群中高于热群，因快群光学厚度较小。

**讨论与结论**
研究表明，LITMM作为一种新的低阶输运方法，通过截断ITM中低于阈值的元素，显著降低了矩阵构建和求解成本。在固定源、屏蔽和临界问题中，LITMM均能产生与高阶解高度一致的近似解，且计算效率随光学厚度增加而提升。动态截断算法实现了近线性的构建成本，使方法适用于大规模问题。然而，LITMM需要更多内存，且误差受截断阈值、散射比和光学厚度影响，需根据问题特征选择适当阈值。未来工作将扩展至高阶空间离散方法（如DGFEM），并开发基于问题特性的阈值选择准则。最终，LITMM可非侵入式地集成到现有输运代码中，作为后处理步骤，为中性粒子输运计算提供一种准确且高效的低阶替代方案。