Advanced Computational Modeling of Multiphase Boiling Flow and Heat Transfer

Author(s):  
Huiying Li ◽  
Sergio A. Vasquez ◽  
Peter Spicka

Numerical simulation of boiling flow and heat transfer presents a number of unique challenges in both theoretical modeling and developing robust numerical methodology. The major difficulty arises due to the heat transfer and phase changes between heated walls and fluid (liquid and vapor). Furthermore, modeling of the liquid-vapor interfacial transfers of momentum, heat and mass proves to be equally challenging. The multiphase boiling modeling approach described in this paper has been found to be capable of addressing these issues and is therefore suitable for inclusion in an advanced general purpose CFD solver. In the present approach, boiling flows are modeled within the framework of the Eulerian multifluid model. The governing equations solved are phase continuity, momentum and energy equations. Turbulence effects can be accounted for using mixture, dispersed or per-phase multiphase turbulence models. Wall boiling phenomena are modeled using the baseline mechanistic RPI model for nucleate boiling, and its extensions to non-equilibrium boiling and critical heat flux regime. A range of sub-models are considered to account for the interfacial momentum, mass and heat transfer, and flow regime transitions. An advanced numerical scheme has been developed for solving the model equations which can handle the heat partition between heated walls and fluid, provide for wall and interfacial mass transfer source terms in phase volume fraction equations, and address the coupling between the phase change rates and the pressure correction equation. The wall boiling models and numerical algorithm have been implemented in an advanced, general-purpose CFD code, FLUENT. Validations have been carried out for a range of 2D and 3D boiling flows, including pressurized water through a vertical pipe with heated walls, R-113 liquid in a vertical annulus with internal heated walls, a 3D BRW core channel geometry with vertical heated rods, and water in a vertical circular pipe under critical heat flux and post dry-out conditions. The results demonstrate that the wall boiling models are able to correctly predict the wall temperature and vapor volume fraction distribution. The predictions in all the cases are in reasonable good agreement with available experiments. Tests also indicate that the present implementation is fast and robust, as compared to previous approaches. All the cases are able to be simulated with the use of the FLUENT steady-state multiphase solver with reasonable numbers of iterations.

Author(s):  
Muhamad Zuhairi Sulaiman ◽  
Masahiro Takamura ◽  
Kazuki Nakahashi ◽  
Tomio Okawa

Boiling heat transfer (BHT) and critical heat flux (CHF) performance were experimentally studied for saturated pool boiling of water-based nanofluids. In present experimental works, copper heaters of 20 mm diameter with titanium-oxide (TiO2) nanocoated surface were produced in pool boiling of nanofluid. Experiments were performed in both upward and downward facing nanofluid coated heater surface. TiO2 nanoparticle was used with concentration ranging from 0.004 until 0.4 kg/m3 and boiling time of tb = 1, 3, 10, 20, 40, and 60 mins. Distilled water was used to observed BHT and CHF performance of different nanofluids boiling time and concentration configurations. Nucleate boiling heat transfer observed to deteriorate in upward facing heater, however; in contrast effect of enhancement for downward. Maximum enhancements of CHF for upward- and downward-facing heater are 2.1 and 1.9 times, respectively. Reduction of mean contact angle demonstrate enhancement on the critical heat flux for both upward-facing and downward-facing heater configuration. However, nucleate boiling heat transfer shows inconsistency in similar concentration with sequence of boiling time. For both downward- and upward-facing nanocoated heater's BHT and CHF, the optimum configuration denotes by C = 400 kg/m3 with tb = 1 min which shows the best increment of boiling curve trend with lowest wall superheat ΔT = 25 K and critical heat flux enhancement of 2.02 times.


1999 ◽  
Author(s):  
Yasuo Koizumi ◽  
Hiroyasu Ohtake ◽  
Manabu Mochizuki

Abstract The effect of solid particle introduction on subcooled-forced flow boiling heat transfer and a critical heat flux was examined experimentally. In the experiment, glass beads of 0.6 mm diameter were mixed in subcooled water. Experiments were conducted in a range of the subcooling of 40 K, a velocity of 0.17–6.7 m/s, a volumetric particle ratio of 0–17%. When particles were introduced, the growth of a superheated liquid layer near a heat trasnsfer surface seemed to be suppressed and the onset of nucleate boiling was delayed. The particles promoted the condensation of bubbles on the heat transfer surface, which shifted the initiation of a net vapor generation to a high heat flux region. Boiling heat trasnfer was augmented by the particle introduction. The suppression of the growth of the superheated liquid layer and the promotion of bubble condensation and dissipation by the particles seemed to contribute that heat transfer augmentation. The wall superheat at the critical heat flux was elevated by the particle introduction and the critical heat flux itself was also enhanced. However, the degree of the critical heat flux improvement was not drastic.


