[论文解读] On the Feynman-alpha formula for fast neutrons

Johan Anderson, Pal Lenard+1 more|arXiv (Cornell University)|May 4, 2011

Statistical Mechanics and Entropy参考文献 1被引用 3

一句话总结

本文基于柯尔莫哥洛夫前向（主方程）的两能群随机模型，推导出适用于快中子与热中子的广义费曼-alpha公式。模型整合了复合泊松源、中子吸收、漏出散射、裂变及探测过程，得到一个具有两个不同衰减速率的方差-均值公式，其相关行为由两个指数衰减项主导，特别适用于核保障与次临界系统诊断。

ABSTRACT

In this contribution, a stochastic theory for a branching process in a neutron population with two energy levels is investigated. In particular, a variance to mean or Feynman-alpha formula is derived in this generalized scenario using the Kolmogorov forward or master equation theory for the probabilities in a system with a compound Poisson source.

研究动机与目标

• 基于柯尔莫哥洛夫前向方程，发展一种包含两个能群（快中子与热中子）的中子种群动力学随机理论。
• 整合真实物理过程：复合泊松源、吸收、漏出散射、裂变及探测。
• 推导适用于快堆类系统脉冲中子实验的方差-均值（费曼-alpha）公式。
• 实现对差分衰减自询问（DDSI）与加速器驱动次临界系统等应用中中子相关性衰减行为的精确建模。
• 评估相关函数中双指数衰减结构在易裂变材料探测中的诊断潜力。

提出的方法

• 使用连续时间马尔可夫过程对中子种群建模，包含两类粒子：快中子（N₁）与热中子（N₂）。
• 建立联合概率 P(N₁, N₂, Z₁, t) 的主方程（柯尔莫哥洛夫前向方程），包含吸收、散射、裂变、探测与源注入引起的跃迁。
• 引入生成函数 G(X, Y, Z, t) 将主方程转化为控制阶乘矩时间演化的偏微分方程。
• 通过在 X=Y=Z=1 处对生成函数求导，推导一阶与二阶矩方程，引入源统计量（r₁, r₂）与裂变产物数（ν₁, ν₂）。
• 应用拉普拉斯变换求解矩方程，得到解中包含两个特征衰减速率 ω₁ 与 ω₂，其分母为 H(s) = s² + (ω₁+ω₂)s + ω₁ω₂。
• 构建最终的费曼-alpha公式：σ²(Z)/⟨Z⟩ = 1 + Y₁(1 - (1−e⁻ω₁T)/(ω₁T)) + Y₂(1 - (1−e⁻ω₂T)/(ω₂T))，其中 Y₁ 与 Y₂ 为涉及系统参数的复杂表达式。

实验结果

研究问题

• RQ1 与单能群情况相比，引入两个能群（快中子与热中子）如何改变费曼-alpha公式？
• RQ2 在复合泊松源与探测条件下，两能群中子系统中方差-均值比的函数形式为何？
• RQ3 由系统反应速率导出的衰减速率 ω₁ 与 ω₂ 如何影响中子相关函数的时间依赖行为？
• RQ4 在实际应用中，相关函数中的双指数成分在多大程度上可被观测到，尤其与罗斯基-alpha方法相比如何？
• RQ5 当测量时间 T → ∞ 时，费曼-alpha公式的渐近行为如何？其对系统参数的依赖关系如何？

主要发现

• 推导出的费曼-alpha公式表现出由特征方程 H(s) = s² + (ω₁+ω₂)s + ω₁ω₂ 所决定的两个独立指数衰减分量。
• 当 T → ∞ 时，费曼-alpha的渐近值由 Y₀ = Y₁ + Y₂ = q₂(λ_d λ₂ λ_R λ₂f)/(ω₁² ω₂²) 决定，其依赖于探测率与反应率。
• 数值结果表明，费曼-alpha曲线的最大值对衰减速率比值及热化率 λ_R 极为敏感，随着 λ_R 增大，渐近值亦随之升高。
• 在典型相关曲线中，双指数结构难以通过肉眼直接识别，表明与罗斯基-alpha方法相比，该方法在易裂变材料探测中可能不够直观。
• 该模型结果与后向柯尔莫哥洛夫方法结果一致，验证了两能群随机公式的自洽性。
• 该公式可实现对实验数据的精确曲线拟合，有望比DDSI方法更精确地提取系统参数（如 λ_R 与 λ₂f），尽管其诊断解释仍有待全面建立。