[发明专利]一种动态提升拟合阶数的搜索堆芯临界棒位计算方法

摘要

本发明公开了一种动态提升拟合阶数的搜索堆芯临界棒位计算方法，包括：1)在第三步判断堆芯在当前棒位是否达到临界状态时，记录了每一次迭代过程中的堆芯有效增殖系数和控制棒组棒位；2)第四步采用最小二乘法拟合得到堆芯有效增殖系数增量与第i组控制棒高度增量的N阶非线性多项式；3)第五步采用牛顿迭代法求解第四步中的N阶非线性多项式，得到第i组控制棒进入第N+1迭代时的预估临界棒位，本方法需要控制棒组初始棒位、调棒顺序和对应的控制棒组调节起止位置，能够计算多种反应堆堆型临界状态的控制棒棒位，能够随着调棒临界搜索计算迭代次数的增加动态提升预估多项式的阶数，提高调棒临界搜索计算临界棒位的效率和精度。

主权项

1.一种动态提升拟合阶数的搜索堆芯临界棒位计算方法，其特征在于，所述方法包括：第一步，判断控制棒位于初始棒位时堆芯是否达到临界；第二步：若第一步堆芯未达到临界状态，则进入迭代搜索堆芯临界棒位的过程；第三步：计算控制棒组在当前棒位处的截面参数，采用堆芯物理程序计算堆芯有效增殖系数keff^N，若则堆芯达到临界，此时的棒位即为临界棒位；反之堆芯未达到临界状态，并记录和keff^N，其中，N为堆芯调棒临界搜索计算的迭代次数，i为控制棒组的编码，为第i组控制棒在第N次临界搜索时距离堆芯底部的高度，Error^N为第N个临界搜索迭代步时的相对偏差，堆芯临界有效增值系数为Keff_critical，堆芯临界准则为eps；第四步，若第三步中堆芯未达到临界，根据记录堆芯临界搜索第0步至第N步的堆芯有效增殖系数和第i组控制棒的棒位，获得有效增殖系数增量ΔKeff与第i组控制棒高度增量ΔH_i的非线性多项式估计函数ΔKeff＝f(ΔH_i)，当调棒临界搜索计算的迭代次数每增加一次，非线性多项式估计函数的阶数则动态增加一阶；其中，f(ΔH_i)为第i组控制棒高度增量ΔH_i的多项式，第n次临界搜索的有效增殖系数增量ΔKeffⁿ＝Keffⁿ‑Keff⁰,n∈[1,N]，第n次临界搜索第i组控制棒高度增量和Keff⁰为堆芯初始棒位时的高度和对应的堆芯有效增殖系数，和Keffⁿ分别为第n次临界搜索时第i组控制棒的高度和堆芯有效增殖系数；第五步，假设第N+1次临界搜索时堆芯达到临界，则求解该非线性多项式估计函数，得到第i组控制棒的预估临界棒位为第N+1次临界搜索时第i组控制棒的棒位；第六步：重复第二步至第五步，直到第三步中堆芯达到临界状态；第七步：进入下一个燃耗步，重复第一步至第六步，直到最后一个燃耗步或控制棒全部提出堆芯。