一种堆芯最小流量估算方法技术

本发明专利技术公开了一种堆芯最小流量估算方法。该方法步骤如下：1、输入参数：棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数。2、计算出通道总功率。3、给定初始堆芯总流量，确定通道最优流量分配方案。4、进行组件出口过冷度计算。5、判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内，若否，则进行堆芯总流量修正，返回步骤4；若是，则进行下一个步骤计算。6、堆芯温度分布计算。7、输出修正后的堆芯最小流量及最小流量下的堆芯温度分布。该方法优点如下：1、采用通道最优流量分配方案。2、对堆芯总流量修正构造修正量的关于迭代步数的函数，提高迭代运算的计算速度。高迭代运算的计算速度。高迭代运算的计算速度。

【技术实现步骤摘要】
一种堆芯最小流量估算方法


[0001]本专利技术属于反应堆热工水力
，具体涉及一种堆芯最小流量估算方法。

技术介绍

[0002]采用带盒燃料组件的紧凑型压水堆，应用流量分区设计能够显著展平堆芯出口温度分布，提高出口过冷度，或降低冷却剂流量，从而更加有效地利用冷却剂，大幅提升堆芯设计性能。寿期内堆芯三维功率分布形状会发生较大变化，为了合理分配流量，达到组件出口冷却剂温度展平、过冷度提高，通过研究流量分区依据并设计考虑多种因素的优化方法，设计闭式通道堆芯流量分区及水密度反馈程序QSubTH，并使程序能够满足以下计算功能：
[0003]快速计算闭式通道燃料组件轴向水密度分布、芯块燃料温度分布、包壳温度分布、组件出口温度及过冷度；
[0004]在给定堆芯热功率及三维分布、最小出口过冷度条件下，估算堆芯最小冷却剂质量流量及最大出口温度。

技术实现思路

[0005]为了解决上述问题，本专利技术提供一种堆芯最小流量估算方法。在给定堆芯热功率及三维分布、最小出口过冷度条件下，估算堆芯最小冷却剂质量流量及最大出口温度。
[0006]为了达到上述目的，本专利技术的采用如下技术方案：
[0007]一种堆芯最小流量估算方法，包括如下步骤：
[0008]一种堆芯最小流量估算方法，包括如下步骤：
[0009]步骤1：输入参数：包括棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数，其中堆芯功率分布参数及组件不均匀系数由堆芯物理计算后提供相应数据，包括组件所在堆芯位置、组件功率分布和组件不均匀系数；步骤2：使用堆芯功率分布参数及组件不均匀系数，对每个通道轴向功率求和，计算出通道总功率；
[0010]步骤3：给定初始堆芯总流量，结合堆芯入口温度以及堆芯总功率，由能量守恒方程计算堆芯出口温度一致时的焓值，由各个通道的总功率以及计算出的焓值确定通道最优流量分配方案；
[0011]步骤4：进行组件出口过冷度计算，在计算结果中搜索组件出口最小过冷度；
[0012]步骤5：判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内，若否，则进行堆芯总流量修正，返回步骤4重新进行组件出口过冷度计算；若是，则进行下一个步骤计算；
[0013]步骤6：堆芯温度分布计算：使用棒状或板状燃料元件堆芯参数、堆芯计算控制变量参数，通过能量守恒方程确定各组件中冷却剂温度，再由换热系数关系式以及材料的导热系数，通过傅里叶定律，牛顿冷却公式由外向内依次计算包壳外层、包壳内层、燃料表面、燃料中心温度的轴向分布情况；
[0014]步骤7：输出修正后的堆芯最小总流量及最小总流量下的堆芯温度分布。
[0015]1、本专利技术与现有技术相比，具有如下优点：步骤1中输入参数包括组件所在堆芯位置、组件功率分布和组件不均匀系数，将堆芯寿期内功率变化考虑入内；输入参数中堆芯流量控制变量提高流量优化效率。
[0016]2、步骤3中采用最优流量分配方案，是对于单一功率分布形状，当各个通道出口处的温度相等时，可以得到组件最大出口过冷度。由用户提供的堆芯总流量，入口温度以及堆芯总功率，由能量守恒方程，可以确定出口温度一致时的焓值，则由各个通道的功率以及计算出的出口焓值确定通道最优流量。
[0017]3、步骤5中的允许范围，由用户提供，程序中也提供了自定义接口，方便用户自由控制，建议值为10
‑7；步骤5中的堆芯总流量修正，是对堆芯总流量给出带符号的修正量，在原有修正量上添加不同量级的比例系数，构造修正量的关于迭代步数的函数，提高迭代运算的计算速度。
附图说明
[0018]图1是本专利技术方法的流程图。
具体实施方式
[0019]下面结合附图和具体实施例对本专利技术作进一步详细说明。
[0020]如图1所示，本专利技术一种堆芯最小流量估算方法，包括如下步骤：
[0021]步骤1：输入参数：包括棒状或板状燃料元件堆芯参数(包括堆芯旁流系数、包壳与燃料芯块之间的间隙导热系数、相邻燃料棒中心距、堆芯活性区高度、堆芯最小流量估算流量迭代系数、组件中导向管数目、堆芯组件数目、轴向划分控制体数目、计算区域中组件数目、划分计算区域数目、组件中燃料棒数、燃耗点数、堆芯压力、堆芯总功率、燃料棒包壳内半径、燃料棒包壳内外径、燃料芯块半径、导向管半径、组件出口最小过冷度设计值、堆芯入口冷却剂温度、通道冷却剂流量、方型组件宽度和堆芯冷却剂净流量，若为板状燃料，还需要输入板元件长度、燃料板数目、包壳厚度、板元件厚度和板元件宽度)、堆芯功率分布参数及组件不均匀系数(由堆芯物理计算后提供相应数据，包括组件所在堆芯位置、组件功率分布和组件不均匀系数)、堆芯计算控制变量参数(包括堆芯流量控制变量、组件功率分布、不同流量分配下堆芯出口温度比较计算控制变量、燃料类型控制变量、计算类型控制变量和参考流量输入控制变量(用于不同流量比较计算))。
[0022]步骤2：使用输入参数，对每个通道轴向功率求和，计算出通道总功率。
[0023]步骤3：给定初始堆芯总流量，结合堆芯入口温度以及堆芯总功率，由能量守恒方程计算堆芯出口温度一致时的焓值，由各个通道的堆芯功率以及计算出的焓值确定通道最优流量分配方案。对计算得到的每套相应燃耗点下最优分配方案进行搜索，确定每个组件在整个寿期中所需的最大流量，为了保证堆芯总流量守恒，以确定的各组件最大流量确定初始流量分配比例因子。
[0024]对于单一功率分布形状，当各个通道出口处的温度相等时，可以得到组件最大出口过冷度。由已知的堆芯总流量，入口温度以及堆芯总功率，由能量守恒方程，可以确定出口温度一致时的组件出口焓值：
[0025][0026]上式中：
[0027]h
out
——流量分配最优时，组件出口焓值；
[0028]h
in
——冷却剂在入口处的焓值；
[0029]P
tot
——堆芯总功率；
[0030]W
tot
——堆芯总冷却剂流量。
[0031]则由各个通道的功率以及计算出的组件出口焓值确定通道最优流量：
[0032][0033]上式中：
[0034]W(i)——通道i的对应的最佳分配流量；
[0035]P(i)——通道i的对应的总功率。
[0036]在整个寿期内，堆芯功率分布会随着燃耗加深而有所改变，为了得到整个寿期内的最优流量分配，需要将功率分布的时间效应考虑进去。在程序中，以组件最大出口温度作为衡量流量分配方案好坏的标准。
[0037]对于每个燃耗点下，都有一套最优流量分配方案，每套分配方案相对独立，没有将功率分布的时间效应考虑进去，在确定初始分配方案时，首先对计算得到的每套相应燃耗点下最优分配方案进行搜索，确定每个组件在整个寿期中所需的最大流量，为了保证堆芯总流量守恒，以确定的各组件最大流量确定初始流量分配比例因子：
[0038][0039]上式中：
...