2019 ◽  
Vol 141 (4) ◽  
Author(s):  
S. R. Darr ◽  
J. W. Hartwig ◽  
J. Dong ◽  
H. Wang ◽  
A. K. Majumdar ◽  
...  

Recently, two-phase cryogenic flow boiling data in liquid nitrogen (LN2) and liquid hydrogen (LH2) were compared to the most popular two-phase correlations, as well as correlations used in two of the most widely used commercially available thermal/fluid design codes in Hartwig et al. (2016, “Assessment of Existing Two Phase Heat Transfer Coefficient and Critical Heat Flux on Cryogenic Flow Boiling Quenching Experiments,” Int. J. Heat Mass Transfer, 93, pp. 441–463). Results uncovered that the correlations performed poorly, with predictions significantly higher than the data. Disparity is primarily due to the fact that most two-phase correlations are based on room temperature fluids, and for the heating configuration, not the quenching configuration. The penalty for such poor predictive tools is higher margin, safety factor, and cost. Before control algorithms for cryogenic transfer systems can be implemented, it is first required to develop a set of low-error, fundamental two-phase heat transfer correlations that match available cryogenic data. This paper presents the background for developing a new set of quenching/chilldown correlations for cryogenic pipe flow on thin, shorter lines, including the results of an exhaustive literature review of 61 sources. New correlations are presented which are based on the consolidated database of 79,915 quenching points for a 1.27 cm diameter line, covering a wide range of inlet subcooling, mass flux, pressure, equilibrium quality, flow direction, and even gravity level. Functional forms are presented for LN2 and LH2 chilldown correlations, including film, transition, and nucleate boiling, critical heat flux, and the Leidenfrost point.


Author(s):  
Matjazˇ Leskovar ◽  
Bosˇtjan Koncˇar

An ex-vessel steam explosion may occur when during a severe reactor accident the reactor vessel fails and the molten core pours into the water in the reactor cavity. A steam explosion is a fuel coolant interaction process where the heat transfer from the melt to water is so intense and rapid that the timescale for heat transfer is shorter than the timescale for pressure relief. This can lead to the formation of shock waves and production of missiles at later times, during the expansion of the highly pressurized water vapor, that may endanger surrounding structures. In contrast to specialized steam explosion CFD codes, where the steam explosion is modeled on micro-scale using fundamental averaged multiphase flow conservation equations, in the presented approach the steam explosion is modeled in a simplified manner as an expanding high-pressure premixture of dispersed molten fuel, liquid water and vapor. Applying the developed steam explosion model, a comprehensive analysis of the ex-vessel steam explosion in a typical PWR reactor cavity was done using the CFD code CFX-10. At four selected locations, which are of importance for the assessment of the vulnerability of cavity structures, the pressure histories were recorded and the corresponding pressure impulses calculated. The pressure impulses determine the destructive potential of the steam explosion and represent the input for the structural mechanical analysis of the cavity structures. The simulation results show that the pressure impulses depend mainly on the steam explosion energy conversion ratio, whereas the influence of the premixture vapor volume fraction, which is a parameter in our model and determines the maximum steam explosion pressure, is not significant.


Micromachines ◽  
2021 ◽  
Vol 12 (8) ◽  
pp. 980
Author(s):  
Kairui Tang ◽  
Jingjing Bai ◽  
Siyu Chen ◽  
Shiwei Zhang ◽  
Jie Li ◽  
...  

With the rapid development of electronics, thermal management has become one of the most crucial issues. Intense research has focused on surface modifications used to enhance heat transfer. In this study, multilayer copper micromeshes (MCMs) are developed for commercial compact electronic cooling. Boiling heat transfer performance, including critical heat flux (CHF), heat transfer coefficients (HTCs), and the onset of nucleate boiling (ONB), are investigated. The effect of micromesh layers on the boiling performance is studied, and the bubbling characteristics are analyzed. In the study, MCM-5 shows the highest critical heat flux (CHF) of 207.5 W/cm2 and an HTC of 16.5 W(cm2·K) because of its abundant micropores serving as nucleate sites, and outstanding capillary wicking capability. In addition, MCMs are compared with other surface structures in the literature and perform with high competitiveness and potential in commercial applications for high-power cooling.


2001 ◽  
Vol 1 (1) ◽  
pp. 32
Author(s):  
P. M. Carrica ◽  
V. Masson

