Numerical Simulation of Critical Flow in a Narrow Rectangular Channel
-
摘要: 冷却剂通过细微裂缝泄漏的主要机理为临界流流动,准确预测临界流流量是实现破前漏(LBB)应用的关键。基于空化模型对窄矩形通道的临界流进行了数值模拟,探讨了表征气泡半径和气化核心密度的常数项(C)对于临界流的影响。结果表明,当C=1.25时模拟结果与实验值符合良好;在实验工况范围内,临界流质量流速的模拟结果与实验值的误差在±15%以内;修正后的空化模型可用于模拟计算上游温度、压力工况下的临界流质量流速以及归纳临界流流量变化规律。Abstract: The main mechanism of coolant leakage through micro-cracks is critical flow, and accurate prediction of critical flow is the key to realize the application of Leak-Before-Break (LBB). Based on the cavitation model, the critical flow in a narrow rectangular channel is numerically simulated, and the influence of the constant term (C) representing the bubble radius and nucleation site density on the critical flow is discussed. The results show that when C=1.25, the simulation results are in good agreement with the experimental values; Within the range of experimental conditions, the error between the simulation results of critical flow and the experimental values is within ± 15%; The modified cavitation model can be used to simulate and calculate the critical flow under upstream temperature and pressure conditions.
-
Key words:
- Critical flow /
- Cavitation model /
- Bubble radius /
- Nucleation site density /
- Mass flow
-
表 1 湍流模型敏感性分析
Table 1. Sensitivity Analysis of Turbulence Models
湍流模型 质量流量/(kg·s−1) 相对误差/% 计算值 实验值 标准 k−ε 0.145 0.167 13.17 RNG k−ε 0.146 12.57 Realizable k−ε 0.147 11.98 标准 k−ω 0.143 14.37 BSL k−ω 0.142 14.97 SST k−ω 0.144 13.77 表 2 计算工况表
Table 2. Calculation Condition Table
工况 rb/mm ${\alpha _{{\text{nuc}}}}$ Fvap C 工况1 0.01 0.0005 50 2.50 工况2 0.02 1.25 工况 3 0.03 0.83 工况 4 0.04 0.63 工况 5 0.05 0.50 工况 6 0.06 0.42 -
[1] 陈听宽,徐进良,罗毓珊. 两相临界流实验研究[J]. 工程热物理学报,2002, 23(5): 623-626. doi: 10.3321/j.issn:0253-231X.2002.05.028 [2] 陈玉宙,赵民富,李伟卿,等. 临界流稳态试验研究[J]. 中国原子能科学研究院年报,2014: 133-134. [3] FAUSKE H K. The discharge of saturated water through tubes[J]. Chemical Engineering Progress, Symposium Series, 1965, 61(59): 210-216. [4] CELATA G P, CUMO M, D'ANNIBALE F, et al. The influence of non-condensible gas on two-phase critical flow[J]. International Journal of Multiphase Flow, 1988, 14(2): 175-187. doi: 10.1016/0301-9322(88)90004-3 [5] 陈玉宙,杨春生,张曙明,等. 基于喷嘴临界流实验对现有物理模型的评价[J]. 原子能科学技术,2009, 43(6): 485-490. [6] 吕玉凤,赵民富,杜开文. 过冷水流经通道的闪蒸临界流模型研究[J]. 原子能科学技术,2019, 53(12): 2375-2381. doi: 10.7538/yzk.2018.youxian.0831 [7] 王晓峰. 窄矩形通道内两相流动实验研究与数值模拟[D]. 哈尔滨: 哈尔滨工程大学, 2015. [8] 赵建福,张征明,张忠海. 窄缝内汽液两相临界流实验研究[J]. 工程热物理学报,2007, 28(S1): 173-176. [9] 黄波,张波涛,孙运达,等. 矩形狭缝流道临界流研究[J]. 西安交通大学学报,2020, 54(5): 26-31. [10] 周源,夏兆阳,汪杨乐,等. 超临界CO2临界流数值模拟[J]. 核科学与技术,2017, 5(3): 177-184. [11] 吕玉凤, 赵民富, 陈玉宙, 等. 两相临界流试验的RELAP5/MOD3.3程序模拟[C]//第十五届全国反应堆热工流体学术会议暨中核核反应堆热工水力技术重点实验室学术年会论文集. 荣成: 中国核学会, 2017. [12] BEN-MANSOUR R, HABIB M A, KHALIFA A, et al. Computational fluid dynamic simulation of small leaks in water pipelines for direct leak pressure transduction[J]. Computers & Fluids, 2012, 57: 110-123. [13] SINGHAL A K, ATHAVALE M M, LI H Y, et al. Mathematical basis and validation of the full cavitation model[J]. Journal of Fluids Engineering, 2002, 124(3): 617-624. doi: 10.1115/1.1486223 [14] BRENNEN C E. Cavitation and Bubble Dynamics[M]. Cambridge: Cambridge University Press, 2013: 30-58. [15] ZWART P J, GERBER A G, BELAMRI T. A two-phase flow model for predicting cavitation dynamics[C]//Proceedings of ICMF 2004 International Conference on Multiphase Flow. Yokohama: 2004.