Flow Patterns and Heat Transfer for Flow Boiling in Small to Micro Diameter Tubes

An overview of the recent developments in the study of flow patterns and boiling heat transfer in small to micro diameter tubes is presented. The latest results of a long-term study of flow boiling of R134a in five vertical stainless-steel tubes of internal diameter 4.26, 2.88, 2.01, 1.1, and 0.52 mm are then discussed. During these experiments, the mass flux was varied from 100 to 700 kg/m²s and the heat flux from as low as 1.6 to 135 kW/m². Five different pressures were studied, namely, 6, 8, 10, 12, and 14 bar. The flow regimes were observed at a glass section located directly at the exit of the heated test section. The range of diameters was chosen to investigate thresholds for macro, small, or micro tube characteristics. The heat transfer coefficients in tubes ranging from 4.26 mm down to 1.1 mm increased with heat flux and system pressure, but did not change with vapor quality for low quality values. At higher quality, the heat transfer coefficients decreased with increasing quality, indicating local transient dry-out, instead of increasing as expected in macro tubes. There was no significant difference between the characteristics and magnitude of the heat transfer coefficients in the 4.26 mm and 2.88 mm tubes but the coefficients in the 2.01 and 1.1 mm tubes were higher. Confined bubble flow was first observed in the 2.01 mm tube, which suggests that this size might be considered as a critical diameter to distinguish small from macro tubes. Further differences have now been observed in the 0.52 mm tube: A transitional wavy flow appeared over a significant range of quality/heat flux and dispersed flow was not observed. The heat transfer characteristics were also different from those in the larger tubes. The data fell into two groups that exhibited different influences of heat flux below and above a heat flux threshold. These differences, in both flow patterns and heat transfer, indicate a possible second change from small to micro behavior at diameters less than 1 mm for R134a.     

换热表面结垢对U型管蒸汽发生器传热性能的影响分析

U型管蒸汽发生器的壳侧沉积了来自二回路系统中的腐蚀产物,结垢导致热量聚积在金属换热管上,容易造成垢下热点腐蚀,危害设备安全。为了明确结垢对蒸汽发生器传热性能的影响,本研究基于仿真平台APROS建立了U型管蒸汽发生器的分布式模型,并根据已公开论文中的数据进行了模型准确性验证;推导了污垢热阻与表面换热系数之间的关系式,分析了不同结垢厚度、位置对U型管蒸汽发生器换热区域的传热管壁面温度、流体温度、传热系数、热流密度等的影响程度。研究结果表明：随着结垢程度的加剧,蒸汽发生器的换热效率不断降低,出口蒸汽品质不断下降;结垢对沸腾段换热效率的影响比对过冷段换热效率的影响更大。     

The effect of liquid film evaporation on flow boiling heat transfer in a micro tube

Flow boiling in micro channels is attracting large attention since it leads to large heat transfer area per unit volume. Generated vapor bubbles in micro channels are elongated due to the restriction of channel wall, and thus slug flow becomes one of the main flow regimes. In slug flow, sequential bubbles are confined by the liquid slugs, and thin liquid film is formed between tube wall and bubble. Liquid film evaporation is one of the main heat transfer mechanisms in micro channels and liquid film thickness is a very important parameter which determines heat transfer coefficient. In the present study, liquid film thickness is measured by laser focus displacement meter under flow boiling condition and compared with the correlation proposed for an adiabatic flow. The relationship between liquid film thickness and heat transfer coefficient is also investigated. Initial liquid film thickness under flow boiling condition can be predicted well by the correlation proposed under adiabatic condition. Under flow boiling condition, liquid film surface fluctuates due to high vapor velocity and shows periodic pattern against time. Frequency of periodic pattern increases with heat flux. At low quality, heat transfer coefficients calculated from measured liquid film thickness show good accordance with heat transfer coefficients obtained directly from wall temperature measurements.     

Correlation for boiling heat transfer of R-134a in horizontal tubes including effect of tube diameter

