液氮过冷沸腾空泡份额分布规律的研究

为了分析低温流体在过冷沸腾中的相分布特性,对液氮在竖直圆管内的过冷沸腾过程进行了CFD数值模拟。通过修正气泡参量模型和改变过程的操作参数,系统地分析了空泡份额在加热上升管内的分布情况。结果表明,气液相间的热质传递决定着空泡份额的分布,截面平均空泡份额沿轴向非线性递增,沿径向呈U形分布。轴向截面平均空泡份额随着热流密度的增大而增大,随着入口质量流率和过冷度的增大而减小。     

加热上升管内过冷流动沸腾数值模拟

采用计算流体动力学（CFD）程序CFX4.4对加热上升管内过冷流动沸腾工况下气水两相流动局部两相流参数（空泡份额和汽泡尺寸）进行了数值模拟。对数值差分方法、相关模型（界面力和气泡诱导的紊流）和汽泡尺寸进行了敏感性分析。空泡份额分布计算结果与实验结果比较表明，在低空泡份额工况下，两者符合较好，在高空泡份额工况下两者存在一定偏差，并且气相速度和汽泡尺寸的计算结果不理想。计算结果与实验结果之间的差异说明程序模型对于加热上升管内过冷流动沸腾模拟并不完善，建立更为合理的汽泡尺寸模型，考虑汽泡的合并和撕裂是必要的。     

摇摆条件下圆管内过冷沸腾局部空泡时空分布特性

过冷流动沸腾在核能、动力、化工等工业领域广泛存在,局部空泡分布特性的准确预测对构建两相流数理模型及两相流动压降和传热特性计算均具有重要意义。实验借助光学探针测量手段研究摇摆条件下圆管内过冷沸腾局部空泡分布特性,并提出多周期同相位叠加方法获得摇摆运动条件下局部瞬时空泡份额。实验结果表明:摇摆条件下圆管内局部空泡分布呈现周期性波动规律,摇摆条件下圆管轴心区及近壁区均会出现局部空泡份额峰值,且流道近壁区局部空泡份额波动幅度最大;与静止条件下不同,摇摆运动周期内过冷沸腾局部空泡份额空间分布形态随时间发生变化;摇摆周期和摇摆角度对过冷沸腾局部空泡分布特性也存在一定影响。     

摇摆条件下圆管内过冷沸腾局部空泡时空分布特性

过冷流动沸腾在核能、动力、化工等工业领域广泛存在,局部空泡分布特性的准确预测对构建两相流数理模型及两相流动压降和传热特性计算均具有重要意义。实验借助光学探针测量手段研究摇摆条件下圆管内过冷沸腾局部空泡分布特性,并提出多周期同相位叠加方法获得摇摆运动条件下局部瞬时空泡份额。实验结果表明:摇摆条件下圆管内局部空泡分布呈现周期性波动规律,摇摆条件下圆管轴心区及近壁区均会出现局部空泡份额峰值,且流道近壁区局部空泡份额波动幅度最大;与静止条件下不同,摇摆运动周期内过冷沸腾局部空泡份额空间分布形态随时间发生变化;摇摆周期和摇摆角度对过冷沸腾局部空泡分布特性也存在一定影响。     

高压流体过冷流动沸腾的数值模拟

在二流体模型基础上,采用RPI模型计算加热壁面的相间质量传输,对竖直圆管内高压流体过冷沸腾流动过程进行数值模拟.研究发现:壁面温度沿流动方向有明显拐点.根据壁温变化趋势可以判断,转折点位置即为过冷沸腾起始点;管内流体速度的增长曲线并不是简单单调递增的过程,而是速度增长由刚开始的缓慢上升到急剧上升,与空泡份额的分布曲线相...     

液氮过冷沸腾壁面换热特性的数值研究

