基于分子动力学理论界面传质的过冷沸腾汽泡生长特性研究

基于分子动力学理论的准平衡态界面处界面蒸发/冷凝因素,以及汽泡底部微液层传热因素建立综合传热传质相变模型,对窄通道内汽泡过冷流动沸腾条件下的生长情况进行模拟。相变模型体现了汽泡底部微液层蒸发、近壁过热液体传热、汽泡顶部主流冷凝等多方面机制对汽泡生长的影响。模拟结果体现了汽泡底部微液层厚度的变化情况,与实验结果相吻合;微液层蒸发机制在汽泡生长初期对汽泡生长有较大影响,流道壁面效应对汽泡生长有显著影响。     

过冷流动沸腾中汽泡脱离直径的理论研究

汽泡脱离直径模型是壁面沸腾计算中的一个重要子模型。为了正确预测过冷流动沸腾中的壁面传热情况,研究结合新改进的汽泡生长模型,采用力平衡方法对过冷流动沸腾中的汽泡脱离直径进行了模拟。汽泡生长模型同时考虑了微液层、过热层和汽泡顶部过冷液体层对汽泡生长所做的贡献,并采用饱和沸腾与过冷沸腾2个实验对其进行了验证,结果表明预测曲线与实验值吻合良好。另外,选取了3个过冷流动沸腾实验来验证汽泡脱离直径模型,模拟结果均在可接受的误差范围内。     

矩形窄通道内汽泡滑移与冷凝前期生长模型研究

矩形窄流道内汽泡生长会直接改变相界面浓度,从而影响流道的传热传质性能。为获得适用于窄流道内不同类型的汽泡生长模型,基于通体可视的实验本体,开展壁面沸腾流动换热实验。基于传热能量方程,研究过冷沸腾中汽泡滑移与冷凝前期两种情况下汽泡生长模型。实验结果表明汽泡呈现两种形式的生长,即汽泡滑移生长以及冷凝前期生长。建立了两种情况下的汽泡生长模型,实验数据验证模型误差在20%以内。因此,本研究能为沸腾两相数值模拟提供更加精细化的汽泡生长模型,从而提高汽泡行为的预测精度。     

低压下水欠热流动沸腾的两相CFD数值模拟研究

采用两流体(汽相和液相)基本数学模型,结合汽相和液相之间的界面传热、传质和动量交换封闭模型、汽泡平均直径模型、汽泡脱离直径模型、汽泡成核模型、汽泡脱离频率模型、欠热沸腾起始点模型和壁面热流密度分配模型,在CFD软件CFX4.4中采用用户自定义函数将相变引起的传热、传质和动量交换作为源项分别添加到汽相和液相的能量、质量和动量守恒方程中,对低压下内管加热外管绝热的环形通道内的欠热沸腾进行了数值研究,得到了欠热流动沸腾下汽相体积份额、液相速度、汽相速度分布等。采用Lee等的环形通道内低压下欠热沸腾体积份额实验数据对计算结果进行了验证,吻合良好。     

微尺度核态沸腾汽泡聚合特性研究

从微观角度深入探究汽泡的生长、聚合及传热规律,采用数值仿真的方法,建立相变物理模型。将建立的模型嵌入计算流体动力学(CFD)软件,利用流体体积函数(VOF)相界面追踪方法,对汽泡的生长、运动、聚合及脱离进行探究。设立1个尺寸为3 mm×2 mm×2 mm的三维微型通道,底部壁面设置5个单元,用于加热并对汽泡底部热流密度进行监测。获得微通道加热条件下汽泡的合并运动形态和底面热流密度分布,并与实验观测值比较。通过对汽泡形态的模拟和底部热流密度的监测,建立汽泡动力学规律与底面热流密度分布的关系,进一步深化对汽泡合并的动力学规律和沸腾传热机理的认识。     

竖直窄流道内过冷流动沸腾的汽泡生长过程流场特性分析

采用VOF模型建立了过冷流动沸腾中的汽泡生长模型,使用UDF接口自编程序,对汽泡生长过程进行了数值模拟。考虑了主流速度、主流过冷度、壁面过热度、汽-液接触角等对汽泡生长过程的影响,获得了过冷流动沸腾条件下汽泡生长曲线,并与相应实验条件下的实验结果进行了对比,同时对汽泡的生长过程中的汽泡内外的流场进行了分析。     