A Chen-type correlation for flow boiling heat transfer of R-134a in horizontal tubes was modified taking into account the effect of tube diameter. The effect of tube diameter on flow boiling heat transfer coefficient was characterized by the Weber number in gas phase. Results showed that this correlation could be applied to a wide range of tube diameters (0.5–11-mm-ID). In addition, the dryout point and the heat transfer characteristics after the dryout point were also investigated based on the annular flow model. The proposed experimental expressions to predict both the dryout quality and the post-dryout heat transfer coefficient could also be applied to a wide range of tube diameter (0.5–11-mm-ID).     

Thermal performance simulation of a solar cavity receiver under windy conditions

Solar cavity receiver plays a dominant role in the light-heat conversion. Its performance can directly affect the efficiency of the whole power generation system. A combined calculation method for evaluating the thermal performance of the solar cavity receiver is raised in this paper. This method couples the Monte-Carlo method, the correlations of the flow boiling heat transfer, and the calculation of air flow field. And this method can ultimately figure out the surface heat flux inside the cavity, the wall temperature of the boiling tubes, and the heat loss of the solar receiver with an iterative solution. With this method, the thermal performance of a solar cavity receiver, a saturated steam receiver, is simulated under different wind environments. The highest wall temperature of the boiling tubes is about 150 °C higher than the water saturation temperature. And it appears in the upper middle parts of the absorbing panels. Changing the wind angle or velocity can obviously affect the air velocity inside the receiver. The air velocity reaches the maximum value when the wind comes from the side of the receiver (flow angle α = 90°). The heat loss of the solar cavity receiver also reaches a maximum for the side-on wind.     

表面强化管外池沸腾传热特性研究

为实现节能降耗,开发了多种强化沸腾传热的高效换热管。以水为工质,在0.1MPa下对垂直光管、烧结多孔管和T槽管进行了池沸腾传热实验研究,并分析了沿管子轴向的温度分布。实验结果表明,烧结多孔管与T槽管能显著降低起始沸腾过热度、强化沸腾传热:烧结多孔管和T槽管的起始沸腾过热度比光管的低1.5K左右;烧结多孔管和T槽管的核态沸腾传热系数分别为光管的2.4～3.2倍和1.6～2.0倍。此外,烧结多孔管和T槽管能降低相同热流密度下的壁面温度,且有利于降低管子轴向的温差。     

Concerning the magnitude of the maximum heat flux and the mechanisms of superintensive bubble boiling

We have developed a stabilizer of the temperature of a thermoresistor wire electric heater based on a PID controller. Using this stabilizer, we investigated heat exchange of subcooled water in pool boiling. We found that on stabilization of the heater temperature up to that of the subcooled water, transition from convection to the regime of bubble boiling and vice versa occurs spontaneously and is accompanied by a jumpwise change in heat transfer. It is shown that in the regime of stable bubble boiling, the law of heat transfer is independent of the liquid temperature and the heater diameter and that the maximum heat loading may attain 50 MW/m², which is much above the values cited earlier in the literature. Based on the results obtained, a mechanism of implementation of bubble boiling for the regimes of a constant heat flux and a constant temperature is suggested. The assumption is made that the regime of heterogeneous vapour generation is possible only in the case of the heater constant temperature. In the regime of a stabilized heat flux on the heater, the spatially inhomogeneous regime of heat transfer is established. This regime represents a spatially distributed combination of three regimes: convective heat transfer, homogeneous boiling, manifesting itself in periodic boiling-up of overheated layers of the liquid near the surface and an unstable regime of heterogeneous vapour generation.     

Effect of tube bundles on nucleate boiling and critical heat flux

In order to elucidate boiling heat transfer characteristics for each tube and the critical heat flux (CHF) for tube bundles, an experimental investigation of pool and flow boiling of Freon-113 at 0.1 MPa was performed using two typical tube arrangements. A total of fifty heating tubes of 14 mm diameter, equipped with thermocouples and cartridge heaters, were arrayed at pitches of 18.2 and 21.0 mm to simulate both square in-line and equilateral staggered bundles. For the flow boiling tests the same bundles as were used in pool boiling were installed in a vertical rectangular channel, to which the fluid was supplied with an approach velocity varying from 0.022 to 0.22 m/s. It was found in this study that the boiling heat transfer coefficient of each tube in a bundle was higher than that for an isolated single tube in pool boiling. This enhancement increases for tubes at higher locations, but decreases as heat flux is increased. At heat fluxes exceeding certain values, the heat transfer coefficient becomes the same as that for an isolated tube. As the heat flux approaches the CHF, flow pulsations occurred in the pool boiling experiments although the heat transfer coefficient was invariant even under this situation. The approach velocity has an appreciable effect on heat transfer up to a certain level of heat flux. In this range of heat flux, the heat transfer coefficient exceeds the values observed for pool boiling. An additive method with two contributions, i.e., single phase convection and boiling, was used to predict the heat transfer coefficient for bundles. The predicted results showed reasonable agreement with the measured results. The critical heat flux in tube bundles tended to increase as more bubbles were rising through the tube clearance. © 1998 Scripta Technica, Heat Trans Jpn Res, 27(4): 312–325, 1998     

