基于CFD-PBM方法熔盐堆换热器内夹带气泡流动规律分析

王子业; 朱贵凤; 陈亮; 梁晓琳; 邹杨

doi:10.3724/j.0253-3219.2026.hjs.49.250426

您当前的位置：

首页 >

文章列表页 >

基于CFD-PBM方法熔盐堆换热器内夹带气泡流动规律分析

中国核学会研究堆与新堆分会第二届学术年会优秀论文专刊 | 更新时间：2026-07-29

- 基于CFD-PBM方法熔盐堆换热器内夹带气泡流动规律分析
- Analysis of the flow pattern of entrained bubbly flow in a molten salt reactor heat exchanger based on CFD-PBM method
- 核技术 2026年49卷第7期文章编号：071122
- 作者机构：
  
  1.中国科学院上海应用物理研究所上海 201800
  2.钍基核裂变能全国重点实验室上海 201800
  3.中国科学院大学北京 100049
- 作者简介：
  
  王子业，男，1998年出生，2024年于上海科技大学获硕士学位，研究领域为熔盐堆夹带气泡CFD计算分析和实验验证
  朱贵凤，E-mail：zhuguifeng@sinap.ac.cn
  邹杨，E-mail：zouyang@sinap.ac.cn
- 基金信息：
  
  上海市自然科学基金(25ZR1401383);钍基熔盐堆核能系统产业基金项目(SINAP-CYJJ-202401;SINAP-CYJJ-202501)
- DOI：10.3724/j.0253-3219.2026.hjs.49.250426
  中图分类号： TL99
- 收稿：2025-12-05，
  
  修回：2026-02-28，
  
  网络首发：2026-07-29，
  
  纸质出版：2026-07-15
- 稿件说明：
移动端阅览
王子业,朱贵凤,陈亮等.基于CFD-PBM方法熔盐堆换热器内夹带气泡流动规律分析[J].核技术,2026,49(07):071122.

WANG Ziye,ZHU Guifeng,CHEN Liang,et al.Analysis of the flow pattern of entrained bubbly flow in a molten salt reactor heat exchanger based on CFD-PBM method[J].NUCLEAR TECHNIQUES,2026,49(07):071122.
王子业,朱贵凤,陈亮等.基于CFD-PBM方法熔盐堆换热器内夹带气泡流动规律分析[J].核技术,2026,49(07):071122. DOI： 10.3724/j.0253-3219.2026.hjs.49.250426. CSTR： 32193.14.hjs.CN31-1342/TL.2026.49.250426.

WANG Ziye,ZHU Guifeng,CHEN Liang,et al.Analysis of the flow pattern of entrained bubbly flow in a molten salt reactor heat exchanger based on CFD-PBM method[J].NUCLEAR TECHNIQUES,2026,49(07):071122. DOI： 10.3724/j.0253-3219.2026.hjs.49.250426. CSTR： 32193.14.hjs.CN31-1342/TL.2026.49.250426.

摘要

夹带气泡的引入对熔盐堆（Molten Salt Reactor，MSR）反应性控制、换热性能以及难溶裂变产物的迁移产生影响。其中换热器中夹带气泡流动十分复杂，且缺少实验数据，边界条件与流动分布规律的关系并不明确，使得气泡对换热器内换热性能和难溶裂变产物迁移行为等研究难以进一步开展。本文运用计算流体力学（Computational Fluid Dynamics，CFD）手段，采用基于群体平衡模型（Population Balance Model，PBM）的非均匀多尺寸组离散形式（inhomogeneous Multiple Size Group Model，iMUSIG），分析了两种换热器在不同入口速度、气体份额和气泡尺寸边界条件下，从入口到出口的速度场、气体份额、气泡尺寸和界面面积浓度（Interfacial Area Concentration，IAC）分布。结果表明，两种换热器中流场具有相同的涡结构而具有类似的气体聚集现象；气泡尺寸增加促进了气体聚集，高流速促使气体向涡心聚集；高气体份额和更低的速度促进气泡聚并，尺寸增大；从入口到出口，界面面积浓度由于单位气体份额内气泡尺寸增大呈下降趋势，低气体份额和大气泡尺寸的界面面积浓度变化平稳。以上结果为熔盐换热器内夹带气泡带来的影响研究提供参考。

Abstract

Background

Entrained bubbles in molten salt reactors affect reactivity control

heat transfer performance

and the migration of insoluble fission products. However

the relationship between inlet boundary conditions and bubbly-flow distribution in molten-salt heat exchangers remains unclear because of complex internal flow structures and limited experimental data.

Purpose

This study aims to investigate the flow patterns of entrained bubbly flow in molten-salt heat exchangers under different inlet velocities

gas fractions

and bubble-size boundary conditions.

Methods

Firstly

an inhomogeneous multiple size group model based on the population balance model was coupled with computational fluid dynamics to simulate polydisperse bubbly flow in two simplified shell-and-tube heat exchangers. Then

horizontal and vertical heat exchanger configurations were compared under varying inlet velocities

inlet gas fractions

and inlet bubble diameters. Finally

the distributions of velocity

gas fraction

bubble size

and interfacial area concentration from the inlet to the outlet were analyzed using consistent mesh and solver settings. Because simplified geometries were used

the results are mainly applicable to compact heat exchangers or similar simplified flow units with comparable symmetry and large-scale periodic flow structures.

Results

The vortex structures of the flow fields remain stable when the inlet velocity is 1.8~5.3 m·s

-1

in the horizontal heat exchanger and 2~5 m·s

-1

in the vertical heat exchanger

leading to similar gas accumulation behavior in both configurations. Under different inlet bubble diameters

the gas fraction ranges from 0.4% to 1.1% and increases with increasing inlet bubble diameter

while a local gas fraction exceeding 3% appears in the vortex region when the inlet gas fraction is only 0.5%. Higher inlet gas fractions and lower velocities promote bubble coalescence and increase the bubble size; the maximum bubble size in the horizontal heat exchanger is 1.9 times that in the vertical heat exchanger under the maximum gas-fraction condition and 2.3 ti

mes that under the minimum-velocity condition. From the inlet to the outlet

the interfacial area concentration generally decreases because the bubble size per unit gas fraction increases

whereas the interfacial area concentration changes steadily at a low inlet gas fraction of 0.25% and a large inlet bubble diameter of 0.5 mm.

Conclusion

This study reveals the effects of heat exchanger orientation and inlet boundary conditions on gas accumulation

bubble coalescence

and interfacial area concentration evolution in molten-salt heat exchangers. The quantitative results

including the 0.4%~1.1% gas-fraction range

the local gas fraction exceeding 3%

and the 1.9-2.3-fold difference in maximum bubble size between the two configurations

provide a reference for evaluating bubble entrainment

mass transfer

and heat exchanger optimization in molten salt reactor systems.

关键词

Keywords

references

Robertson R . MSRE design & operations report Part 1 description of reactor design [R ] . Oak Ridge, Tennessee : Oak Ridge Laboratory , 1965 . DOI: 10.2172/4654707 http://dx.doi.org/10.2172/4654707 .

陈玉爽 , 田健 , 丁梦婷 , 等 . 基于流热固耦合的U型管式熔盐换热器温度场与应力场分析 [J ] . 核技术 , 2023 , 46 ( 1 ): 010604 . DOI: 10.11889/j.0253-3219.2023.hjs.46.010604 http://dx.doi.org/10.11889/j.0253-3219.2023.hjs.46.010604 .

CHEN Yushuang , TIAN Jian , DING Mengting , et al . Temperature field and stress field analysis of molten salt U-tube heat exchanger based on fluid-thermal-solid coupling method [J ] . Nuclear Techniques , 2023 , 46 ( 1 ): 010604 . DOI: 10.11889/j.0253-3219.2023.hjs.46.010604 http://dx.doi.org/10.11889/j.0253-3219.2023.hjs.46.010604 .

Lee A J H , Kozlowski T . Effects of sparging gas entrainment on void and temperature coefficients in the Transatomic Power molten salt reactor core [J ] . Annals of Nuclear Energy , 2023 , 192 : 110013 . DOI: 10.1016/j.anucene.2023.110013 http://dx.doi.org/10.1016/j.anucene.2023.110013 .

Bajpai P , Lorenzi S , Cammi A . A multiphysics model for analysis of inert gas bubbles in Molten Salt Fast Reactor [J ] . The European Physical Journal Plus , 2020 , 135 ( 5 ): 409 . DOI: 10.1140/epjp/s13360-020-00428-2 http://dx.doi.org/10.1140/epjp/s13360-020-00428-2 .

Mochizuki H . Study on start-up of molten salt fast modular reactor from shutdown based on two-phase flow CFD analysis [J ] . Nuclear Engineering and Design , 2023 , 414 : 112560 . DOI: 10.1016/j.nucengdes.2023.112560 http://dx.doi.org/10.1016/j.nucengdes.2023.112560 .

Caruggi F , Cammi A , Cervi E , et al . Multiphysics modelling of gaseous fission products in the Molten Salt Fast Reactor [J ] . Nuclear Engineering and Design , 2022 , 392 : 111762 . DOI: 10.1016/j.nucengdes.2022.111762 http://dx.doi.org/10.1016/j.nucengdes.2022.111762 .

