液态金属内单个气泡上升行为的MPS法数值模拟

液态金属冷却核反应堆采用气泡泵的概念设计来提升堆芯自然循环能力。液态金属内气液两相流动特征将直接影响核反应系统一回路的自然循环能力及堆芯安全。本研究通过采用移动粒子半隐式（MPS）方法,对液态金属中单个上升气泡的气泡动力学行为进行数值模拟。分析了铅铋合金中3种初始直径不同的单个氮气泡在上升过程中的气泡形状和速度的变化趋势;对比了初始直径相同的单个氮气泡在液钾、液钠、铅铋合金、钾钠合金和锂铅合金5种液态金属中的上升行为;同时将模拟得到的气泡形状与Grace经验关系图进行了对比,验证了MPS方法数值模拟结果的正确性。     

基于扩散界面法的液态LBE中单个弹状气泡上升行为CFD研究

利用扩散界面法对液态铅铋合金中弹状气泡在垂直管中的上升过程进行数值模拟研究,用来验证扩散界面法在模拟液态铅铋合金(LBE)中弹状气泡上升行为的准确性和可行性,并对模拟结果进行进一步分析,来解释弹状气泡在LBE中的上升行为。数值模拟得到不同液体流速情况下弹状气泡在上升过程中的形态和速度变化,数值模拟结果与文献中的经验关系式以及实验中的结果都吻合良好,证明了扩散界面法在模拟液态铅铋合金中弹状气泡上升行为的准确性和可行性,为液态金属和弹状气泡的两相流提供了一种新方法。对数值模拟结果进行进一步分析,发现气泡尾部形态变化较大,会在尾部两端分裂出子气泡,尾部附近流场产生旋涡,受到强烈干扰而出现湍流,对于提高换热效率起到积极作用。气泡上升对尾部区域产生影响的最远距离约80 mm,为连续气泡的注入提供了气泡间距参考值。     

气泡在液态铅铋合金内上升行为的数值模拟

为利用气泡提升泵的概念设计来增加铅铋合金非等温回路的自然循环能力,需对充入回路的气泡大小进行选择。气泡的初始直径能显著影响气泡提升泵提升自然循环的能力。利用VOF模型,对不同初始直径的氦气泡在液态铅铋合金中的提升作用进行了数值模拟研究。结果表明,对于单个气泡,随着初始直径增大,提升能力呈先增加后趋于不变的趋势;初始直径较大的气泡易分裂生成较小的子气泡附着在管壁上,有影响传热效率的潜在可能。研究了在静止液态铅铋合金中,不同初始直径的氦气泡上升速度和上升高度与时间的关系及气泡分裂的原因。     

基于MPS方法的液态铅铋合金内气泡上升流数值模拟

移动粒子半隐式(MPS)方法捕捉多相流体动力学相界面的能力比传统网格方法具有明显优势.本研究采用MPS方法,使用FORTRAN语言自编程序,对单个氩气气泡在液态铅铋合金内从静止到充分发展整个过程中的瞬态动力学行为进行二维数值模拟,得到气泡变形特性与上升速度的关系.结果表明:在气泡上升过程中,气泡由球形先变成酒窝形状,最...     

基于MPS方法模拟液态铅铋合金中弹状气泡上升

基于移动粒子半隐式(MPS)方法对液态铅铋合金中单个弹状气泡的垂直上升行为进行数值模拟,得到弹状气泡上升过程的形态变化、气泡最终上升速度拟合直线及下降液膜厚度与下降液膜中的轴向速度分布。将部分数值模拟结果与文献中的实验结果进行比较,揭示了弹状气泡最终上升速度、下降液膜厚度及其中轴向速度分布所满足的规律。比较结果证明了所选取计算模拟模型的正确性与合理性,以及MPS-MAFL方法用于模拟弹状气泡上升行为的准确性与可靠性。     

环形通道内液态铅铋合金流动换热特性实验研究

通过对环形通道内液态铅铋合金的流动换热特性进行实验研究,得到了气泡泵注气对液态金属流动的影响,并拟合出环形通道内液态铅铋合金的摩擦系数关系式和换热特性关系式。结果表明：采用气泡泵注气能有效提升铅铋合金的质量流速;相同Reynolds数下环形通道内液态铅铋合金的摩擦系数大于由布拉休斯公式计算得到的摩擦系数;液态铅铋合金对流换热过程中,导热项占主导地位,并且Nusselt数随Peclet数的增大而增大。     

环形通道内液态铅铋合金流动换热特性实验研究

通过对环形通道内液态铅铋合金的流动换热特性进行实验研究,得到了气泡泵注气对液态金属流动的影响,并拟合出环形通道内液态铅铋合金的摩擦系数关系式和换热特性关系式。结果表明:采用气泡泵注气能有效提升铅铋合金的质量流速;相同Reynolds数下环形通道内液态铅铋合金的摩擦系数大于由布拉休斯公式计算得到的摩擦系数;液态铅铋合金对流换热过程中,导热项占主导地位,并且Nusselt数随Peclet数的增大而增大。     

基于扩散界面法的较大气泡上升过程数值模拟

基于扩散界面法和有限元法,对较大气泡在上升阶段的形态和速度进行了模拟,结果与实验吻合较好,说明该方法能准确地模拟气泡的运动特性。利用该模型,对初始直径不同的较大气泡上升过程中的形态、速度和振荡随时间变化的规律进行了分析。并分析了14mm直径的气泡在不同尺寸通道中上升过程的形态、速度的变化规律。结果表明:气泡的稳定形态随着气泡初始直径的增大由椭球形变为球帽形,且达到稳定形状的时间更长。气泡初始直径越大,气泡的顶端速度越快,并稍有波动。而气泡的底端速度开始快速增大使气泡向内凹陷,随后回落并在气泡顶端速度上下振荡。气泡上升通道越窄,气泡达到稳定形态的时间越长,顶端速度越小,气泡的高宽比越大。     

液态锂铅合金的氚释放行为

为了完成聚变堆液态锂铅包层鼓泡提氚系统的工程设计和建造,以金属与氢的作用理论为基础,建立了氚从液态锂铅中的动力学释放行为的数学模型.计算和分析了温度、饱和器氚分压、氦流量对解吸器顶部气相中的氚分压的影响以及氚在液态锂铅中的传质系数、解吸率和吸附率.结果表明:在633～723 K的解吸温度范围内,氚从液态锂铅到气相的整个释放过程虽然包含了氚在熔融合金气泡中的扩散与对流、氚通过与气-液界面相连合金层的扩散、在界而发生的氚原子重组多相反应、氚通过气相边界层的扩散和气相中氚的扩散与对流5个子过程,但起决定作用的是氚在合金内的扩散和气.液界面的多相反应重组,其他子过程意义不大.     

铅水反应中铅铋合金凝固的数值模拟

为了研究铅铋合金在蒸汽发生器传热管破裂（SGTR）事故所引发的铅铋合金与水反应过程中的凝固机理,通过耦合VOF模型、Realizable k-ε湍流模型、凝固传热模型,利用FLUENT软件建立了铅铋合金与水反应过程的二维仿真模型,并将该模型与现有反应实验的结果进行对比验证。随后基于热焓法建立可以直观描述铅铋合金凝固现象的凝固传热特性热焓方程,通过控制模型变量研究影响铅铋合金凝固发生的因素及条件,最后将该模型应用于复杂结构场景中。结果表明,铅铋合金与水的温差、水流喷射初始速度、注水管径是影响铅铋合金凝固的主导因素,本文提出的模型具有较高可靠性,能够模拟实际工况中铅铋合金的凝固现象。本研究所得到的机理性结论与现象学结论能够为铅基快堆安全分析提供理论支撑。      

Two-dimensional numerical simulation of single bubble rising behavior in liquid metal using moving particle semi-implicit method

