CN111553022A - 瞬变外力对动态自反馈条件下堆芯dnbr限值修正系数的整体式求解方法 - Google Patents

瞬变外力对动态自反馈条件下堆芯dnbr限值修正系数的整体式求解方法 Download PDF

Info

Publication number
CN111553022A
CN111553022A CN202010362512.8A CN202010362512A CN111553022A CN 111553022 A CN111553022 A CN 111553022A CN 202010362512 A CN202010362512 A CN 202010362512A CN 111553022 A CN111553022 A CN 111553022A
Authority
CN
China
Prior art keywords
chf
under
transient
external force
correction coefficient
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
CN202010362512.8A
Other languages
English (en)
Other versions
CN111553022B (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 CN202010362512.8A priority Critical patent/CN111553022B/zh
Publication of CN111553022A publication Critical patent/CN111553022A/zh
Application granted granted Critical
Publication of CN111553022B publication Critical patent/CN111553022B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明公开了瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,通过乘积法求解耦合瞬变外力对反应堆动态自反馈条件下通道内CHF的影响因子,最终获得瞬变外力对动态自反馈条件下DNBR限值修正系数;所述CHF的影响因子为相应瞬变外力下的CHF值与静止条件下的CHF值的比值;所述瞬变外力包括倾斜、起伏、摇摆和耦合运动。本发明提供的方法可简单、直观、准确的获得瞬变外力场对动态自反馈条件下堆芯通道内DNBR限值修正系数,该修正系数可直接应用于核动力舰艇反应堆的设计。

Description

瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体 式求解方法
技术领域
本发明涉及核反应堆技术领域,具体涉及一种瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法。
背景技术
为了提高核动力舰船的隐蔽性,自然循环运行工况是目前核动力舰船研究的重点。另外,核动力反应堆堆芯通道内的临界热流密度(Critical Heat Flux,CHF)关系到反应堆的安全,通过通道内偏离泡核沸腾比(Departure from Nucleate Boiling Ratio,DNBR)限值来表示通道内发生CHF的裕量,在路基反应堆中DNBR限值是反应堆堆芯设计的准则,当DNBR限值应用到核动力舰船反应堆中时由于瞬变外力的影响需要对其进行修正,其修正值为瞬变外力对CHF的影响因子的最小值。核动力舰船在工作过程中会受海浪等影响发生倾斜、起伏和摇摆等运动,而倾斜条件对CHF的影响尤为重要。
目前国内DNBR限值修正系数是参考日本isshi的研究以及国内相关研究选取的,此值是在强迫循环条件下获得,强迫循环条件下瞬变外力对CHF的影响远小于自然循环的影响,因此,此值不能直接应用到动态自反馈条件下的堆芯设计。
经专利查新,目前国内外没有针对瞬变外力对动态自反馈条件下CHF影响因子的分析方法,因此,有必要研究一种瞬变外力对动态自反馈条件下CHF影响因子的分析方法,清晰、准确的表达瞬变外力对动态自反馈条件下通道内CHF的影响,能够应用于瞬变外力动态自反馈条件下DNBR限值的修正,从而在此设计准则条件下的对象能够运行于动态自反馈条件下,从而提高舰艇的隐蔽性和作战能力。
发明内容
针对上述技术问题,本发明提供了解决上述问题的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,通过考虑各瞬变外力对堆芯CHF的影响,对瞬变外力动态自反馈条件下的DNBR限值进行修正,通过此分析方法得到的DNBR限值的修正系数能够直接应用于核动力舰艇反应堆的设计。
本发明通过下述技术方案实现:
一种瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,通过乘积法求解耦合瞬变外力对反应堆动态自反馈条件下通道内CHF的影响因子,最终获得瞬变外力对动态自反馈条件下DNBR限值修正系数;所述CHF的影响因子为相应瞬变外力下的CHF值与静止条件下的CHF值的比值;所述瞬变外力包括倾斜、起伏、摇摆和耦合运动。
优选地,所述瞬变外力对动态自反馈条件下DNBR限值修正系数的计算公式如下式所示:
K=KQX×KSQ×KYB
其中,K表示瞬变外力对动态自反馈条件下DNBR限值修正系数;
KQX表示倾斜瞬变外力下CHF的最小影响因子;
KSQ表示瞬变外力下CHF的最小影响因子;
KYB表示摇摆瞬变外力下CHF的最小影响因子。
优选地,所述倾斜瞬变外力下CHF的影响因子通过以下方法获得:
通过考察倾斜角度、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内倾斜瞬变外力对CHF的最小影响因子;CHF的影响因子kQX的求解方法为:
Figure BDA0002475597530000021
其中,qCHF,QX表示倾斜瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
优选地,所述起伏瞬变外力下CHF的影响因子通过以下方法获得:
通过考察起伏周期、起伏幅值、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内起伏瞬变外力对CHF的最小影响因子;CHF的影响因子kSQ的求解方法为:
Figure BDA0002475597530000022
其中,qCHF,SQ表示起伏瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
优选地,所述摇摆瞬变外力下CHF的影响因子通过以下方法获得:
通过考察摇摆周期、摇摆角度、摇摆幅值、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内摇摆瞬变外力对CHF的最小影响因子;CHF的影响因子kYB的求解方法为:
Figure BDA0002475597530000023
其中,qCHF,YB表示倾斜瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
优选地,上述瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,是基于以下试验系统进行计算:建造自然循环试验回路,将自然循环试验回路系统搭建在模拟海洋运动的平台上,然后开展静止条件和瞬变外力下的沸腾临界实验。
本发明具有如下的优点和有益效果:
1、本发明提供的方法可简单、直观、准确的获得瞬变外力场对动态自反馈条件下堆芯通道内DNBR限值修正系数,该修正系数可直接应用于核动力舰艇反应堆的设计,利于提高舰艇的隐蔽性和作战能力。
2、本发明提供了一种适用于自然循环条件下的瞬变外力对CHF影响的整体式分析方法,填补瞬变外力对动态自反馈条件下通道内CHF影响因子的分析方法的空白,并直观、准确的获得瞬变外力对动态自反馈条件下CHF的定量影响,使得通过此分析方法得到的DNBR限值的修正系数能够应用于核动力舰艇反应堆堆芯设计所需的准则中。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例
本实施例提供了一种瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,具体步骤如下所示
步骤1,设计自然循环实验回路:
保证自然循环实验回路的运行参数能够覆盖原型装置的动态自反馈条件下的热工参数,同时保证在相同的还原条件下自然循环实验回路流量的波动幅值大于原型的流量波动,且原型的稳定性优于自然循环实验回路。
步骤2,搭建实验台架:
将自然循环实验回路搭建在模拟海洋运动平台上,开展沸腾临界实验研究;
步骤3,研究倾斜瞬变外力的影响:
首先研究倾斜条件对沸腾临界的影响规律,获得倾斜角度、通道入口温度、压力对CHF的影响规律,最终获得实验参数范围内倾斜对CHF的最小影响因子。倾斜瞬变外力对CHF的影响因子的求解方法为:瞬变外力下的CHF值与对应静止条件下的CHF的比值,如下式所示;其中对应静止条件是指回路上边界条件不做改变只对运动平台进行操作。倾斜条件是在静止条件下保持回路的温度、压力等边界条件不变,只改变瞬变外力参数。
Figure BDA0002475597530000041
其中,qCHF,QX表示倾斜瞬变外力下的CHF值;qCHF,JZ表示静止条件下的CHF值。
步骤4,研究起伏条件的影响:
然后研究起伏条件对沸腾临界的影响规律,获得起伏周期、起伏幅值、通道入口温度、压力对CHF的影响规律,最终获得实验参数范围内起伏对CHF的最小影响因子。起伏瞬变外力对CHF的影响因子的求解方法为:起伏条件下的CHF值与对应静止条件下的CHF的比值,如下式所示:
Figure BDA0002475597530000042
其中,qCHF,SQ表示起伏瞬变外力下的CHF值;qCHF,JZ表示静止条件下的CHF值。
步骤5,研究摇摆条件的影响:
最后研究摇摆条件对沸腾临界的影响规律,获得摇摆周期、摇摆角度、摇摆幅值、通道入口温度、压力对CHF的影响规律,获得实验参数范围内摇摆对CHF的最小影响因子。摇摆瞬变外力对CHF的影响因子的求解方法为:摇摆条件下的CHF值与对应静止条件下的CHF的比值,如下式所示:
Figure BDA0002475597530000043
其中,qCHF,YB表示倾斜瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
步骤6,计算DNBR限值修正系数:
通过乘积法求解耦合瞬变外力对动态自反馈条件下通道内CHF的最小影响因子,即瞬变外力动态自反馈条件下DNBR限值修正系数,如下式所示:
K=KQX×KSQ×KYB
其中,K表示瞬变外力对动态自反馈条件下DNBR限值修正系数;
KQX表示倾斜瞬变外力下CHF的最小影响因子;
KSQ表示瞬变外力下CHF的最小影响因子;
KYB表示摇摆瞬变外力下CHF的最小影响因子。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,通过乘积法求解耦合瞬变外力对反应堆动态自反馈条件下通道内CHF的影响因子,最终获得瞬变外力对动态自反馈条件下DNBR限值修正系数;
所述CHF的影响因子为相应瞬变外力下的CHF值与静止条件下的CHF值的比值;
所述瞬变外力包括倾斜、起伏、摇摆和耦合运动。
2.根据权利要求1所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,所述瞬变外力对动态自反馈条件下DNBR限值修正系数的计算公式如下式所示:
K=KQX×KSQ×KYB
其中,K表示瞬变外力对动态自反馈条件下DNBR限值修正系数;
KQX表示倾斜瞬变外力下CHF的最小影响因子;
KSQ表示瞬变外力下CHF的最小影响因子;
KYB表示摇摆瞬变外力下CHF的最小影响因子。
3.根据权利要求2所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,所述倾斜瞬变外力下CHF的影响因子通过以下方法获得:
通过考察倾斜角度、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内倾斜瞬变外力对CHF的最小影响因子;CHF的影响因子kQX的求解方法为:
Figure FDA0002475597520000011
其中,qCHF,QX表示倾斜瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
4.根据权利要求2所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,所述起伏瞬变外力下CHF的影响因子通过以下方法获得:
通过考察起伏周期、起伏幅值、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内起伏瞬变外力对CHF的最小影响因子;CHF的影响因子kSQ的求解方法为:
Figure FDA0002475597520000012
其中,qCHF,SQ表示起伏瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
5.根据权利要求2所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,所述摇摆瞬变外力下CHF的影响因子通过以下方法获得:
通过考察摇摆周期、摇摆角度、摇摆幅值、通道入口温度、压力对CHF值的影响规律,获得实验参数范围内摇摆瞬变外力对CHF的最小影响因子;CHF的影响因子kYB的求解方法为:
Figure FDA0002475597520000021
其中,qCHF,YB表示倾斜瞬变外力下的CHF值;
qCHF,JZ表示静止条件下的CHF值。
6.根据权利要求1至5任一项所述的瞬变外力对动态自反馈条件下堆芯DNBR限值修正系数的整体式求解方法,其特征在于,基于以下试验系统进行计算:建造自然循环试验回路,将自然循环试验回路系统搭建在模拟海洋运动的平台上,然后开展静止条件和瞬变外力下的沸腾临界实验。
CN202010362512.8A 2020-04-30 2020-04-30 动态自反馈条件下堆芯dnbr限值修正系数的求解方法 Active CN111553022B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010362512.8A CN111553022B (zh) 2020-04-30 2020-04-30 动态自反馈条件下堆芯dnbr限值修正系数的求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010362512.8A CN111553022B (zh) 2020-04-30 2020-04-30 动态自反馈条件下堆芯dnbr限值修正系数的求解方法

Publications (2)

Publication Number Publication Date
CN111553022A true CN111553022A (zh) 2020-08-18
CN111553022B CN111553022B (zh) 2020-12-08

Family

ID=72007880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010362512.8A Active CN111553022B (zh) 2020-04-30 2020-04-30 动态自反馈条件下堆芯dnbr限值修正系数的求解方法

Country Status (1)

Country Link
CN (1) CN111553022B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112507631A (zh) * 2020-11-19 2021-03-16 中国核动力研究设计院 一种窄通道流动失稳出口界限含气率限值测试方法及系统
CN115688488A (zh) * 2022-12-30 2023-02-03 中国核动力研究设计院 准则确定方法、装置、设备、计算机存储介质及程序产品
CN115862909A (zh) * 2023-02-14 2023-03-28 中国核动力研究设计院 堆芯热工安全准则应用方法、装置、设备、介质及产品
CN116884655A (zh) * 2023-09-08 2023-10-13 中国核动力研究设计院 外力场对热工安全影响确定方法、装置、核反应堆及设备
CN117252045A (zh) * 2023-11-17 2023-12-19 中国核动力研究设计院 反应堆堆芯参数的设计方法、设备、介质及核能动力系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4318778A (en) * 1973-05-22 1982-03-09 Combustion Engineering, Inc. Method and apparatus for controlling a nuclear reactor
CN109192342A (zh) * 2018-09-12 2019-01-11 中国核动力研究设计院 一种用于低流速chf实验的临界热流密度近似测量方法
CN110134906A (zh) * 2019-03-06 2019-08-16 中国核电工程有限公司 一种非对称功率分布下堆芯dnbr计算方法
CN110289111A (zh) * 2019-07-09 2019-09-27 上海交通大学 自然循环回路的流量确定方法、装置及服务器
CN110322976A (zh) * 2019-08-06 2019-10-11 中国核动力研究设计院 一种用于反应堆lpd和dnbr在线保护和监测的实现方法
CN110633454A (zh) * 2019-09-19 2019-12-31 中国核动力研究设计院 一种基于修正法的chf关系式dnbr限值统计学确定方法
CN110727920A (zh) * 2019-09-19 2020-01-24 中国核动力研究设计院 一种基于分组法的chf关系式dnbr限值统计学确定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4318778A (en) * 1973-05-22 1982-03-09 Combustion Engineering, Inc. Method and apparatus for controlling a nuclear reactor
CN109192342A (zh) * 2018-09-12 2019-01-11 中国核动力研究设计院 一种用于低流速chf实验的临界热流密度近似测量方法
CN110134906A (zh) * 2019-03-06 2019-08-16 中国核电工程有限公司 一种非对称功率分布下堆芯dnbr计算方法
CN110289111A (zh) * 2019-07-09 2019-09-27 上海交通大学 自然循环回路的流量确定方法、装置及服务器
CN110322976A (zh) * 2019-08-06 2019-10-11 中国核动力研究设计院 一种用于反应堆lpd和dnbr在线保护和监测的实现方法
CN110633454A (zh) * 2019-09-19 2019-12-31 中国核动力研究设计院 一种基于修正法的chf关系式dnbr限值统计学确定方法
CN110727920A (zh) * 2019-09-19 2020-01-24 中国核动力研究设计院 一种基于分组法的chf关系式dnbr限值统计学确定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
EVANS DAMENORTEY KITCHER: "PERFORMANCE AND SAFETY ANALYSIS OF A GENERIC SMALL MODULAR REACTOR", 《CORE.AC.UK》 *
RONG CAI 等: "Development of a thermal-hydraulic subchannel analysis code for motion conditions", 《PROGRESS IN NUCLEAR ENERGY》 *
张玉相 等: "CHF关系式开发与DNBR限值确定方法研究", 《核动力工程》 *
陈立新: "脉冲反应堆过冷沸腾传热特性与堆芯空泡份额计算方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112507631A (zh) * 2020-11-19 2021-03-16 中国核动力研究设计院 一种窄通道流动失稳出口界限含气率限值测试方法及系统
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 中国核动力研究设计院 外力场对热工安全影响确定方法、装置、核反应堆及设备
CN117252045A (zh) * 2023-11-17 2023-12-19 中国核动力研究设计院 反应堆堆芯参数的设计方法、设备、介质及核能动力系统
CN117252045B (zh) * 2023-11-17 2024-01-26 中国核动力研究设计院 反应堆堆芯参数的设计方法、设备、介质及核能动力系统

Also Published As

Publication number Publication date
CN111553022B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN111553022B (zh) 动态自反馈条件下堆芯dnbr限值修正系数的求解方法
Korde et al. Hydrodynamic control of wave energy devices
Faraggiana et al. A review of numerical modelling and optimisation of the floating support structure for offshore wind turbines
Wang et al. Stochastic dynamic response analysis of a floating vertical‐axis wind turbine with a semi‐submersible floater
CN111581806B (zh) 瞬变外力对动态自反馈条件下通道内chf的影响分析方法
Antonutti et al. An investigation of the effects of wind-induced inclination on floating wind turbine dynamics: heave plate excursion
Spanos et al. Efficient dynamic analysis of a nonlinear wave energy harvester model
Lupton Frequency-domain modelling of floating wind turbines
CN111507016B (zh) 一种运动条件动态下并联窄通道流动失稳边界的确定方法
Matha et al. Efficient critical design load case identification for floating offshore wind turbines with a reduced nonlinear model
Peiffer et al. Design of an oscillating wave surge converter on the windfloat structure
Zhang et al. Novel method for designing and optimising the floating platforms of offshore wind turbines
Alkhabbaz et al. The aero-hydrodynamic interference impact on the NREL 5-MW floating wind turbine experiencing surge motion
Zuo et al. Computationally inexpensive approach for pitch control of offshore wind turbine on barge floating platform
Deng et al. Frequency-dependent aerodynamic damping and its effects on dynamic responses of floating offshore wind turbines
Li Coupled computational fluid dynamics/multibody dynamics method with application to wind turbine simulations
CN108153976B (zh) 一种海浪中舰船横摇运动的仿真方法及上位机
CN115859523A (zh) 一种作动器力控加载浮式风机混合实验系统及方法
Yu et al. Model test and validation of the CROWN floating offshore wind turbine
CN113779698A (zh) 一种水流作用下的船舶系泊系统简化设计方法
CN111368474A (zh) 海流流速预报方法、装置及电子设备
Chen et al. Systematic comparisons among OpenFAST, Charm3D-FAST simulations and DeepCWind model test for 5 MW OC4 semisubmersible offshore wind turbine
Forbush et al. Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code
Zhao et al. Numerical investigation on the effects of current headings on vortex induced motions of a semi-submersible
Ding¹ et al. Influence of Mooring Clump Weight on the Dynamic Responses of a Floating Offshore

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