Convective boiling heat transfer and two-phase flow characteristics inside a small horizontal helically coiled tubing once-through steam generator

The pressure drop and boiling heat transfer characteristics of steam-water two-phase flow were studied in a small horizontal helically coiled tubing once-through steam generator. The generator was constructed of a 9-mm ID 1Cr18Ni9Ti stainless steel tube with 292-mm coil diameter and 30-mm pitch. Experiments were performed in a range of steam qualities up to 0.95, system pressure 0.5-3.5 MPa, mass flux 236-943 kg/m²s and heat flux 0-900 kW/m². A new two-phase frictional pressure drop correlation was obtained from the experimental data using Chisholm’s B-coefficient method. The boiling heat transfer was found to be dependent on both of mass flux and heat flux. This implies that both the nucleation mechanism and the convection mechanism have the same importance to forced convective boiling heat transfer in a small horizontal helically coiled tube over the full range of steam qualities (pre-critical heat flux qualities of 0.1-0.9), which is different from the situations in larger helically coiled tube where the convection mechanism dominates at qualities typically >0.1. Traditional single parameter Lockhart-Martinelli type correlations failed to satisfactorily correlate present experimental data, and in this paper a new flow boiling heat transfer correlation was proposed to better correlate the experimental data.     

Flow boiling in horizontal flattened tubes: Part II – Flow boiling heat transfer results and model

Experiments of flow boiling heat transfer were conducted in four horizontal flattened smooth copper tubes of two different heights of 2 and 3 mm. The equivalent diameters of the flattened tubes are 8.6, 7.17, 6.25, and 5.3 mm. The working fluids were R22 and R410A. The test conditions were: mass velocities from 150 to 500 kg/m² s, heat fluxes from 6 to 40 kW/m² and saturation temperature of 5 °C. The experimental heat transfer results are presented and the effects of mass flux, heat flux, and tube diameter on heat transfer are analyzed. Furthermore, the flow pattern based flow boiling heat transfer model of Wojtan et al. [L. Wojtan, T. Ursenbacher, J.R. Thome, Investigation of flow boiling in horizontal tubes: Part I – A new diabatic two-phase flow pattern map, Int. J. Heat Mass Transfer 48 (2005) 2955–2969; L. Wojtan, T. Ursenbacker, J.R. Thome, Investigation of flow boiling in horizontal tubes: Part II – Development of a new heat transfer model for stratified-wavy, dryout and mist flow regimes, Int. J. Heat Mass Transfer 48 (2005) 2970–2985], using the equivalent diameters, were compared to the experimental data. The model predicts 71% of the entire database of R22 and R410A ±30% overall. The model predicts well the flattened tube heat transfer coefficients for R22 while it does not predicts well those for R410A. Based on several physical considerations, a modified flow boiling heat transfer model was proposed for the flattened tubes on the basis of the Wojtan et al. model and it predicts the flattened tube heat transfer database of R22 and R410A by 85.8% within ±30%. The modified model is applied to the reduced pressures up to 0.19.     

一种避免重沸器在不稳定过渡沸腾区操作的控制方法