应用计算流体力学(CFD)研究手段,对液氮过冷沸腾过程进行了模拟。系统分析了过冷沸腾壁面换热过程中存在的各传热模式,得到3部分热流沿管长的分布规律。计算结果表明,在过冷沸腾后期,蒸发热流成为壁面换热的主要模式。同时,就壁面热流密度的变化对3部分热流的影响进行了分析。随着壁面热流的增加,过冷沸腾出口处淬火热流和对流热流将减小,只有蒸发热流增大,且所占总热流份额绝对占优,壁面热流的增加促进了沸腾换热。     

高压高过冷度下过冷流动沸腾数值模拟

采用双流体模型对高压高过冷度下垂直圆管中水的过冷流动沸腾进行了数值模拟。通过对比不同气泡直径模型,揭示了气泡直径对于壁面传热方式的影响,确定了适合高压工况的气泡直径模型。考察了压力及壁面热通量对流动及传热特性的影响。计算结果表明,压力增加气泡脱离直径减小,单相对流传热所占比例增加,表面传热系数减小。高压高过冷度特征决定了气泡相分布极不均匀,随着热通量的增加,壁面附近容易形成气泡的密集,对过冷流动沸腾中的传热特性有重要影响。     

竖直圆管内液氮过冷流动沸腾数值模拟研究

在现有数学模型的基础上,从沸腾换热的机理入手,对过冷流动沸腾模型中的气泡参量模型进行了修正,分析确定了壁面热流密度的拆分方法,构建出适应于液氮的过冷沸腾计算模型。将新模型应用于CFX4.4中,对液氮在三维竖直圆管内的过冷流动沸腾进行了数值模拟。研究发现,在过冷流动沸腾形成之后,蒸发热流成为壁面换热的主要部分,沸腾换热占据换热的主导地位。数值模拟结果与已有的实验数据吻合较好,证明了新建模型的正确性。     

不同宽度窄缝通道过冷沸腾

以去离子水为实验工质,在窄缝宽度δ=3、4 mm,质量流速G=143、300 kg·m^-2·s^-1,主流过冷度ΔT_sub=17、25℃,热通量q=1～20 W·cm^-2的参数范围内,对常压下竖直窄缝通道内向上流动过冷沸腾的换热规律进行了实验研究。对不同宽度窄缝通道内的同一区域过冷沸腾气泡演变过程进行了可视化实验分析,发现窄缝宽度因素对过冷流动沸腾的流动换热特性和壁面核化特性影响显著,其中包括沸腾起始点ONB,压降ΔP,传热系数h,汽化核心密度N_a,气泡脱离直径D_d,气泡脱离频率f等。     

竖直环形通道内液氮流动沸腾的数值模拟

在多尺寸组模型的基础上,从加热壁面上脱离汽泡的受力分析人手,对液氮过冷流动沸腾模型进行了修正.将新模型应用于环形通道内液氮过冷流动沸腾的数值模拟,同时为了比较,采用基于Kirichenko,Fritz汽泡脱离直径公式的多尺寸组模型对同一管道内液氮过冷流动进行了数值模拟.结果表明:结合脱离汽泡受力分析模型的多尺寸组模型可...     

基于MUSIG模型的竖直圆管内液氮过冷沸腾数值模拟研究(英文)

Multiple size group （MUSIG） model combined with a threedimensional twofluid model were em ployed to predict subcooled boiling flow of liquid nitrogen in a vertical upward tube. Based on the mechanism of boiling heat transfer, some important bubble model parameters were amended to be applicable to the modeling of liquid nitrogen. The distribution of different discrete bubble classes was demonstrated numerically and the distribu tion patterns of void fraction in the wallheated tube were analyzed. It was found that the average void fraction in creases nonlinearly along the axial direction with wall heat flux and it decreases with inlet mass flow rate and sub cooled temperature. The local void fraction exhibited a Ushape distribution in the radial direction. The partition of the wall heat flux along the tube was obtained. The results showed that heat flux consumed on evaporation is the leading part of surface heat transfer at the rear region of subcooled boiling. The turning point in the pressure drop curve reflects the instability of bubbly flow. Good agreement was achieved on the local heat transfer coefficient aalnst experimental measurements, which demonstrated the accuracy of the numerical model.     