【技术保护点】

【技术特征摘要】
1.一种堆芯最小流量估算方法，其特征在于：包括如下步骤：步骤1：输入参数：包括棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数，其中堆芯功率分布参数及组件不均匀系数由堆芯物理计算后提供相应数据，包括组件所在堆芯位置、组件功率分布和组件不均匀系数；步骤2：使用堆芯功率分布参数及组件不均匀系数，对每个通道轴向功率求和，计算出通道总功率；步骤3：给定初始堆芯总流量，结合堆芯入口温度以及堆芯总功率，由能量守恒方程计算堆芯出口温度一致时的焓值，由各个通道的总功率以及计算出的焓值确定通道最优流量分配方案；步骤4：进行组件出口过冷度计算，在计算结果中搜索组件出口最小过冷度；步骤5：判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内，若否，则进行堆芯总流量修正，返回步骤4重新进行组件出口过冷度计算；若是，则进行下一个步骤计算；步骤6：堆芯温度分布计算：使用棒状或板状燃料元件堆芯参数、堆芯计算控制变量参数，通过能量守恒方程确定各组件中冷却剂温度，再由换热系数关系式以及材料的导热系数，通过傅里叶定律，牛顿冷却公式由外向内依次计算包壳外层、包壳内层、燃料表面、燃料中心温度的轴向分布情况；步骤7：输出修正后的堆芯最小总流量及最小总流量下的堆芯温度分布。2.根据权利要求1所述的一种堆芯最小流量估算方法，其特征在于：步骤1中所述的棒状或板状燃料元件堆芯参数包括堆芯旁流系数、包壳与燃料芯块之间的间隙导热系数、相邻燃料棒中心距、堆芯活性区高度、堆...

【专利技术属性】
技术研发人员：张魁，蔡佳仪，田文喜，苏光辉，秋穗正，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：