蒸汽发生器传热管一、二次侧耦合换热及管外过冷沸腾数值研究

采用两流体欧拉数学模型,结合气相和液相之间的界面传热、传质和动量交换封闭模型以及RPI壁面沸腾模型,利用ANSYS CFX 12.0对蒸汽发生器局部传热管束二次侧的过冷沸腾进行数值研究。数值研究结果与单管内过冷沸腾实验数据对比验证符合良好。结果表明,采用壁面沸腾模型能准确预测沸腾起始点的位置,同时梅花孔板的存在对二次侧流动换热特性影响显著。     

多孔球层内沸腾现象与传热特性研究

采用池式沸腾实验系统,在常压底部加热条件下分别对由直径4、6、8mm玻璃球构建的多孔结构内沸腾过程进行了可视化研究.结果表明,过冷沸腾时,加热壁面上产生孤立汽泡,小汽泡可聚合为主汽泡,主汽泡脱离频率较低,汽相以分散的小汽泡为主;饱和沸腾初期,汽泡生长变快,主汽泡体积变大,连续汽相范围广阔;主汽泡形成频率随热流密度增加而增加;膜态沸腾时,底面被汽膜包围,液相占据球层空间.球体直径越大,产生同类现象需要的热流密度越大,传热系数的极值越大.饱和沸腾存在传热强化区和抑制区.直径4、8mm玻璃球构建的多孔介质传热系数随热流密度的增加而增加,6mm多孔介质则相反.     

管内过冷流动沸腾CFD模型参数敏感性研究

采用CFD方法对燃料组件进行过冷流动沸腾数值模拟研究是反应堆热工水力分析的一项重要内容。本研究使用STAR CCM+基于欧拉双流体模型结合壁面沸腾模型对管内过冷流动沸腾进行数值模拟,得到了壁面温度、主流温度及空泡份额的分布。基于实验结果对网格模型、湍流模型、壁面沸腾模型及相间作用力模型的参数设置进行了敏感性分析。研究结果表明,对于欧拉双流体模型,并非网格量越多结果越准确,加热面第1层网格的高度对结果影响显著。湍流模型和曳力模型对计算结果影响较小,非曳力中的湍流耗散力及升力对结果影响较大。Li Quan或Hibiki Ishii汽化核心密度模型与Kocamustafaogullari气泡脱离直径模型组合对壁面温度及空泡份额的计算较准确。本研究可为反应堆燃料组件内过冷流动沸腾数值模拟提供参考依据。     

窄通道内的汽泡核化以及滑移汽泡的影响

针对不同压力下窄通道内的过冷流动沸腾汽泡核化进行了实验研究,通过高速摄影技术观察了汽泡的核化和滑移过程。实验结果表明,与常规通道有所不同,窄通道内的汽泡一般不会离开加热壁面而产生浮升现象,汽泡主要沿加热壁面进行滑移运动。不同压力条件下的汽泡核化有较大的差别,较低压力条件下汽泡核化点沿加热壁面分布比较均匀,而压力升高后的汽泡核化点主要集中在沸腾起始点附近,下游核化点数目则相对较少。核化点分布形式的不同主要是由不同压力下汽泡滑移特性的不同所导致的。     

单气泡生长的高速激光干涉可视化研究

为研究过冷沸腾中气泡的生成机理并为壁面沸腾模型的建立提供数据支撑,本研究采用透光的氧化铟锡（ITO）薄膜作为加热材料,使用激光干涉法和高速相机对池沸腾中单气泡的生长过程进行同步测量,获取了气泡底部微液层的结构和变化过程,同时得到气泡干斑直径、气泡直径等相关参数。研究结果表明,气泡生长过程可分为底部存在微液层的生长阶段和底部完全蒸干后的脱离阶段两部分,两阶段的时长大致相当。在生长阶段,干斑直径和接触直径逐渐增大;在脱离阶段,接触直径逐渐减小。在气泡稳定生长时期,接触直径与气泡直径的比值约为0.6。实验结果验证了微液层边缘存在弯曲结构。本研究获得的气泡微液层演变数据可为气泡生长模型提供实验支撑。     