当精馏塔釜温度较高、且蒸汽温位过低而不能用做重沸器热源时,炼油厂通常采用高温导热油、工艺液体等不发生相变的介质作为精馏塔重沸器热源.当这些加热介质温度高至使重沸器传热处于过渡沸腾区时,重沸器传热速率则会随传热温差的增大而下降,造成精馏控制不稳定及产品质量不合格.本文结合重沸器传热过程原理,提供了一种可以避免重沸器在不稳定的过渡沸腾区传热的控制方法,即“加热介质部分循环、部分补充”的控制方法.该方法将来自重沸器出口的、温度较低的加热介质,与来自系统的、温度较高的加热介质按一定比例混合,以降低重沸器入口的加热介质温度,确保重沸器在稳定的泡核沸腾区操作,并使精馏塔在平稳的状态下运行.应用实例表明,应用本方法,系统总能耗增加很少,甚至会有所降低,精馏塔的操作性能得到明显改善,产品质量大幅提高.     

Modeling and numerical investigation of hydrogen nucleate flow boiling heat transfer

Liquid hydrogen flow boiling heat transfer in tubes is of great importance in the hydrogen applications such as superconductor cooling, hydrogen fueling. In the present study, a numerical model for hydrogen nucleate flow boiling based on the wall partition heat flux model is established. The key parameters in the model such as active nucleation site density, bubble departure diameter and frequency are carefully discussed and determined to facilitate the modeling and simulation of hydrogen flow boiling. Simulation results of the numerical model show reasonably well agreement with experimental data from different research groups in a wide operation condition range with the means absolute error (MAE) of 10.6% for saturated and 5.3% for subcooled flow boiling. Based on the model, wall heat flux components and void fraction distribution of hydrogen flow boiling are studied. Effects of mass flow rate and wall heat flux on the flow boiling heat transfer performance are investigated. It is found that in the hydrogen nucleate flow boiling, the predominated factor is the Boiling number, rather than the vapor quality. A new simple correlation is proposed for predicting hydrogen saturated nucleate flow boiling Nusselt number. The MAE between the correlation predicted and experimentally measured Nusselt number is 13.6% for circular tubes and 12.5% for rectangular tubes. The new correlation is applicable in the range of channel diameter 4–6.35 mm, Reynolds number 64000–660,000, saturation temperature 22–29 K, Boiling number 8.37 × 10^?5–2.33 × 10^?3.     

The Effect of Tube Flattening on Flow Boiling Heat Transfer Enhancement

The present study investigated the effect of smooth tube flattening on heat transfer enhancement in an evaporator. The tubes with internal diameter of 8.7 mm were flattened into an oblong shape with different inside heights. The test setup was basically a vapor compression refrigeration system equipped with all necessary measuring instruments. Refrigerant R-134a flowing inside the tube was heated by an electrical coil heater wrapped around it. The ranges of mass velocities were from 74 to 106 kg/m²-s and vapor quality varied from 25% to 95%. Analysis of the collected data indicated that the heat transfer coefficient elevates by increasing the mass velocity and vapor quality in flattened tubes just like the round tube. The flow boiling heat transfer coefficient increases when the flattened tube is used instead of the round tube. The highest heat transfer coefficient enhancement of 172% was achieved for the tube with the lowest inside height at mass velocity of 106 kg/m²-s and vapor quality of 85%. Finally, based on the present experimental results, a correlation was developed to predict the heat transfer coefficient in flattened tubes.     

New prediction methods for CO2 evaporation inside tubes: Part II—An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO₂ using the Cheng–Ribatski–Wojtan–Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082–4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: “New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes” [Heat Mass Transfer 49 (21–22) (2006) 4082–4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO₂ was developed based on the CO₂ data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO₂ covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m² s, heat fluxes from 1.8 to 46 kW/m² and saturation temperatures from ?28 to 25 °C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within ±30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes.     

Experimental and theoretical studies on subcooled flow boiling of pure liquids and multicomponent mixtures

To improve the design of modern industrial reboilers, accurate knowledge of boiling heat transfer coefficients is essential. In this study flow boiling heat transfer coefficients for binary and ternary mixtures of acetone, isopropanol and water were measured over a wide range of heat flux, subcooling, flow velocity and composition. The measurements cover the regimes of convective heat transfer, transitional boiling and fully developed subcooled flow boiling. Two models are presented for the prediction of flow boiling heat transfer coefficients. The first model is the combination of the Chen model with the Gorenflo correlation and the Schlünder model for single and multicomponent boiling, respectively. This model predicts flow boiling heat transfer coefficients with acceptable accuracy, but fails to predict the nucleate boiling fraction NBF reasonably well. The second model is based on the asymptotic addition of forced convective and nucleate boiling heat transfer coefficients. The benefit of this model is a further improvement in the accuracy of flow boiling heat transfer coefficient over the Chen type model, simplicity and the more realistic prediction of the nucleate boiling fraction NBF.     

Condensation in a vertical tube bundle passive condenser – Part 2: Complete condensation

An experimental study of the tube bundle effect on heat removal capabilities in complete condensation mode of a passive condenser was performed. A full scale test section, with four condenser tubes, was designed and constructed to simulate operating conditions of a passive containment cooling system. For complete condensation analysis, pure steam was supplied to the test section and heat transfer properties were measured for pressure from 100 to 280 kPa. The condensation heat transfer results were similar to the findings from single tubes, except for a slightly higher condensate mass flux. This was determined to be a result of turbulent mixing in the secondary boiling water caused by the tube bundle.     

Cooling Augmentation with Microchanneled Structures

Experiments were conducted to investigate the heat transfer characteristics and cooling performance of subcooled liquid, water, flowing through rectangular cross-section microchanneled structures machined on a stainless steel plate. Heat transfer or flow mode transition was observed when the heating rate or wall temperature was increased. This transition was found to be suggestively induced by the variation in liquid thermophysical properties due to the significant rise of liquid temperature in the microstructures. The influence of such parameters as liquid velocity, subcooling, property variation, and microchannel geometric configuration on the heat transfer behavior, cooling performance and the heat transfer and liquid flow mode transition were also investigated. The experiments indicated that both slngle-phase forced convection and flow boiling characteristics were quite different from those in normal-sized tubes and the heat transfer was obviously intensified.     

Heat transfer and pressure drop during flow boiling of pure refrigerants and refrigerant/oil mixtures in tube with porous coating

The experimental stand and procedure for flow boiling investigations are described. Experimental data for pure R22, R134a, R407C and their mixtures with polyester oil FUCHS Reniso/Triton SEZ 32 in a tube with porous coating and smooth, stainless steel reference tube are presented. Mass fraction of oil was equal to 1% or 5%. During the tests inlet vapour quality was set at 0 and outlet quality at 0.7. Mass velocity varied from about 250 to 500 kg/m²s. The experiments have been conducted for average saturation temperature 0 °C. In the case of flow boiling of pure refrigerants, the application of a porous coating on inner surface of a tube results in higher average heat transfer coefficient and simultaneously in lower pressure drop in comparison with the flow boiling in a smooth tube for the same mass velocity. Correlation equation for heat transfer coefficient calculation during the flow boiling of pure refrigerants inside a tube with porous coating has been proposed.     

缩放管管外流动沸腾换热的数值模拟与场协同分析

本文利用FLUENT软件对制冷剂R134a在光管和缩放管水平管外沸腾传热进行了三维数值模拟,得到了其饱和泡状沸腾过程中体积含汽率的分布规律,并比较了它们的换热系数,结果表明:缩放管外侧能够很好地强化沸腾传热.此外,通过改变边界条件分析了质量流量、热流密度的变化对缩放管管外沸腾换热系数的影响.最后应用场协同理论,从局部换热角度分析其强化机理.研究表明:缩放管水平管外沸腾换热得到强化的原因是其凹槽前后的速度场与温度梯度场之间夹角更小,协同程度更好.     

Steam condensing – liquid CO2 boiling heat transfer in a steam condenser for a new heat recovery system

In a new waste heat recovery system, waste heat is recovered from steam condensers through cooling by liquid CO₂ instead of seawater, taking advantage of effective boiling heat transfer performance; the heat is subsequently used for local heat supply. The steam condensing – liquid CO₂ boiling heat transfer performance in a steam condenser with a shell and a helical coil non-fin tube was studied both numerically and experimentally. A heat transfer numerical model was constructed from two models developed for steam condensation and for liquid CO₂ boiling. Experiments were performed to verify the model at a steam pressure range of 3.2–5 kPa and a CO₂ saturation pressure range of 5–6 MPa. Overall heat transfer coefficients obtained from the numerical model agree with the experimental data within ±5%. The numerical estimations show that the boiling local heat transfer coefficient reaches a maximum value of 26 kW/m² K. This value is almost one order higher than that of a conventional water-cooled condenser.