CN112231990B - 一种核动力系统一二回路耦合计算方法 - Google Patents
一种核动力系统一二回路耦合计算方法 Download PDFInfo
- Publication number
- CN112231990B CN112231990B CN202011035522.7A CN202011035522A CN112231990B CN 112231990 B CN112231990 B CN 112231990B CN 202011035522 A CN202011035522 A CN 202011035522A CN 112231990 B CN112231990 B CN 112231990B
- Authority
- CN
- China
- Prior art keywords
- steam generator
- loop
- steam
- model
- inlet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种核动力系统一二回路耦合计算方法,步骤如下:1、建立反应堆蒸汽发生器的三维几何模型并划分网格,建立核反应堆一二回路一维简化模型;2、使用计算流体动力学软件设置蒸汽发生器几何模型的初始边界条件和计算域初始状态;3、计算流体动力学软件对蒸汽发生器模型进行计算,得到计算结果,将计算结果传递给系统分析程序,作为系统分析程序此时间步长的边界条件;4、系统分析程序对一二回路模型进行计算,得到计算结果,将计算结果传递给计算流体动力学软件,作为计算流体动力学软件下一时间步长计算时的边界条件;5、重复步骤3、步骤4,直至反应堆系统模型完成全时段的计算。完成对反应堆系统多维度、全系统高保真的耦合分析。
Description
技术领域
本发明属于核反应堆系统技术领域,具体涉及到一种核动力系统一二回路耦合计算方法。
背景技术
蒸汽发生器是核反应堆系统中连接一二回路的重要设备,具有在正常运行和某些事故工况下将一回路热量传递给二回路的重要功能,是一二回路耦合计算实现的重点。对蒸汽发生器的分析方法有使用系统分析程序和计算流体动力学程序两种。系统分析程序一般对一回路进行整体一维建模,蒸汽发生器使用集总参数模型或分布参数模型,根据蒸汽流量、给水流量等因素设置与二回路相连部分的边界条件,其计算结果能够反映出在稳态和事故工况下蒸汽发生器的流动传热特性,但对蒸汽发生器内部自由体积较大的空间,如入口腔室、出口腔室、蒸气室等结构,不能得出比较好的模拟结果。计算流体动力学软件受到计算条件的限制一般只对蒸汽发生器进行整体三维建模,使用雷诺平均NS模型、大涡模拟、分离涡等湍流求解模型,边界条件使用系统分析程序的拟合结果,被动接受一二回路工况变化,无法将计算结果传递回系统分析程序。
综上所述,耦合程序结合了系统分析程序和计算流体动力学软件的优点,能够得出蒸汽发生器在一二回路复杂工况下的详细热工水力信息,同时一二回路的系统分析程序接受蒸汽发生器进出口的准确信息,能够得出系统行为的最佳估算,计算结果与实验值符合得更好。故本发明针对蒸汽发生器在核反应堆一二回路中的热工水力分析具有重要的参考价值和指导意义。
发明内容
本发明的目的在于提供一种核动力系统一二回路耦合计算方法,该方法结合了系统分析程序计算效率高和计算流体动力学软件计算准确的特点,能够在同时考虑一二回路影响因素下对蒸汽发生器进行耦合计算,完成对反应堆系统多维度、全系统高保真的耦合分析。
为达到上述目的,本发明通过以下技术方案实现:
一种核动力系统一二回路耦合计算方法,在使用计算流体动力学软件对蒸汽发生器的流动换热特性分析时,耦合程序使用区域重叠法进行建模,其中计算流体动力学软件通过与系统分析程序在蒸汽发生器与一二回路的两个进口与两个出口处进行数据传递,实现跨纬度耦合,能够考虑一二回路系统工况变化对蒸汽发生器计算结果的影响;包括以下步骤:
步骤1:建立蒸汽发生器的几何模型,包含入口腔室、U型管、出口腔室、给水管、上升段、汽水分离器、蒸汽室和下降段,得到几何模型后进行网格划分;建立一二回路系统的节点图,并按节点图填写系统分析程序的输入卡片,一回路包括压力容器、热腿、蒸汽发生器、稳压器、蒸汽发生器一次侧部分、阀门、冷却水泵和冷腿,二回路包括蒸汽发生器二次侧部分、汽水分离再热器、汽轮机、凝汽器、凝结水泵、低压给水加热器、除氧器、给水泵和高压给水加热器;
步骤2:根据设计工况使用计算流体动力学软件选择蒸汽发生器求解计算时的湍流模型、传热模型、多相流模型和组分输运模型,设置T0时刻(即初始时刻)的边界条件和计算域的初始状态;
步骤3:计算流体动力学软件根据Tn(n=0、1、2、3……)时刻的边界条件对蒸汽发生器模型进行一个时间步长的计算,将计算结果传递给系统分析程序,作为系统分析程序Tn时刻的边界条件,包括分别将蒸汽发生器入口腔室入口处和出口腔室出口处的温度、压力、冷却剂流速和硼浓度传递给一回路系统热腿出口处和冷腿入口处,分别将蒸汽发生器给水管入口处和蒸汽室蒸汽出口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给二回路系统给水管道出口处和主蒸汽管道入口处;
步骤4:系统级分析程序根据Tn时刻的边界条件对一二回路模型进行一个时间步长的计算,将计算结果传递给计算流体动力学软件,作为计算流体动力学软件Tn+1时刻的边界条件,包括分别将一回路系统热腿出口处和冷腿入口处的温度、压力、冷却剂流速和硼浓度传递给蒸汽发生器入口腔室入口处和出口腔室出口处,分别将二回路系统给水管道出口处和主蒸汽管道入口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给蒸汽发生器给水管入口处和蒸汽室蒸汽出口处;
步骤5:重复步骤3、步骤4,直至核动力系统的一二回路模型和蒸汽发生器模型完成全时段的计算,完成蒸汽发生器在核反应堆系统一二回路的跨维度耦合。
本发明提出的耦合分析方法适用于同时考虑一二回路工况对蒸汽发生器的影响,与现有直接对蒸汽发生器的分析方法相比,本发明具有以下有益效果:
1.本发明使用系统分析程序同时对反应堆一二回路进行建模计算,在保障计算效率的同时为蒸汽发生器两处入口、两处出口提供更加精准的边界条件,有利于准确分析一二回路工况变化对蒸汽发生器的影响;
2.本发明使用区域重叠法的区域划分方法进行系统模型建立,系统分析程序拥有完整回路,在瞬态工况的计算过程中收敛性表现较好;
3.本发明使用计算流体动力学软件对蒸汽发生器进行蒸汽三维建模,能够重点分析U型管内冷却剂和二回路给水的流动换热特性,提高反应堆安全分析的准确性;
4.本发明中数据传递发生在每一个时间步长计算结束后,能够在每一时间步长内均进行耦合,重新设置边界条件,保障计算精度;
5.本发明在维持不同维度程序完整性的前提下完成耦合,适用于现有的大多数热工水力分析程序,便于耦合的实现与计算模型与模拟条件的修改。
附图说明
图1为本发明耦合计算方法的一种具体耦合过程示意图。
图2为本发明耦合计算方法对一二回路系统和蒸汽发生器建模示意图。
具体实施方式
下面以计算流体动力学软件通过用户自定义函数功能调用系统分析程序实现耦合分析为例,选取显式耦合方法与区域重叠的区域划分方法,结合图1的耦合过程示意图与图2的建模示意图进行进一步描述:
本发明提供了一种核动力系统一二回路耦合计算方法,图1给出耦合过程的示意图,图中实线框中为系统程序计算域,虚线框中为计算流体动力学软件的计算域,耦合具体方法如下:
步骤1:使用几何模型建立软件对蒸汽发生器的内部构件进行三维建模,忽略细小结构,只对主要结构进行建模分析,如图2中虚线框中内容所示,主要结构包含入口腔室、U型管、出口腔室、给水管、上升段、汽水分离器、蒸汽室和下降段,得到几何模型后使用网格划分工具以六面体网格进行划分,生成非结构网格;对一二回路系统进行一维简化,建立节点图,并按节点图与实验工况填写输入卡片,如图2所示,一回路包括压力容器、热腿、蒸汽发生器、稳压器、蒸汽发生器一次侧、阀门、冷却水泵、冷腿等设备,二回路包括蒸汽发生器二次侧、汽水分离再热器、汽轮机、凝汽器、凝结水泵、低压给水加热器、除氧器、给水泵、高压给水加热器等设备;
步骤2:启动计算流体动力学软件,根据设计工况选择蒸汽发生器求解计算时的湍流模型、传热模型、多相流模型和组分输运模型,设置T0时刻(即初始时刻)的边界条件和计算域的初始状态,编译用户自定义函数实现数据传递与在每一时间步长后调用系统分析程序的功能;
步骤3:计算流体动力学软件根据Tn(n=0、1、2、3……)时刻的边界条件对蒸汽发生器模型进行一个时间步长的计算,将计算结果传递给系统分析程序,作为系统分析程序Tn时刻的边界条件,包括分别将蒸汽发生器入口腔室入口处和出口腔室出口处的温度、压力、冷却剂流速和硼浓度传递给一回路系统热腿出口处和冷腿入口处,分别将蒸汽发生器给水管入口处和蒸汽室蒸汽出口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给二回路系统给水管道出口处和主蒸汽管道入口处;
步骤4:系统级分析程序根据Tn时刻的边界条件对一二回路模型进行一个时间步长的计算,将计算结果传递给计算流体动力学软件,作为计算流体动力学软件Tn+1时刻的边界条件,包括分别将一回路系统热腿出口处和冷腿入口处的温度、压力、冷却剂流速和硼浓度传递给蒸汽发生器入口腔室入口处和出口腔室出口处,分别将二回路系统给水管道出口处和主蒸汽管道入口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给蒸汽发生器给水管入口处和蒸汽室蒸汽出口处;
步骤5:重复步骤3和步骤4,直至核动力系统一二回路模型和蒸汽发生器模型完成全时段的计算,完成蒸汽发生器在核反应堆系统一二回路中的跨维度耦合。
以上内容是针对具体方案对耦合方法的进一步阐述,未详细说明部分属本领域技术人员公知常识。
Claims (1)
1.一种核动力系统一二回路耦合计算方法,其特征在于:在使用计算流体动力学软件对蒸汽发生器的流动换热特性分析时,耦合程序使用区域重叠法进行建模,其中计算流体动力学软件通过与系统分析程序在蒸汽发生器与一二回路的两个进口与两个出口处进行数据传递,实现跨纬度耦合,能够考虑一二回路系统工况变化对蒸汽发生器计算结果的影响;
该方法包含以下步骤:
步骤1:建立蒸汽发生器的几何模型,包含入口腔室、U型管、出口腔室、给水管、上升段、汽水分离器、蒸汽室和下降段,得到几何模型后进行网格划分;建立一二回路系统的节点图,并按节点图填写系统分析程序的输入卡片,一回路包括压力容器、热腿、蒸汽发生器、稳压器、蒸汽发生器一次侧部分、阀门、冷却水泵和冷腿,二回路包括蒸汽发生器二次侧部分、汽水分离再热器、汽轮机、凝汽器、凝结水泵、低压给水加热器、除氧器、给水泵和高压给水加热器;
步骤2:根据设计工况使用计算流体动力学软件选择蒸汽发生器求解计算时的湍流模型、传热模型、多相流模型和组分输运模型,设置T0时刻即初始时刻的边界条件和计算域的初始状态;
步骤3:计算流体动力学软件根据Tn时刻的边界条件对蒸汽发生器模型进行一个时间步长的计算,其中n=0、1、2、3……,将计算结果传递给系统分析程序,作为系统分析程序Tn时刻的边界条件,包括分别将蒸汽发生器入口腔室入口处和出口腔室出口处的温度、压力、冷却剂流速和硼浓度传递给一回路系统热腿出口处和冷腿入口处,分别将蒸汽发生器给水管入口处和蒸汽室蒸汽出口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给二回路系统给水管道出口处和主蒸汽管道入口处;
步骤4:系统级分析程序根据Tn时刻的边界条件对一二回路模型进行一个时间步长的计算,将计算结果传递给计算流体动力学软件,作为计算流体动力学软件Tn+1时刻的边界条件,包括分别将一回路系统热腿出口处和冷腿入口处的温度、压力、冷却剂流速和硼浓度传递给蒸汽发生器入口腔室入口处和出口腔室出口处,分别将二回路系统给水管道出口处和主蒸汽管道入口处的温度、压力、流体流速、蒸汽流速和空泡份额传递给蒸汽发生器给水管入口处和蒸汽室蒸汽出口处;
步骤5:重复步骤3、步骤4,直至核动力系统的一二回路模型和蒸汽发生器模型完成全时段的计算,完成蒸汽发生器在核反应堆系统一二回路的跨维度耦合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011035522.7A CN112231990B (zh) | 2020-09-27 | 2020-09-27 | 一种核动力系统一二回路耦合计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011035522.7A CN112231990B (zh) | 2020-09-27 | 2020-09-27 | 一种核动力系统一二回路耦合计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112231990A CN112231990A (zh) | 2021-01-15 |
CN112231990B true CN112231990B (zh) | 2021-12-28 |
Family
ID=74119381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011035522.7A Active CN112231990B (zh) | 2020-09-27 | 2020-09-27 | 一种核动力系统一二回路耦合计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112231990B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113191093B (zh) * | 2021-03-22 | 2024-02-06 | 国核电力规划设计研究院有限公司 | 压水堆二回路汽水流动模型建立方法及多工况模拟方法 |
CN113312727B (zh) * | 2021-06-21 | 2022-04-12 | 哈尔滨工程大学 | 一种高效的核电站阀门在线设计优化方法 |
CN113486608B (zh) * | 2021-07-05 | 2022-07-26 | 西安交通大学 | 一种核动力系统二回路模块化建模仿真方法 |
CN116362155B (zh) * | 2023-03-22 | 2024-01-30 | 西安交通大学 | 一种液态金属直流蒸汽发生器腔室换热系数计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015061641A1 (en) * | 2013-10-24 | 2015-04-30 | Holtec International | Steam generator for nuclear steam supply system |
CN109002741A (zh) * | 2018-06-14 | 2018-12-14 | 中国电力科学研究院有限公司 | 一种压水堆核电机组一、二回路系统传递功率模拟方法及系统 |
CN109903870A (zh) * | 2019-03-15 | 2019-06-18 | 西安交通大学 | 一种核动力系统跨维度耦合模拟方法 |
CN110020476A (zh) * | 2019-04-08 | 2019-07-16 | 西安交通大学 | 一种反应堆u型管式蒸汽发生器全三维耦合模型建立方法 |
CN111261232A (zh) * | 2020-01-14 | 2020-06-09 | 西安热工研究院有限公司 | 反应堆一回路冷却剂流场、温度场和应力场耦合计算方法 |
-
2020
- 2020-09-27 CN CN202011035522.7A patent/CN112231990B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015061641A1 (en) * | 2013-10-24 | 2015-04-30 | Holtec International | Steam generator for nuclear steam supply system |
CN109002741A (zh) * | 2018-06-14 | 2018-12-14 | 中国电力科学研究院有限公司 | 一种压水堆核电机组一、二回路系统传递功率模拟方法及系统 |
CN109903870A (zh) * | 2019-03-15 | 2019-06-18 | 西安交通大学 | 一种核动力系统跨维度耦合模拟方法 |
CN110020476A (zh) * | 2019-04-08 | 2019-07-16 | 西安交通大学 | 一种反应堆u型管式蒸汽发生器全三维耦合模型建立方法 |
CN111261232A (zh) * | 2020-01-14 | 2020-06-09 | 西安热工研究院有限公司 | 反应堆一回路冷却剂流场、温度场和应力场耦合计算方法 |
Non-Patent Citations (4)
Title |
---|
THREE DIMENSIONAL THERMAL HYDRAULIC COUPLING METHOD ON PWR STEAM GENERATOR USING CFD METHOD;Xiaohan Zhao, Mingjun Wang, Wenxi Tian, G.H. Su, Suizheng QIU;《•The Proceedings of the International Conference on Nuclear Engineering (ICONE) 》;20191231;第1-5页 * |
丛腾龙 ; 田文喜 ; 秋穗正 ; 苏光辉.耦合一、二次侧换热的蒸汽发生器二次侧流场分析.《原子能科学技术》.2014, * |
王坤 ; 董秀臣 ; 张鑫 ; 袁江涛.小型压水堆偏环路运行堆芯入口流量分配耦合分析.《舰船科学技术》.2019, * |
田文喜 ; 王明军 ; 秋穗正 ; 苏光辉.基于CFD方法的核动力系统热工安全特性研究进展.《原子能科学技术》.2019, * |
Also Published As
Publication number | Publication date |
---|---|
CN112231990A (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112231990B (zh) | 一种核动力系统一二回路耦合计算方法 | |
CN108595752A (zh) | 一种面向核动力装置仿真应用的单相水流体网络仿真方法 | |
CN104484528A (zh) | 一种基于pscad的核电机组动力、电气混合仿真平台 | |
CN116362155B (zh) | 一种液态金属直流蒸汽发生器腔室换热系数计算方法 | |
CN114239432B (zh) | 基于Modelica的船舶二回路系统换热器动态仿真模型开发方法 | |
Papukchiev et al. | Multiscale analysis of forced and natural convection including heat transfer phenomena in the TALL-3D experimental facility | |
CN112182905B (zh) | 一种用于综合能源系统的供热管网仿真方法和装置 | |
Allison et al. | The development of RELAP/SCDAPSIM/MOD4. 0 for advanced fluid systems design analysis | |
Kim et al. | Preliminary coupling of MATRA code for multi-physics analysis | |
Toti et al. | Development and preliminary validation of a STH-CFD coupling method for multiscale thermal-hydraulic simulations of the MYRRHA reactor | |
Soler-Martinez | Semi-implicit thermal-hydraulic coupling of advanced subchannel and system codes for pressurized water reactor transient applications | |
CN113779713B (zh) | 考虑壳侧横向交混的钠-水直流蒸汽发生器精细建模方法 | |
CN114528783A (zh) | 一种基于信号滤波核电系统程序的两相流耦合方法 | |
Crump et al. | Nonlinear transient analysis of light water reactor steam generators using an implicit Eulerian method | |
CN116484764A (zh) | 一种钠冷快堆多维度耦合计算方法 | |
Saphier | DSNP models used in the pebble-bed HTGR dynamic simulation. V. 2 | |
Bai et al. | Long-Term Simulation of Sodium Dynamics During a Large Leakage Sodium-Water Reaction | |
Yang et al. | Boundary condition coupling methods and its application to BOP-integrated transient simulation of SMART | |
CN117828830A (zh) | 一种核电厂多回路热工水力系统混合阶次一体化仿真方法 | |
Fu et al. | Qualitative and Quantitative Evaluation of Coupling Approaches for Coupling of RELAP5 and LabVIEW | |
Nassersharif et al. | Comparison of three-dimensional fine and coarse nodalization of the reactor vessel for a small-break loss-of-coolant accident | |
Williams et al. | Thermal Extraction Modeling within HYBRID: Nominal and Extraction Conditions in Advanced Reactor Systems | |
Bin et al. | Modeling of main steam and two-phase heat exchanger for nuclear power unit | |
Martelli et al. | Validation of the coupled calculation between RELAP5 STH code and Ansys FLUENT CFD code | |
Agrawal et al. | Simulation of transients in liquid-metal fast breeder reactor systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |