一种反应堆小破口多维度耦合分析方法技术

本发明专利技术公开了一种反应堆小破口多维度耦合分析方法，步骤如下：1、建立带有小破口的压力壳流体域几何模型，并进行节点划分；2、建立安全壳几何模型，并进行网格划分；3、设置边界条件以及初始条件，并计算T0时刻压力壳流体域参数；4、将T

【技术实现步骤摘要】
一种反应堆小破口多维度耦合分析方法


[0001]本专利技术属于方法专利技术
，具体涉及一种反应堆小破口多维度耦合分析方法。

技术介绍

[0002]反应堆小破口事故指由于反应堆中冷却剂管道(部件)出现破口，导致冷却剂丧失速率大于补给速率时的冷却剂丧失事故。事故严重程度与破口位置、破口尺寸等影响因素有关。反应堆冷却系统由于冷却剂减少会引起的后果，包括压力下降、堆芯冷却恶化、冷却剂泄漏到安全壳和潜在的放射性向工厂外的泄漏。因此开展小破口事故分析将具有重要意义。
[0003]对此事故当前通常采用系统程序和安全壳分析程序开展独立分析，此方法可以获得相对保守的计算结果。但对于先紧压水堆而言，其采用大量非能动装置，事故过程中需耦合压力壳与安全壳内部过程，以了解事故下二者的相互影响。
[0004]采用一维系统程序与三维计算流体力学程序耦合方法对反应堆小破口事故开展耦合分析，将有助于对反应堆小破口事故下压力壳内的三维热工水力特性开展分析。开展实时耦合数据更新迭代将更好地反应真实物理过程。本文这提出的分析方法可以为反应堆小破口事故多维度耦合分析提供方法参考。

技术实现思路

