Three-Dimensional Numerical Simulation of Burnout on Horizontal Surface in Pool Boiling

2003 ◽  
Author(s):  
Zoran V. Stosic ◽  
Vladimir D. Stevanovic
Keyword(s):  
Heat Flux ◽  
Residence Time ◽  
Pool Boiling ◽  
Three Dimensional ◽  
Heating Surface ◽  
Heat Fluxes ◽  
Heating Wall ◽  
Two Phase ◽  
Nucleation Sites ◽  
Phase Mixture

Prediction of nucleate boiling mechanisms and burnout conditions, when heat transfer coefficient sharply drops and the heating surface destruction could occur are one of the crucial topics in thermal design and safety analyses of various thermal equipment. Although these phenomena have been intensively investigated for decades, various influencing factors and complexity of coupled thermal and fluid dynamic processes have not yet been fully understood. The integral approach towards prediction of nucleate boiling and burnout conditions requires modelling and numerical simulation of micro level phenomena of bubble rise and departure at a numbers of nucleation sites, as well as macroscopic two-phase mixture behaviour on the heating surface. In this paper multidimensional numerical simulation of the atmospheric saturated pool boiling is performed under high heat fluxes, near to and at the occurrence of burnout conditions. Micro level phenomena on the heating surface are modelled with the key parameters of vapour generation on the heating wall, such as bubble nucleation site density, bubble residence time on the heating wall and certain level of randomness in the location of bubble nucleation. Heat flux is non-uniformly distributed on the heating surface with peaks at the nucleation sites. The nucleation sites are determined by a random function, while the mean number of nucleation sites is prescribed according to the material characteristics and roughness of the heating surface. The applied numerical grids are able to represent the nucleation sites on the heating wall for both fresh (polished) and aged (rough) heaters at the atmospheric pool boiling conditions. The macro level phenomena are modelled with two-fluid model of liquid-vapour flow. The interfacial drag is modelled with appropriate closure laws. The applied modelling and numerical methods enable full representation of the two-phase mixture behaviour on the heating surface with inclusion of the swell level prediction. In this way the integral conditions of nucleate pool boiling with the possibility of burnout are simulated and the critical heat flux conditions are predicted. The result of the three-dimensional numerical simulations and analyses are presented as the extension of the previously published two-dimensional numerical results. Here presented three-dimensional investigation is performed in order to take into account more realistically spatial effects of vapour generation and two-phase flow, such as phase dispersion within the two-phase mixture, than it was able with previously performed two-dimensional investigation. Results are presented for short time period after the initiation of heat supply and vapour generation on the heating surface, as well as for quasi steady-state conditions after several seconds from pool boiling initiation. A replenishment of the heating surface with water and partial surface wetting for lower heat fluxes is shown, while heating surface dry-out is observed for high heat fluxes. The influence of the density of nucleation sites and the bubble residence time on the wall on the pool boiling dynamics is investigated. Also, the influence of the heat flux intensity on the pool boiling dynamics is analysed. Numerical simulations show that decrease of the density of nucleation sites and increase of bubble residence time on the heating surface (characteristics pertinent to fresh-polished heaters) lead to the reduction of critical heat flux values. Obtained results are in excellent agreement with the recent experimental investigations of the upward facing burnout conditions on the horizontal heated plate. Details of the developed numerical procedure are presented. The introduced method of random spatial and temporal generation of the vapour at the heated wall is a new approach. It enables the macroscopic representation of the population of microscopic vapour bubbles, which are generated at nucleation sites on the heater wall, and which burst through liquid micro-layer in thermal-hydraulic conditions close to the burnout. The applied numerical and modelling method has shown robustness by allowing stable calculations for wide ranges of applied modelling boiling parameters (density of nucleation sites and bubble residence time).

Multi-Dimensional Numerical Simulation of Burnout on Horizontal Surface in Pool Boiling

2002 ◽  
Author(s):  
Zoran V. Stosic ◽  
Vladimir D. Stevanovic
Keyword(s):  
Numerical Simulation ◽  
Heat Flux ◽  
Residence Time ◽  
Pool Boiling ◽  
Heating Surface ◽  
High Heat ◽  
Heat Fluxes ◽  
Heating Wall ◽  
Nucleation Sites ◽  
Vapour Generation

