› 反应堆动特性模拟器 返回
物理模拟器

反应堆动特性·点堆模型模拟器

使用RK4对点堆动特性方程进行数值积分。实时可视化反应度投入、停堆、延迟中性子6群的行为。具有即发临界警告功能。

参数设置
预设
反应度 ρ [美元]
$
美元单位：ρ/$ (β=0.0065)
反应度投入类型
即发中性子寿命 Λ [μs]
μs
模拟时间 [s]
s
启用停堆
实时读数（点堆动力学）
1.00
功率比 n/n₀
0.00
反应度 ρ [$]
∞
反应堆周期 T [s]
—
倍增时间 [s]
次临界
瞬发/缓发状态
1.00
峰值功率比
实时点堆动力学动画
n(t)/n₀（中子密度） 先驱核 C/C₀ 无缓发中子（对照） 反应度 ρ(t)
理论·主要公式
$$\frac{dn}{dt}= \frac{\rho - \beta}{\Lambda}\,n + \sum_{i=1}^{6}\lambda_i c_i + S$$ $$\frac{dc_i}{dt}= \frac{\beta_i}{\Lambda}\,n - \lambda_i c_i \quad (i=1,\ldots,6)$$

反时间方程（in-hour equation）：$\dfrac{\rho}{\beta}= \dfrac{T}{\ell}+ \sum_i \dfrac{a_i}{1+\lambda_i T}$

即发临界条件：$\rho \geq \beta$（$1美元以上）→ 急速输出上升

反应堆动特性·点堆模型介绍

🙋
我看到这个模拟器可以显示反应堆的输出如何随时间变化。"点堆"是什么意思？我们不需要考虑空间上的变化吗？
🎓
完全正确。简单来说，点堆模型把整个反应堆当作一个"中性子密度均匀的单点"来处理。在实际应用中，当我们要快速预测调整控制棒时反应堆输出如何随时间变化时，这个模型是第一步工具。你可以用模拟器中的"反应度ρ"滑块改变反应堆状态，就能实时看到"这一点"的输出如何响应。
🙋
我看到一个警告说"即发临界"。这是很危险的状态吗？"ρ≧β"是什么意思？
🎓
这是一个极其重要的概念。当反应度ρ超过延迟中性子割合β时，不需要等待延迟中性子就能维持连锁反应，这叫"即发临界"。比如说，在这个工具里把反应度设为1美元（ρ=β）以上，再点击运行，你就会看到输出一下子飙升。在实际操作中，我们必须非常谨慎地控制，绝不能让反应度进入这个范围。
🙋
还有个"停堆"功能。那是做什么的？
🎓
那是紧急停止。在紧急情况下，我们会把所有控制棒一下子插入反应堆，投入一个很大的负反应度来立刻停止反应。在这个工具里，你可以启用"停堆"，设置停堆时刻和停堆反应度，就能看到输出如何急速下降。图表会显示中性子密度和延迟中性子先驱核浓度是怎样变化的。在实际的安全设计中，对停堆过程要进行非常详细的模拟。

常见问题

是的。投入正反应度会导致中性子密度增加，特别是当达到即发临界（ρ>β）时，会出现无法控制的急速上升（即发跳跃）。本模拟器会显示即发临界警告，现实中的反应堆也会出现同样现象，所以反应度投入必须谨慎进行。
本模拟器中，延迟中性子6群的衰变常数λᵢ和割合βᵢ是固定的（U-235标准值），无法修改。每一群对应具有不同半衰期的先驱核，这些会影响反应堆的响应速度和可控性。延迟中性子是支撑反应堆稳定性的关键因素。
把反应度滑块快速设为负值（比如-0.01）就能模拟停堆。即发中性子寿命Λ越短，中性子密度衰减就越快。屏幕上的图表能实时显示中性子密度和延迟中性子先驱核浓度的变化。
Λ越小，中性子密度的响应就越快，对反应度变化的灵敏度越高。比如Λ=10⁻⁴秒时响应很陡峭，Λ=10⁻³秒时响应就比较缓和。在实际反应堆设计中，Λ的值直接关系到可控性和安全性。用模拟器改变它，你能亲身感受到其影响。

实际应用

