结合阻力和能量分布包含格架搅混效应的子通道分析方法

摘要

结合阻力和能量分布包含格架搅混效应的子通道分析方法，通过收集包含待分析反应堆堆芯搅混格架在内的多种搅混格架在不同工况下的实验数据，拟合能够反映搅混格架关键部件搅混性能的动量源项关系式和能量源项关系式，并将该动量源项关系式添加至对应的动量守恒方程中，将该能量源项关系式添加至对应的能量守恒方程中；通过求解质量守恒方程、动量守恒方程和能量守恒方程，获得反应堆堆芯内更准确的局部热工水力参数，从而使得临界热流密度(CHF)值和CHF位置的预测更加准确。本发明不受子通道程序结构和解法的限制，可以广泛应用在各种类型的子通道程序计算中，主要应用于两流体模型的子通道程序。

说明书

技术领域

本发明属于反应堆燃料组件流体动力学模拟和分析技术领域，具体涉及一种结合阻力和能量分布包含格架搅混效应的子通道分析方法。

背景技术

在核电站的设计和运行中，临界热流密度(CHF)、临界功率是限制核电站运行的最重要的参数之一。在临界功率以下，反应堆可以最大限度安全地将核能转换成电能。燃料棒束CHF和堆芯的局部热工水力参数密切相关，鉴于反应堆堆芯子通道复杂的几何结构以及宽参数范围的运行工况，并不能在理论上对CHF完全给出解析解。

传统核燃料棒束的CHF关系式是基于棒束实验数据和子通道程序计算出的局部热工水力参数开发得到。而为了获得这些反应堆工况下的实验数据需耗费相当大的财力和较长的时间周期。基于此开发的CHF关系式用来进行核电站运行安全分析并指导核燃料设计。子通道分析方法是目前安审必备的工具，用计算反应堆堆芯的局部热工水力参数并藉以做出安全分析和评估其CHF限值，这是反应堆安全分析的关键步骤。目前除了因为集总参数平均效应产生的不准确性以外，子通道分析主要的弊端之一是缺少对绝大多数反应堆堆芯中所使用格架的搅混性能的模拟和计算。因此，目前所有子通道程序在计算局部热工水力参数时普遍存在很大的不确定性和误差。这种不确定性和误差，在CHF关系式开发拟定时，经常需借助大量棒束CHF实验数据来弥补，并导致其拟定的CHF关系式应用范围的局限性。同时，在安全分析时，这种子通道程序在计算局部热工水力参数时的不确定性也经常导致CHF计算中很大的误差，并因而产生很大的设计限值及设计偏差，使得反应堆运行效率偏低或安全系数降低。

由于对搅混格架功能缺乏认知，目前所有的子通道程序均将搅混格架的模拟简化为普通的形阻系数或总体性的所谓“搅混系数”(β或TDC)，完全无法对格架上各个关键部件(如搅混翼、弹簧、刚凸和焊点等)的不同搅混性能进行模拟。此外，有人提出了在子通道分析方法中，采用基于雷诺数开发的形阻系数关系式来代替原来的形阻系数常量。而采用上述方法，因为所用的形阻系数关系式与格架几何结构仍然无关，依旧不能精确地反映格架的搅混特性，无法消除局部热工水力参数计算的不准确性。

发明内容

本发明目的是提供一种结合阻力和能量分布包含格架搅混效应的子通道分析方法，该方法通过拟合能够反映搅混格架关键部件搅混性能的动量源项关系式和能量源项关系式，并将该动量源项关系式添加至对应的动量守恒方程中，将该能量源项关系式添加至对应的能量守恒方程中，求解新的守恒方程，从而获得反应堆堆芯内更准确的局部热工水力参数，使得CHF值和CHF位置的预测更加准确。

本发明的技术解决方案是：

一种结合阻力和能量分布包含格架搅混效应的子通道分析方法，包括以下步骤：

1)预设待分析的反应堆堆芯的CHF值和CHF位置的预测值的准确度；

2)收集搅混格架在不同工况下的实验数据；

所述搅混格架为待分析的反应堆堆芯的搅混格架，或包含待分析的反应堆堆芯搅混格架在内的多种搅混格架；

所述实验数据包括至少四类：第一类和第二类实验数据为棒束轴向和径向非均匀加热工况下的CHF实验获得的CHF值和CHF位置，其余实验数据是在棒束通道压力分布、棒束通道流速分布、搅混实验获得的出口温度分布以及棒束通道流场温度中任选至少两类；

3)采取以下任一方法建立待分析的反应堆堆芯的搅混格架的动量源项关系式和能量源项关系式：

方法A)利用待分析的反应堆堆芯的搅混格架在不同工况下的至少四类实验数据直接拟合动量源项关系式和能量源项关系式；