Multidimensional numerical simulation of the atmospheric saturated pool boiling is performed under high heat fluxes, near to and at the occurrence of burnout conditions. Heat flux through the vessel bottom wall is varied and its influence on the pool boiling dynamics is analysed. Dynamics of vapour generation on the heating wall is modelled through the density of nucleation sites and the bubble residence time on the wall. The nucleation sites are determined by a random function. The applied numerical grid is able to represent the nucleation sites on the heating wall for both fresh (polished) and aged (rough) heaters at the atmospheric pool boiling conditions. Results are presented for short time period after the initiation of heat supply and vapour generation on the heating surface, as well as for quasi steady-state conditions after two seconds from pool boiling initiation. The results show a replenishment of the heating surface with water and partial surface wetting for lower heat fluxes, while heating surface dry-out is predicted for high heat fluxes. The influence of the density of nucleation sites and the bubble residence time on the wall on the pool boiling dynamics is investigated. Numerical simulations show that decrease of the density of nucleation sites and increase of bubble residence time on the heating surface (characteristics pertinent to fresh-polished heaters) lead to the reduction of critical heat flux values. Obtained results are in excellent agreement with the recent experimental investigations of the upward facing burnout conditions on the horizontal heated plate. Details of the developed numerical procedure are presented. The introduced method of random spatial and temporal generation of the vapour at the heated wall is a new approach. It enables the macroscopic representation of the population of microscopic vapour bubbles, which are generated at nucleation sites on the heater wall, and which burst through liquid micro-layer in thermal-hydraulic conditions close to the burnout. The applied numerical and modelling method has shown robustness by allowing stable calculations for wide ranges of applied modelling boiling parameters (density of nucleation sites and bubble residence time).

A Review of Two-Phase Forced Cooling in Three-Dimensional Stacked Electronics: Technology Integration

2015 ◽  
Vol 137 (4) ◽  
Author(s):  
Craig Green ◽  
Peter Kottke ◽  
Xuefei Han ◽  
Casey Woodrum ◽  
Thomas Sarvey ◽  
...  
Keyword(s):  
Heat Flux ◽  
Heat Sink ◽  
Three Dimensional ◽  
High Heat ◽  
Heat Fluxes ◽  
Thermal Design ◽  
High Heat Flux ◽  
Two Phase ◽  
Forced Cooling ◽  
Fluidic Routing

Three-dimensional (3D) stacked electronics present significant advantages from an electrical design perspective, ranging from shorter interconnect lengths to enabling heterogeneous integration. However, multitier stacking exacerbates an already difficult thermal problem. Localized hotspots within individual tiers can provide an additional challenge when the high heat flux region is buried within the stack. Numerous investigations have been launched in the previous decade seeking to develop cooling solutions that can be integrated within the 3D stack, allowing the cooling to scale with the number of tiers in the system. Two-phase cooling is of particular interest, because the associated reduced flow rates may allow reduction in pumping power, and the saturated temperature condition of the coolant may offer enhanced device temperature uniformity. This paper presents a review of the advances in two-phase forced cooling in the past decade, with a focus on the challenges of integrating the technology in high heat flux 3D systems. A holistic approach is applied, considering not only the thermal performance of standalone cooling strategies but also coolant selection, fluidic routing, packaging, and system reliability. Finally, a cohesive approach to thermal design of an evaporative cooling based heat sink developed by the authors is presented, taking into account all of the integration considerations discussed previously. The thermal design seeks to achieve the dissipation of very large (in excess of 500 W/cm2) background heat fluxes over a large 1 cm × 1 cm chip area, as well as extreme (in excess of 2 kW/cm2) hotspot heat fluxes over small 200 μm × 200 μm areas, employing a hybrid design strategy that combines a micropin–fin heat sink for background cooling as well as localized, ultrathin microgaps for hotspot cooling.

scholarly journals Numerical prediction of nucleate pool boiling heat transfer coefficient under high heat fluxes