Gas-lift pump in liquid metal cooling fast reactor (LMFR) is an innovative conceptual design to enhance the natural circulation ability of reactor core. The two phase flow characteristics of gas–liquid metal make significant improvement of the natural circulation capacity and reactor safety. It is important to study bubble flow in liquid metal. In present study, the rising behaviors of a single nitrogen bubble in 5 kinds of common stagnant liquid metals (lead bismuth alloy (LBE), liquid kalium (K), sodium (Na), potassium sodium alloy (Na–K) and lithium lead alloy (Li–Pb)) and in flowing lead bismuth alloy have been numerically simulated using two-dimensional moving particle semi-implicit (MPS) method. The whole bubble rising process in liquid was captured. The bubble shape, rising velocity and aspect ratio during rising process of single nitrogen bubble were studied. The computational results show that, in the stagnant liquid metals, the bubble rising shape can be described by the Grace's diagram, the terminal velocity is not beyond 0.3 m/s, the terminal aspect ratio is between 0.5 and 0.6. In the flowing lead bismuth alloy, as the liquid velocity increases, both the bubble aspect ratio and terminal velocity increase as well. This work is the fundamental research of two phase flow and will be important to the study of the natural circulation capability of Accelerator Driven System (ADS) by using gas-lift pump.     

Numerical simulation of single bubble condensation in subcooled flow using OpenFOAM

The single condensing bubble behavior in subcooled flow has been numerical investigated using the open source code OpenFOAM. A coupled Level Set (LS) and Volume of Fluid (VOF) method (CLSVOF) model with a phase change model for condensation was developed and implemented in the code. The simulated results were firstly compared with the experimental results, they were in great agreements, and thus the simulation model was validated. The validated numerical model was then used to analyze the condensing bubble deformation, bubble lifetime, bubble size history, condensate Nusselt number and other interesting parameters with different variables in subcooled flow. The numerical results indicated that the initial bubble size, subcooling of liquid and system pressure play an important role to influence the condensing bubble behaviors significantly and bubble will be pierced when the subcooling and initial diameter reach a certain value at the later condensing stage. The bubble diameter history and condensate Nusselt number were found in good agreement with the empirical correlation. The drag force coefficient was predicted well by introducing a reduced drag coefficient.     

静水中单气泡运动特性实验研究

气泡运动特性是两相流研究的重要领域之一。本文基于图像识别技术对静水中单气泡的上升过程及其动力学特性进行研究，分析了不同等效体积直径下气泡上升过程中轨迹、变形程度、瞬时速度等参数。研究结果表明：椭球型单气泡上升过程中横向振荡幅度峰值与横向速度峰值呈正比，瞬时纵横比与瞬时纵向速度呈反比；单气泡横向运动位置、纵横比、上升速度具有周期性变化的现象。通过定量分析气泡末速度与等效体积直径的关系，对已有预测关系式进行了比较和评价，并基于实验数据拟合得到了一个更为精确的显式末速度关系式。     

Experimental Investigation on Characteristic of Single Bubble Motion in Stagnant Water

The bubble motion characteristic is one of the most important fields in two-phase flow research. Based on the image recognition method, the characteristic and dynamics of single bubble rising process in stagnant water were analyzed. The trajectory, deformation extent and instantaneous terminal velocity of single bubble rising process with different equivalent volume diameters were analyzed. The research results show that the peak value of lateral oscillation amplitude is proportional to the peak value of lateral velocity. The instantaneous aspect ratio is inversely proportional to instantaneous axial velocity. The phenomenon that the lateral motion position, aspect ratio and rising speed of a single bubble has periodical change. Based on the quantitative analysis of the relationship between bubble terminal velocity and bubble equivalent volume diameter, the existing prediction formulas were compared and evaluated. A more accurate explicit terminal velocity correlation based on experimental data was obtained.     

上升气泡内气溶胶沉降机理研究

