[论文解读] Simplified P$_N$ Equations for Nonclassical Transport with Isotropic Scattering

Richard Vasques, Rachel Slaybaugh|arXiv (Cornell University)|Oct 14, 2016

Nuclear reactor physics and engineering参考文献 6被引用 5

一句话总结

本文通过渐近分析，为具有各向同性散射的非经典输运方程推导出一组新的扩散近似——非经典 SP_N 方程。该方法通过引入自由程分布的矩，将经典的简化 P_N (SP_N) 方程推广至非经典介质。数值结果表明，当系统趋近于扩散极限时，非经典 SP_N 解收敛于基准结果，且随着 N 和系统尺寸 M 的增加，精度提高。

ABSTRACT

An asymptotic analysis is used to derive a set of diffusion approximations to the nonclassical transport equation with isotropic scattering. These approximations are shown to reduce to the simplified P$_N$ equations under the assumption of classical transport, and therefore are labeled nonclassical SP$_N$ equations. In addition, the nonclassical SP$_N$ equations can be manipulated into a classical form with modified parameters, which can be implemented in existing SP$_N$ codes. Numerical results are presented for an one-dimensional random periodic system, validating the theoretical predictions.

研究动机与目标

• 为具有各向同性散射的非经典输运方程开发一种高阶渐近近似，适用于具有空间相关散射中心的随机介质。
• 通过引入自由程分布的矩，将经典的简化 P_N (SP_N) 框架扩展至非经典输运。
• 通过将非经典 SP_N 方程重写为具有修正参数的经典形式，确保与现有 SP_N 代码的兼容性。
• 在一维随机周期性系统中，通过集合平均基准解对理论框架进行数值验证。
• 为未来扩展至多维系统、角度相关自由程分布及非各向同性散射奠定基础。

提出的方法

• 对非经典输运方程进行渐近展开，使用小参数 ε，对总截面、散射比和源项进行缩放，以推导高阶近似。
• 通过在渐近级数展开中匹配至 O(ε^2) 阶项，推导出非经典 SP_1（扩散）、SP_2 和 SP_3 方程的显式形式。
• 为使非经典 SP_N 方程在 N = M 时有效，要求自由程分布 p(s) 的前 2M 阶原始矩有限。
• 将非经典 SP_N 方程重写为具有修正截面和吸收系数的经典 SP_N 形式，从而可使用标准 Marshak 真空边界条件。
• 采用一维随机周期性系统，包含交替的固体层和空腔层，层宽 ℓ = 0.5 以确保平均自由程尺度，M = 1/ε 控制扩散极限。
• 通过在多个随机系统实现上进行蒙特卡罗式平均生成基准解，使用 S_16 离散纵坐标和钻石差分法数值求解输运方程。

实验结果

研究问题

• RQ1 对具有各向同性散射的非经典输运方程进行高阶渐近展开，能否得到一组一致的扩散近似，从而推广经典 SP_N 方程？
• RQ2 当自由程分布退化为指数分布（经典输运情形）时，非经典 SP_N 方程的行为如何？
• RQ3 非经典 SP_N 方程能否被重写为具有修正参数的经典 SP_N 形式，以实现对现有 SP_N 代码的集成？
• RQ4 如渐近理论预测，随着系统尺寸增大并趋近于扩散极限，非经典 SP_N 近似的精度是否提高？
• RQ5 在 1D 随机周期性系统中，非经典 SP_N 解相对于集合平均基准结果的收敛行为如何？

主要发现

• 非经典 SP_N 方程作为具有各向同性散射的非经典输运方程的高阶渐近近似被推导出来，其有效性依赖于自由程分布的前 2N 阶原始矩有限。
• 当自由程分布为指数分布（经典输运情形）时，非经典 SP_N 方程精确退化为经典 SP_N 方程，证实与已有理论的一致性。
• 非经典 SP_N 方程可转化为具有修正总截面和吸收系数的经典 SP_N 形式，从而可使用标准 Marshak 真空边界条件，并可集成至现有 SP_N 代码中。
• 针对一维随机周期性系统的数值结果表明，x=0 处标量通量估计的相对误差随 M 增加而减小（即 ε→0），证实其收敛至扩散极限。
• 对于每个 N，误差随 M 增加而减小，且高阶 SP_N 近似（如 SP_3）的误差低于低阶近似（如 SP_1），验证了渐近理论中关于 N 增大时精度提升的预测。
• 基准解通过在多个随机系统实现上进行集合平均获得，采用 S_16 离散纵坐标和钻石差分法，作为误差评估的参考标准。