2016 ◽  
Vol 20 (suppl. 1) ◽  
pp. 113-123 ◽  
Author(s):  
Milada Pezo ◽  
Vladimir Stevanovic
Keyword(s):  
Heat Transfer ◽  
Heat Transfer Coefficient ◽  
Transfer Coefficient ◽  
Pool Boiling ◽  
Heating Surface ◽  
High Heat ◽  
Bubble Nucleation ◽  
Heat Fluxes ◽  
Two Phase ◽  
Boiling Heat Transfer Coefficient

This paper presents CFD (Computational Fluid Dynamics) approach to prediction of the heat transfer coefficient for nucleate pool boiling under high heat fluxes. Three-dimensional numerical simulations of the atmospheric saturated pool boiling are performed. Mathematical modelling of pool boiling requires a treatment of vapor-liquid two-phase mixture on the macro level, as well as on the micro level, such as bubble growth and departure from the heating surface. Two-phase flow is modelled by the two-fluid model, which consists of the mass, momentum and energy conservation equations for each phase. Interface transfer processes are calculated by the closure laws. Micro level phenomena on the heating surface are modelled with the bubble nucleation site density, the bubble resistance time on the heating wall and with the certain level of randomness in the location of bubble nucleation sites. The developed model was used to determine the heat transfer coefficient and results of numerical simulations are compared with available experimental results and several empirical correlations. A considerable scattering of the predictions of the pool boiling heat transfer coefficient by experimental correlations is observed, while the numerically predicted values are within the range of results calculated by well-known Kutateladze, Mostinski, Kruzhilin and Rohsenow correlations. The presented numerical modeling approach is original regarding both the application of the two-fluid two-phase model for the determination of heat transfer coefficient in pool boiling and the defined boundary conditions at the heated wall surface.

scholarly journals Heat transfer and critical heat flux for two-phase immersion cooling of an inverter power module

2021 ◽  
Vol 2116 (1) ◽  
pp. 012007
Author(s):  
I T’ Jollyn ◽  
J Nonneman ◽  
M De Paepe
Keyword(s):  
Heat Transfer ◽  
Heat Flux ◽  
Critical Heat Flux ◽  
Pool Boiling ◽  
Flux Measurement ◽  
Base Plate ◽  
Heat Fluxes ◽  
Power Module ◽  
Two Phase ◽  
Heat Flux Measurement

Abstract Heat transfer and critical heat flux measurement are reported for pool boiling cooling of the base plate of an inverter power module. Novec 649 is used as refrigerant. Heat fluxes up to 14.6 W/cm2 were applied with refrigerant saturation temperatures of 36 °C, 41 °C and 46 °C. The measured boiling curves are comparable to those reported for similar refrigerants. The critical heat fluxes range from 12.1 W/cm2 to 14.6 W/cm2, which corresponds within 10% to the correlation of Zuber. The critical heat flux is significantly lower than the highest heat fluxes expected from the power module, indicating that methods to increase the critical heat flux are needed to enable two-phase power module cooling.

scholarly journals Numerical investigation of nucleate pool boiling heat transfer

2016 ◽  
Vol 20 (suppl. 5) ◽  
pp. 1301-1312
Author(s):  
Andrijana Stojanovic ◽  
Vladimir Stevanovic ◽  
Milan Petrovic ◽  
Dragoljub Zivkovic
Keyword(s):  
Heat Flux ◽  
Pool Boiling ◽  
Nucleation Site ◽  
Bubble Nucleation ◽  
Heated Wall ◽  
Wall Surface ◽  
Two Phase ◽  
Site Density ◽  
Nucleation Sites ◽  
Nucleation Site Density

Multidimensional numerical simulation of the atmospheric saturated pool boiling is performed. The applied modelling and numerical methods enable a full representation of the liquid and vapour two-phase mixture behaviour on the heated surface, with included prediction of the swell level and heated wall temperature field. In this way the integral behaviour of nucleate pool boiling is simulated. The micro conditions of bubble generation at the heated wall surface are modelled by the bubble nucleation site density, the liquid wetting contact angle and the bubble grow time. The bubble nucleation sites are randomly located within zones of equal size, where the number of zones equals the nucleation site density. The conjugate heat transfer from the heated wall to the liquid is taken into account in wetted heated wall areas around bubble nucleation sites. The boiling curve relation between the heat flux and the heated wall surface temperature in excess of the saturation temperature is predicted for the pool boiling conditions reported in the literature and a good agreement is achieved with experimentally measured data. The influence of the nucleation site density on the boiling curve characteristic is confirmed. In addition, the influence of the heat flux intensity on the spatial effects of vapour generation and two-phase flow are shown, such as the increase of the swell level position and the reduced wetting of the heated wall surface by the heat flux increase.