方法B)利用待分析的反应堆堆芯的搅混格架在不同工况下的至少四类实验数据标定计算流体动力学软件，再用该软件在标定的范围内计算局部热工水力参数，拟合动量源项关系式和能量源项关系式；

方法C)利用所收集的多种搅混格架在不同工况下的至少四类实验数据拟合与所收集的多种搅混格架相匹配的动量源项通用关系式和能量源项通用关系式，再由待分析的反应堆堆芯的搅混格架在不同工况下的实验数据求得动量源项通用关系式和能量源项通用关系式的系数，最终获得反映该搅混格架关键部件搅混性能的动量源项关系式和能量源项关系式；所述关键部件包括搅混翼、条带、刚凸、弹簧和焊点；

4)在子通道程序中添加步骤3)获得的动量源项关系式和能量源项关系式，再用该子通道程序计算待分析的反应堆堆芯的热工水力参数，从中得到该反应堆堆芯不同工况下CHF值和CHF位置的预测值；

5)分别判断该反应堆堆芯不同工况下CHF值的预测值与步骤2)中所收集的CHF值之间的比值，以及不同工况下CHF位置的预测值与步骤2)中所收集的CHF位置之间的差值是否满足步骤1)中所预设的准确度；若任意一项不满足，则对子通道程序进行优化和/或对CHF关系式进行改进，再重复步骤4)和5)；

6)输出待分析反应堆堆芯的CHF值和CHF位置的预测值。

上述步骤3)的方法C)中拟合动量源项关系式和能量源项关系式的方法具体如下：

C1)将收集的同一类实验数据进行对比，找出与格架搅混性能相关的几何结构参数和局部热工水力参数；

C2)选择搅混格架的多个几何结构参数和局部热工水力参数分别作为动量源项通用关系式和能量源项通用关系式的变量；

C3)结合收集的实验数据和选择的变量，采用最小二乘法拟合动量源项通用关系式和能量源项通用关系式；

C4)针对待分析的反应堆堆芯搅混格架的实验数据确定动量源项通用关系式和能量源项通用关系式中的系数，最终得到动量源项关系式和能量源项关系式。

上述步骤C3)中所拟合的动量源项通用关系式包括轴向动量源项通用关系式Φ_A和横向动量源项通用关系式Φ_L，具体如下：

式中：

A″_W是包含格架时的润湿周长，单位为m；

ρ是流体密度，单位为kg/m³；

V_total是流场局部速度，单位为m/s；

A_R是不包含格架时的润湿周长，单位为m；

θ是流场局部速度与棒束方向的夹角，单位为度；它是流场的局部轴向速度和局部横向速度的函数，其计算式可表达如下：

D″_V是包含格架时的子通道当量直径，单位为m；

P是子通道的湿周，单位为m；

u是棒束轴向分速度，单位为m/s；

v是棒束横向分速度，单位为m/s；

u₁是搅混翼切向分速度，单位为m/s；

v₁是搅混翼法向分速度，单位为m/s；

α是搅混翼与棒束方向的夹角，单位为度；

K_u，K_v，K_u1和K_v1是动量源项通用关系式中对应待分析的反应堆堆芯的搅混格架的系数；K_u，K_v，K_u1和K_v1的值跟该搅混格架的几何结构密切相关，由该搅混格架的实验数据来确定；

所述步骤C3)中所拟合的能量源项通用关系式跟气液两相的界面面积、液相比焓、气相比焓、液相速度、气相速度、压力密切相关，具体如下：

式中：

下标l表示液相，下标g表示气相；

ρ_l和ρ_g分别是液相密度和气相密度，单位为kg/m³；

h_l和h_g分别是液相比焓和气相比焓，单位为kJ/kg；

V_l和V_g分别是液相速度和气相速度，单位为m/s；

L是搅混翼长度，单位为m；

E₁，E₂和E₃是能量源项通用关系式中对应待分析的反应堆堆芯的搅混格架的系数；E₁，E₂和E₃的值跟该搅混格架的几何结构和流场的流型相关，包括泡状流、弹状流、弥散流和环状流，在不同的流型工况下，E₁，E₂和E₃的值不同，需要根据不同流型工况下的实验数据来确定。

上述步骤4)具体如下：

4.1)在子通道程序中添加步骤3)获得的动量源项关系式和能量源项关系式；同时，确定待分析的反应堆堆芯的几何条件和计算所要满足的残差要求；所述几何条件包括燃料棒个数、燃料棒直径、燃料棒中心间距、棒束长度和棒束流通面积、格架间距；所述残差要求包括求解质量守恒方程、动量守恒方程或能量守恒方程时所要满足的残差要求；