反应堆控制系统设计·验证：在调整控制棒时，用点堆模型快速评估反应堆输出的响应速度和稳定性。模拟器中改变"反应度投入类型"为阶跃、斜坡或正弦波，就是控制系统设计基础培训的一部分。

安全分析（RIA：反应度起始事故）：评估假想的反应度事故（比如控制棒误拉出）是否会导致即发临界，以及输出会上升多少，用点堆模型做初步筛选。工具的警告功能有助于直观理解这条"安全边界"。

起动试验时的堆周期测量：反应堆起动时，测量输出上升的时间周期（堆周期），通过反时间方程（inhour equation）推算反应度。这个模拟器能让你数值上体会这个物理过程。

控制棒价值的校准：逐步移动控制棒，测量反应堆输出变化，把数据拟合到点堆模型，由此确定单根控制棒的反应度价值（美元值）。模拟器为这项工作提供理论基础。

常见误解和注意事项

首先有个误解："点堆模型太简单，实际用不了"。虽然它确实给不出空间上的输出分布，但对理解控制棒操作导致的反应堆全体时间响应是最优工具。实务中，在用更复杂的三维代码前，我们经常用它快速检查参数的敏感性。再次，参数设置有陷阱。比如不小心把"即发中性子寿命Λ"从10⁻⁴秒改成10⁻⁵秒，反应堆的动特性会变成极快，控制会极其困难。这其实就是快中子堆和热中子堆根本区别的体现。还有，"1美元（ρ=β）是安全分界线"的观念是危险的。在这个模拟器中，β被当成常数，但实际反应堆中，随着燃烧核燃料成分会逐渐变化，β的值也会缓缓改变。所以模拟器中的"1美元"不一定对应实际的安全裕度。这点需要注意。

相关工具

热分解反应速率模拟器 免费在线热分解（热解）反应速率模拟器。输入活化能和频率因子，实时可视化转化率、TGA升温速率曲线、产物分布。适用于材料、化学工程教学和研究... 化学反应器设计（CSTR·PFR） 化学反应器设计工具"NovaSolver"，快速计算CSTR和PFR所需体积·比较。考虑反应级数和温度校正，实时可视化转化率和Levenspiel曲线。化工专业人士... 反应堆中性子扩散模拟器 反应堆中性子扩散模拟器"NovaSolver"基于单群扩散理论，实时计算并可视化反应堆临界条件和中性子通量分布。改变炉心设计参数，观察... 双金属温度计 模拟器 — 双层带的热变形 双金属温度计模拟器，用Timoshenko公式实时计算热膨胀系数不同的双层带在温度变化下的弯曲程度。改变温度、带长、厚度、膨胀系数差，观察曲率... 布雷顿循环 模拟器 — 燃气轮机热效率 布雷顿循环模拟器，以空气标准循环，从压力比、比热比、进气温度、涡轮进口温度，实时计算热效率、压缩机出口温度、涡轮出口温度、净功。燃气轮机... 窒息流 模拟器 — 收敛喷嘴的临界条件和质量流量 窒息流模拟器，从上游全压 P_0、全温 T_0、喉部面积 A*、背压 P_b，实时计算收敛喷嘴的临界压力比、最大质量流量、喉部音速。可视化窒息状态...

使用指南

  1. 设置反应度ρ（-0.005～+0.010 [Δk/k]），通过与延迟中性子割合β（0.0065）的大小关系判断临界性；
  2. 调整平均延迟中性子寿命λ（20～80 [s]），控制6群延迟中性子源导致的输出响应时间常数；
  3. 设置模拟结束时刻tEnd（10～300 [s]）和停堆投入时刻，点击运行按钮，使用Runge-Kutta 4阶法对点堆动方程积分；
  4. 实时监视输出比n/n₀、渐近堆周期T、倍增时间、临界判定（ρ≧β警告）。

具体计算例

原子反应堆输出100MW，投入反应度+0.003Δk/k的情况：由于β=0.0065，反应度未达即发临界，会显示稳定的输出上升。在λ_avg=35[s]条件下，初期30秒内由于延迟中性子的作用，渐近堆周期T≈45[s]，输出约增加15倍。设置停堆时刻为15秒，反应度骤然变为-0.008Δk/k，输出倍增时间转负，开始衰减，60秒后降至初始值的10%以下，显示减速特性。

实务中的注意事项