逆时方程(in-hour equation):$\dfrac{\rho}{\beta}= \dfrac{T}{\ell}+ \sum_i \dfrac{a_i}{1+\lambda_i T}$
瞬发临界条件:$\rho \geq \beta$(超过1美元)→ 功率急速上升
用RK4数值积分点堆动力学方程。实时可视化反应性引入、紧急停堆及缓发中子六群的行为。含瞬发临界预警。
逆时方程(in-hour equation):$\dfrac{\rho}{\beta}= \dfrac{T}{\ell}+ \sum_i \dfrac{a_i}{1+\lambda_i T}$
瞬发临界条件:$\rho \geq \beta$(超过1美元)→ 功率急速上升
点堆模型的核心是描述中子密度变化率的方程。它考虑了瞬发中子、六组缓发中子以及外中子源的影响。
$$\frac{dn}{dt}= \frac{\rho - \beta}{\Lambda}\,n + \sum_{i=1}^{6}\lambda_i c_i + S$$n:中子密度(与功率成正比);ρ:反应性(“油门/刹车”);β:总的缓发中子份额(安全缓冲);Λ:瞬发中子寿命(极短,约10⁻⁴秒);λᵢ, cᵢ:第i组缓发中子的衰变常数和先驱核浓度;S:外中子源。
这个方程描述了六组缓发中子先驱核浓度的变化。它们像一个个“延时胶囊”,在中子产生后一段时间才释放出缓发中子,这对控制至关重要。
$$\frac{dc_i}{dt}= \frac{\beta_i}{\Lambda}\,n - \lambda_i c_i \quad (i=1,\ldots,6)$$cᵢ:第i组先驱核浓度;βᵢ:第i组的缓发中子份额;λᵢ:该组先驱核的衰变常数。先驱核由裂变产生(正比于n),然后按各自的寿命(1/λᵢ)衰变,释放出缓发中子。
反应堆控制系统设计与验证:工程师在设计控制棒提棒程序时,会用点堆模型仿真各种操作下的功率响应,确保不会意外引入过大的正反应性导致瞬变过快。比如,模拟控制棒分组提升时,功率是否能平稳上升。
反应性引入事故(RIA)安全分析:分析如控制棒意外弹出等事故场景。通过仿真快速引入大正反应性的后果,评估是否可能触及瞬发临界,从而制定安全限值和保护系统动作逻辑。
反应堆启动与堆期测量:在反应堆启动的低功率阶段,操作员会进行“堆期”测量,即测量功率翻倍所需的时间。点堆模型是解释和预测堆期的基础,通过仿真可以理解不同反应性引入速率下的堆期变化。
控制棒价值刻度与标定:通过有目的地移动某根控制棒并测量引起的功率变化(即反应性变化),再利用点堆模型进行反推,可以精确标定每根控制棒的价值(能引入多少反应性),这是运行规程制定的关键依据。
首先,存在一种“点堆模型过于简化而不实用”的误解。虽然它确实无法给出空间输出分布,但对于快速、本质地理解控制棒操作引起的全堆时间响应而言,它是最佳工具。在实际工作中,我们会在运行更复杂的三维程序之前,先在此检查参数敏感性。其次是参数设置的陷阱。例如,若随意缩短“瞬发中子寿命Λ”(比如从10^-7秒改为10^-5秒),动态特性会急剧加快,导致控制变得极其困难——这是亲身体验快堆与热中子堆根本差异的典型案例。另外,认为“反应性1美元(ρ=β)是安全边界”的固有观念也很危险。本工具将β设为常数,但实际反应堆中核燃料成分会随燃耗变化,β值也会逐渐改变。需注意:模拟中的“1美元”并不总是与实际安全裕度完全一致。