4.2)在添加了动量源项关系式和能量源项关系式的子通道程序中，将反应堆堆芯的流场按照其轴向功率分布、横向功率分布以及几何条件划分成以子通道为流通截面的多个轴向控制体；所述子通道如图2所示，是由相邻的燃料棒表面和相邻燃料棒的中心连线围成的最小流通通道A，或者由相邻的燃料棒表面、壁面和相邻燃料棒的中心连线的延伸线围成的最小流通通道B或C；

4.3)求解该反应堆堆芯的燃料模型获得热边界条件，或根据该反应堆堆芯的热流密度获得热边界条件；

4.4)求解能量守恒方程、动量守恒方程和质量守恒方程，得到满足残差要求的所有轴向控制体内的热工水力参数；其中动量守恒方程的求解需要将该搅混格架的动量源项关系式添加至对应的流体气相、液相和/或液滴相的动量守恒方程中，能量守恒方程的求解需要将该搅混格架的能量源项关系式添加至对应的流体气相、液相和/或液滴相的能量守恒方程中；

4.5)从求解得到的所有轴向控制体内的热工水力参数中提取局部热工水力参数，从中得到该反应堆堆芯不同工况下CHF值和CHF位置的预测值。

上述步骤4.4)第一种解法具体如下：

4.4.1)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的能量守恒方程，并将该反应堆堆芯中的搅混格架的能量源项关系式添加至对应的流体气相、液相和/或液滴相的能量守恒方程中；再根据热边界条件求解该能量守恒方程，获得每个轴向控制体的比焓、燃料棒表面温度、流体密度和流体温度；

4.4.2)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的动量守恒方程，并将该反应堆堆芯中的搅混格架的动量源项关系式添加至对应的流体气相、液相和/或液滴相的动量守恒方程中；再根据每个轴向控制体的比焓、燃料棒表面温度、流体密度和流体温度，求解该动量守恒方程，获得横向流速和压力；

4.4.3)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的质量守恒方程；根据每个轴向控制体内的横向流速和压力，求解该质量守恒方程，获得轴向流速；

4.4.4)判断质量守恒方程、动量守恒方程和能量守恒方程是否满足残差要求，若不满足，则按照步骤4.4.1)至步骤4.4.3)的顺序进行反复迭代，若满足，则获得每个轴向控制体内的热工水力参数；所述热工水力参数包括燃料棒表面温度、轴向流速、横向流速、压力、流体密度、流体温度、比焓。

上述步骤4.4)第二种解法具体如下：

4.4.1)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的动量守恒方程，并将该反应堆堆芯中的搅混格架的动量源项关系式添加至对应的流体气相、液相和/或液滴相的动量守恒方程中；再根据热边界条件求解该动量守恒方程，获得初始的横向流速和轴向流速；

4.4.2)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的质量守恒方程，根据每个轴向控制体的横向流速和轴向流速，求解该质量守恒方程，获得更新后的横向流速、轴向流速和压力；

4.4.3)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的能量守恒方程，并将该反应堆堆芯中的搅混格架的能量源项关系式添加至对应的流体气相、液相和/或液滴相的能量守恒方程中；再根据每个轴向控制体内更新后的横向流速、轴向流速和压力，求解该能量守恒方程，获得比焓；

4.4.4)判断动量守恒方程、质量守恒方程和能量守恒方程是否满足残差要求，若不满足，则按照步骤4.4.1)至步骤4.4.3)的顺序进行反复迭代，若满足，则获得每个轴向控制体内的热工水力参数；所述热工水力参数包括轴向流速、横向流速、压力、流体密度、流体温度、比焓。

上述步骤4)中添加了动量源项关系式的子通道程序的动量守恒方程具体如下：

轴向动量守恒方程：

式中：

下标i和j是子通道编号，其中i表示所要求解的子通道，j表示与子通道i相邻的子通道；

m_i是子通道i内的轴向流量，单位为kg/s；

t是时间，单位为s；

Z是单个轴向控制体的轴向高度，单位为m；

ρ_i是子通道i内的流体密度，单位为kg/m³；

A_i是子通道i的流通面积，单位为m²；

w_ij是从子通道i到子通道j的横向流量，单位为kg/(m*s)；w_ij为正值时表示横向流量从子通道i流到子通道j，w_ij为负值时表示横向流量从子通道j流到子通道i；

u^*是施主单元控制体内的轴向流速，单位为m/s；其意义表达如下式：

f_T是横向湍流搅混摩擦系数，单位为1；

w′_ij是从子通道i到子通道j的湍流搅混量，单位为kg/(m*s)；

w′_ji是从子通道j到子通道i的湍流搅混量，单位为kg/(m*s)；

u′_i是子通道i的脉动速度，单位为m/s；

u′_j是子通道j的脉动速度，单位为m/s；

p_i是子通道i内的压力，单位为Pa；

g是重力加速度，单位为m/s²；

f是轴向摩擦系数，单位为1；

D_h是子通道的当量直径，单位为m；

