Marangoni heat transfer in subcooled nucleate pool boiling

The liquid motion induced by surface tension variation, termed the Marangoni effect, and its contribution to boiling heat transfer has been an issue of much controversy. Boiling heat transfer theory, although acknowledging its existence, considers its contribution to heat transfer to be insignificant in comparison with buoyancy induced convection. However, recent microgravity experiments have shown that although the boiling mechanism in a reduced gravity environment is different, the corresponding heat transfer rates are similar to those obtained under normal gravity conditions, raising questions about the validity of the assumption. An experimental investigation was performed in which distilled water was gradually heated to boiling conditions on a copper heater surface at four different levels of subcooling. Photographic investigation of the bubbles appearing on the surface was carried out in support of the measurements. The results obtained indicate that Marangoni convection associated with the bubbles formed by the air dissolved in the water which emerged from solution when the water was heated sufficiently, significantly influenced the heat transfer rate in subcooled nucleate pool boiling. A heat transfer model was developed in order to explain the phenomena observed.     

Investigation of Enhanced Boiling Heat Transfer from Porous Surfaces

ＳｕｒｆａｃｅｓＩｎｖｅｓｔｉｇａｔｉｏｎｏｆＥｎｈａｎｃｅｄＢｏｉｌｉｎｇＨｅａｔＴｒａｎｓｆｅｒｆｒｏｍＰｏｒｏｕｓＳｕｒｆａｃｅｓ￥ＬｉｎＺｈｉｐｉｎｇ；ＭａＴｏｎｇｚｅ；ＺｈａｎｇＺｈｅｎｇｆａｎｇ（ＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎ...     

Capillary column enhanced CHF microgravity flow boiling

Flow boiling heat transfer under microgravity conditions can be extended and enhanced by means of using porous stacks, or capillary columns, arranged on top of a flat heated surface. Under these conditions, body forces are negligible to remove the generated vapor away from the hot surface, which eventually hinders liquid from reaching it. It is possible to increase the critical heat flux (CHF) by having porous stacks symmetrically arranged on this surface; which draws the liquid phase towards it by means of capillary forces. Various flow regimes in the capillary enhanced surface flow boiling can be identified. These include: the regime where the liquid is supplied between the columns, the regime where the liquid flow is controlled by liquid capturing and the viscous drag-capillarity in the columns, and the critical heat flux. For the theoretical model, the expression for the interfacial lift-off model critical heat flux was interpreted based on customizable parameters instead of those imposed by the physics of the flow. This study indicates a potential improvement in CHF by having an inter-column spacing smaller than the critical wavelength for a plain surface. There is also a potential benefit of having the wetting contact to wavelength ratio to be larger than the constant of 0.2 found in experimental studies. The CHF regime can occur by a limitation of the stacks to have access to the liquid phase, as it happens when they are completely submerged in a vapor phase, or by reaching the maximum capillary pressure drop in the stack (as per the Darcy–Ergun momentum equation), or by reaching an entrainment limit of the vapor flow passed the capillary columns. Therefore the critical heat flux can also be extended as long as the capillary columns protrude over the vapor layer and their viscous capillary and entrainment limits are not reached.     

Effect of microgravity on flow boiling heat transfer of liquid hydrogen in transportation pipes

The flow boiling phenomenon of liquid hydrogen (LH2) during transportation in microgravity is very different from that under terrestrial condition. In this study, a saturated flow boiling of LH2 in a horizontal tube has been simulated under microgravity condition using coupled level-set and volume of fluid method. The validation of the developed model shows good agreement with the experimental data from the literature. The changes of heat fluxes and pressure drops under different gravitational accelerations were analyzed. And, the variation of heat fluxes with different wall superheat and contact angle were compared between microgravity (10⁻⁴g) and normal gravity (1g) condition. Also, the influence of surface tension were studied under microgravity. The numerical results indicate that the heat flux decrease with the decrement of gravitational acceleration. And the heat transfer ratio decrease with the increment of wall superheat in the nucleate boiling regime. The heat transfer slightly reduce when considering surface tension. In addition, the changes of contact angle have a more significant impact on heat transfer under microgravity condition.     

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

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

An experimental study of void fraction and heat transfer under upward and downward flow boiling conditions

Heat transfer coefficient and actual void fraction have been measured during upflow and downflow boiling of water in an annular channel. At the same values of pressure, mass flux, heat flux and flow quality significant difference of void fraction has been established in upflow and downflow. The upflow and downflow heat transfer coefficients did not deviate significantly from each other, if compared at identical values of pressure, mass flux, heat flux and actual void fraction.     

