根据临界热流密度Benchmark本体试验，研究定位格架间距对临界热流密度大小和发生点以及发生临界通道的热工水力参数的影响。分析结果表明，随着定位格架间距的增大，临界点向下一个定位格架上游移动，临界热流密度增大，发生临界通道的出口质量流速减小。 Based on the ontology test of Benchmark, it was studied that how the positioning lattice spacing influenced on the critical heat flux and its occurrence point as well as the thermal hydraulic pa-rameters of the critical passage. The results show that as the distance between the positioning lattice increases, the critical point moves upstream of the next positioning lattice, the critical heat flux increases, and the exit mass velocity of the critical passage decreases.
赖诚，白金川，金聪，徐樟楠，袁亮
苏州热工研究院有限公司，广东 深圳
收稿日期：2018年9月28日；录用日期：2018年10月13日；发布日期：2018年10月20日
根据临界热流密度Benchmark本体试验，研究定位格架间距对临界热流密度大小和发生点以及发生临界通道的热工水力参数的影响。分析结果表明，随着定位格架间距的增大，临界点向下一个定位格架上游移动，临界热流密度增大，发生临界通道的出口质量流速减小。
关键词 :Benchmark试验，定位格架，临界热流密度
Copyright © 2018 by authors and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
定位格架是压水堆燃料组件的关键部件之一，它除了在结构上起支撑和定位之外，还能改善反应堆内流体的热工水力特性，提高反应堆的安全性。燃料组件的临界热流密度(CHF)综合反映了反应堆燃料组件内冷却剂的流动和传热特性，直接影响到组件的安全裕量和最小临界热流密度比(MDNBR)。由于定位格架交混长度 [
本文在临界热流密度Benchmark 5 × 5试验本体的基础上，对沿轴向方向等间距设置13个定位格架的工况进行分析。通过改变定位格架间距，计算分析定位格架间距对MDNBR点的位置以及发生临界的通道热工水力参数的影响。
在压水堆燃料组件中，定位格架沿轴向长度等间距排列。定位格架的存在提高了冷却剂的传热系数，但同时增加了当地压降。如图1、表1所示，基于临界热流密度BENCHMARK 5 × 5实验本体结构 [
图1. 5 × 5组件子通道划分
| 组件类型 |
|
|||||
|---|---|---|---|---|---|---|
| 加热棒数 | 25 | |||||
| 棒间距(mm) | 12.6 | |||||
| 轴向加热长度(mm) | 3658 | |||||
| 流道内宽(mm) | 64.9 | |||||
| 径向功率分布 |
|
|||||
| 轴向功率分布 | 余弦分布 | |||||
| 定位格架数 | 13 | |||||
| 定位格架间距(mm) | 254 | 264 | 274 | 284 | 294 | 304 |
表1. 5 × 5格架参数
将图1所示的燃料组件截面示意图划分为四种典型类型的子通道：角子通道、边子通道和两种中心通道，如图2所示。沿轴向将燃料组件划分为100个节点，不考虑流体所受定位格架的流动压降，只考虑其形阻压降。燃料组件的轴向功率分布采用余弦分布，径向功率分布如图3所示，其中心燃料棒相对功率份额与边通道燃料棒相对功率份额之比C/P为1:0.85。
为了便于采用子通道概念分析燃料组件中的冷却剂流动机理，本文作以下假设：
1) 沸腾时，各子通道内的流体为一维两相可分离的滑移流模型。
图2. 子通道分布图
图3. 径向功率分布示意图
2) 空泡份额为焓，压力，质量流率，轴向位置和时间的函数。
3) 相邻子通道间的紊流交混不会引起流量的再分配。
4) 紊流交混可以叠加在能引起流量再分配的横流上。
5) 忽略声速传递的影响。
6) 在子通道内，横流速度远小于轴向速度。
基于以上假设和流动的三大定律，可以得到燃料组件的冷却剂流动的守恒方程 [
连续性方程：
A i ∂ ρ m , i ∂ t + ∂ m i ∂ x = − ∑ j ( w i j + w ′ i j )
能量守恒方程：
m i u ″ ∂ h i ∂ t + m i ∂ h i ∂ x = q ′ − ∑ j w ′ i j ( h i − h j ) − ∑ j c ( t i − t j ) + h i ∑ j w i j − ∑ j h * w i j
轴向动量守恒方程：
1 A i ∂ m i ∂ t − 2 u i ∂ ρ i ∂ t + ∂ ρ ∂ x = a ′ + 2 u i ∑ j w i j − ∑ j u * w i j A i
a ′ = ( m i A i ) 2 ( v i f ϕ 2 D + K v ′ i 2 Δ x + A i ∂ ( v ′ i / A i ) ∂ x + ρ i g cos θ ) − f T ∑ j w ′ i j ( u i − u j ) A i
横向动量守恒方程：
∂ w i j ∂ t + ∂ w i j u i ∂ x = s l c w i j + s l ( P i − P j )
子通道程序的计算流程图如图4所示。
本文选用了两个参数范围较广且广泛应用于压水堆安全分析的临界热流密关系式B&W-2和W-3关系式来预测临界热流密度Benchmark 5 × 5试验，并比较两关系式与试验的结果，然后分析定位格架间距对临界发生点位置以及热工水力参数的影响。表2列出了B&W-2 [
图4. 子通道程序计算流程图
| 关系式 | 表达式 | 数据基础 |
|---|---|---|
| B&W-2关系式 | q ″ C H F , E U = 3.155 × 10 − 3 ( 1.155 − 16.025 D e ) × [ 3.702 × 10 7 ( 4.3604 × 10 − 4 G ) B − 48.21 G H f g ] / [ 12.71 ( 2.252 × 10 − 3 G ) A ] A = 0.712 + 3.006 × 10 − 5 ( p − 13790 ) B = 0.834 + 9.93 × 10 − 5 ( p − 13790 ) | L h = 1.83 m D e = 0.0051 ~ 0.0127 m p = 13800 ~ 16550 kPa G = 1020 ~ 5425 kg / ( m 2 ⋅ s ) x = − 0.03 ~ 0.2 |
| F c = q ″ C H F , E U q ″ C H F , N U = 1.025 C ∫ 0 z D N B , N q ″ ( z ) exp [ − C ( l C H F − z ) ] d z q ( z D N B , N ) [ 1 − exp ( − C z D N B , E U ) ] C = 9.8 ( 1 − X D N B ) 7.82 ( G / 10 6 ) 0.457 | ||
| W-3关系式 | q ″ C H F , E U = 3.14 × 10 6 { ( 2.022 − 6.238 × 10 − 8 p ) + ( 0.1722 − 1.43 × 10 − 8 p ) × exp [ ( 18.177 − 5.987 × 10 − 7 p ) x e ] } × [ ( 0.1484 − 1.596 x e + 0.1729 x e | x e | ) × 0.2049 G 10 6 + 1.037 ] × ( 1.157 − 0.869 x e ) [ 0.2664 + 0.8357 exp ( − 124 D e ) ] × [ 0.8258 + 0.341 × 10 − 6 ( h f s − h f , i n ) ] | L h = 0.254 ~ 3.668 m D e = 0.00508 ~ 0.0178 m p = 6895 ~ 15860 kPa G = ( 4.9 ~ 24.5 ) × 10 6 kg / ( m 2 ⋅ h ) x = − 0.03 ~ 0.2 |
| F c = q ″ C H F , E U q ″ C H F , N U = C ∫ 0 z D N B , N q ″ ( z ) exp [ − C ( l C H F − z ) ] d z q ( z D N B , N ) [ 1 − exp ( − C z D N B , E U ) ] C = 12.64 ( 1 − X D N B ) 4.31 ( G / 10 6 ) 0.478 |
表2. 临界热流密度公式
本文通过临界热流密度关系式的比较，选用合适的计算关系式分析定位格架间距对5 × 5燃料组件的临界发生点位置和该点热工水力参数的影响。
图5给出了临界热流密度Benchmark 5 × 5试验的数据与B&W-2和W-3公式所计算的CHF值的对比。从图3可以看到在临界热流密度Benchmark 5 × 5的试验中W-3公式与试验结果符合较好，而B&W-2在某些工况下的计算结果不可靠。
在临界热流密度Benchmark 5 × 5试验的条件下，由图5确定采用W-3关系式分析沿轴向长度分布13个定位格架的工况进行定位格架间距对组件临界热流密度发生点和发生临界的通道热工水力参数的影响。
由图6给出了沿轴向分布13个定位格架工况在改变定位格架间距时，组件临界热流密度发生点和临界热流密度大小的模拟结果。图6显示当间距比较小时，临界发生点位于定位格架下游处且临界热流密度较小，随着定位格架间距的变大，轴向流体交混越来越充分临界发生点向下一个定位格架上游靠近且临界热流密度大致保持不变。
图7是发生临界的通道质量流速随定位格架位置的变化结果。图8可以看出临界通道的质量流速沿轴向高度逐渐减少，随着格架间距的增大，出口质量流速逐渐变小。结合图9，临界通道空泡份额沿轴向的变化可以知道，由于空泡产生前，通道中冷却剂由单相变为两相，摩擦系数将增加，这时出现一定的交混，因此在空泡起始处质量流密度发生突变。
图5. 临界热流密度关系式计算值与试验值对比图
图6. 格架间距变化下的MDNBR分布
图7. 格架间距变化下的临界通道轴向质量流速
图8. 格架间距变化下临界通道出口质量流速
图9. 格架间距变化下临界通道空泡份额轴向分布
本文以临界热流密度Benchmark试验本体为研究对象，首先比较临界热流密度关系式B&W-2和W-3，并选用W-3关系式分析定位格架间距对5 × 5燃料组件临界热流密度的影响，结论如下：
1) 定位格架的存在加强了流体的交混。随着定位格架间距的增大，交混越充分，使得临界热流密度变大，同时临界点向下一个定位格架的上游移动。
2) 当定位格架的间距大于定位格架的交混长度时，组件临界热流密度和临界点的位置基本保持不变。
3) 随着定位格架间距的增大，发生临界的通道的质量流速沿轴线高度方向逐渐减小，且出口质量流速也相应变小。
赖 诚,白金川,金 聪,徐樟楠,袁 亮. 定位格架对燃料棒热工水力特性影响的研究Study on the Influence of Positioning Lattice on Thermal Hydraulic Characteristics of Fuel Rod[J]. 核科学与技术, 2018, 06(04): 99-107. https://doi.org/10.12677/NST.2018.64012