Φ_A是反映该搅混格架的搅混性能的轴向动量源项关系式；

横向动量守恒方程：

式中：

l是横向动量控制体的等效长度，单位为m；

v^*是施主单元控制体内的横向流速，单位为m/s；其意义表达如下式：

K_G是横向流动阻力系数，单位为1；

S_ij是子通道i与子通道j之间的间隙宽度，单位为m；

p_i和p_j分别是子通道i与子通道j内的压力，单位为Pa；

Φ_L是反映该搅混格架的搅混性能的横向动量源项关系式；

上述步骤4)中添加了能量源项关系式的子通道程序的能量守恒方程具体如下：

式中：

h_i和h_j分别是子通道i和子通道j内的比焓，单位为kJ/kg；

h^*是施主单元控制体内的比焓，单位为kJ/kg；其意义表达如下式：

λ_i是子通道i内的流体导热系数，单位为W/(m*K)；

T_i和T_j分别是子通道i和子通道j内的流体温度，单位为℃；

L_C是几何修正常数，单位为1；

n是与子通道i相邻的加热壁面个数，单位为1；

是与子通道i相邻的加热壁面面积，单位为m²；

q_n是与子通道i相邻的加热壁面的热流密度，单位为J/(m²*s)；

Θ是反映该搅混格架搅混性能的能量源项关系式；当能量守恒方程代表气相的能量变化规律时，Θ前面的符号为正，当能量守恒方程代表液相的能量变化规律时，Θ前面的符号为负。

本发明的优点：

1、本发明通过引入代表搅混格架影响的参数，包括搅混翼长度、搅混翼面积、搅混翼角度以及代表了条带、刚凸、弹簧和焊点对流场影响的系数，从而更准确地反映了搅混格架的几何特征以及其对流场带来的影响，提高了子通道方法预测反应堆堆芯局部热工水力参数的准确性；

2、本发明通过引入代表搅混格架几何结构的参数，包括搅混翼长度、搅混翼面积、搅混翼角度以及代表了条带、刚凸、弹簧和焊点对流场影响的系数，并且收集多个搅混格架的实验数据进行动量源项通用关系式和能量源项通用关系式的开发，因此能广泛应用于不同类型搅混格架对流场影响的分析，能够反映不同类型搅混格架对流场的影响，提高了子通道方法的适用性；

3、本发明计算所得到的包含格架搅混效应的局部热工水力参数，比传统方法计算所得到的局部热工水力参数更加准确，而局部热工水力参数是CHF关系式开发的关键因素之一，因此本发明计算所得到的包含格架搅混效应的局部热工水力参数可用于更准确的反映格架效应的CHF关系式的开发；

4、本发明计算所得到的包含格架搅混效应的局部热工水力参数，比传统方法计算所得到的局部热工水力参数更加准确，而局部热工水力参数是CHF值和位置预测的关键因素之一，因此本发明计算所得到的包含格架搅混效应的局部热工水力参数可用于更准确地预测CHF值和CHF位置；

5、本发明能够反映特定类型搅混格架对流场的影响，可提供更可靠的安全分析，提高安全系数，增加经济效益；

6、本发明能够反映不同类型搅混格架对流场的影响，可作为格架设计开发的工具。

附图说明

图1是权利要求3中所述轴向动量源项通用关系式和横向动量源项通用关系式中的部分变量示意图。

图2是本发明中代表核电厂典型反应堆堆芯的小范围(5x5)测试用子通道示意图。

图3和图4是阻力分布(DRM)包含格架搅混效应的子通道分析的两种流程图。

图5是5x5棒束子通道划分和编号示意图。

图6是子通道(2，8)之间的横向流速曲线图。

图7是子通道(9，10)之间的横向流速曲线图。

具体实施方式

本发明是结合阻力和能量分布包含格架搅混效应的子通道分析方法，包括以下步骤：

1)预设待分析的反应堆堆芯的CHF值和CHF位置的预测值的准确度；

所述CHF值的预测值的准确度可用下式表示：

式中CHF_m是实验获得的CHF值，CHF_p是计算得到的CHF值的预测值，ε是预设的CHF值的预测值的准确度；

所述CHF位置的预测值的准确度可用下式表示：

|L_m-L_p|＜δ>

式中L_m是实验获得的CHF位置，L_p是计算得到的CHF位置的预测值，δ是预设的CHF位置的预测值的准确度；

2)收集搅混格架在不同工况下的实验数据；

所述搅混格架为待分析的反应堆堆芯的搅混格架，或包含待分析的反应堆堆芯的搅混格架在内的多种搅混格架；

