核反应堆中子扩散

稳态下，中子扩散过程由以下方程控制，它本质是一个平衡关系：中子泄漏出去的量 + 被吸收的量 = 中子源产生的量。

其中，$\Phi$是中子通量（中子数/cm²·s），$D$是扩散系数（cm），$\Sigma_a$是宏观吸收截面（1/cm），表示中子被吸收的概率，$S$是外中子源强度。$\nabla^2\Phi$描述了中子从高浓度区域向低浓度区域的扩散（泄漏）。

对于一个有限大小的裸堆（没有反射层），可以求解得到反应堆达到临界状态的条件，由此引出至关重要的有效增殖系数$k_{eff}$。

这里，$L = \sqrt{D/\Sigma_a}$称为扩散长度，表示中子从产生到被吸收所穿行的平均直线距离。$B^2$是几何曲率，只与反应堆的形状和尺寸有关，比如对于球形堆，$B^2 = (\pi / R)^2$，$R$是堆芯半径。当$k_{eff}= 1$时，反应堆临界。

反应堆堆芯设计：工程师利用扩散理论计算不同燃料富集度、布置方式下的中子通量分布和$k_{eff}$，以确定最经济、安全的堆芯装载方案，确保反应堆在整个运行周期内都能保持临界。

屏蔽设计：通过计算中子从堆芯泄漏出来的通量，可以确定生物屏蔽层（如混凝土、水）需要多厚，才能将辐射剂量降低到安全水平，保护工作人员和环境。

反应性控制：控制棒、硼酸溶液等反应性控制手段的设计，都依赖于对$\Sigma_a$（吸收截面）的精确计算。通过插入控制棒改变局部$\Sigma_a$，从而快速、精确地调整$k_{eff}$，实现反应堆的启动、功率调节和停堆。

燃料管理：在核电站换料时，需要计算燃烧过的旧燃料组件和新燃料组件混合布置后的中子学特性。扩散方程的解可以帮助预测新燃料布置下的功率分布，避免出现局部功率过高的“热点”。

首先，请务必牢记“单群模型”并非万能。这终究只是处理“平均中子”的超简化模型。例如，快中子和热中子的扩散方式与吸收方式完全不同，但模型却将其合并为一类。因此，它虽适用于定性趋势和灵敏度分析，却无法用于实际反应堆的详细设计。下一步需要使用多群扩散代码。

其次，参数设置中容易陷入的误区是混淆“扩散长度 L”与“迁徙面积 M²”。模拟器中两者均有滑块控制，但L²表示中子“从减速结束到被吸收前”的平均移动距离平方。而M²表示“从核裂变产生到被吸收前（包含减速过程）”的距离平方。多数情况下存在关系式 M² ≒ L² + τ（τ为费米年龄）。例如，轻水堆中L为厘米量级，M为数十厘米量级。若忽视这个区别，可能会误解参数变更的物理意义。

最后，要警惕“只要增大k∞就必然能临界”的误解。虽然k∞确实表征材料特性，但当堆芯尺寸过小时，泄漏效应占主导地位，即使无限增大k∞也可能无法使keff达到1。反之，当尺寸超过某个阈值（临界尺寸）后，keff会逐渐趋近于k∞。这个工具的优势，正是能让使用者切身感受“泄漏”与“材料”之间的博弈关系。