一种高温气冷堆螺旋管蒸汽发生器热工水力特性三维计算方法与流程

1.本发明属于核反应堆热工水力计算技术领域，具体涉及到一种高温气冷堆螺旋管蒸汽发生器热工水力特性三维计算方法。


背景技术：

2.蒸汽发生器是核反应堆动力系统中至关重要的组成部分，是反应堆一二回路热量传递的枢纽。螺旋管蒸汽发生器是一种结构紧凑、换热效率高的蒸汽发生器，其主要由：螺旋管束、内筒、壳体外壁面组成。螺旋管束绕套筒按照不同的径向位置分层布置，每层内均匀布置多根螺旋管。螺旋管蒸汽发生器因其独特结构，较常规换热器增加了两侧流体的换热面积与湍流强度，换热能力得到增强，得益于其独特优势，螺旋管蒸汽发生器被广泛应用于高温气冷堆、船用压水堆等。
3.高温气冷堆螺旋管蒸汽发生器结构特殊，管内流体除了受到重力、摩擦力的影响，还受到离心力的影响；管外流体为外掠管束流动，受到管束的搅浑作用，管外流体的湍流程度大大提升，强化了两侧流体的换热，提高了蒸汽发生器的换热效率。为了高温气冷堆结构设计与进一步优化、运行控制以及安全性能分析，有必要对高温气冷堆螺旋管蒸汽发生器换热性能进行分析。在分析方法上，有实验模拟与计算流体动力学计算两种方法，实验模拟成本较高，使用计算流体动力学手段进行计算是有效且方便的。已有的方法多关注单侧流体流动换热，无法真实模拟螺旋管蒸汽发生器内部流动过程，或是开发系统程序进行模拟，但是缺乏对细节现象的描述能力。使用计算流体动力学手段对高温气冷堆螺旋管蒸汽发生器两侧流体耦合换热进行研究，可以对全流域进行精细模拟，同时当前还没有针对高温气冷堆螺旋管蒸汽发生器换热特性进行三维耦合的研究，因此该方法将对高温气冷堆螺旋管蒸汽发生器的进一步发展具有重要意义。


技术实现要素：