所述实验数据包括至少四类：第一类和第二类实验数据为棒束轴向和径向非均匀加热工况下的CHF实验获得的CHF值和CHF位置，其余实验数据在棒束通道压力分布、棒束通道流速分布、搅混实验获得的出口温度分布以及棒束通道流场温度中任选至少两类；要求收集的每种搅混格架的实验数据至少包含200点以上的实验数据。

3)采取以下任一方法建立待分析的反应堆堆芯的搅混格架的动量源项关系式和能量源项关系式：

方法A)利用待分析的反应堆堆芯的搅混格架在不同工况下的至少四类实验数据直接拟合动量源项关系式和能量源项关系式；

方法B)利用待分析的反应堆堆芯的搅混格架在不同工况下的至少四类实验数据标定计算流体动力学软件，再用该软件在标定范围内计算局部热工水力参数，拟合动量源项关系式和能量源项关系式；标定范围是指实验数据对应工况所包含的压力、温度和流速范围，计算流体动力学软件在标定范围内的计算才是准确可靠的；

方法C)利用所收集的多种搅混格架在不同工况下的至少四类实验数据拟合与所收集的多种搅混格架相匹配的动量源项通用关系式和能量源项通用关系式，再由待分析的反应堆堆芯的搅混格架在不同工况下的实验数据求得动量源项通用关系式和能量源项通用关系式的系数，最终获得反映该搅混格架关键部件搅混性能的动量源项关系式和能量源项关系式；所述关键部件包括搅混翼、条带、刚凸、弹簧和焊点；

方法C)中拟合动量源项关系式和能量源项关系式的方法具体如下：

C1)将收集的同一类实验数据进行对比，根据物理现象找出与格架搅混性能相关的几何结构参数和局部热工水力参数，对找出的参数进行敏感性分析，或根据经验列出重要度排序表；所述与搅混性能相关的格架几何结构参数包括但不限于搅混翼面积、搅混翼弯角、搅混翼扭角、搅混翼长度、搅混翼宽度；与搅混性能相关的局部热工水力参数包括但不限于轴向流速、横向流速、温度、压力、CHF值、CHF位置；

C2)依据敏感性分析结果或者重要度排序表，选择搅混格架的多个几何结构参数和局部热工水力参数作为动量源项通用关系式和能量源项通用关系式的变量；动量源项通用关系式和能量源项通用关系式的变量不完全相同；所述变量个数取决于所收集的实验数据的多少，所收集的实验数据越多，变量个数越多；变量个数越多，动量源项通用关系式和能量源项通用关系式越能够更精确地反映实际的搅混性能；

C3)结合收集的实验数据和选择的变量，采用最小二乘法拟合动量源项通用关系式和能量源项通用关系式；

所拟合的动量源项通用关系式包括轴向动量源项通用关系式Φ_A和横向动量源项通用关系式Φ_L，它们都跟选择的变量密切相关，具体如下：

式中：

A″_W是包含格架时的润湿周长，单位为m；

ρ是流体密度，单位为kg/m³；

V_total是流场局部速度，单位为m/s；

A_R是不包含格架时的润湿周长，单位为m；

θ是流场局部速度与棒束方向的夹角，单位为度；它是流场的局部轴向速度和局部横向速度的函数，其计算式可表达如下：

D″_V是包含格架时的子通道当量直径，单位为m；

P是子通道的湿周，单位为m；

u是棒束轴向分速度，单位为m/s；

v是棒束横向分速度，单位为m/s；

u₁是搅混翼切向分速度，单位为m/s；

v₁是搅混翼法向分速度，单位为m/s；

α是搅混翼与棒束方向的夹角，单位为度；

式中K_u，K_v，K_u1和K_v1是动量源项通用关系式中对应待分析的反应堆堆芯的搅混格架的系数。上述轴向动量源项通用关系式和横向动量源项通用关系式适用于气相、液相和/或液滴相，对于具体的各相，式中各项的表达式可作进一步细化。

结合收集的实验数据和选择的变量，采用跟拟合动量源项通用关系式相同的方法拟合能量源项通用关系式；所拟合的能量源项通用关系式跟气液两相的界面面积、液相比焓、气相比焓、液相速度、气相速度、压力密切相关，具体如下：

式中：

下标l表示液相，下标g表示气相；

ρ_l和ρ_g分别是液相密度和气相密度，单位为kg/m³；

h_l和h_g分别是液相比焓和气相比焓，单位为kJ/kg；

V_l和V_g分别是液相速度和气相速度，单位为m/s；

L是搅混翼长度，单位为m；

E₁，E₂和E₃是能量源项通用关系式中对应待分析的反应堆堆芯的搅混格架的系数。

E₁，E₂和E₃的值跟待分析的反应堆堆芯的搅混格架的几何结构和流场的流型相关，包括泡状流、弹状流、弥散流和环状流。

