CN110826247A - 一种核反应堆重力型非能动安注设备特性的模拟方法 - Google Patents
一种核反应堆重力型非能动安注设备特性的模拟方法 Download PDFInfo
- Publication number
- CN110826247A CN110826247A CN201911126756.XA CN201911126756A CN110826247A CN 110826247 A CN110826247 A CN 110826247A CN 201911126756 A CN201911126756 A CN 201911126756A CN 110826247 A CN110826247 A CN 110826247A
- Authority
- CN
- China
- Prior art keywords
- safety injection
- simulation
- dimensionless
- passive safety
- equipment
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种核反应堆重力型非能动安注设备特性的模拟方法,包括以下步骤:1)、获得模拟准则数:通过对重力驱动型安注建立质量和动量守恒方程,对方程进行合理简化和无量纲后获得模拟原型非能动安注设备特性的模拟准则数;2)、确保试验装置与原型设备的模拟准则数相同,获得试验装置主要参数的模拟比例;3)、在满足压力容器、非能动安注设备和安注管线等比例缩小的前提下,根据模拟比例和等高模拟的要求设置试验装置和参数控制。本发明能够以小规模的实验规模比例,实时、准确地模拟核反应堆发生失水事故后重力型非能动安注设备的安注特性以及对堆芯的冷却能力。
Description
技术领域
本发明涉及核反应堆领域,具体涉及一种核反应堆重力型非能动安注设备特性的模拟方法。
背景技术
当核反应堆发生失水事故时,一回路压力急剧下降,随后控制棒落棒实施紧急堆芯,主泵短暂惰转后停止运转。随着一回路冷却剂不断从破口喷放,反应堆堆芯内的液位不断下降,堆芯出现裸露后燃料元件的衰变热功率可能烧毁燃料元件包壳,导致放射性物质外泄。因此当反应堆发生失水事故后,往反应堆压力容器注入冷却剂对堆芯进行有效冷却是防止核反应堆发生放射性物质外泄的重要措施。
目前核反应堆采用安注泵对失水事故下的堆芯进行补水。日本福岛核事故后,全世界范围内对反应堆的安全提出了更高的要求。采用非能动安注方式往堆芯注入冷流体,对失水事故条件下的核反应堆堆芯进行有效冷却是当今先进核反应堆设计和研发的一个重要方向。
由于非能动安注方式产生的安注驱动力较小,且受诸多因素的影响,如安注位差和回路阻力等。因此新型核反应堆失水事故下非能动安注系统的安注和冷却堆芯能力需要经过实验验证。
发明内容
本发明的目的在于提供一种核反应堆重力型非能动安注设备特性的模拟方法,能够以小规模的实验规模比例,实时、准确地模拟核反应堆发生失水事故后重力型非能动安注设备的安注特性以及对堆芯的冷却能力。
本发明通过下述技术方案实现:
一种核反应堆重力型非能动安注设备特性的模拟方法,包括以下步骤:
1)、获得模拟准则数:
a1)、根据质量守恒方程和能量守恒方程获得对安注过程进行描述的公式:
式中,AE为安注管线出口流通截面积,AT为安注设备流通截面积,ΠTE为摩擦数,ΔP为静压差,ρl为流体密度,L为液位,g为重力加速度,M为流体质量,t为时间;
b1)、考虑重力型非能动安注设备的特性对步骤1)获得的公式进行简化,令ΔP=0,M=ρlATaL,(AE/AT)2≈0,获得简化公式:
c1)、引入无量纲参数将步骤2)获得的公式进行无量纲化得到无量纲方程:
式中,L+为无量纲液位,t+为无量纲时间,M+为无量纲质量,L0为无量纲初始液位,M0为无量纲初始质量,W0为无量纲初始质量流速;
将无量纲方程解析得到:
通过步骤1)获得的公式得到重力驱动型非能动设备初始安注流速的表达式:
d1)、获得非能动安注模拟需要满足的模拟准则数:
式中,LR为液位比,d为直径,K为局部形阻系数,f为摩擦系数,l为管道长度;
2)、根据试验装置与原型设备相似性的原则,确保试验装置与原型设备的模拟准则数相同,根据模拟准则数的要求,获得试验装置主要参数的模拟比例,进而获得所有设备尺寸和结构参数的模拟比例;
3)、在满足压力容器、非能动安注设备和安注管线等比例缩小的前提下,根据模拟比例和等高模拟的要求设置试验装置和参数控制,依次满足试验装置与原型密度比准则数相似、液位比准则数相似、时间比准则数相似、摩擦数相似和初始速度相似。
由于核反应堆的功率以及几何尺寸,使得大多数情况下不能进行全尺寸的实验研究和验证。因此进行缩小功率规模,减小系统尺寸的方法的模拟研究是非常必要的。
所谓反应堆非能动安注是指,反应堆发生失水事故时,反应堆压力容器内大量冷却剂通过破口流出,堆芯燃料冷却冷却能力骤降。为了保护反应堆堆芯的安全,非能动安注设备布置在反应堆压力容器上部,非能动安注设备内流体在重力的作用下,沿非能动系统安注管线注入反应堆压力容器,对堆芯进行应急冷却,避免堆芯融化放射性物质外泄。为了研究反应堆失水事故下的安全特性,试验装置需要对非能动安注设备特性进行准确模拟。
本发明通过对重力驱动型安注建立质量和动量守恒方程,对方程进行合理简化和无量纲后获得模拟原型非能动安注设备特性的模拟准则数,通过保证试验装置的模拟准则数与原型一一相等,即可真实获得原型反应堆系统非能动系统的安注特性。如此,本发明能够以小规模的实验规模比例,实时、准确地模拟核反应堆发生失水事故后重力型非能动安注设备的安注特性以及对堆芯的冷却能力。
本发明提供了一种可模拟核反应堆重力型非能动安注设备特性的方法,为设计建造大型非能动安注系统热工水力整体性能试验装置提供技术参考。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明通过对重力驱动型安注建立质量和动量守恒方程,对方程进行合理简化和无量纲后获得模拟原型非能动安注设备特性的模拟准则数,通过保证试验装置的模拟准则数与原型一一相等,即可真实获得原型反应堆系统非能动系统的安注特性。如此,本发明能够以小规模的实验规模比例,实时、准确地模拟核反应堆发生失水事故后重力型非能动安注设备的安注特性以及对堆芯的冷却能力。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明所述试验装置的示意图;。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1:
如图1所示,本发明的试验装置主要包括反应堆压力容器模拟体、非能动安注设备模拟体、阀门、流量计和回路管道,试验装置为根据原型等比缩小,根据试验装置的缩比因子λ可确定试验装置所有设备和管道的面积缩比因子为λ,进而可以得到试验装置所有设备和管道的流道面积。
一种核反应堆重力型非能动安注设备特性的模拟方法,包括以下步骤:
1)、获得模拟准则数:
a1)、根据质量守恒方程和能量守恒方程获得对安注过程进行描述的公式:
式中,AE为安注管线出口流通截面积,AT为安注设备流通截面积,ΠTE为摩擦数,ΔP为静压差,ρl为流体密度,L为液位,g为重力加速度,M为流体质量,t为时间;
具体地,质量守恒方程:
式中,质量流量W=ρluTAT=ρluEAE,ρ为流体密度,u为流速,A为流通截面积,脚注T(ank)为安注设备,E(xit)为安注管线出口。
动量守恒方程(Bernoulli方程):
联立式(1)和式(2)可以得到:
b1)、考虑重力型非能动安注设备的特性对步骤1)获得的公式进行简化,令ΔP=0,M=ρlATaL,(AE/AT)2≈0,获得简化公式:
具体地,
重力型非能动安注设备的安注驱动力主要来源于高度差,静压差可以忽略,即ΔP=0;不妨假设安注箱为柱状筒结构,则有M=ρlATaL;若安注箱截面积远远大于安注管线的截面积,即AT>>AE,则(AE/AT)2≈0;根据上述条件及假设,式(3)可转化为:
c1)、引入无量纲参数将步骤2)获得的公式进行无量纲化得到无量纲方程:
式中,L+为无量纲液位,t+为无量纲时间,M+为无量纲质量,L0为无量纲初始液位,M0为无量纲初始质量,W0为无量纲初始质量流速;
式(6)就是“重力驱动型”非能动安注设备内冷却剂无量纲水位变化微分方程。根据式(6)可以得到解析解:
式中,0≤L+≤1,0≤t+≤2。
根据式(3)可以得到重力驱动型非能动设备初始安注流速的表达式:
d1)、获得非能动安注模拟需要满足的模拟准则数:
式中,LR为液位比,d为直径,K为局部形阻系数,f为摩擦系数,l为管道长度;
2)、根据试验装置与原型设备相似性的原则,确保试验装置与原型设备的模拟准则数相同,根据模拟准则数的要求,获得试验装置主要参数的模拟比例,进而获得所有设备尺寸和结构参数的模拟比例;
3)、在满足压力容器、非能动安注设备和安注管线等比例缩小的前提下,根据模拟比例和等高模拟的要求设置试验装置和参数控制,依次满足试验装置与原型密度比准则数相似、液位比准则数相似、时间比准则数相似、摩擦数相似和初始速度相似:
具体地:
根据等高模拟准则的要求,试验装置的反应堆压力容器、非能动安注设备的设备高度与原型核反应堆的反应堆压力容器、非能动安注设备的设备高度相同。试验装置反应堆压力容器与非能动安注设备的相对高度与原型核反应堆中反应堆压力容器与非能动安注设备的相对高度相同。试验装置非能动安注系统管道高度和长度均与原型反应堆非能动安注系统管道和长度一致。由于原型反应堆非能动安注系统的管道直径大,相同自然循环流速下沿程摩擦阻力系数下,为了保证试验装置非能动安注系统的总阻力系数与原型反应堆一回路系统的总阻力系数相同,增加了阻力系数非常小的流量计用于测量非能动安注流量,省略了部分阀门的模拟(如阻力系数较大的隔离阀)(根据摩擦数相同的原则,通过减小阻力系数的方法保证试验装置的总阻力系数与原型反应堆的总阻力系数相同)。根据相同流体工质模拟的需求,采用与原型反应堆冷却剂相同的工质材料——轻水作为实验用冷却剂。根据等温等压模拟的要求,原型非能动安注设备初始流体压力为15.0MPa,温度50℃,试验装置采用相同的运行压力和温度:非能动安注设备运行压力15.0MPa,温度50℃。假设试验装置与原型反应堆的缩比因子为100,原型反应堆非能动安注设备的内径为900mm,为保证非能动安注设备内的冷却剂总质量为原型设备的1/100,那么试验装置非能动安注设备的内径为90mm。原型反应堆非能动安注管线的内径为200mm,长度为7.6m,那么试验装置非能动安注管线的内径为20mm,长度为7.6m。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种核反应堆重力型非能动安注设备特性的模拟方法,其特征在于,包括以下步骤:
1)、获得模拟准则数:
a1)、根据质量守恒方程和能量守恒方程获得对安注过程进行描述的公式:
式中,AE为安注管线出口流通截面积,AT为安注设备流通截面积,ΠTE为摩擦数,ΔP为静压差,ρl为流体密度,L为液位,g为重力加速度,M为流体质量,t为时间;
b1)、考虑重力型非能动安注设备的特性对步骤1)获得的公式进行简化,令ΔP=0,M=ρlATaL,(AE/AT)2≈0,获得简化公式:
c1)、引入无量纲参数将步骤2)获得的公式进行无量纲化得到无量纲方程:
式中,L+为无量纲液位,t+为无量纲时间,M+为无量纲质量,L0为无量纲初始液位,M0为无量纲初始质量,W0为无量纲初始质量流速;
将无量纲方程解析得到:
通过步骤1)获得的公式得到重力驱动型非能动设备初始安注流速的表达式:
d1)、获得非能动安注模拟需要满足的模拟准则数:
式中,LR为液位比,d为直径,K为局部形阻系数,f为摩擦系数,l为管道长度;
2)、根据试验装置与原型设备相似性的原则,确保试验装置与原型设备的模拟准则数相同,根据模拟准则数的要求,获得试验装置主要参数的模拟比例,进而获得所有设备尺寸和结构参数的模拟比例;
3)、在满足压力容器、非能动安注设备和安注管线等比例缩小的前提下,根据模拟比例和等高模拟的要求设置试验装置和参数控制,依次满足试验装置与原型密度比准则数相似、液位比准则数相似、时间比准则数相似、摩擦数相似和初始速度相似。
2.根据权利要求1所述的一种核反应堆重力型非能动安注设备特性的模拟方法,其特征在于,试验装置采用与原型采用相同的运行压力、温度和相同流体工质实现试验装置采的物性可以与原型一致,以满足密度比准则数相似:
6.根据权利要求1所述的一种核反应堆重力型非能动安注设备特性的模拟方法,其特征在于,满足L0=1、后,带入重力驱动型非能动设备初始安注流速的表达式,获得初始速度相似:u0=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911126756.XA CN110826247B (zh) | 2019-11-18 | 2019-11-18 | 一种核反应堆重力型非能动安注设备特性的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911126756.XA CN110826247B (zh) | 2019-11-18 | 2019-11-18 | 一种核反应堆重力型非能动安注设备特性的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110826247A true CN110826247A (zh) | 2020-02-21 |
CN110826247B CN110826247B (zh) | 2022-07-29 |
Family
ID=69556224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911126756.XA Active CN110826247B (zh) | 2019-11-18 | 2019-11-18 | 一种核反应堆重力型非能动安注设备特性的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110826247B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016218A (zh) * | 2020-09-24 | 2020-12-01 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977529A (en) * | 1973-02-23 | 1990-12-11 | Westinghouse Electric Corp. | Training simulator for a nuclear power plant |
US20050279867A1 (en) * | 2003-03-14 | 2005-12-22 | Ismailov Murad M | Systems and methods for operating an electromagnetic actuator |
CN107527664A (zh) * | 2017-07-31 | 2017-12-29 | 中国核动力研究设计院 | 一种模拟核反应堆流体自然循环流动特性的方法 |
CN109765067A (zh) * | 2017-11-10 | 2019-05-17 | 国核华清(北京)核电技术研发中心有限公司 | 大型先进压水堆核电站非能动堆芯冷却系统整体性能试验平台 |
-
2019
- 2019-11-18 CN CN201911126756.XA patent/CN110826247B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4977529A (en) * | 1973-02-23 | 1990-12-11 | Westinghouse Electric Corp. | Training simulator for a nuclear power plant |
US20050279867A1 (en) * | 2003-03-14 | 2005-12-22 | Ismailov Murad M | Systems and methods for operating an electromagnetic actuator |
CN107527664A (zh) * | 2017-07-31 | 2017-12-29 | 中国核动力研究设计院 | 一种模拟核反应堆流体自然循环流动特性的方法 |
CN109765067A (zh) * | 2017-11-10 | 2019-05-17 | 国核华清(北京)核电技术研发中心有限公司 | 大型先进压水堆核电站非能动堆芯冷却系统整体性能试验平台 |
Non-Patent Citations (4)
Title |
---|
PENG CHUANXIN等: ""Research and Evaluation for Passive Safety System in Low Pressure Reactor"", 《SCIENCE AND TECHNOLOGY OF NUCLEAR INSTALLATIONS》 * |
XIANGBIN LI等: ""Application of dynamical system scaling method on simple gravity-driven draining process"", 《JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY》 * |
卢冬华等: "压水堆自然循环比例模化基本方程及相似准则数的研究", 《核动力工程》 * |
张盼等: "核电厂全厂断电工况下蒸汽发生器自然循环特性的比例分析", 《核科学与工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016218A (zh) * | 2020-09-24 | 2020-12-01 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
CN112016218B (zh) * | 2020-09-24 | 2022-07-22 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN110826247B (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Saadati et al. | Safety margin and fuel cycle period enhancements of VVER-1000 nuclear reactor using water/silver nanofluid | |
CN107527664A (zh) | 一种模拟核反应堆流体自然循环流动特性的方法 | |
Li et al. | Experimental investigation on thermal stratification induced by steam direct contact condensation with non-condensable gas | |
Li et al. | Experimental investigation on thermal stratification induced by steam-air mixture vertical injection with shallow submergence depth | |
Sun et al. | Scaling analysis of AP1000 ADS-4 entrainment and depressurization | |
CN110826247B (zh) | 一种核反应堆重力型非能动安注设备特性的模拟方法 | |
Yanhua et al. | Analysis of diffusion process and influence factors in the air ingress accident of the HTR-PM | |
Avery et al. | Critical heat flux in TRIGA-fueled reactors cooled by natural convection | |
Sun et al. | Experimental investigation of upper plenum entrainment in AP1000 | |
Tong et al. | Numerical investigation on erosion of hydrogen stratification by steam jet within a local compartment | |
Chikhi et al. | Effect of water entrainment on the coolability of a debris bed surrounded by a by-pass: Integral reflood experiments and modelling | |
Chun et al. | Critical heat flux in uniformly heated vertical annulus under a wide range of pressures 0.57 to 15.0 MPa | |
Moon et al. | An experimental study on the critical heat flux for low flow of water in a non-uniformly heated vertical rod bundle over a wide range of pressure conditions | |
Torres et al. | A MTR fuel element flow distribution measurement preliminary results | |
Cheng et al. | Suppression pool mixing and condensation tests in PUMA facility | |
Rashid | Coolability of volumetrically heated particle beds | |
Repetto et al. | Core coolability in loss of coolant accident: the PERFROI project | |
Modak et al. | Experimental investigation of quenching phenomena in high temperature conical debris bed–Part 1: MONET tests | |
Sarchami et al. | Comparison between surface heating and volumetric heating methods inside CANDU reactor moderator test facility (MTF) using 3D numerical simulation | |
Chavez et al. | Further development of a thermal-hydraulics two-phase flow tool | |
Yang et al. | Numerical investigation on heat transfer to supercritical water flowing upward in a 4-m long bare vertical circular tube | |
Jeltsov | Application of CFD to safety and thermal-hydraulic analysis of lead-cooled systems | |
Zhu et al. | DILUTION AND MIXING CHARACTERISTICS OF BORIC ACID IN REACTOR PRESSURE VESSEL WITH DIRECT VESSEL INJECTION | |
Wang et al. | CFD Simulation of Natural Ventilation Performance of the Interim Region in Spent Fuel Dry-Storage System of HTR-PM | |
Nomura et al. | Study on control method of natural circulation by injection of helium gas |
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 |