Pool boiling heat transfer of functionalized nanofluid under sub-atmospheric pressures

A water-based functionalized nanofluid was made by surface functionalizing the ordinary silica nanoparticles. The functionalized nanoparticles were water-soluble and could still keep dispersing well even at the mass concentration of 10% and no sedimentation was observed. An experimental study was carried out to investigate the pool boiling heat transfer characteristics of functionalized nanofluid at atmospheric and sub-atmospheric pressures. The same work was also performed for DI water and traditional nanofluid consisted of water and ordinary silica nanoparticles for the comparison. Experimental results show that there exist great differences between pool boiling heat transfer characteristics of functionalized and traditional nanofluid. The differences mainly result from the changes of surface characteristics of the heated surface during the boiling. A porous deposition layer exists on the heated surface during the boiling of traditional nanofluid; however, no layer exists for functionalized nanofluid. Functionalized nanofluid can slightly increase the heat transfer coefficient comparing with the water case, but has nearly no effects on the critical heat flux. It is mainly due to the changes of the thermoproperties of nanofluids. Traditional nanofluid can significantly enhance the critical heat flux, but conversely deteriorates the heat transfer coefficient. It is mainly due to effect of surface characteristics of the heated surface during the boiling. Therefore, the pool boiling heat transfer of nanofluids is governed by both the thermoproperties of nanofluids and the surface characteristics of the heated surface.     

Pressure drop across micro-pin heat sinks under unstable boiling conditions

Flow boiling was investigated under unstable boiling conditions in three different micro-pin fin heat sinks using water and R-123 as working fluids. Once boiling was initiated severe temperature fluctuations were recorded for all the tested (three) micro-pin fin heat sinks.Flow images and fast-Fourier transform (FFT) of pressure signals during flow boiling were used to explain experimental results. The boiling instability mechanisms behind unstable boiling were discussed for both water and R-123. Accordingly, no significant pressure fluctuations with respect to time averaged pressure drop were evident for the tested micro-pin fin heat sinks before and after flow boiling instability initiates. However, a step change in the pressure signals were recorded with the inception of unstable boiling, and a sharp increase in the magnitude peaks of the FFT profiles was observed in the device operated with R-123, while there was no significant change in the FFT profiles in the devices operated with water. According to complementary flow visualization studies, the oscillation frequency of the periodic flow patterns for the device operated with R-123 was higher (f～80 Hz) than that of the devices operated with water (f～20 Hz).     

CFD simulation of heat transfer and phase change characteristics of the cryogenic liquid hydrogen tank under microgravity conditions

The heat transfer and phase change processes of cryogenic liquid hydrogen (LH₂) in the tank have an important influence on the working performance of the liquid hydrogen-liquid oxygen storage and supply system of rockets and spacecrafts. In this study, we use the RANS method coupled with Lee model and VOF (volume of fraction) method to solve Navier-stokes equations. The Lee model is adopted to describe the phase change process of liquid hydrogen, and the VOF method is utilized to calculate free surface by solving the advection equation of volume fraction. The model is used to simulate the heat transfer and phase change processes of the cryogenic liquid hydrogen in the storage tank with the different gravitational accelerations, initial temperature, and liquid fill ratios of liquid hydrogen. Numerical results indicate greater gravitational acceleration enhances buoyancy and convection, enhancing convective heat transfer and evaporation processes in the tank. When the acceleration of gravity increases from 10^?2 g0 to 10^?5 g0, gaseous hydrogen mass increases from 0.0157 kg to 0.0244 kg at 200s. With the increase of initial liquid hydrogen temperature, the heat required to raise the liquid hydrogen to saturation temperature decreases and causes more liquid hydrogen to evaporate and cools the gas hydrogen temperature. More cryogenic liquid hydrogen (i.e., larger the fill ratio) makes the average fluid temperature in the tank lower. A 12.5% reduction in the fill ratio resulted in a decrease in fluid temperature from 20.35 K to 20.15 K (a reduction of about 0.1%, at 200s).     

Experimental study on Marangoni effect induced by heat and mass transfer