C4)K_u，K_v，K_u1和K_v1需要结合该搅混格架的实验数据来确定，最终获得该搅混格架的动量源项关系式。K_u，K_v，K_u1和K_v1的值跟搅混格架的几何结构密切相关，包括搅混翼形状和大小，弹簧、刚凸、条带和焊点的大小和排布方式。

E₁，E₂和E₃需要结合该搅混格架的实验数据来确定，最终获得该搅混格架的能量源项关系式。E₁，E₂和E₃的值除了跟几何结构相关外，还跟流型有关，包括泡状流、弹状流、弥散流和环状流，在不同的流型工况下，E₁，E₂和E₃的值不同，需要根据不同流型工况下的实验数据来确定。

动量源项通用关系式和能量源项通用关系式中系数的确定需要根据该搅混格架的实验数据确定其压力范围、流场温度范围、流场速度范围和加热功率范围，并利用最小二乘法确定各个待定系数。确定后的系数只适用于该搅混格架，且只适用于该搅混格架的实验数据所确定的参数范围内。

4)在子通道程序中添加步骤3)获得的动量源项关系式和能量源项关系式，再用该子通道程序计算待分析的反应堆堆芯的热工水力参数，从中得到该反应堆堆芯不同工况下CHF值和CHF位置的预测值；

具体步骤为：

4.1)在子通道程序中添加步骤3)获得的动量源项关系式和能量源项关系式；同时，确定待分析的反应堆堆芯的几何条件和计算所要满足的残差要求；所述几何条件包括燃料棒个数、燃料棒直径、燃料棒中心间距、棒束长度和棒束流通面积、格架间距；所述残差要求包括求解质量守恒方程、动量守恒方程或能量守恒方程时所要满足的残差要求；

4.2)在添加了动量源项关系式和能量源项关系式的子通道程序中，将反应堆堆芯的流场按照其轴向功率分布、横向功率分布以及几何条件划分成以子通道为流通截面的多个轴向控制体；所述子通道如附图2所示，是由相邻的燃料棒表面和相邻燃料棒的中心连线围成的最小流通通道A，或者由相邻的燃料棒表面、壁面和相邻燃料棒的中心连线的延伸线围成的最小流通通道B或C；

4.3)求解该反应堆堆芯的燃料模型获得热边界条件，或根据该反应堆堆芯的热流密度获得热边界条件；反应堆堆芯的燃料模型是工业上应用比较成熟的模块，不是本发明的内容，因此不在此详述；热流密度是由用户根据实际情况自定义的参数；

4.4)求解能量守恒方程、动量守恒方程和质量守恒方程，得到满足残差要求的所有轴向控制体内的热工水力参数；其中动量守恒方程的求解需要将该搅混格架的动量源项关系式添加至对应的流体气相、液相和/或液滴相的动量守恒方程中，能量守恒方程的求解需要将该搅混格架的能量源项关系式添加至对应的流体气相、液相和/或液滴相的能量守恒方程中；

4.5)从求解得到的所有轴向控制体内的热工水力参数中提取局部热工水力参数，从中得到该反应堆堆芯不同工况下CHF值和CHF位置的预测值。所有轴向控制体内的热工水力参数中包含热流密度的预测值，热边界条件中包含所有轴向控制体内加热的热流密度；当某个轴向控制体内热流密度的预测值与对应的该轴向控制体内加热的热流密度的比值为1时，对应的该轴向控制体的热流密度的预测值就是CHF值的预测值，对应的该轴向控制体的位置就是CHF位置的预测值。

步骤4.4)具体有多种解法；虽然不同的子通道分析的具体解法各不相同，但是本发明根据实验数据拟合了动量源项通用关系式和能量源项通用关系式，并且根据该搅混格架的实验数据，给不同形状的搅混格架的动量源项通用关系式和能量源项通用关系式赋予不同的系数，因此可以广泛应用于各种类型的子通道分析中；本发明具体给出了两种解法。

第一种解法：

4.4.1)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的能量守恒方程，并将反应堆堆芯中该搅混格架的能量源项关系式添加至该能量守恒方程中，示例如下：

式中：

i和j是子通道编号，其中i表示所要求解的子通道，j表示与子通道i相邻的子通道；

A_i是子通道i的流通面积，单位为m²；

ρ_i是子通道i内的流体密度，单位为kg/m³；

h_i和h_j分别是子通道i和子通道j内的比焓，单位为kJ/kg；

t是时间，单位为s；

m_i是子通道i内的轴向流量，单位为kg/s；

Z是单个轴向控制体的轴向高度，单位为m；

w_ij是从子通道i到子通道j的横向流量，单位为kg/(m*s)；w_ij为正值时表示横向流量从子通道i流到子通道j，w_ij为负值时表示横向流量从子通道j流到子通道i；

h^*是施主单元控制体内的比焓，单位为kJ/kg；其意义表达如下式：