气溶胶池洗过滤是反应堆严重事故中去除放射性源项的重要手段。本文以严重事故条件下上升气泡中气溶胶的滞留过程为背景，设计搭建了可视化单气泡鼓泡实验装置。通过该装置研究了气溶胶在上升气泡中的沉降效率，并与MELCOR中的气溶胶沉降模型计算结果进行了对比。结果表明，气溶胶沉降效率对气泡尺寸的变化较为敏感，当气体流量大于0.1 L/min时，气泡等效直径迅速增加，相应的气溶胶沉降效率快速降低；与MELCOR模型计算结果的对比表明，两者在总体趋势上呈现出较好的一致性，但计算结果低估了液相对气溶胶的实际去除能力，导致这种偏差的主要原因是气泡在上升过程中存在无规则的晃动以及气液界面的波动。     

Analysis of Control Rod Effect on Neutronics in FBR by Three-Dimensional Calculation

Experiments were carried out to study the mass transfer of iodine from a xenon-iodine mixed gas bubble to a liquid sodium pool. The bubble was generated in the pool by cracking a quartz glass ball which contains the xenon-iodine mixed gas. The mixed gas released into the argon cover gas was collected to determine the transferred iodine mass into the sodium pool. A rising velocity of the bubble was measured by Chen-type sensors arranged vertically in the pool. From the measured rising velocity and another observation of bubble behavior in simulated water experiments, it is found that the generated bubble breaks up into several smaller bubbles of spherical cap type during the rising period. Transferred iodine mass per the unit initial bubble volume from the bubbles to the sodium pool shows increases with increasing time and the initial iodine concentration. A mass transfer rate obtained by differentiating the transferred iodine mass with respect to the time indicates a rapid decrease just after the bubble generation and a slow decrease for the successive period. Measured decontamination factors are described as a function of the time and the iodine concentration.     

Numerical investigation of bubble wake properties in the moving liquid with LES model

Wake flow caused by the relative motions between bubble and liquid phase influences bubble motions and enhances turbulent properties in the liquid phase. This phenomenon has been stressed for a better understanding of bubbly flow. In this paper, large eddy simulation (LES) is performed to simulate a single bubble rising in the moving liquid, with volume of fluid (VOF) method to capture the interface movements between bubble and liquid phase. The simulation results are firstly compared with the numerical and experimental data from the literature. A good agreement demonstrated the capability of the employed computational fluid dynamics (CFD) approach to predict turbulent properties in the liquid phase and capture interface movement as well as its induced wake flow. Consequently, the dynamic behaviors of a single bubble rising in the moving liquid were investigated. An ensemble averaged has been employed to evaluate the velocity distribution composed by wake velocity and liquid velocity quantitatively as well as the velocity fluctuations enhanced by bubble motion. Their dependency was also evaluated based on a systematic CFD simulation which covers a wide range of liquid velocity. With comparisons of the single phase flow, the influence from the existence of bubble on turbulent properties was determined.     

Testing and qualification of CIRCE instrumentation based on bubble tubes

The paper is focused on the development and the qualification of the instrumentation suitable for the level and the differential pressure measurements in the experimental facility CIRCE. This large-scale facility is installed at the ENEA Brasimone Centre for studying the fluid-dynamic and operating behavior of ADS reactor plants. This a rather challenging objective since the facility adopts the molten lead bismuth eutectic (LBE) alloy as a coolant and, at present, instrumentation qualified for operating under such conditions is scarce or does not exist. Bubble tubes have been installed in CIRCE to transfer pressure signals from the LBE to differential pressure cells operating with gas at room temperature. The bubble tube is a simple measuring device, but its use in LBE must be carefully assessed. Therefore, preliminary tests of bubble tubes in representative conditions have been carried out in a smaller test section. Experimental tests were performed at several temperatures, with LBE in stagnant conditions. The results obtained in these tests, aiming at checking the performance of the bubble tubes adopted in measuring pressure, differential pressure and level in the CIRCE facility, are discussed here. The obtained information will allow to calibrate the related measuring systems and to verify the accuracy and repeatability of the measurements, as a function of the injected gas flow rate, the tube diameter and the geometry of the tube exit section.