4.本发明的目的在于提供一种高温气冷堆螺旋管蒸汽发生器热工水力特性三维计算方法，该方法能够运用计算流体动力学手段对高温气冷堆螺旋管蒸汽发生器管内的蒸汽蒸发以及螺旋管束外的高温气体流动进行耦合计算。
5.为了达到上述目的，本发明采用如下技术方案：
6.本发明方法采用多孔介质方法以及计算流体动力学软件中的组分输运模型以及质量扩散理论来模拟高温气冷堆螺旋管蒸汽发生器的蒸汽蒸发现象；通过编写udf导入该软件的方法模拟螺旋管蒸汽发生器外部的高温气体流动；同时根据傅里叶导热定律计算螺旋管管壁的温度分布。
7.一种高温气冷堆螺旋管蒸汽发生器热工水力特性三维计算方法，包括如下步骤：
8.步骤1：建立高温气冷堆螺旋管蒸汽发生器计算域简化模型，包括螺旋管蒸汽发生器外部封壳模型和螺旋管蒸汽发生器中心支撑柱模型，具体步骤如下：
31.——汽相相对于混合物质量加权平均速度的相对速度，m/s
32.s
e
——能量源项，w
·
m
‑333.通过组分输运方程，得到高温气冷堆螺旋管蒸汽发生器管侧流体组分浓度场分布，其中组分输运方程为：
[0034][0035]
其中：
[0036]
t——时间，s
[0037]
ρ
m
——混合物密度，kg/m3[0038]
——混合物速度，m/s
[0039]
y
i
——混合物中的第i相组分浓度即蒸汽和水的质量份额，％
[0040]
d
i,m
——混合物中第i相的层流扩散系数，m2·
s
‑1[0041]
d
i,t
——混合物中第i相的湍流扩散系数，m2·
s
‑1[0042]
步骤3
‑
4：求解二次侧能量守恒方程，获得高温气冷堆螺旋管蒸汽发生器二次侧流体温度场分布，其中二次侧能量守恒方程为:
[0043][0044]
其中：
[0045]
t——时间，s
[0046]
β——孔隙率
[0047]
α
g
——汽相体积份额
[0048]
ρ
g
——汽相密度，kg/m3[0049]
h
g
——汽相焓，j
·
kg
‑1[0050]
α
l
——液相体积份额
[0051]
ρ
l
——液相密度，kg/m3[0052]
h
l
——液相焓，j
·
kg
‑1[0053]
——汽相速度，m/s
[0054]
p——压力，pa
[0055]
——液相速度，m/s
[0056]
k
m
——混合物导热系数，w
·
m
‑1·
k
‑1[0057]
s
g
——汽相质量源项，kg
·
m
‑3·
s
‑1[0058]
步骤3
‑
5：在求得的二次侧流体速度场、温度场与组分浓度场的基础上，通过得到的蒸汽质量份额w
s
能够计算出螺旋管蒸汽发生器管侧流体域与螺旋管束管壁固体域交界面处的蒸汽蒸发率为：
[0059][0060]
其中：
[0061]
——蒸汽蒸发率，kg/s
[0062]
ρ
m
——混合物密度，kg/m3[0063]
d
m
——层流扩散系数，m2·
s
‑1[0064]
w
s
——蒸汽质量份额，％
[0065]
n——面法向量
[0066]
步骤3
‑
6：由步骤3
‑
5中得到的蒸汽蒸发率计算蒸汽蒸发吸收螺旋管管壁结构的固体域热量，其中蒸汽蒸发吸收螺旋管壁热量的吸热功率为：
[0067][0068]
其中：
[0069]
——蒸汽蒸发率，kg/s
[0070]
h
fg
——蒸汽汽化潜热，j/kg
[0071]
q
e
——蒸汽吸热功率，w
[0072]
步骤4：计算螺旋管管壁结构固体域内的温度分布，并得到螺旋管管壁与螺旋管蒸汽发生器壳侧气体流体域交界面处的网格分布，具体的，在螺旋管管壁固体域内，热量主要通过热传导的方式传递，根据傅里叶导热定律：
[0073][0074]
其中：
[0075]
q
e
——蒸汽吸热功率，w
[0076]
λ——螺旋管管壁结构固体域的导热系数，w/m2·
k
[0077]
a——螺旋管蒸汽发生器管侧流体域和螺旋管管壁结构的固体域的交界面面积，m2[0078]
t——螺旋管管壁结构的固体域温度，k
[0079]
x——螺旋管管壁结构的固体域由内指向外的位置坐标，m
[0080]
由(7)式得到螺旋管管壁结构的固体域与螺旋管蒸汽发生器壳侧气体流体域交界面处的温度分布；
[0081]
步骤5：以步骤4中得到的螺旋管管壁结构的固体域与螺旋管蒸汽发生器壳侧气体流体域交界面处的温度分布，采用高温气体的计算方法计算螺旋管蒸汽发生器壳侧气体流体域内的速度场与温度场，具体步骤如下：
[0082]
步骤5
‑
1：将螺旋管蒸汽发生器壳侧高温气体流体域下部面设置为高温气体的出口，上部面设置为高温气体进口；
[0083]
步骤5
‑
2：求解高温气体的质量、动量和能量守恒方程，得到螺旋管蒸汽发生器壳侧气体流体域内的速度场与温度场分布；考虑到高温气体的特殊物理性质，采用真实气体模拟，采用式(8)作为气体物性关系式:
[0084]
pv＝rt
g
b(t)p
ꢀꢀꢀꢀ
(8)
[0085]
其中：
[0086]
p——气体压力，pa；
[0087]
v——气体比热容，j
·
kg
‑1·
k
‑1；
[0088]
r——气体常数；
[0089]
t
g
——气体温度，k；
[0090]
b(t)——维里系数；
[0091]
并将其写成脚本形式，导入计算流体动力学软件；对高温气体的计算需对计算流体动力学软件进行修正，具体的，选取适配的湍流普朗特关系式，并将其写成脚本形式并导入计算流体动力学软件，实现对湍流模型的修正；
[0092]
步骤5
‑
3：按照实际情况，对不同区域内的组分浓度进行初始化；同时对入口边界进行组分设置，保证流场与实际情况符合。
[0093]
本发明具有以下有益效果：
[0094]
1)能够实现高温气冷堆螺旋管蒸汽发生器管侧蒸发过程与壳侧高温气体冷却过程的耦合计算；
[0095]
2)模型独立，方法通用性强，可以适用于不同的计算流体动力学计算分析程序；
[0096]
3)该计算模型能够实现高温气体、液态水、蒸汽的同时模拟，可以给高温气冷堆中的螺旋管蒸汽发生器工程设计提供更准确地计算数据。
附图说明
[0097]
图1为典型高温气冷堆螺旋管式蒸汽发生器三维管束模型。
[0098]
图2为典型高温气冷堆螺旋管式蒸汽发生器流场示意图。
[0099]
图3为本发明计算方法流程图。
具体实施方式
[0100]
以下结合图3所示流程图，以典型高温气冷堆螺旋管蒸汽发生器计算过程为例，对本发明作进一步的详细描述，典型高温气冷堆螺旋管蒸汽发生器管束结构如图1所示，其壳侧、管侧流体流动方向如图2所示。
[0101]
本发明高温气冷堆螺旋管蒸汽发生器热工水力三维计算方法，包括如下步骤：
[0102]
步骤1：建立高温气冷堆螺旋管蒸汽发生器计算域简化模型，包括螺旋管蒸汽发生器外部封壳模型，螺旋管蒸汽发生器中心支撑柱模型，具体步骤如下：
[0103]
步骤1
‑
1：运用几何模型建立软件solidworks建立高温气冷堆螺旋管蒸汽发生器外部封壳模型，模型为圆柱型结构；
[0104]
步骤1
‑
2：运用几何模型建立软件solidworks建立高温气冷堆螺旋管蒸汽发生器内部支撑柱模型，模型为与步骤1
‑
1建立的高温气冷堆螺旋管蒸汽发生器外部封壳模型同心等高的圆柱体；
[0105]
步骤1
‑
3：将步骤1
‑
1与步骤1
‑
2建立的模型进行装配，得到高温气冷堆螺旋管蒸汽发生器简化环形通道模型；
[0106]
步骤2：在步骤1得到的高温气冷堆螺旋管蒸汽发生器简化环形通道模型的基础上
运用ansys
‑
icem进行网格划分，具体过程如下：利用“define_property”设置管束区多孔介质模型，定义二次侧多孔介质孔隙率；在高温气冷堆螺旋管蒸汽发生器其他部位内沿轴向、周向、径向划分不同层数的四面体结构化网格，通过网格独立性分析判定网格数量。
[0107]
步骤3：在高温气冷堆螺旋管蒸汽发生器管侧计算域网格模型范围内进行流动换热计算，具体步骤如下：
[0108]
步骤3
‑
1：将高温气冷堆螺旋管蒸汽发生器管侧计算域下部表面设置为流体域入口，并根据真实条件设置入口处混合物中气相体积份额、混合物速度、混合物密度、混合物温度和压力；
[0109]
步骤3
‑
2：将高温气冷堆螺旋管蒸汽发生器管侧计算域上部表面设置成流体域出口，并设置出口处的压力和温度；
[0110]
步骤3
‑
3：求解二次侧阻力源项，耦合求解二次侧质量、动量守恒方程，获得高温气冷堆螺旋管蒸汽发生器二次侧流体速度场分布；其中二次侧质量守恒方程为:
[0111][0112]
二次侧动量守恒方程为：
[0113][0114]
其中：
[0115]
t——时间，s
[0116]
β——孔隙率
[0117]
ρ
m
——混合物密度，kg/m3[0118]
——混合物速度，m/s
[0119]
μ
m,eff
——混合物有效粘度，pa
·
s
[0120]
——重力加速度，m
·
s
‑2[0121]
α
g
——汽相体积份额
[0122]
ρ
g
——汽相密度，kg/m3[0123]
α
l
——液相体积份额
[0124]
ρ
l
——液相密度，kg/m3[0125]
——汽相相对于混合物质量加权平均速度的相对速度，m/s
[0126]
s
e
——能量源项，w
·
m
‑3[0127]
通过组分输运方程，得到高温气冷堆螺旋管蒸汽发生器管侧流体组分浓度场分布，其中组分输运方程为：
[0128][0129]
其中：
[0130]
t——时间，s
[0131]
ρ
m
——混合物密度，kg/m3[0132]
——混合物速度，m/s
[0133]
y
i
——混合物中的第i相组分浓度即蒸汽和水的质量份额，％
[0134]
d
i,m
——混合物中第i相的层流扩散系数，m2·
s
‑1[0135]
d
i,t
——混合物中第i相的湍流扩散系数，m2·
s
‑1[0136]
步骤3
‑
4：求解二次侧能量守恒方程，获得高温气冷堆螺旋管蒸汽发生器二次侧流体温度场分布，其中二次侧能量守恒方程为:
[0137][0138]
其中：
[0139]
t——时间，s
[0140]
β——孔隙率
[0141]
α
g
——汽相体积份额
[0142]
ρ
g
——汽相密度，kg/m3[0143]
h
g
——汽相焓，j
·
kg
‑1[0144]
α
l
——液相体积份额
[0145]
ρ
l
——液相密度，kg/m3[0146]
h
l
——液相焓，j
·
kg
‑1[0147]
——汽相速度，m/s
[0148]
p——压力，pa
[0149]
——液相速度，m/s
[0150]
k
m
——混合物导热系数，w
·
m
‑1·
k
‑1[0151]
s
g
——汽相质量源项，kg
·
m
‑3·
s
‑1[0152]
步骤3
‑
5：在求得的二次侧流体速度场、温度场与组分浓度场的基础上，通过得到的蒸汽质量份额w
s
能够计算出螺旋管蒸汽发生器管侧流体域与螺旋管束管壁固体域交界面处的蒸汽蒸发率为：
[0153][0154]
其中：
[0155]
——蒸汽蒸发率，kg/s
[0156]
ρ
m
——混合物密度，kg/m3[0157]
d
m
——层流扩散系数，m2·
s
‑1[0158]
w
s
——蒸汽质量份额，％
[0159]
n——面法向量
[0160]
步骤3
‑
6：由步骤3
‑
5中得到的蒸汽蒸发率计算蒸汽蒸发吸收螺旋管管壁结构的固体域热量，其中蒸汽蒸发吸收螺旋管壁热量的吸热功率为：
[0161][0162]
其中：
[0163]
——蒸汽蒸发率，kg/s
[0164]
h
fg
——蒸汽汽化潜热，j/kg
[0165]
q
e
——蒸汽吸热功率，w
[0166]
步骤4：计算螺旋管管壁结构固体域内的温度分布，并得到螺旋管管壁与螺旋管蒸汽发生器壳侧气体流体域交界面处的网格分布，具体的，在螺旋管管壁固体域内，热量主要通过热传导的方式传递，根据傅里叶导热定律：
[0167][0168]
其中：
[0169]
q
e
——蒸汽吸热功率，w
[0170]
λ——螺旋管管壁结构固体域的导热系数，w/m2·
k
[0171]
a——螺旋管蒸汽发生器管侧流体域和螺旋管管壁结构的固体域的交界面面积，m2[0172]
t——螺旋管管壁结构的固体域温度，k
[0173]
x——螺旋管管壁结构的固体域由内指向外的位置坐标，m
[0174]
由(7)式得到螺旋管管壁结构的固体域与螺旋管蒸汽发生器壳侧气体流体域交界面处的温度分布；
[0175]
步骤5：以步骤4中得到的螺旋管管壁结构的固体域与螺旋管蒸汽发生器壳侧气体流体域交界面处的温度分布，采用高温气体的计算方法计算螺旋管蒸汽发生器壳侧气体流体域内的速度场与温度场，具体步骤如下：
[0176]
步骤5
‑
1：将螺旋管蒸汽发生器壳侧高温气体流体域下部面设置为高温气体的出口，上部面设置为高温气体进口；
[0177]
步骤5
‑
2：求解高温气体的质量、动量和能量守恒方程，得到螺旋管蒸汽发生器壳侧气体流体域内的速度场与温度场分布；考虑到高温气体的特殊物理性质，采用真实气体模拟，采用式(8)作为气体物性关系式:
[0178]
pv＝rt
g
b(t)p
ꢀꢀꢀꢀ
(8)
[0179]
其中：
[0180]
p——气体压力，pa；
[0181]
v——气体比热容，j
·
kg
‑1·
k
‑1；
[0182]
r——气体常数；
[0183]
t
g
——气体温度，k；
[0184]
b(t)——维里系数；
[0185]
利用“define_property”宏将此关系式撰写成udf文件格式，最后导入fluent软件；对高温气体的计算需对计算流体动力学软件进行修正，具体的，选取适配的湍流普朗特
关系式，并利用“define_property”宏将此关系式撰写成udf文件格式，最后导入fluent软件，实现对湍流模型的修正；
[0186]
步骤5
‑
3：按照实际情况，对不同区域内的组分浓度进行初始化。将螺旋管蒸汽发生器壳侧高温气体流域初始化为氦相，螺旋管蒸汽发生器管侧两相流体域初始化为单相水；同时对于入口边界进行组分设置，将螺旋管蒸汽发生器壳侧高温气体流体域入口处氦气质量分数设置为1，将螺旋管蒸汽发生器管侧两相流体域入口处液态水质量分数设置为1，保证流场与实际情况符合。