[论文解读] On the Accuracy of the Non-Classical Transport Equation in 1-D Random Periodic Media

Richard Vasques, Kai Krycki|arXiv (Cornell University)|Dec 9, 2014

Advanced Mathematical Modeling in Engineering参考文献 10被引用 3

一句话总结

本论文首次对一维随机周期介质中的非经典输运方程进行了数值验证，其中由于存在空间相关性，粒子碰撞距离服从非指数分布。通过解析推导路径长度分布函数 $ p(s) $，计算相应的 $ \Sigma_t(s) $，并求解非经典输运方程，作者证明了在预测系综平均标量通量方面，该方法相较于原子混合模型具有更高的精度，尤其在低散射区域和源附近区域表现更优。

ABSTRACT

We present a first numerical investigation of the accuracy of the recently proposed {\em non-classical transport equation}. This equation contains an extra independent variable (the path-length $s$), and models particle transport taking place in random media in which a particle's distance-to-collision is {\em not} exponentially distributed. To solve the non-classical equation, one needs to know the $s$-dependent ensemble-averaged total cross section $Σ_t(s)$, or its corresponding path-length distribution function $p(s)$. We consider a 1-D spatially periodic system consisting of alternating solid and void layers, randomly placed in the infinite line. In this preliminary work, we assume transport in rod geometry: particles can move only in the directions $μ=\pm 1$. We obtain an analytical expression for $p(s)$, and use this result to compute the corresponding $Σ_t(s)$. Then, we proceed to solve the non-classical equation for different test problems. To assess the accuracy of these solutions, we produce "benchmark" results obtained by (i) generating a large number of physical realizations of the system, (ii) numerically solving the transport equation in each realization, and (iii) ensemble-averaging the solutions over all physical realizations. We show that the results obtained with the non-classical equation accurately model the ensemble-averaged scalar flux in this 1-D random system, generally outperforming the widely-used atomic mix model. We conclude by discussing plans to extend the present work to slab geometry, as well as to more general random mixtures.

研究动机与目标

• 评估非经典输运方程在模拟具有空间相关散射中心的一维随机周期介质中粒子输运时的准确性。
• 推导一维交替固态与空洞层结构中路径长度分布函数 $ p(s) $ 的解析表达式。
• 利用所推导的 $ \Sigma_t(s) $ 求解非经典输运方程，并将其解与通过大量物理实现的系综平均结果进行对比。
• 评估非经典模型在预测系综平均标量通量方面相对于广泛使用的原子混合模型的性能。
• 为未来向 slab 几何及更一般随机介质的扩展奠定基础。

提出的方法

• 建立一个在一维无限直线上随机排列的交替固态与空洞层构成的周期性空间系统。
• 假设为棒状几何，粒子仅沿 $ \mu = \pm 1 $ 方向运动，从而简化角度依赖性。
• 利用更新理论和更新-奖励定理，解析推导路径长度分布函数 $ p(s) $。
• 利用恒等式 $ \Sigma_t(s) = p(s) / \left(1 - \int_0^s p(s') ds' \right) $，从 $ p(s) $ 计算依赖于 $ s $ 的总截面 $ \Sigma_t(s) $。
• 采用带有 HLL Riemann 求解器的有限体积格式，对非经典输运方程进行数值求解。
• 通过模拟大量物理实现，分别求解每个实现中的经典输运方程，并对标量通量进行系综平均，生成基准解。

实验结果

研究问题

• RQ1 非经典输运方程能否准确预测具有空间相关散射的 1-D 随机周期介质中系综平均标量通量？
• RQ2 在不同散射比和源配置下，非经典模型与原子混合模型的准确性相比如何？
• RQ3 与原子混合模型相比，非经典模型是否能更好地保持源附近标量通量的振荡结构？
• RQ4 数值扩散对非经典模型与原子混合模型相对性能的影响如何？
• RQ5 能否将 $ p(s) $ 的解析推导方法推广至更复杂的非周期随机介质？

主要发现

• 非经典输运方程能够准确模拟 1-D 随机周期介质中的系综平均标量通量，其解与通过系综平均物理实现获得的基准结果高度一致。
• 在低散射比情况下（例如 $ c = 0.1 $），非经典模型优于原子混合模型，尤其在捕捉源附近最大通量和保持解的振荡形态方面表现更优。
• 如图 4 所示，在 $ c = 0.1 $ 的问题集 A 中，非经典模型相比原子混合模型将相对误差降低了最多达 10%（绝对值）。
• 原子混合模型在散射比 $ c $ 增大时精度有所提升，而非经典模型的精度在此条件下略有下降，可能归因于 HLL 格式中的数值扩散。
• 在远离源的区域，非经典模型与基准结果保持更好的定性一致性，尤其是在低散射情况下。
• 结果验证了非经典输运理论框架的正确性，并支持其在临界性计算中的应用，特别是在需要更高精度的中子增殖因子和本征函数估计时。