On Correlating the Peak and Minimum Boiling Heat Fluxes With Pressure and Heater Configuration

1966 ◽  
Vol 88 (1) ◽  
pp. 94-99 ◽  
Author(s):  
John H. Lienhard ◽  
Kiyokazu Watanabe
Keyword(s):  
Heat Flux ◽  
Pressure Range ◽  
Pool Boiling ◽  
Three Dimensional ◽  
Nucleate Boiling ◽  
Heat Fluxes ◽  
Reduced Pressure ◽  
Minimum Heat ◽  
Dimensional Surface ◽  
Horizontal Wire

The peak nucleate boiling heat flux has been observed for five fluids during saturated pool boiling on horizontal wire heaters, ranging in radius from 0.0025 in. to 0.0254 in., over a reduced pressure range from 0.0010 to 0.0197. A scheme for correlating the peak and minimum heat fluxes is developed heuristically and successfully applied to these data. The result is a single three-dimensional surface which represents all of the data. The surface can be represented as the product of a function of geometric scale, and a function of pressure. The function of pressure appears to be the same in any configuration.

THE MECHANISM OF RAISING AND QUANTIFICATION OF SPECIFIC HEAT FLUX AT BOILING OF NANOFLUIDS IN FREE CONVECTION CONDITIONS

2017 ◽  
pp. 25-34
Author(s):  
V.N. Moraru
Keyword(s):  
Heat Transfer ◽  
Heat Flux ◽  
Heat Capacity ◽  
Specific Heat ◽  
Chemical Properties ◽  
Heating Surface ◽  
Physical Models ◽  
Superheated Liquid ◽  
Two Phase ◽  
Boiling Process

The results of our work and a number of foreign studies indicate that the sharp increase in the heat transfer parameters (specific heat flux q and heat transfer coefficient _) at the boiling of nanofluids as compared to the base liquid (water) is due not only and not so much to the increase of the thermal conductivity of the nanofluids, but an intensification of the boiling process caused by a change in the state of the heating surface, its topological and chemical properties (porosity, roughness, wettability). The latter leads to a change in the internal characteristics of the boiling process and the average temperature of the superheated liquid layer. This circumstance makes it possible, on the basis of physical models of the liquids boiling and taking into account the parameters of the surface state (temperature, pressure) and properties of the coolant (the density and heat capacity of the liquid, the specific heat of vaporization and the heat capacity of the vapor), and also the internal characteristics of the boiling of liquids, to calculate the value of specific heat flux q. In this paper, the difference in the mechanisms of heat transfer during the boiling of single-phase (water) and two-phase nanofluids has been studied and a quantitative estimate of the q values for the boiling of the nanofluid is carried out based on the internal characteristics of the boiling process. The satisfactory agreement of the calculated values with the experimental data is a confirmation that the key factor in the growth of the heat transfer intensity at the boiling of nanofluids is indeed a change in the nature and microrelief of the heating surface. Bibl. 20, Fig. 9, Tab. 2.

The Mechanism of and Stability Criterion for Nucleate Pool Boiling of Sodium

1969 ◽  
Vol 91 (3) ◽  
pp. 315-328 ◽  
Author(s):  
I. Shai ◽  
W. M. Rohsenow
Keyword(s):  
Heat Flux ◽  
Liquid Metals ◽  
Bubble Formation ◽  
Pool Boiling ◽  
Heat Fluxes ◽  
Nucleate Pool Boiling ◽  
Thermal Layer ◽  
Minimum Heat ◽  
Bubble Growth And Departure ◽  
Recorded Temperature

