CN115935625A - 一种堆芯最小流量估算方法 - Google Patents
一种堆芯最小流量估算方法 Download PDFInfo
- Publication number
- CN115935625A CN115935625A CN202211489328.5A CN202211489328A CN115935625A CN 115935625 A CN115935625 A CN 115935625A CN 202211489328 A CN202211489328 A CN 202211489328A CN 115935625 A CN115935625 A CN 115935625A
- Authority
- CN
- China
- Prior art keywords
- core
- flow
- reactor core
- outlet
- total
- 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.)
- Pending
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种堆芯最小流量估算方法。该方法步骤如下:1、输入参数:棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数。2、计算出通道总功率。3、给定初始堆芯总流量,确定通道最优流量分配方案。4、进行组件出口过冷度计算。5、判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内,若否,则进行堆芯总流量修正,返回步骤4;若是,则进行下一个步骤计算。6、堆芯温度分布计算。7、输出修正后的堆芯最小流量及最小流量下的堆芯温度分布。该方法优点如下:1、采用通道最优流量分配方案。2、对堆芯总流量修正构造修正量的关于迭代步数的函数,提高迭代运算的计算速度。
Description
技术领域
本发明属于反应堆热工水力技术领域,具体涉及一种堆芯最小流量估算方法。
背景技术
采用带盒燃料组件的紧凑型压水堆,应用流量分区设计能够显著展平堆芯出口温度分布,提高出口过冷度,或降低冷却剂流量,从而更加有效地利用冷却剂,大幅提升堆芯设计性能。寿期内堆芯三维功率分布形状会发生较大变化,为了合理分配流量,达到组件出口冷却剂温度展平、过冷度提高,通过研究流量分区依据并设计考虑多种因素的优化方法,设计闭式通道堆芯流量分区及水密度反馈程序QSubTH,并使程序能够满足以下计算功能:
快速计算闭式通道燃料组件轴向水密度分布、芯块燃料温度分布、包壳温度分布、组件出口温度及过冷度;
在给定堆芯热功率及三维分布、最小出口过冷度条件下,估算堆芯最小冷却剂质量流量及最大出口温度。
发明内容
为了解决上述问题,本发明提供一种堆芯最小流量估算方法。在给定堆芯热功率及三维分布、最小出口过冷度条件下,估算堆芯最小冷却剂质量流量及最大出口温度。
为了达到上述目的,本发明的采用如下技术方案:
一种堆芯最小流量估算方法,包括如下步骤:
一种堆芯最小流量估算方法,包括如下步骤:
步骤1:输入参数:包括棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数,其中堆芯功率分布参数及组件不均匀系数由堆芯物理计算后提供相应数据,包括组件所在堆芯位置、组件功率分布和组件不均匀系数;步骤2:使用堆芯功率分布参数及组件不均匀系数,对每个通道轴向功率求和,计算出通道总功率;
步骤3:给定初始堆芯总流量,结合堆芯入口温度以及堆芯总功率,由能量守恒方程计算堆芯出口温度一致时的焓值,由各个通道的总功率以及计算出的焓值确定通道最优流量分配方案;
步骤4:进行组件出口过冷度计算,在计算结果中搜索组件出口最小过冷度;
步骤5:判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内,若否,则进行堆芯总流量修正,返回步骤4重新进行组件出口过冷度计算;若是,则进行下一个步骤计算;
步骤6:堆芯温度分布计算:使用棒状或板状燃料元件堆芯参数、堆芯计算控制变量参数,通过能量守恒方程确定各组件中冷却剂温度,再由换热系数关系式以及材料的导热系数,通过傅里叶定律,牛顿冷却公式由外向内依次计算包壳外层、包壳内层、燃料表面、燃料中心温度的轴向分布情况;
步骤7:输出修正后的堆芯最小总流量及最小总流量下的堆芯温度分布。
1、本发明与现有技术相比,具有如下优点:步骤1中输入参数包括组件所在堆芯位置、组件功率分布和组件不均匀系数,将堆芯寿期内功率变化考虑入内;输入参数中堆芯流量控制变量提高流量优化效率。
2、步骤3中采用最优流量分配方案,是对于单一功率分布形状,当各个通道出口处的温度相等时,可以得到组件最大出口过冷度。由用户提供的堆芯总流量,入口温度以及堆芯总功率,由能量守恒方程,可以确定出口温度一致时的焓值,则由各个通道的功率以及计算出的出口焓值确定通道最优流量。
3、步骤5中的允许范围,由用户提供,程序中也提供了自定义接口,方便用户自由控制,建议值为10-7;步骤5中的堆芯总流量修正,是对堆芯总流量给出带符号的修正量,在原有修正量上添加不同量级的比例系数,构造修正量的关于迭代步数的函数,提高迭代运算的计算速度。
附图说明
图1是本发明方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
如图1所示,本发明一种堆芯最小流量估算方法,包括如下步骤:
步骤1:输入参数:包括棒状或板状燃料元件堆芯参数(包括堆芯旁流系数、包壳与燃料芯块之间的间隙导热系数、相邻燃料棒中心距、堆芯活性区高度、堆芯最小流量估算流量迭代系数、组件中导向管数目、堆芯组件数目、轴向划分控制体数目、计算区域中组件数目、划分计算区域数目、组件中燃料棒数、燃耗点数、堆芯压力、堆芯总功率、燃料棒包壳内半径、燃料棒包壳内外径、燃料芯块半径、导向管半径、组件出口最小过冷度设计值、堆芯入口冷却剂温度、通道冷却剂流量、方型组件宽度和堆芯冷却剂净流量,若为板状燃料,还需要输入板元件长度、燃料板数目、包壳厚度、板元件厚度和板元件宽度)、堆芯功率分布参数及组件不均匀系数(由堆芯物理计算后提供相应数据,包括组件所在堆芯位置、组件功率分布和组件不均匀系数)、堆芯计算控制变量参数(包括堆芯流量控制变量、组件功率分布、不同流量分配下堆芯出口温度比较计算控制变量、燃料类型控制变量、计算类型控制变量和参考流量输入控制变量(用于不同流量比较计算))。
步骤2:使用输入参数,对每个通道轴向功率求和,计算出通道总功率。
步骤3:给定初始堆芯总流量,结合堆芯入口温度以及堆芯总功率,由能量守恒方程计算堆芯出口温度一致时的焓值,由各个通道的堆芯功率以及计算出的焓值确定通道最优流量分配方案。对计算得到的每套相应燃耗点下最优分配方案进行搜索,确定每个组件在整个寿期中所需的最大流量,为了保证堆芯总流量守恒,以确定的各组件最大流量确定初始流量分配比例因子。
对于单一功率分布形状,当各个通道出口处的温度相等时,可以得到组件最大出口过冷度。由已知的堆芯总流量,入口温度以及堆芯总功率,由能量守恒方程,可以确定出口温度一致时的组件出口焓值:
上式中:
hout——流量分配最优时,组件出口焓值;
hin——冷却剂在入口处的焓值;
Ptot——堆芯总功率;
Wtot——堆芯总冷却剂流量。
则由各个通道的功率以及计算出的组件出口焓值确定通道最优流量:
上式中:
W(i)——通道i的对应的最佳分配流量;
P(i)——通道i的对应的总功率。
在整个寿期内,堆芯功率分布会随着燃耗加深而有所改变,为了得到整个寿期内的最优流量分配,需要将功率分布的时间效应考虑进去。在程序中,以组件最大出口温度作为衡量流量分配方案好坏的标准。
对于每个燃耗点下,都有一套最优流量分配方案,每套分配方案相对独立,没有将功率分布的时间效应考虑进去,在确定初始分配方案时,首先对计算得到的每套相应燃耗点下最优分配方案进行搜索,确定每个组件在整个寿期中所需的最大流量,为了保证堆芯总流量守恒,以确定的各组件最大流量确定初始流量分配比例因子:
上式中:
f——流量分配比例因子;
Wk(i)——组件i在寿期内为确保最小出口温度所需的最大流量。
再用总量按比例因子分配各通道的初始流量。这样可以初步将功率分布的时间效应考虑在流量分配中。
步骤4:进行组件出口过冷度计算,在计算结果中搜索组件出口最小过冷度。
冷却剂换热,通过采用合适的经验关系式来进行计算。由于要求堆芯出口存在一定的过冷度,冷却剂处于单相液态,其换热模型采用单相液体换热关系式,根据流量的不同分为高、中、低流量率部分:
(1)高流量率(Re>2500)
Dittus-Boelter关系式
Sieder-Tate关系式
Mihaiyeph关系式
(2)低流量率(Re<1800)
Collier关系式
以上关系式中:
h——换热系数,W/(m2·K);
λc——导热率,W/(m·K);
De——流道当量直径,m;
Re——雷诺数;
Pr——普朗特数;
Prst-water——饱和水普朗特数;
Gr——格拉晓夫数;
g——重力加速度,m·s-2;
μ——动力黏性系数,N·s/m2;
μf——流体动力黏性系数,N·s/m2;
μwall——壁面动力黏性系数,N·s/m2;
ρ——密度,kg/m3;
Twall——壁面温度,K;
Tf——流体温度,K;
β——体积变化系数。
(3)中流量率(1800≤Re≤2500)
在上述两种情况的端点值下进行线性插值,得到对应雷诺数下的换热系数。
包壳导热相关计算,根据燃料几何特点,选用不同的稳态导热计算模型:
以上关系式中:
ql——单位长度上热流量,W/m;
λ——材料导热系数,W/(m·K);
t1-t2——圆筒壁内外温差,K;
d2——圆筒外径,m;
d1——圆筒内径,m;
Lwi——板宽度,m;
Lth——板厚度,m。
步骤5:判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内,若否,则进行堆芯总流量修正,返回步骤4重新进行组件出口过冷度计算;若是,则进行下一个步骤计算。
步骤6:堆芯温度分布计算:通过能量守恒方程确定各组件中冷却剂温度,再由换热系数关系式以及材料的导热系数,通过傅里叶定律,牛顿冷却公式由外向内依次计算:包壳外层、包壳内层、燃料表面、燃料中心温度的轴向分布情况。
堆芯中涉及到的传热问题主要包括:燃料芯块的导热,燃料棒间隙的导热,包壳导热以及包壳外层与冷却剂之间的换热,由于堆芯内部相邻部件的温度差距不是很大,这里忽略辐射换热。将同一组件中的燃料棒均匀化处理,认为同一组件中的燃料棒温度分布一致。
燃料棒包壳导热系数,根据Zr-4合金物性关系式确定:
以上关系式中:
λclad——燃料棒包壳导热系数,W/(m·K);
T——包壳温度,℃。
间隙导热率由外部文件输入经验值,在根据柱坐标系下一维导热公式计算间隙热流密度;燃料中心温度由积分导热率进行计算:
上式中:
tuo——半径为ruo处芯块温度,K;
tui——半径为rui处芯块温度,K;
ku——燃料芯块热导率,W/(m·K)。
步骤7:输出修正后的堆芯最小流量及最小流量下的堆芯温度分布。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,但不能认定本发明的具体实施方式仅限于此,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所提交的权利要求书确定专利保护范围。
Claims (5)
1.一种堆芯最小流量估算方法,其特征在于:包括如下步骤:
步骤1:输入参数:包括棒状或板状燃料元件堆芯参数、堆芯功率分布参数及组件不均匀系数、堆芯计算控制变量参数,其中堆芯功率分布参数及组件不均匀系数由堆芯物理计算后提供相应数据,包括组件所在堆芯位置、组件功率分布和组件不均匀系数;
步骤2:使用堆芯功率分布参数及组件不均匀系数,对每个通道轴向功率求和,计算出通道总功率;
步骤3:给定初始堆芯总流量,结合堆芯入口温度以及堆芯总功率,由能量守恒方程计算堆芯出口温度一致时的焓值,由各个通道的总功率以及计算出的焓值确定通道最优流量分配方案;
步骤4:进行组件出口过冷度计算,在计算结果中搜索组件出口最小过冷度;
步骤5:判断出口最小过冷度计算值与设计值的相对偏差是否在允许范围内,若否,则进行堆芯总流量修正,返回步骤4重新进行组件出口过冷度计算;若是,则进行下一个步骤计算;
步骤6:堆芯温度分布计算:使用棒状或板状燃料元件堆芯参数、堆芯计算控制变量参数,通过能量守恒方程确定各组件中冷却剂温度,再由换热系数关系式以及材料的导热系数,通过傅里叶定律,牛顿冷却公式由外向内依次计算包壳外层、包壳内层、燃料表面、燃料中心温度的轴向分布情况;
步骤7:输出修正后的堆芯最小总流量及最小总流量下的堆芯温度分布。
2.根据权利要求1所述的一种堆芯最小流量估算方法,其特征在于:步骤1中所述的棒状或板状燃料元件堆芯参数包括堆芯旁流系数、包壳与燃料芯块之间的间隙导热系数、相邻燃料棒中心距、堆芯活性区高度、堆芯最小流量估算流量迭代系数、组件中导向管数目、堆芯组件数目、轴向划分控制体数目、计算区域中组件数目、划分计算区域数目、组件中燃料棒数、燃耗点数、堆芯压力、堆芯总功率、燃料棒包壳内半径、燃料棒包壳内外径、燃料芯块半径、导向管半径、组件出口最小过冷度设计值、堆芯入口冷却剂温度、通道冷却剂流量、方型组件宽度和堆芯冷却剂净流量,若为板状燃料,还需要输入板元件长度、燃料板数目、包壳厚度、板元件厚度和板元件宽度;
所述堆芯计算控制变量参数包括堆芯流量控制变量、组件功率分布、不同流量分配下堆芯出口温度比较计算控制变量、燃料类型控制变量、计算类型控制变量和用于不同流量比较计算的参考流量输入控制变量。
3.根据权利要求1所述的一种堆芯最小流量估算方法,其特征在于:步骤3中所述的确定通道最优流量分配方案,是对于单一功率分布形状,当各个通道出口处的温度相等时,得到组件最大出口过冷度;由用户提供的堆芯总流量,堆芯入口温度以及堆芯总功率,由能量守恒方程,确定堆芯出口温度一致时的组件出口焓值,则由各个通道的功率以及计算出的组件出口焓值确定通道最优流量。
4.根据权利要求1所述的一种堆芯最小流量估算方法,其特征在于:步骤5中所述的允许范围为10-7。
5.根据权利要求1所述的一种堆芯最小流量估算方法,其特征在于:步骤5中所述的堆芯总流量修正,是对堆芯总流量给出带符号的修正量,在原有修正量上添加不同量级的比例系数,构造修正量的关于迭代步数的函数,提高迭代运算的计算速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211489328.5A CN115935625A (zh) | 2022-11-25 | 2022-11-25 | 一种堆芯最小流量估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211489328.5A CN115935625A (zh) | 2022-11-25 | 2022-11-25 | 一种堆芯最小流量估算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115935625A true CN115935625A (zh) | 2023-04-07 |
Family
ID=86655341
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211489328.5A Pending CN115935625A (zh) | 2022-11-25 | 2022-11-25 | 一种堆芯最小流量估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115935625A (zh) |
-
2022
- 2022-11-25 CN CN202211489328.5A patent/CN115935625A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113094947B (zh) | 一种核反应堆堆芯核热耦合分析方法 | |
CN111261232B (zh) | 反应堆一回路冷却剂流场、温度场和应力场耦合计算方法 | |
CN109657372B (zh) | 一种新型的管壳式换热器多尺度耦合协同换热模拟方法 | |
Li et al. | Similarity and generalized analysis of efficiencies of thermal energy storage systems | |
CN114282460B (zh) | 一种铅铋快堆堆芯热工水力特性分析方法 | |
CN111274748B (zh) | 池式钠冷快堆非能动余热排出系统跨维度耦合计算方法 | |
CN112699620B (zh) | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 | |
CN109783936B (zh) | 压水堆核电站核岛变工况蓄热增量的计算方法 | |
CN111027112B (zh) | 一种针对快堆棒束组件耦合传热模型的多孔介质模拟方法 | |
US5638413A (en) | Apparatus and method for estimating core performance | |
CN113255249B (zh) | 针对热管固态堆的多物理耦合瞬态计算方法及装置 | |
CN110310710B (zh) | 获得燃料元件矩形窄缝通道内干涸型临界热流密度的方法 | |
CN114444413B (zh) | 一种板状燃料堆芯亚通道级三维热工水力分析方法 | |
You et al. | A numerical study on the turbulent heat transfer enhancement of Rodbaffle heat exchanger with staggered tubes supported by round rods with arc cuts | |
Wei et al. | Fluid flow distribution optimization for minimizing the peak temperature of a tubular solar receiver | |
CN111259596B (zh) | 基于有限体积理论的管壳式热交换器全三维耦合模拟方法 | |
CN115935625A (zh) | 一种堆芯最小流量估算方法 | |
Zima et al. | Mathematical modelling of heat transfer in liquid flat-plate solar collector tubes | |
CN115862901A (zh) | 一种给定总流量下堆芯寿期内最优流量分配方法 | |
CN115112323A (zh) | 一种液态铅铋冲刷绕丝定位棒束流致振动实验装置及方法 | |
CN113065241B (zh) | 一种预测超临界二氧化碳冷却堆燃料元件主要参数的方法 | |
Hernandez et al. | Development of a CFD-based model to simulate loss of flow transients in a small lead-cooled reactor | |
Mao et al. | Numerical and experimental investigation on heat transfer performance of a solar single storage tank | |
CN106156460B (zh) | 获得核燃料棒内部燃料芯块温度分布的方法 | |
Chen et al. | Configurational explorations and optimizations of a phase change material based thermal energy storage with tube bundle |
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 |