Surface tension gradients along a fluid–fluid interface provoke strong convective activity, such effect is called Marangoni effect. The Marangoni effect might be induced by heat and mass transfer. In both cases, instabilities develop at the interface. The importance of Marangoni effect has led to many investigations over several years [1], [2], [3], [4], [5], [6] [L. Scriven and C. Sternling, Nature, 187, pp. 186–188 (1960), H. Lu, Y. Yang and J. Maa, Ind. Eng. Chem. Res., 35, pp. 1921–1928 (1996), V. Kaminsky, A. Vyaz'min, N. Kulov and V. Dil'man, Chem. Eng. Sci., 53, No. 19 (1998), P. Lyford, H. Pratt, F. Greiser and D. Shallcross, Can. J. Chem. Eng., 76 (1998a), P. Lyford, H. Pratt, F. Greiser and D. Shallcross, Can. J. Chem. Eng., 76 (1998b), H. Lu, Y. Yang and J. Maa, Ind. Eng. Chem. Res., 36, pp. 474–482 (1997)]. The aim of this work is to study experimentally the Marangoni effect induced by heat and mass transfer, and to determine the Marangoni number (Ma) in both cases. Four aliphatic alcohol were studied, from C1 to C4. In the experiments, interfacial instabilities were observed and Ma was calculated for different temperatures and pressures. It was observed that Ma decreases as the carbon number atoms (n) increase, this behavior can be explained through the changes of some properties (surface tension, diffusivity) as heat and mass transfer occurred. In heat transfer experiments, the highest temperature gradients were reached for methanol, the Bond number (Bd) was also calculated and it was found that the natural convection predominated. In mass transfer experiments, CO₂ was used as gas phase, and it was observed that as pressure increases the Ma increases, it might be explained by the decrease of CO₂ diffusion in hydrocarbons as pressure increases. In conclusion, the Marangoni effect was observed for aliphatic alcohol, the influence of temperature and pressure was also observed, and finally, the Marangoni number decreases as n increases.     

Correlation for flow boiling heat transfer in mini-channels

In view of practical significance of a correlation of heat transfer coefficient in the aspects of such applications as engineering design and prediction, some efforts towards correlating flow boiling heat transfer coefficients for mini-channels have been made in this study. Based on analyses of existing experimental investigations of flow boiling, it was found that liquid-laminar and gas-turbulent flow is a common feature in many applications of mini-channels. Traditional heat transfer correlations for saturated flow boiling were developed for liquid-turbulent and gas-turbulent flow conditions and thus may not be suitable in principle to be used to predict heat transfer coefficients in mini-channels when flow conditions are liquid-laminar and gas-turbulent. By considering flow conditions (laminar or turbulent) in the Reynolds number factor F and single-phase heat transfer coefficient h_sp, the Chen correlation has been modified to be used for four flow conditions such as liquid-laminar and gas-turbulent one often occurring in mini-channels. A comparison of the newly developed correlation with various existing data for mini-channels shows a satisfactory agreement. In addition, an extensive comparison of existing general correlations with databases for mini-channels has also been made.     

Modeling of nucleate boiling heat transfer under an impinging free jet

A model using an analytical/empirical approach has been developed to predict the rate of heat transfer in the stagnation region of a planar jet impinging on a horizontal flat surface. The model has been developed based on the hypothesis that bubble-induced mixing would result in enhanced or additional diffusivity. The additional diffusivity has been included in the diffusion term of the conservation equations. The value of the effective diffusivity has been correlated with jet parameters (velocity and temperature) and surface temperature using experimental data. The important aspects of the bubble dynamics (generation frequency and average bubble diameter) have been acquired using high-speed imaging and an intrusive optical probe. The applicability of the proposed model has been investigated under conditions of partial and fully-developed nucleate boiling. Experiments have been carried out using water at atmospheric pressure, mass flux in the range of 388–1649 kg/m² s, degree of sub-cooling in the range of 10–28 °C, and surface temperature in the range of 75–120 °C. Results showed that the proposed model is able to predict the surface heat flux with reasonable accuracy (+30% and ?15%).     

沸腾传热的分形学分析

分形论是当代兴起的非线性科学的一个重要组成部分。简要地阐述了分形论形的特征，解释了几何学中的维和物理学中的量纲的不同涵义。把分形论分数维的概念扩展应用到沸腾传热中的非线性问题的分析中，结合沸腾传热的操作特征和基本原理，以沸腾传热中的推动力△Tsat作为一个动力学维，进行了q-△Tsat^Df的标绘；对Df的物理意义作了分析和探讨。     

Models for gaseous radiative heat transfer applied to oxy-fuel conditions in boilers