Experimental data for sodium boiling on horizontal surfaces containing artificial cavities at heat fluxes of 20,000 to 300,000 Btu/ft2 hr and pressures between 40 to 106 mm Hg were obtained. Observations are made for stable boiling, unstable boiling and “bumping.” Some recorded temperature variations in the solid close to the nucleating cavity are presented. It is suggested that for liquid metals the time for bubble growth and departure is a very small fraction of the total bubble cycle, hence the delay time during which a thermal layer grows is the most significant part of the process. On this basis the transient conduction heat transfer is solved for a periodic process, and the period time is found to be a function of the degree of superheat, the heat flux and the liquid thermal properties. A simplified model for stability of nucleate pool boiling of liquid metals is postulated from which the minimum heat flux for stable boiling can be found as a function of liquid-solid properties, liquid pressure, the degree of superheat, and the cavity radius and depth. At relatively low heat fluxes, convection currents have significant effects on the period time of bubble formation. An empirical correlation is proposed, which takes into account the convection effects, to match the experimental results.

Forced Convective Boiling of Refrigerant HCFC123 in a Mini-Tube

2010 ◽  
Author(s):  
Koichi Araga ◽  
Keisuke Okamoto ◽  
Keiji Murata
Keyword(s):  
Heat Transfer ◽  
Heat Flux ◽  
Mass Flux ◽  
Pool Boiling ◽  
Constant Heat Flux ◽  
Transfer Coefficients ◽  
Two Phase ◽  
Convective Boiling ◽  
Frictional Pressure Drop ◽  
Heat Transfer Coefficients

This paper presents an experimental investigation of the forced convective boiling of refrigerant HCFC123 in a mini-tube. The inner diameters of the test tubes, D, were 0.51 mm and 0.30 mm. First, two-phase frictional pressure drops were measured under adiabatic conditions and compared with the correlations for conventional tubes. The frictional pressure drop data were lower than the correlation for conventional tubes. However, the data were qualitatively in accord with those for conventional tubes and were correlated in the form φL2−1/Xtt. Next, heat transfer coefficients were measured under the conditions of constant heat flux and compared with those for conventional tubes and for pool boiling. The heat transfer characteristics for mini-tubes were different from those for conventional tubes and quite complicated. The heat transfer coefficients for D = 0.51 mm increased with heat flux but were almost independent of mass flux. Although the heat transfer coefficients were higher than those for a conventional tube with D = 10.3 mm and for pool boiling in the low quality region, they decreased gradually with increasing quality. The heat transfer coefficients for D = 0.30 mm were higher than those for D = 0.51 mm and were almost independent of both mass flux and heat flux.

scholarly journals Determination of Pool Boiling Critical Heat Flux Enhancement in Nanofluids

2007 ◽  
Author(s):  
Bao H. Truong
Keyword(s):  
Heat Flux ◽  
Critical Heat Flux ◽  
Pool Boiling ◽  
Nucleate Boiling ◽  
Porous Coating ◽  
Transfer Coefficients ◽  
Two Phase ◽  
Sem Images ◽  
Heat Transfer Coefficients ◽  
Nucleate Boiling Heat Transfer

Nanofluids are engineered colloids composed of nano-size particles dispersed in common fluids such as water or refrigerants. Using an electrically controlled wire heater, pool boiling Critical Heat Flux (CHF) of Alumina and Silica water-based nanofluids of concentration less than or equal to 0.1 percent by volume were measured. Silica nanofluids showed a CHF enhancement up to 68% and there seems to be a monotonic relationship between the nanoparticle concentration and the magnitude of enhancement. Alumina nanofluids had a CHF enhancement up to 56% but the peak occurred at the intermediate concentration. The boiling curves in nanofluid were found to shift to the left of that of water and correspond to higher nucleate boiling heat transfer coefficients in the two-phase flow regime. Scanning Electron Microscopy (SEM) images show a porous coating layer of nanoparticles on wires subjected to nanofluid CHF tests. These coating layers change the morphology of the heater’s surface, and are responsible for the CHF enhancement. The thickness of the coating was estimated using SEM and was found ranging from 3.0 to 6.0 micrometers for Alumina, and 3.0 to 15.0 micrometers for Silica.

Sign in / Sign up

Export Citation Format

Share Document