Experiments on microgravity boiling heat transfer by using transparent heaters

To clarify the relation between the liquid–vapor behavior and the heat transfer characteristics in the boiling phenomena, the structures of transparent heaters were developed for both flow boiling and pool boiling experiments and were applied to the microgravity environment realized by the parabolic flight of aircraft. In the flow boiling experiment, a transparent heated tube makes the heating, the observation of liquid–vapor behavior and the measurement of heat transfer data simultaneously possible. The heat transfer coefficient in the annular flow regime at moderate quality has distinct dependence on gravity provided that the mass velocity is not so high, while no noticeable gravity effect is seen at high quality and in the bubbly flow regime. The measured gravity effect was directly related to the behavior of annular liquid film observed through the transparent tube wall. In the pool boiling experiment, a structure of transparent heating surface realizes both the observation of the macrolayer or microlayer behavior from underneath and the measurements of local surface temperatures and the layer thickness. It was clarified in the microgravity experiments that no vapor stem exists but tiny bubbles are observed in the macrolayer underneath a large coalesced bubble at high heat flux. The heat flux evaluated by the heat conduction across the layer assumes less than 30% of the total to be transferred. The evaporation of the microlayers underneath primary bubbles just after the generation dominates the heat transfer in the microgravity, not only in the isolated bubble region but also in the coalesced bubble region.     

Study on characteristics of vapor–liquid two-phase flow in mini-channels

To explore the mechanism of boiling bubble dynamics in narrow channels, two types of channels are investigated which have I- and Z-shaped with width of 2 mm. Using VOF model and self-programming, the whole flow field is simulated with two different kinds of media, namely, water and ethanol. The influence of wall contact angle on the process of bubble generating and growth is studied, and the relationship between different channel shapes and the pressure drop is also investigated taking into account the effects of gravity, viscosity, surface tension and wall adhesion. The bubble generation, growth and departure processes are analyzed through numerical simulation and self-programming, and the influence of interface movements and changes on internal pressure difference and average surface heat transfer coefficient is investigated by using geometry reconstruction and interface tracking. It is found that wall contact angle has a great influence on the morphology of bubble. The smaller the wall contact angles are, the more round the bubbles are, and the less time the bubbles take to depart from the wall. The variation of contact angle also has effect upon the heat transfer coefficient. The greater the wall contact angle is, the larger the bubble-covered area is, thus the wall thermal resistance gets higher, and bubble nucleation is suppressed, and the heat transfer coefficient becomes lower. The role of surface tension in the process of boiling heat transfer is much more important than the gravity in narrow channels. The generation of bubbles dramatically disturbs the boundary layer, and the bubble bottom micro-layer can enhance the heat transfer. The heat transfer coefficient of Z-shaped channels is larger than that of I-shaped ones, while the pressure drop of the former is obviously higher. In addition, surface tension and viscosity significantly impact the pressure drop of boiling system, and different specific heat and boiling point values result in different heat transfer coefficients. The simulation results in this paper match well with the experimental data revealed in other sources, both show that the heat transfer coefficient of water is higher than that of ethanol and Z-shaped channels have better heat transfer capability.     

An investigation of transition boiling mechanisms of subcooled water under forced convective conditions

A mechanistic model for forced convective transition boiling has been developed to investigate transition boiling mechanisms and to predict transition boiling heat flux realistically. This model is based on a postulated multi-stage boiling process occurring during the passage time of the elongated vapor blanket specified at a critical heat flux (CHF) condition. Between the departure from nucleate boiling (DNB) and the departure from film boiling (DFB) points, the boiling heat transfer is established through three boiling stages, namely, the macrolayer evaporation and dryout governed by nucleate boiling in a thin liquid film and the unstable film boiling characterized by the frequent touches of the interface and the heated wall. The total heat transfer rates after the DNB is weighted by the time fractions of each stage, which are defined as the ratio of each stage duration to the vapor blanket passage time. The model predictions are compared with some available experimental transition boiling data. The parametric effects of pressure, mass flux, inlet subcooling on the transition boiling heat transfer are also investigated. From these comparisons, it can be seen that this model can identify the crucial mechanisms of forced convective transition boiling, and that the transition boiling heat fluxes including the maximum heat flux and the minimum film boiling heat flux are well predicted at low qualities/high pressures near 10 bar. In future, this model will be improved in the unstable film boiling stage and generalized for high quality and low pressure situations.     