Models of gas radiation properties have been evaluated for conditions relevant to oxy-fired boilers, characterized by larger pressure path-lengths and possibly different ratios of H₂O/CO₂ compared to air-fired boilers. Statistical narrow band (SNB) models serve as reference. The other radiation models tested are the weighted-sum-of-grey-gases model, the spectral line-based weighted-sum-of-grey-gases model and two grey-gas approximations. The range of validity of the existing coefficients of the weighted-sum-of-grey-gases model is limited, and new coefficients have therefore been determined to cover the conditions of interest. Several assumed test cases, involving both uniform and non-uniform paths, have been studied to evaluate the accuracy of the models. Comparisons with experimental data are also included. The results show that a grey approximation can give accurate wall fluxes, but at the expense of errors in the radiative source term. The weighted-sum-of-grey-gases model with the new coefficients yields predictions within 20% of those of the reference model in most cases, while the spectral line-based weighted-sum-of-grey-gases model usually gives results within 10%. There are, however, discrepancies between the SNB models at high temperatures. The weighted-sum-of-grey-gases model with its low computational cost is recommended for computationally demanding applications where predictions of both wall fluxes and the radiative source term are important.     

Surface wettability effects on critical heat flux of boiling heat transfer using nanoparticle coatings

This study investigates the effects of surface wettability on pool boiling heat transfer. Nano-silica particle coatings were used to vary the wettability of the copper surface from superhydrophilic to superhydrophobic by modifying surface topography and chemistry. Experimental results show that critical heat flux (CHF) values are higher in the hydrophilic region. Conversely, CHF values are lower in the hydrophobic region. The experimental CHF data of the modified surface do not fit the classical models. Therefore, this study proposes a simple model to build the nexus between the surface wettability and the growth of bubbles on the heating surface.     

Enhanced boiling heat transfer in mesochannels

Heat transfer characteristics are studied for a hybrid boiling case that combine features of spray cooling and flow boiling. In such a hybrid system, a liquid is atomized and the surrounding vapor is entrained into the droplet cone to provide an initial quality for enhanced boiling. An in-house experimental setup was developed to obtain surface temperature and heat flux measurements in a series of converged mesochannels for hybrid boiling. To compare the heat transfer performance of this hybrid technique, a flow boiling module was also developed using the same series of converged mesochannels. The inlet and exit hydraulic diameter of the mesochannels was 1.55 and 1.17 mm, respectively. The heat flux was in the range of 15–45 kW/m² and the estimated mass flux varied from 45 kg/m²s at the channel inlet to 110 kg/m²s at the channel outlet. Moreover, a model was presented to predict surface temperatures and heat transfer coefficients for flow boiling and hybrid boiling in mesochannels. This model was developed based on Chen’s formulation (1966) [21] but with two essential modifications. First, the laminar entry length effect was taken into consideration for heat transfer coefficient calculation. Second, the boiling enhancement factor was calculated based on the fluid properties. The model was compared to the experimental data and several other correlations for both cases. This model shows good agreement with the experimental data (mean deviations on the order of 12–16%).     

Marangoni凝结表面传热系数的影响因素分析

对水和酒精混合蒸气在竖直平板上的凝结传热进行了研究。利用液膜覆盖率,分别按考虑和忽略浓度边界层扩散热阻计算了凝结表面传热系数,并将计算结果与实验值进行了比较。结果发现同时考虑液膜导热热阻及浓度边界层扩散热阻的计算值在极小酒精质量分数时高于实验值,在较小酒精质量分数时和实验值接近,但是在高酒精质量分数时低于实验值。由此可见计算凝结表面传热系数时,在小酒精质量分数条件下扩散热阻可以忽略,但在高酒精质量分数条件下扩散热阻对整个热阻的贡献较大,必须考虑其对传热的影响。     

A wall heat transfer model for subcooled boiling flow

High compactness, low weight and little space requirement are gaining attention as prominent design criteria in the development of modern cooling systems in many applications. The resulting demand for highest possible heat transfer rates has lead to the very promising concept of providing for a controlled transition from pure single-phase convection to subcooled boiling flow in thermally highly loaded regions. For its application in modern engineering design this approach requires a realistic modeling of the complex phenomena associated with the two-phase flow heat transfer. The present work proposes for the computation of the specific wall heat transfer rate a modified superposition model, where the total heat flux is assumed to be additively composed of a forced convective and a nucleate boiling component. Since the present model requires only local input quantities, it is well suited to CFD of geometrically very complex coolant flows, where the definition of global length or velocity scales would be impractical. The wall heat fluxes predicted by the present model were compared against experimental data which were obtained by in-house measurements with water being the working fluid. The overall agreement is very good, particularly, in the partially nucleate boiling regime, where the effect of the bulk flow rate on the heat transfer is significant. Deviations are primarily observed at higher wall superheats, where a strong two-way coupling between the motion of the liquid and the motion of the bubbles as well as considerable bubble–bubble interactions typically occur.