w′_ij是从子通道i到子通道j的湍流搅混量，单位为kg/(m*s)；

w′_ji是从子通道j到子通道i的湍流搅混量，单位为kg/(m*s)；

λ_i是子通道i内的流体导热系数，单位为W/(m*K)；

T_i和T_j分别是子通道i和子通道j内的流体温度，单位为℃；

S_ij是子通道i与子通道j之间的间隙宽度，单位为m；

L_C是几何修正常数，单位为1；

n是与子通道i相邻的加热壁面个数，单位为1；

是与子通道i相邻的加热壁面面积，单位为m²；

q_n是与子通道i相邻的加热壁面的热流密度，单位为J/(m²*s)；

Θ是反映该搅混格架搅混性能的能量源项关系式；当能量守恒方程代表气相的能量变化规律时，Θ前面的符号为正，当能量守恒方程代表液相的能量变化规律时，Θ前面的符号为负。

式中左边第一项是每单位轴向长度上子通道i内能量随时间的变化率；第二项是单位轴向长度上子通道i内能量的空间变化；第三项是单位轴向长度上子通道i内全部连接间隙上的横向能量之和，即单位轴向长度上从所有相邻子通道j流至子通道i的横向能量之和；第四项是单位轴向长度上由于湍流搅混引起的子通道i与所有相邻子通道j之间的能量交换；式中右边表示单位轴向长度上子通道i内的总能量变化。

根据热边界条件求解流体气相、液相和/或液滴相的能量守恒方程，获得每个轴向控制体的比焓、燃料棒表面温度、流体密度和流体温度；

4.4.2)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的动量守恒方程，并将反应堆堆芯中该搅混格架的动量源项关系式添加至该动量守恒方程中；因为动量是标量，有方向性，因此子通道分析方法中一般将动量分解为轴向动量和横向动量两部分，对应的动量守恒方程包括轴向动量守恒方程和横向动量守恒方程，对应的动量源项关系式也包括轴向动量源项关系式和横向动量源项关系式；所述搅混格架的动量源项关系式是由步骤3)中任一方法建立的；所述轴向动量守恒方程和横向动量守恒方程分别示例如下：

轴向动量守恒方程：

式中：

u^*是施主单元控制体内的轴向流速，单位为m/s；其意义表达如下式：

f_T是横向湍流搅混摩擦系数，单位为1；

u′_i和u′_j分别是子通道i和子通道道j的脉动速度，单位为m/s；

p_i是子通道i内的压力，单位为Pa；

g是重力加速度，单位为m/s²；

f是轴向摩擦系数，单位为1；

D_h是子通道的当量直径，单位为m；

Φ_A是反映该搅混格架的搅混性能的轴向动量源项关系式。

式中左边第一项是每单位轴向长度上轴向动量随时间的变化，第二项是单位轴向长度上轴向动量的空间变化，第三项是单位轴向长度上全部连接间隙上对轴向动量的贡献，第四项是湍流搅混对轴向动量的贡献；右边第一项至第三项分别是轴向压力梯度、重力和摩擦阻力对轴向动量的影响，最后一项是反映该搅混格架的搅混性能的轴向动量源项对轴向动量的影响。上述轴向动量守恒方程适用于气相、液相和/或液滴相，对于具体的各相，式中各项的表达式可作进一步细化。

横向动量守恒方程：

式中：

l是横向动量控制体的等效长度，单位为m；

v^*是施主单元控制体内的横向流速，单位为m/s；其意义表达如下式：

K_G是横向流动阻力系数，单位为1；

Φ_L是反映该搅混格架的搅混性能的横向动量源项关系式。

式中左边第一项是每单位轴向长度上横向动量随时间的变化，第二项是单位轴向长度上横向动量的空间变化，第三项是单位轴向长度上全部连接间隙上对横向动量的贡献；右边第一项和第二项分别是横向阻力对横向阻力和横向压力梯度对横向动量的影响，第三项是反映该搅混格架的搅混性能的横向动量源项对横向动量的影响。上述横向动量守恒方程适用于气相、液相和/或液滴相，对于具体的各相，式中各项的表达式可作进一步细化。

根据每个轴向控制体的比焓、燃料棒表面温度、流体密度和流体温度，同时联立求解该轴向动量守恒方程和横向动量守恒方程，获得横向流速和压力；

4.4.3)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的质量守恒方程，示例如下：

式中第一项是每单位轴向长度上质量随时间的变化，第二项是单位轴向长度上轴向质量流量的空间变化，第三项是单位轴向长度上湍流搅混引起的质量交换。上述能量守恒方程适用于气相、液相和/或液滴相，对于具体的各相，式中各项的表达式可作进一步细化。

根据每个轴向控制体内的横向流速和压力，求解该质量守恒方程，获得轴向流速；