Scioscioli F , Cammi A , Lorenzi S . Analysis of the operational and safety features of the in-core bubbling system of the molten salt fast reactor [J ] . Nuclear Science and Engineering , 2024 , 198 ( 6 ): 1288 – 1307 . DOI: 10.1080/00295639.2023.2250144 http://dx.doi.org/10.1080/00295639.2023.2250144 .

Frederix E M A , Komen E M J . Simulation of noble metal particle growth and removal in the molten salt fast reactor [J ] . Nuclear Engineering and Design , 2023 , 415 : 112690 . DOI: 10.1016/j.nucengdes.2023.112690 http://dx.doi.org/10.1016/j.nucengdes.2023.112690 .

Tas-Koehler S , Liao Y X , Hampel U . A critical analysis of drag force modelling for disperse gas-liquid flow in a pipe with an obstacle [J ] . Chemical Engineering Science , 2021 , 246 : 117007 . DOI: 10.1016/j.ces.2021.117007 http://dx.doi.org/10.1016/j.ces.2021.117007 .

Li C X , Cui Y Z , Zhai L X , et al . Investigation into drag coefficient for sub-millimeter bubbles in gas–liquid bubbly flow: Experiments and CFD simulations [J ] . Chemical Engineering Journal , 2023 , 478 : 147236 . DOI: 10.1016/j.cej.2023.147236 http://dx.doi.org/10.1016/j.cej.2023.147236 .

Nishi Y , Kogawa H , Utsumi Y , et al . Improvement of bubble distribution characteristics through multi-objective optimization of flow characteristics of a swirling flow type microbubble generator with fixed blades [J ] . Chemical Engineering Research and Design , 2024 , 212 : 457 – 471 . DOI: 10.1016/j.cherd.2024.11.022 http://dx.doi.org/10.1016/j.cherd.2024.11.022 .

Catrawedarma I G N B , Ton S , Pranowo D D , et al . Hydrodynamic characteristics of the microbubble dissolution in water using an ejector-type bubble generator [J ] . Case Studies in Chemical and Environmental Engineering , 2025 , 11 : 101043 . DOI: 10.1016/j.cscee. 2024.101043 http://dx.doi.org/10.1016/j.cscee.2024.101043 .

Lopez J , Ratkovich N , Pereyra E , et al . Computational fluid dynamics modeling of gas bubble separation efficiency in downhole separators and vertical/inclined annulus: an experimental validation and analysis [J ] . Geoenergy Science and Engineering , 2024 , 243 : 213281 . DOI: 10.1016/j.geoen.2024.213281 http://dx.doi.org/10.1016/j.geoen.2024.213281 .

Zeng X B , Fan G M , Wang M , et al . Experimental and numerical study on the length of LZVV in the bubble separator for molten salt reactor [J ] . Nuclear Engineering and Design , 2023 , 401 : 112100 . DOI: 10.1016/j.nucengdes.2022.112100 http://dx.doi.org/10.1016/j.nucengdes.2022.112100 .

Wang Z Y , Zhu G F , Zou Y , et al . Numerical simulation of entrained bubbles flow in the shell-tube heat exchanger of MSRs based on population balance model [J ] . International Journal of Energy Research , 2025 , 2025 ( 1 ): 6907471 . DOI: 10.1155/er/6907471 http://dx.doi.org/10.1155/er/6907471 .

Frank T , Zwart P J , Krepper E , et al . Validation of CFD models for mono- and polydisperse air-water two-phase flows in pipes [J ] . Nuclear Engineering and Design , 2008 , 238 ( 3 ): 647 – 659 . DOI: 10.1016/j.nucengdes.2007. 02.056 http://dx.doi.org/10.1016/j.nucengdes.2007.02.056 .

Krepper E , Lucas D , Frank T , et al . The inhomogeneous MUSIG model for the simulation of polydispersed flows [J ] . Nuclear Engineering and Design , 2008 , 238 ( 7 ): 1690 – 1702 . DOI: 10.1016/j.nucengdes.2008.01.004 http://dx.doi.org/10.1016/j.nucengdes.2008.01.004 .

Wang T F , Wang J F , Jin Y . A CFD-PBM coupled model for gas-liquid flows [J ] . AIChE Journal , 2006 , 52 ( 1 ): 125 – 140 . DOI: 10.1002/aic.10611 http://dx.doi.org/10.1002/aic.10611 .

Ansys Inc . Ansys fluent user’s guide [EB/OL ] . Pennsylvania : ANSYS , Inc, 2024 . https://ansyshelp.ansys.com/public/account/secured?returnurl=/Views/Secured/corp/v242/en/flu_ug/flu_ug.html https://ansyshelp.ansys.com/public/account/secured?returnurl=/Views/Secured/corp/v242/en/flu_ug/flu_ug.html .