[0005]为了解决上述问题，本专利技术提供了一种反应堆小破口多维度耦合分析方法，该方法将一维系统分析程序(经过二次开发，修改子程序)、三维流体力学分析程序(经过二次开发，用户自定义文件)通过数据接口程序耦合起来，完整而准确地模拟反应堆小破口事故下安全壳与压力壳耦合热工水力过程。
[0006]为达到上述目的，本专利技术采用了如下技术方案：
[0007]一种反应堆小破口多维度耦合分析方法，包括以下步骤：
[0008]步骤1：在反应堆热工水力系统分析程序中根据压力壳几何结构参数建立几何模型，得到带有小破口的压力壳高压水
‑
空气流体域几何模型，并将小破口位置确定为数据交互的耦合面；具体地，对压力壳高压水
‑
空气流体域几何模型进行节点划分，并设置边界条件和初始条件，基于简化目的，在计算中压力壳高压水
‑
空气流体域几何模型被处理为边界存在小破口的圆柱体，并将破口位置的控制体以及接管设置为耦合边界；
[0009]步骤2：获得安全壳的几何结构参数，在三维几何建模软件中对安全壳空气
‑
水流体域进行建模，得到安全壳空气
‑
水流体域几何模型；对安全壳空气
‑
水流体域几何模型进行网格划分，并设置边界条件和初始条件，具体步骤如下：
[0010]步骤2
‑
1：运用三维几何建模软件根据安全壳几何结构参数建立安全壳空气
‑
水流体域几何模型，基于简化目的，在计算中安全壳空气
‑
水流体域被处理为内部存在小破口的圆柱体；
[0011]步骤2
‑
2：在步骤2
‑
1获得的安全壳空气
‑
水流体域几何模型的基础上进行网格划分，得到安全壳空气
‑
水流体域网格模型；
[0012]步骤3：对反应堆热工水力系统分析程序以及计算流体动力学计算程序进行边界条件以及初始条件的设置，并根据T0时刻的初始条件进行计算，得到压力壳高压水
‑
空气流体域中小破口处出口流体温度、流速、组分浓度以及物性参数；
[0013]步骤4：将T
n
时刻反应堆热工水力系统分析程序计算得到的压力壳高压水
‑
空气流体域小破口处出口流体温度、流速、组分浓度以及物性参数作为安全壳空气
‑
水流体域计算时的输入参数，作为计算流体动力学计算程序的入口条件，计算得到T
n+1
时刻安全壳空气
‑
水流体域内部热工水力参数分布，包括温度、压力、组分浓度以及物性参数，其中n＝0、1、2
……
，具体步骤如下：
[0014]步骤4
‑
1：将T
n
时刻反应堆热工水力系统分析程序计算得到的压力壳高压水
‑
空气流体域小破口处出口流体温度、流速、组分浓度以及物性参数作为安全壳空气
‑
水流体域输入参数，作为计算流体力学软件的入口条件；
[0015]步骤4
‑
2：将安全壳空气
‑
水流体域网格模型的侧壁面以及上下底面设置为壁面边界，依据真实事故条件将其设置为恒定壁面温度边界；
[0016]步骤4
‑
3：通过求解混合物的质量、动量、能量及组分输运方程，得到安全壳内部空气
‑
水流体域的混合物速度场、温度场、压力场以及组分浓度场，求解组分浓度场的组分输运方程为：
[0017][0018]方程(1)为组分输运方程的具体形式，其中：
[0019]ρ——混合物的密度，kg/m3；
[0020]Y
i
——混合物中的第i相质量份额，％；
[0021]——混合物的速度，m/s；
[0022]t——时间,s；
[0023]D
i,m
——混合物中第i相的层流扩散系数m2·
s
‑1；
[0024]D
i,t
——混合物中第i相的湍流扩散系数m2·
s
‑1；
[0025]步骤4
‑
4：以步骤4
‑
3中获得的混合物速度场、温度场、压力场以及组分浓度场为基础，得到高压水由于压力突降造成的蒸发质量，以及蒸汽受到低温壁面冷却的冷凝质量，并进一步更新得到的安全壳内部空气
‑
水流体域的混合物速度场、温度场、压力场以及组分浓度场；具体地，蒸发质量和冷凝质量计算如下：
[0026][0027][0028]当流体温度大于饱和温度时，有：
[0029][0030]当流体温度小于饱和温度时，有：
[0031][0032]式中：
[0033]α
v
——蒸汽份额；
[0034]ρ
v
——蒸汽的密度，kg/m3；
[0035]——蒸汽速度，m/s；
[0036]——蒸汽速度，m/s；
[0037]α
l
——液态水份额；
[0038]ρ
v
——水的密度，kg/m3；
[0039]——蒸发质量，kg；
[0040]——冷凝质量，kg；
[0041]T
sat
——饱和温度，K；
[0042]T
l
——液态水温度，K；
[0043]T
g
——蒸汽温度，K；
[0044]coeff——蒸发冷凝系数；
[0045]步骤5：将T
n+1
时刻安全壳空气
‑
水流体域内部温度、压力作为压力壳高压水
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种反应堆小破口多维度耦合分析方法，其特征在于：针对反应堆小破口中的热工水力现象，由不同的反应堆热工水力分析程序包括计算流体动力学计算程序、反应堆热工水力系统分析程序通过压力壳与安全壳之间实时双向的数据传递，实现耦合计算；该方法包括以下步骤：步骤1：在反应堆热工水力系统分析程序中根据压力壳几何结构参数建立几何模型，得到带有小破口的压力壳高压水
‑
空气流体域几何模型，并将小破口位置确定为数据交互的耦合面；具体地，对压力壳高压水
‑
空气流体域几何模型进行节点划分，并设置边界条件和初始条件，基于简化目的，在计算中压力壳高压水
‑
空气流体域几何模型被处理为边界存在小破口的圆柱体，并将破口位置的控制体以及接管设置为耦合边界；步骤2：获得安全壳的几何结构参数，在三维几何建模软件中对安全壳空气
‑
水流体域进行建模，得到安全壳空气
‑
水流体域几何模型；对安全壳空气
‑
水流体域几何模型进行网格划分，并设置边界条件和初始条件，具体步骤如下：步骤2
‑
1：运用三维几何建模软件根据安全壳几何结构参数建立安全壳空气
‑
水流体域几何模型，基于简化目的，在计算中安全壳空气
‑
水流体域被处理为内部存在小破口的圆柱体；步骤2
‑
2：在步骤2
‑
1获得的安全壳空气
‑
水流体域几何模型的基础上进行网格划分，得到安全壳空气
‑
水流体域网格模型；步骤3：对反应堆热工水力系统分析程序以及计算流体动力学计算程序进行边界条件以及初始条件的设置，并根据T0时刻的初始条件进行计算，得到压力壳高压水
‑
空气流体域中小破口处出口流体温度、流速、组分浓度以及物性参数；步骤4：将T
n
时刻反应堆热工水力系统分析程序计算得到的压力壳高压水
‑
空气流体域小破口处出口流体温度、流速、组分浓度以及物性参数作为安全壳空气
‑
水流体域计算时的输入参数，作为计算流体动力学计算程序的入口条件，计算得到T
n+1
时刻安全壳空气
‑
水流体域内部热工水力参数分布，包括温度、压力、组分浓度以及物性参数，其中n＝0、1、2
……
，具体步骤如下：步骤4
‑
1：将T
n
时刻反应堆热工水力系统分析程序计算得到的压力壳高压水
‑
空气流体域小破口...

【专利技术属性】
技术研发人员：王成龙，杨宇鹏，张大林，秋穗正，田文喜，苏光辉，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：