参数设置
中子扩散方程
$$D\nabla^2\Phi - \Sigma_a\Phi + S = 0$$$k_{eff}= \frac{1}{1+L^2 B^2}$
结果 1
—
结果 2
—
计算结果
从单群和双群中子扩散方程的解析解计算堆芯中子通量分布和有效增殖系数keff。
$k_{eff}= \frac{1}{1+L^2 B^2}$
计算结果
稳态下，中子扩散过程由以下方程控制，它本质是一个平衡关系：中子泄漏出去的量 + 被吸收的量 = 中子源产生的量。
$$D\nabla^2\Phi - \Sigma_a\Phi + S = 0$$其中，$\Phi$是中子通量（中子数/cm²·s），$D$是扩散系数（cm），$\Sigma_a$是宏观吸收截面（1/cm），表示中子被吸收的概率，$S$是外中子源强度。$\nabla^2\Phi$描述了中子从高浓度区域向低浓度区域的扩散（泄漏）。
对于一个有限大小的裸堆（没有反射层），可以求解得到反应堆达到临界状态的条件，由此引出至关重要的有效增殖系数$k_{eff}$。
$$k_{eff}= \frac{1}{1+L^2 B^2}$$这里，$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∞。这个工具的优势，正是能让使用者切身感受“泄漏”与“材料”之间的博弈关系。
构成此模拟器核心的“扩散方程”，实际上也出现在许多原子堆物理之外的输运现象中。例如半导体工艺中的杂质扩散。在硅晶圆注入掺杂剂后，通过热处理进行扩散的过程，正是用扩散方程描述的。尽管存在无源项（$\nabla^2 C = (1/D) \partial C/\partial t$）等差异，但从数学角度看它们如同兄弟。
此外，求解中子通量分布钟形曲线的计算，在数学上与结构力学的特征值分析高度相似。将扩散方程变形可得 $\nabla^2\Phi + B_g^2\Phi = 0$ 这种亥姆霍兹方程，这与求解固定边界薄膜的振动模态（固有频率与振型）的方程形式完全相同。堆芯的“临界”对应方程的“特征值”，中子通量分布则对应“特征函数”。CAE工程师应当能在此处领悟到有限元求解器所处理问题的本质关联。
进一步拓展应用范围，光子扩散理论（生物光学成像）与地下水流及污染物迁移建模中同样使用类似方程。这意味着通过此工具学习的，正是“产生·扩散·消失的平衡”这一普遍物理概念的第一步。
下一步建议从学习“双群扩散理论”开始。将中子分为快群与热群，并显式引入两者间的“减速”过程。仅此一步，就能理解实际堆芯中关键的“能谱效应”和“慢化剂作用”。在体验NovaSolver后，可以尝试思考：“快群与热群的扩散长度哪个更大？”“为何轻水堆要将燃料布置成栅格结构？”等问题。
若希望深化数学背景，请务必扎实复习“分离变量法”与“边界条件处理”。本模拟器求解的球几何方程通过分离变量设 $\Phi(r) = R(r)$，转化为关于 $R(r)$ 的常微分方程求解。其解具有 $\sin(B_g r)/(B_g r)$ 的形式，并利用边界条件 $\Phi(a)=0$ 求得 $B_g = \pi/a$。若能逐步推导理解这个过程，就能掌握将其应用于圆柱堆等其他几何构型的方法。
最终推荐尝试阅读实际工程代码（如MONJU、SRAC、MVP等）的输入文件。其中定义了数十至数百个能群、详细的燃料成分、复杂的堆芯布置。你既能深切体会到NovaSolver“单群”“均匀堆芯”假设对现实进行了何等程度的简化，同时正因为学习了这种“简化的本质”，才能理解复杂输入数据中每个参数的意义。