Ishii M , Zuber N . Drag coefficient and relative velocity in bubbly, droplet or particulate flows [J ] . AIChE Journal , 1979 , 25 ( 5 ): 843 – 855 . DOI: 10.1002/aic.690250513 http://dx.doi.org/10.1002/aic.690250513 .

Tomiyama A , Tamai H , Zun I , et al . Transverse migration of single bubbles in simple shear flows [J ] . Chemical Engineering Science , 2002 , 57 ( 11 ): 1849 – 1858 . DOI: 10.1016/S0009-2509(02)00085-4 http://dx.doi.org/10.1016/S0009-2509(02)00085-4 .

Hosokawa S , Tomiyama A , Misaki S , et al . Lateral migration of single bubbles due to the presence of wall [C ] . Volume 1: Fora, Parts A and B. Montreal, Quebec, Canada : ASMEDC , 2002 : 855 – 860 . DOI: 10.1115/FEDSM2002-31148 http://dx.doi.org/10.1115/FEDSM2002-31148 .

Burns A , Frank T , Hamill I , et al . The favre averaged drag model for turbulent dispersion in eulerian multi-phase flows [R ] . Yokohama, Japan : Conference on Multiphase Flow, ICMF 2004 : 392 , 2004.

Milne-Thomson L M . Theoretical hydrodynamics, 5th ed [M ] . London : Macmillan Education , 1968 : 246 – 249 . DOI: 10.1115/1.3601335 http://dx.doi.org/10.1115/1.3601335 .

Menter F R . Two-equation eddy-viscosity turbulence models for engineering applications [J ] . AIAA Journal , 1994 , 32 ( 8 ): 1598 – 1605 . DOI: 10.2514/3.12149 http://dx.doi.org/10.2514/3.12149 .

Sato Y , Sekoguchi K . Liquid velocity distribution in two-phase bubble flow [J ] . International Journal of Multiphase Flow , 1975 , 2 ( 1 ): 79 – 95 . DOI: 10.1016/0301-9322(75)90030-0 http://dx.doi.org/10.1016/0301-9322(75)90030-0 .

Liao Y X , Rzehak R , Lucas D , et al . Baseline closure model for dispersed bubbly flow: bubble coalescence and breakup [J ] . Chemical Engineering Science , 2015 , 122 : 336 – 349 . DOI: 10.1016/j.ces.2014.09.042 http://dx.doi.org/10.1016/j.ces.2014.09.042 .

Kedl R J , McGlothlan C K . Tube vibration in msre primary heat exchanger [R ] . Oak Ridge, Tennessee : Oak Ridge Laboratory , 1968 . DOI: 10.2172/4553122 http://dx.doi.org/10.2172/4553122 .

崔蕾 , 陈玉爽 , 李启明 , 等 . 熔盐堆用U型管壳式换热器的传热与流动特性模拟分析 [J ] . 核技术 , 2024 , 47 ( 7 ): 070606 . DOI: 10.11889/j.0253-3219.2024.hjs.47.070606 http://dx.doi.org/10.11889/j.0253-3219.2024.hjs.47.070606 .

CUI Lei , CHEN Yushuang , LI Qiming , et al . Simulation analysis of heat transfer and flow characteristics of a U-tube heat exchanger in a molten salt reactor [J ] . Nuclear Techniques , 2024 , 47 ( 7 ): 070606 . DOI: 10.11889/j.0253-3219.2024.hjs.47.070606 http://dx.doi.org/10.11889/j.0253-3219.2024.hjs.47.070606 .

Janz G J , Gardner G L , Krebs U , et al . Molten salts: volume 4, Part 1, fluorides and mixtures electrical conductance, density, viscosity, and surface tension data [J ] . Journal of Physical and Chemical Reference Data , 1974 , 3 ( 1 ): 1 – 115 . DOI: 10.1063/1.3253134 http://dx.doi.org/10.1063/1.3253134 .

Beall S E , Haubenrich P N , et al . MSRE design and operations report Part 5 Reactor safety analysis report(ORNL-TM-0732) [R ] . Oak Ridge, Tennessee : Oak Ridge Laboratory , 1964 .

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

熔盐堆降功率过程中温度负反馈效应的影响研究

基于先驱核迁移与蒙特卡罗粒子输运耦合的熔盐堆缓发中子共轭通量计算

超压工况下FLiBe熔盐浸渗对T220石墨组件中子学性能的影响

熔盐堆应急排盐冷却系统换热管内层屏障破损传热影响

基于LSTM的熔盐空气换热器关键参数预测模型