4.4.4)判断质量守恒方程、动量守恒方程和能量守恒方程是否满足残差要求，若不满足，则按照步骤4.4.1)至步骤4.4.3)的顺序进行反复迭代，若满足，则获得每个轴向控制体内的热工水力参数；所述热工水力参数包括燃料棒表面温度、轴向流速、横向流速、压力、流体密度、流体温度、比焓。

第二种解法：

4.4.1)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的动量守恒方程，并将该反应堆堆芯中的搅混格架的动量源项关系式添加至对应的流体气相、液相和/或液滴相的动量守恒方程中；因为动量是标量，有方向性，因此子通道分析方法中一般将动量分解为轴向动量和横向动量两部分，对应的动量守恒方程包括轴向动量守恒方程和横向动量守恒方程，对应的动量源项关系式也包括轴向动量源项关系式和横向动量源项关系式；所述搅混格架的动量源项关系式是由步骤3)中任一方法建立的；所述轴向动量守恒方程和横向动量守恒方程分别示例如下：

轴向动量守恒方程：

横向动量守恒方程：

根据热边界条件求解轴向动量守恒方程和横向动量守恒方程，获得初始的横向流速和轴向流速；

4.4.2)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的质量守恒方程，示例如下：

根据每个轴向控制体的横向流速和轴向流速，求解该质量守恒方程，获得更新后的横向流速、轴向流速和压力；

4.4.3)以每个轴向控制体为单元，列出对应的流体气相、液相和/或液滴相的能量守恒方程，并将反应堆堆芯中该搅混格架的能量源项关系式添加至该能量守恒方程中，示例如下：

根据每个轴向控制体内更新后的横向流速、轴向流速和压力，求解该能量守恒方程，获得比焓；

4.4.4)判断动量守恒方程、质量守恒方程和能量守恒方程是否满足残差要求，若不满足，则按照步骤4.4.1)至步骤4.4.3)的顺序进行反复迭代，若满足，则获得每个轴向控制体内的热工水力参数；所述热工水力参数包括燃料棒表面温度、轴向流速、横向流速、压力、流体密度、流体温度、比焓。

5)分别判断该反应堆堆芯不同工况下CHF值的预测值与步骤2)中所收集的CHF值之间的比值，以及不同工况下CHF位置的预测值与步骤2)中所收集的CHF位置之间的差值是否满足步骤1)中所预设的准确度；若任意一项不满足，则对子通道程序进行优化和/或对CHF关系式进行改进，再重复步骤4)和5)；

6)输出待分析反应堆堆芯的CHF值和CHF位置的预测值。

设计思路：

本发明提供一种结合阻力和能量分布包含格架搅混效应的子通道分析方法。鉴于目前的子通道程序未能真正反映搅混格架对流场产生的影响，本发明将搅混格架对流场的影响通过动量源项关系式和能量源项关系式的形式分别添加至对应的动量守恒方程和能量守恒方程中，通过求解新的守恒方程，从而获得搅混翼对子通道内压降、横流和温度分布的影响效应，获得更准确的局部热工水力参数。

本发明的技术路线图如图3和图4所示。首先收集包含待分析反应堆堆芯搅混格架在内的多种搅混格架在不同工况下的实验数据，根据收集到的实验数据拟合与所收集的多种搅混格架相匹配的动量源项通用关系式和能量源项通用关系式，然后对待分析的反应堆堆芯的搅混格架的实验数据，获得动量源项通用关系式和能量源项通用关系式中的系数，最终获得该搅混格架的动量源项关系式和能量源项关系式，再列出质量守恒方程、动量守恒方程和能量守恒方程，将动量源项关系式添加至对应的动量守恒方程中，将能量源项关系式添加至对应的能量守恒方程中，按一定的顺序迭代求解质量、动量和能量守恒方程直到满足残差要求，从而得到更准确的局部热工水力参数，用于更准确的CHF预测或者CHF关系式的开发。

由于本发明通过添加动量源项关系式和能量源项关系式来反映搅混格架带来的影响，不改变程序的求解流程，因此本发明不受子通道程序结构和解法的限制，可以广泛应用在各种类型的子通道程序计算中，主要应用于两流体模型的子通道程序。

图6和图7是带搅混格架的5x5棒束子通道间隙的横向流速的计算结果。其中黑色曲线表示由标定后的流体动力学软件计算的结果，作为对比分析的基准值，蓝色曲线表示由原子通道程序计算的结果，红色曲线表示由本发明的方法改进后的子通道程序计算的结果。5x5棒束的子通道划分和编号如图5所示，图6是子通道(2，8)之间的间隙的横流值，图7是子通道(9，10)之间的间隙的横流值。由图可知，采用本发明的方法改进后的子通道程序显著改善了局部热工水力参数的预测能力，可用于更准确的CHF的预测或者CHF关系式的开发。