Numerical simulation of bubble behaviors in subcooled flow boiling under swing motion

A numerical investigation of bubble behaviors in subcooled flow boiling of water under the effect of additional inertial forces has been performed considering energy and mass transfer during phase change based on the VOF (volume-of-fluid) method. The pressure ranges from 0.1 to 1.0 MPa, and heat flux from 200 to 500 kW/m². The mass flow rate and inlet subcooling are specified at 320 kg/m² s and 10 K, respectively. The liquid-vapor interface is captured using the piecewise linearity interpolation calculation (PLIC) geometry restructuring method. The simulations are carried out on upward water flow in a vertical, rectangular duct with single side heating surface. The pressure, velocity vector and temperature distribution around two isolated bubbles are studied firstly. The behaviors of bubble coalescence, sliding, detachment from the heated wall, and the bubble shape variation during lifetime are further examined. The bubble behaviors in the different pressure and heat flux are investigated. The simulated results of bubble growth rate and wall temperature are agreed well with the correlations in the literatures. The additional inertial forces caused by swing are negligible, but the fluctuation of mass flow rate caused by swing motion influences the forces acting on bubble significantly. Compared with the motionless condition, the pressure drop is increased and the fluctuation becomes acute as heat flux increases under the swing condition.     

The nature of bubble growth under different system pressures in a narrow channel

Bubble growth in a vertical narrow channel under different system pressures (1-10 bar) was photographically studied. It was found that bubble growth rates and bubble sizes decrease with increasing system pressures. It was also found that bubbles grow at nucleate sites and begin to shrink due to condensation after departing from nucleate sites under lower system pressures (1-3 bar). However, bubbles keep growing along the heating wall under higher system pressures (6-10 bar). All bubble growth curves under different system pressures can be predicted by power curve model in the present study. A theoretical analysis on the effect of system pressure on bubble growth shows that the latent heat needed for a bubble with unit volume and the heat transfer on the heating wall are quite different under different system pressures, which results in distinct difference in bubble growth under different system pressure.     

基于CFD方法的矩形窄缝通道内过冷流动沸腾模型研究

板状燃料元件中的矩形窄缝通道具有宽高比大的几何特征,高度方向速度梯度大、分布陡峭,发生过冷沸腾时,近壁面汽泡运动行为将受其影响而改变,其中汽泡滑移现象对沸腾换热影响较大。本文针对矩形窄缝通道中的汽泡滑移行为,构建了包含滑移热流的壁面热流分配模型,并建立机理性的汽泡受力模型和滑移模型计算汽泡脱离直径、浮升直径和滑移距离等辅助参数,开发了一套适用于矩形窄缝通道内向上流动沸腾的壁面沸腾模型。选用Nuthel窄缝通道沸腾实验进行数值模拟验证,结果表明：本文模型可以较好地预测1～4 MPa中低压工况窄缝通道向上流动沸腾的壁面过热度,最大误差相比RPI模型由80%降低至17%;蒸发热流份额和近壁面空泡份额相比RPI模型更低。     

过冷沸腾自然对流两相CFD模拟及应用

采用计算流体力学（CFD）方法,开展过冷沸腾自然对流两相模拟与应用研究。对侧壁加热圆柱水箱过冷沸腾自然对流实验采用两相CFD瞬态模拟,模拟时间为1 500 s,通过模型设置与模拟方法研究,再现了过冷沸腾发生后实验的温度阶跃,得到与实验较一致的温度分布、气泡产生时间与产生位置,确保了数值计算的合理性与准确性。在此基础上,对以欧洲ESBWR（经济简化沸水堆）非能动安全壳冷却系统（PCCS）为原型的ISP-42实验进行了两相CFD模拟,获得与实验一致的温度分布,确定采用两相CFD数值模拟对非能动安全壳冷却系统及非能动余热排出系统进行应用研究可行,为下一步计算传热系数、构建自然对流传热模型建立了良好基础。该项研究对工程应用中探寻非能动安全壳冷却系统及非能动余热排出系统的两相自然循环传热特性具有较大价值。