We present the results of an experimental study of the effects of externally imposed electric fields on boiling heat transfer and critical heat flux (CHF) in dielectric fluids. The study comprises the analysis of geometries that, under the effects of electric fields, cause the bubbles either to be pushed toward the heater or away from it. A local phase detection probe was used to measure the void fraction and the interfacial impact rate near the heater. It was found that the critical heat flux can be either augmented or reduced with the application of an electric field, depending on the direction of . In addition, the heat transfer can be slightly enhanced or degraded depending on the heat flux. The study of the two-phase flow in nucleate boiling, only for the case of favorable dielectrophoretic forces, reveals that the application of an electric field reduces the bubble detection time and increases the detachment frequency. It also shows that the two-phase flow characteristics of the second film boiling regime resemble more a nucleate boiling regime than a film boiling regime.


2018 ◽  
Vol 22 (Suppl. 2) ◽  
pp. 617-627
Author(s):  
Jie Chen ◽  
Weihua Cai ◽  
Shulei Li ◽  
Yan Ren ◽  
Hongqiang Ma ◽  
...  

Plate-fin heat exchanger with rectangular minichannels, as a type of high-perfor- mance compact heat exchangers, has been widely used in liquefied natural gas field. However, the studies on saturated boiling flow and heat transfer for mixture refrigerant in plate-fin heat exchanger have been scarcely explored, which are helpful for designing more effective plate-fin heat exchanger using in liquefied natural gas field. Therefore, in this paper, the characteristics of saturated boiling flow and heat transfer for mixture refrigerant in rectangular minichannels of plate-fin heat exchanger were studied numerally based on validated model. Then, the effect of different parameters (vapor quality, mass flux, and heat flux) on heat transfer coefficient and frictional pressure drop were discussed. The results indicated that the boiling heat transfer coefficient and pressure drop are mainly influenced by quality and mass flux while heat flux has little influence on them. This is due to the fact that the main boiling mechanisms were forced convective boiling and the evaporation of dispersed liquid phase while nucleate boiling is slight.


1993 ◽  
Vol 115 (1) ◽  
pp. 78-88 ◽  
Author(s):  
C. O. Gersey ◽  
I. Mudawar

The effects of chip protrusion on the forced-convection boiling and critical heat flux (CHF) of a dielectric coolant (FC-72) were investigated. The multi-chip module used in the present study featured a linear array of nine, 10 mm x 10 mm, simulated microelectronic chips which protruded 1 mm into a 20-mm wide side of a rectangular flow channel. Experiments were performed in vertical up flow with 5-mm and 2-mm channel gap thicknesses. For each configuration, the velocity and subcooling of the liquid were varied from 13 to 400 cm/s and 3 to 36° C, respectively. The nucleate boiling regime was not affected by changes in velocity and subcooling, and critical heat flux generally increased with increases in either velocity or subcooling. Higher single-phase heat transfer coefficients and higher CHF values were measured for the protruded chips compared to similar flush-mounted chips. However, adjusting the data for the increased surface area and the increased liquid velocity above the chip caused by the protruding chips yielded a closer agreement between the protruded and flush-mounted results. Even with the velocity and area adjustments, the most upstream protruded chip had higher single-phase heat transfer coefficients and CHF values for high velocity and/or highly-subcooled flow as compared the downstream protruded chips. The results show that, except for the most upstream chip, the performances of protruded chips are very similar to those of flush-mounted chips.


Author(s):  
Wei Tong ◽  
Alireza Ganjali ◽  
Omidreza Ghaffari ◽  
Chady Alsayed ◽  
Luc Frechette ◽  
...  

Abstract In a two-phase immersion cooling system, boiling on the spreader surface has been experimentally found to be non-uniform, and it is highly related to the surface temperature and the heat transfer coefficient. An experimentally obtained temperature-dependent boiling heat transfer coefficient has been applied to a numerical model to investigate the spreader's cooling performance. It is found that the surface temperature distribution becomes less uniform with higher input power. But it is more uniform when the thickness is increased. By defining the characteristic temperatures that represent different boiling regimes on the surface, the fraction of the surface area that has reached the critical heat flux has been numerically calculated, showing that increasing the thickness from 1 mm to 6 mm decreases the critical heat flux reached area by 23% at saturation liquid temperatures. Therefore, on the thicker spreader, more of the surface is utilized for nucleate boiling while localized hot regions that lead to surface dry-out are avoided. At a base temperature of 90 oC, the optimal thickness is found to be 4 mm, beyond which no significant improvement in heat removal can be obtained. Lower coolant temperatures can further increase the heat removal; it is reduced from an 18% improvement in the input power for the 1 mm case to only 3% in the 6 mm case for a coolant temperature drop of 24 oC. Therefore, a trade-off exists between the cost of maintaining the low liquid temperature and the increased heat removal capacity.


Sign in / Sign up

Export Citation Format

Share Document