垂直圆管内液氮流动沸腾的理论模型及数值模拟

分析了液氮流动沸腾过程中气液两相间动量、能量以及质量的传输规律,建立了相应的理论模型,新模型重点修正了界面面积浓度和气泡挣脱直径的计算式;采用新建立的理论模型作为封闭方程对CFX-4.3中内建的双流体模型进行了修正,并采用修正后的双流体模型模拟了液氮在垂直圆管内的流动沸腾过程.数值模拟的结果与文献中的实验数据吻合较好,证明了本文所建模型的合理性.通过数值模拟发现,两相流参数分布的不均匀性对液氮流动沸腾过程中的热质传输特性有重要影响.     

Numerical investigation on effects of fuel rod with different shapes of bow deformation on the subcooled boiling of coolant in a fuel assembly

The Eulerian two-phase boiling model of the subcooled boiling of coolant in a 3 × 3 fuel assembly is established and compared with the experimental data to verify its accuracy. The effects of four shapes of bow deformation on subcooled boiling flow and heat transfer characteristics are obtained by comparing and analyzing the distributions of thermal–hydraulic parameters, including the axial wall superheat, liquid phase temperature, axial void fraction, various heat fluxes, heat transfer coefficient, and turbulent kinetic energy. All shapes of bow deformation will lead to the redistribution of coolant among different subchannels, and the distributions of thermal–hydraulic parameters would be greatly affected. The bow deformation of fuel rod easily causes local boiling, which results in local high temperature of liquid phase and bubble accumulation, and a decrease in the area of high liquid phase heat transfer coefficient on the surface of bowing fuel rod. Additionally, the non-uniform distribution of turbulent kinetic energy caused by bow deformation in different axial sections not only affects the heat transfer performance of coolant, but also causes the increment in pressure drop, which has negative effects on the safe operation of the nuclear reactor. This paper can provide data and theoretical support for engineering design.     

三维竖直圆管内流动沸腾的数值模拟

沸腾传热及气液两相流动都是常见的物理现象,其流动形式及传热机理复杂.文章借助于CFX流体计算平台,模拟了水在三维竖直圆管内的过冷流动沸腾过程,给出管内流体状态参数的变化规律.得到沿管长的气相体积分数的变化规律,沿管长沸腾传热系数的变化趋势,明显地看到沸腾对传热的强化作用,得到径向液体温度分布情况.对过冷流动沸腾的内在机...     

Population balance modelling for bubbly flows with heat and mass transfer

Population balance equations combined with a three-dimensional two-fluid model are employed to predict bubbly flows with the presence of heat and mass transfer processes. Subcooled boiling flow belongs to this specific category of bubbly flows is considered. The MUSIG (MUltiple-SIze-Group) model implemented in CFX4.4 is further developed to account for the wall nucleation and condensation in the subcooled boiling regime. Comparison of model predictions against local measurements near the test channel exit is made for the radial distribution of the bubble Sauter diameter, void fraction, interfacial area concentration and gas and liquid velocities covering a range of different mass and heat fluxes and inlet subcooling temperatures. Additional comparison was also performed against existing boiling model in CFX4.4 and the modified model developed in our previous work (Int. J. Heat Mass Transfer 45 (2002) 1197). Good agreement is better achieved with the local radial bubble Sauter diameter, void fraction, interfacial area concentration and liquid velocity profiles against measurements using the newly formulated MUSIG boiling model over the simpler boiling models. However, significant weakness of the model is still evidenced in the prediction of the vapour velocity. Work is in progress to circumvent the deficiency of the model by the consideration of additional momentum equations or an algebraic slip model to account for bubble separation.