CN111581806A - 瞬变外力对动态自反馈条件下通道内chf的影响分析方法 - Google Patents

瞬变外力对动态自反馈条件下通道内chf的影响分析方法 Download PDF

Info

Publication number
CN111581806A
CN111581806A CN202010363721.4A CN202010363721A CN111581806A CN 111581806 A CN111581806 A CN 111581806A CN 202010363721 A CN202010363721 A CN 202010363721A CN 111581806 A CN111581806 A CN 111581806A
Authority
CN
China
Prior art keywords
chf
under
external force
transient
influence
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
Application number
CN202010363721.4A
Other languages
English (en)
Other versions
CN111581806B (zh
Inventor
岳倪娜
黄彦平
昝元峰
曹念
袁德文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202010363721.4A priority Critical patent/CN111581806B/zh
Publication of CN111581806A publication Critical patent/CN111581806A/zh
Application granted granted Critical
Publication of CN111581806B publication Critical patent/CN111581806B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明公开了一种瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，瞬变外力对自然循环工况下通道内CHF的影响因子包括宏观影响因子和微观影响因子；宏观影响因通过瞬变外力下的流量对应的静止条件下的CHF值与静止条件的流量对应的静止条件下的CHF值地比值获得；微观影响因子通过瞬变外力下的通道内的CHF值与进口参数相同的静止条件下的CHF值的比值获得；通过乘积法求解耦合瞬变外力对反应堆自然循环工况下通道内CHF的影响因子，最终获得瞬变外力对自然循环工况下DNBR限值修正系数。本发明可用于分析不同装置中瞬变外力对自然循环工况下通道内DNBR限值修正系数，且具有一定的可推广性。

Description

瞬变外力对动态自反馈条件下通道内CHF的影响分析方法
技术领域
本发明涉及核动力技术领域，具体涉及一种瞬变外力对动态自反馈条件下通道内CHF的影响分析方法。
背景技术
为了提高核动力舰船的隐蔽性，自然循环运行工况是目前核动力舰船研究的重点。另外，核动力反应堆堆芯通道内的临界热流密度(Critical Heat Flux,CHF)关系到反应堆的安全，通过通道内偏离泡核沸腾比(Departure from Nucleate Boiling Ratio,DNBR)限值来表示通道内发生CHF的裕量，在陆基反应堆中DNBR限值是反应堆堆芯设计的准则，当DNBR限值应用到核动力舰船反应堆中时由于瞬变外力的影响需要对其进行修正，其修正值为瞬变外力对CHF的影响因子的最小值。核动力舰船在工作过程中会受海浪等影响发生倾斜、起伏和摇摆等，而倾斜条件对CHF的影响尤为重要。
目前国内DNBR限值修正系数是参考日本isshi的研究以及国内相关研究选取的，此值是在强迫循环条件下获得，强迫循环条件下瞬变外力对CHF的影响远小于自然循环的影响，因此，此值不能直接应用到自然循环工况下的堆芯设计。
经专利查新，目前国内外没有针对瞬变外力对自然循环工况下CHF影响因子的分析方法，因此，有必要研究一种瞬变外力对自然循环工况下CHF影响因子的分析方法，清晰、准确的表达瞬变外力对自然循环工况下通道内CHF的影响，能够应用于瞬变外力自然循环工况下DNBR限值的修正，从而在此设计准则条件下的对象能够运行于自然循环工况下，从而提高舰艇的隐蔽性和作战能力。
发明内容
针对上述问题，本发明提供了一种瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，本发明清晰、准确地表达了瞬变外力对自然循环工况下CHF各方面的影响，使得通过此分析方法得到的DNBR限值的修正系数能够应用于核动力舰艇反应堆的设计。
本发明通过下述技术方案实现：
一种瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，瞬变外力对自然循环工况下通道内CHF的影响因子包括宏观影响因子和微观影响因子：
其中，宏观影响因子kf如式(1)所示：
Figure BDA0002475979500000011
qc(GJZ-ΔG)表示瞬变外力下的流量对应的静止条件下的CHF值；qc(GJZ)表示静止条件的流量对应的静止条件下的CHF值；
其中微观影响因子kL如式(2)所示：
Figure BDA0002475979500000021
qHY,CHF(GJZ-ΔG)表示瞬变外力下的通道内的CHF值；qc(GJZ-ΔG)表示进口参数相同的静止条件下的CHF值；所述进口参数包括温度、流速、压力；
上述GJZ表示静止条件的流量，ΔG表示瞬变外力场带来的流速变化量。
进一步地，静止条件下沸腾临界关系式如式(3)所示：
qc＝f(G,x,P) (3)；
G表示静止条件下通道内流量；x表示出口含汽率，P表示出口处压力。
进一步地，所述瞬变外力包括倾斜、起伏和摇摆。
进一步地，瞬变外力对自然循环工况下通道内CHF的影响因子值k为宏观影响因子和微观影响因子的乘积：k＝kf×kL。
瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的分离式求解方法，通过乘积法求解耦合瞬变外力对反应堆自然循环工况下通道内CHF的影响因子，最终获得瞬变外力对自然循环工况下DNBR限值修正系数；
其中，瞬变外力对反应堆自然循环工况下通道内CHF的影响因子基于上述的瞬变外力对动态自反馈条件下通道内CHF的影响分析方法获得。
进一步地，所述瞬变外力对自然循环工况下DNBR限值修正系数的计算公式如下式所示：
K＝KQX×KSQ×KYB；
其中，K表示瞬变外力对自然循环工况下DNBR限值修正系数；
KQX表示倾斜瞬变外力下CHF的最小影响因子；
KSQ表示瞬变外力下CHF的最小影响因子；
KYB表示摇摆瞬变外力下CHF的最小影响因子。
本发明具有如下的优点和有益效果：
本发明提出的瞬变外力对自然循环工况下CHF影响因子的分解式分析方法，未见国内外存在相关报道，填补了国内空白，填补瞬变外力对自然循环工况下通道内CHF影响因子的分析方法的空白，并清晰、准确的表达瞬变外力对自然循环工况下CHF的各方面的影响；解决了瞬变外力对自然循环工况下CHF的影响因子不具有拓展性的问题，通过此分析方法得到的DNBR限值的修正值能够应用于核动力舰艇反应堆堆芯设计所需的准则中。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明的瞬变外力对自然循环工况下CHF的影响分解分析结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1
本实施例提供了一种瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，本发明首先通过实验对瞬变外力自然循环工况下通道内的CHF影响因素开展了研究工作，并获得了静止条件下的CHF的关系式。通过实验研究发现瞬变外力对自然循环工况下CHF的影响分两个方面，宏观方面和微观方面，CHF影响因子的分离式求解方法是将两方面的影响分开求解，如图1所示。具体步骤如下：
步骤1，根据实验数据拟合得到静止条件下沸腾临界关系式：qc＝f(G,x,P)。
步骤2，求解典型瞬变外力对自然循环工况下通道内CHF的宏观方面的影响，宏观方面的影响主要是由于瞬变外力改变了系统的流量从而引起局部参数的改变。通过实验研究瞬变外力带来的流量的变化量ΔG。
倾斜条件：ΔG＝GQX-GJZ；
瞬变外力场条件：ΔG＝G波谷-GJZ；
根据实验数据拟合得到静止条件下沸腾临界关系式qc＝f(G,x,P)，求解静止条件的流量与瞬变外力下的流量对应的静止条件下沸腾临界值(CHF值)，瞬变外力下流量对应的CHF值与静止条件下流量对应的CHF值的比值为宏观影响因子kf：
Figure BDA0002475979500000031
qc(GJZ-ΔG)表示瞬变外力下的流量对应的静止条件下的CHF值；qc(GJZ)表示静止条件的流量对应的静止条件下的CHF值。
步骤3，求解典型瞬变外力对对自然循环工况下通道内CHF的微观方面的影响，微观方面的影响主要是由于瞬变外力改变了通道内气泡、液膜的分布等局部特性，从而引起CHF的改变。此部分影响通过瞬变外力下对应的通道内的CHF值与进口参数相同的静止条件下的CHF的比值求得微观影响因子kL：
Figure BDA0002475979500000041
qHY,CHF(GJZ-ΔG)表示瞬变外力下的通道内的CHF值；qc(GJZ-ΔG)表示进口参数相同的静止条件下的CHF值；所述进口参数包括温度、压力和流速。
步骤4，求解某典型瞬变外力对自然循环工况下通道内CHF的影响因子，瞬变外力对自然循环工况下通道内CHF的影响因子值k为宏观影响因子和微观影响因子的乘积：k＝kf×kL。
步骤5，依据步骤1至4，依次计算倾斜条件、起伏条件和摇摆条件对自然循环工况下通道内CHF的影响因子，并获得倾斜条件、起伏条件和摇摆条件对自然循环工况下通道内CHF的最小影响因子KQX、KSQ、KYB。
步骤6，通过乘积法求解耦合瞬变外力对反应堆自然循环工况下通道内CHF的影响因子，最终获得瞬变外力对自然循环工况下DNBR限值修正系数。瞬变外力对自然循环工况下DNBR限值修正系数的计算公式如下式所示：
K＝KQX×KSQ×KYB；
其中，K表示瞬变外力对自然循环工况下DNBR限值修正系数；
KQX表示倾斜瞬变外力下CHF的最小影响因子；
KSQ表示瞬变外力下CHF的最小影响因子；
KYB表示摇摆瞬变外力下CHF的最小影响因子。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，其特征在于，瞬变外力对自然循环工况下通道内CHF的影响因子包括宏观影响因子和微观影响因子：
其中，宏观影响因子kf如式(1)所示：
Figure FDA0002475979490000011
qc(GJZ-ΔG)表示瞬变外力下的流量对应的静止条件下的CHF值；qc(GJZ)表示静止条件的流量对应的静止条件下的CHF值；
其中，微观影响因子kL如式(2)所示：
Figure FDA0002475979490000012
qHY,CHF(GJZ-ΔG)表示瞬变外力下的通道内的CHF值；qc(GJZ-ΔG)表示进口参数相同的静止条件下的CHF值；所述进口参数包括温度、流速、压力；
上述GJZ表示静止条件的流量，ΔG表示瞬变外力场带来的流速变化量。
2.根据权利要求1所述的瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，其特征在于，静止条件下沸腾临界关系式如式(3)所示：
qc＝f(G,x,P) (3)；
G表示静止条件下通道内流量；x表示出口含汽率，P表示出口处压力。
3.根据权利要求1所述的瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，其特征在于，所述瞬变外力包括倾斜、起伏和摇摆。
4.根据权利要求1至3任一项所述的瞬变外力对动态自反馈条件下通道内CHF的影响分析方法，其特征在于，瞬变外力对自然循环工况下通道内CHF的影响因子值k为宏观影响因子和微观影响因子的乘积：k＝kf×kL。
5.瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的分离式求解方法，其特征在于，通过乘积法求解耦合瞬变外力对反应堆自然循环工况下通道内CHF的影响因子，最终获得瞬变外力对自然循环工况下DNBR限值修正系数；
其中，瞬变外力对反应堆自然循环工况下通道内CHF的影响因子基于权利要求1至4任一项所述的瞬变外力对动态自反馈条件下通道内CHF的影响分析方法获得。
6.根据权利要求5所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的分离式求解方法，其特征在于，所述瞬变外力对自然循环工况下DNBR限值修正系数的计算公式如下式所示：
K＝KQX×KSQ×KYB；
其中，K表示瞬变外力对自然循环工况下DNBR限值修正系数；
KQX表示倾斜瞬变外力下CHF的最小影响因子；
KSQ表示瞬变外力下CHF的最小影响因子；
KYB表示摇摆瞬变外力下CHF的最小影响因子。
CN202010363721.4A 2020-04-30 2020-04-30 瞬变外力对动态自反馈条件下通道内chf的影响分析方法 Active CN111581806B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010363721.4A CN111581806B (zh) 2020-04-30 2020-04-30 瞬变外力对动态自反馈条件下通道内chf的影响分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010363721.4A CN111581806B (zh) 2020-04-30 2020-04-30 瞬变外力对动态自反馈条件下通道内chf的影响分析方法

Publications (2)

Publication Number Publication Date
CN111581806A true CN111581806A (zh) 2020-08-25
CN111581806B CN111581806B (zh) 2020-12-04

Family

ID=72113305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010363721.4A Active CN111581806B (zh) 2020-04-30 2020-04-30 瞬变外力对动态自反馈条件下通道内chf的影响分析方法

Country Status (1)

Country Link
CN (1) CN111581806B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112420228A (zh) * 2020-11-19 2021-02-26 中国核动力研究设计院 瞬变外力场作用下棒束通道chf的流道选取方法及系统
CN115688488A (zh) * 2022-12-30 2023-02-03 中国核动力研究设计院 准则确定方法、装置、设备、计算机存储介质及程序产品
CN115862909A (zh) * 2023-02-14 2023-03-28 中国核动力研究设计院 堆芯热工安全准则应用方法、装置、设备、介质及产品
CN116884655A (zh) * 2023-09-08 2023-10-13 中国核动力研究设计院 外力场对热工安全影响确定方法、装置、核反应堆及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106531246A (zh) * 2016-11-25 2017-03-22 中国核动力研究设计院 运动条件下非能动余热排出系统实验初始工况建立方法及系统
CN106886686A (zh) * 2017-03-03 2017-06-23 西安交通大学 一种压水堆少群常数历史效应的复合修正方法
CN106897478A (zh) * 2016-12-19 2017-06-27 杨保文 基于阻力分布包含格架搅混效应的子通道分析方法
CN107122546A (zh) * 2017-04-27 2017-09-01 西安交通大学 一种压水堆稳态计算的多物理耦合方法
CN107239606A (zh) * 2017-05-27 2017-10-10 国网福建省电力有限公司 一种用于压水堆系统动态模型参数评价的灵敏度分析方法
CN109448874A (zh) * 2018-09-29 2019-03-08 国核自仪系统工程有限公司 核动力反应堆保护系统的停堆功能的测试装置及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106531246A (zh) * 2016-11-25 2017-03-22 中国核动力研究设计院 运动条件下非能动余热排出系统实验初始工况建立方法及系统
CN106897478A (zh) * 2016-12-19 2017-06-27 杨保文 基于阻力分布包含格架搅混效应的子通道分析方法
CN106886686A (zh) * 2017-03-03 2017-06-23 西安交通大学 一种压水堆少群常数历史效应的复合修正方法
CN107122546A (zh) * 2017-04-27 2017-09-01 西安交通大学 一种压水堆稳态计算的多物理耦合方法
CN107239606A (zh) * 2017-05-27 2017-10-10 国网福建省电力有限公司 一种用于压水堆系统动态模型参数评价的灵敏度分析方法
CN109448874A (zh) * 2018-09-29 2019-03-08 国核自仪系统工程有限公司 核动力反应堆保护系统的停堆功能的测试装置及方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112420228A (zh) * 2020-11-19 2021-02-26 中国核动力研究设计院 瞬变外力场作用下棒束通道chf的流道选取方法及系统
CN112420228B (zh) * 2020-11-19 2022-02-11 中国核动力研究设计院 瞬变外力场作用下棒束通道chf的流道选取方法及系统
CN115688488A (zh) * 2022-12-30 2023-02-03 中国核动力研究设计院 准则确定方法、装置、设备、计算机存储介质及程序产品
CN115862909A (zh) * 2023-02-14 2023-03-28 中国核动力研究设计院 堆芯热工安全准则应用方法、装置、设备、介质及产品
CN116884655A (zh) * 2023-09-08 2023-10-13 中国核动力研究设计院 外力场对热工安全影响确定方法、装置、核反应堆及设备
CN116884655B (zh) * 2023-09-08 2023-11-10 中国核动力研究设计院 外力场对热工安全影响确定方法、装置、核反应堆及设备

Also Published As

Publication number Publication date
CN111581806B (zh) 2020-12-04

Similar Documents

Publication Publication Date Title
CN111581806B (zh) 瞬变外力对动态自反馈条件下通道内chf的影响分析方法
CN111553022B (zh) 动态自反馈条件下堆芯dnbr限值修正系数的求解方法
CN102067241B (zh) 热管裂变燃料元件
Zhang et al. Multiobjective optimization design and experimental investigation on the axial flow pump with orthogonal test approach
CN111507016B (zh) 一种运动条件动态下并联窄通道流动失稳边界的确定方法
CN110516321B (zh) 一种变速水轮机转轮选型计算方法
Kollerup et al. A mathematical model for ethanol production by extractive fermentation in a continuous stirred tank fermentor
CN110898454A (zh) 精馏塔控制方法、存储介质及电子设备
CN109409019A (zh) 一种风力机叶片气动外形的预弯优化方法
Lemos et al. In situ extractive ethanol fermentation in a drop column bioreactor
CN106222438A (zh) 一种从含锇物料中提取锇的方法
CN116752189A (zh) 一种碱性水电解制氢系统及其压力、液位调节阀控制方法
CN104817429A (zh) 从环己醇精馏塔塔底料液中回收环己醇的方法
CN108569739B (zh) 用于回收聚四氟乙烯浓缩分散液废液中op-10的系统及方法
CN105182752B (zh) 一种稀丙酮精馏工业动态优化控制层输出约束的快速设计方法
CN2804059Y (zh) 高温浓硫酸稀释器
CN105679384B (zh) 小型核电站
CN104897364A (zh) 一种水平及微倾斜管内气液两相水动力段塞流的判别方法
CN107179275B (zh) 一种注气井选材界限确定方法
CN111079367B (zh) 适用于不同进水锥管的轴流泵辐条优化设计方法
Zeng et al. A transfer-learning fault diagnosis method considering nearest neighbor feature constraints
CN111680405A (zh) 一种自然循环能力水力特性计算方法
CN106884193A (zh) 合金‑塑胶复合壳体的制备方法
CN115600518B (zh) 基于相分布特征的燃料组件热工性能分析方法及系统
WO2024068650A3 (de) Hydrodynamisches elektrolyseverfahren

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