Numerical Modeling of the Conjugate Heat Transfer Problem for Annular Laminar Film Condensation in Microchannels

2012 ◽  
Vol 134 (5) ◽  
Author(s):  
Stefano Nebuloni ◽  
John R. Thome

This paper presents numerical simulations of annular laminar film condensation heat transfer in microchannels of different internal shapes. The model, which is based on a finite volume formulation of the Navier–Stokes and energy equations for the liquid phase only, importantly accounts for the effects of axial and peripheral wall conduction and nonuniform heat flux not included in other models so far in the literature. The contributions of the surface tension, axial shear stresses, and gravitational forces are included. This model has so far been validated versus various benchmark cases and versus experimental data available in literature, predicting microchannel heat transfer data with an average error of 20% or better. It is well known that the thinning of the condensate film induced by surface tension due to gravity forces and shape of the surface, also known as the “Gregorig” effect, has a strong consequence on the local heat transfer coefficient in condensation. Thus, the present model accounts for these effects on the heat transfer and pressure drop for a wide variety of geometrical shapes, sizes, wall materials, and working fluid properties. In this paper, the conjugate heat transfer problem arising from the coupling between the thin film fluid dynamics, the heat transfer in the condensing fluid, and the heat conduction in the channel wall has been studied. In particular, the work has focused on three external channel wall boundary conditions: a uniform wall temperature, a nonuniform wall heat flux, and single-phase convective cooling are presented. As the scale of the problem is reduced, i.e., when moving from mini- to microchannels, the results show that the axial conduction effects can become very important in the prediction of the wall temperature profile and the magnitude of the heat transfer coefficient and its distribution along the channel.

Author(s):  
Stefano Nebuloni ◽  
John R. Thome

This paper presents numerical simulations of annular laminar film condensation heat transfer in micro-channels of different internal shapes. The model, which is based on a finite volume formulation of the Navier-Stokes and energy equations for the liquid phase only, importantly accounts for the effects of axial and peripheral wall conduction and non-uniform heat flux not included in other models so far in the literature. The contributions of the surface tension, axial shear stresses and gravitational forces are included. This model has so far been validated versus various benchmark cases and versus experimental data available in literature, predicting microchannel heat transfer data with an average error of 20% or better. It is well-known that the thinning of the condensate film induced by surface tension due to gravity forces and shape of the surface, also known as the ‘Grigorig’ effect, has a strong consequence on the local heat transfer coefficient in condensation. Thus, the present model accounts for these effects on the heat transfer and pressure drop for a wide variety of geometrical shapes, sizes, wall materials and working fluid properties. In this paper, the conjugate heat transfer problem arising from the coupling between the thin film fluid dynamics, the heat transfer in the condensing fluid and the heat conduction in the channel wall has been studied. In particular, the work has focused on three external channel wall boundary conditions: a uniform wall temperature, a non uniform wall heat flux and single-phase convective cooling is presented. As the scale of the problem is reduced, i.e. when moving from mini to micro channels, the results shows that the axial conduction effects can become very important in the prediction of the wall temperature profile and the magnitude of the heat transfer coefficient and its distribution along the channel.


Author(s):  
Lei Zhao ◽  
Ting Wang

In film cooling heat transfer analysis, one of the core concepts is to deem film cooled adiabatic wall temperature (Taw) as the driving potential for the actual heat flux over the film-cooled surface. Theoretically, the concept of treating Taw as the driving temperature potential is drawn from compressible flow theory when viscous dissipation becomes the heat source near the wall and creates higher wall temperature than in the flowing gas. But in conditions where viscous dissipation is negligible, which is common in experiments under laboratory conditions, the heat source is not from near the wall but from the main hot gas stream; therefore, the concept of treating the adiabatic wall temperature as the driving potential is subjected to examination. To help investigate the role that Taw plays, a series of computational simulations are conducted under typical film cooling conditions over a conjugate wall with internal flow cooling. The result and analysis support the validity of this concept to be used in the film cooling by showing that Taw is indeed the driving temperature potential on the hypothetical zero wall thickness condition, ie. Taw is always higher than Tw with underneath (or internal) cooling and the adiabatic film heat transfer coefficient (haf) is always positive. However, in the conjugate wall cases, Taw is not always higher than wall temperature (Tw), and therefore, Taw does not always play the role as the driving potential. Reversed heat transfer through the airfoil wall from downstream to upstream is possible, and this reversed heat flow will make Tw > Taw in the near injection hole region. Yet evidence supports that Taw can be used to correctly predict the heat flux direction and always result in a positive adiabatic heat transfer coefficient (haf). The results further suggest that two different test walls are recommended for conducting film cooling experiments: a low thermal conductivity material should be used for obtaining accurate Taw and a relative high thermal conductivity material be used for conjugate cooling experiment. Insulating a high-conductivity wall will result in Taw distribution that will not provide correct heat flux or haf values near the injection hole.


Author(s):  
Zhan Yin ◽  
Jianjun Wen ◽  
Min Zeng ◽  
Qiuwang Wang

A steady three-dimensional numerical simulation of laminar film condensation of vapor in the presence of air inside a 1 mm horizontal tube is presented. The volume of fluid (VOF) method is used to capture the liquid-vapor interface with a phase change model. According to a generally accepted flow regime map, annular flow pattern is to be expected. Uniform wall temperature and interface temperature are assumed to be boundary condition. The influence of gravity is obvious while the effect of surface tension is neglected. At inlet, the liquid film is thin and evenly distributed around tube wall. Moving downstream the tube, film at the bottom half becomes thicker under the influence of gravity, while film on upper half remains almost constant. Correspondingly, local heat transfer coefficient on bottom half declines gradually and global average heat transfer coefficient shows little difference along axial direction. Existence of air makes heat transfer coefficient decrease sharply compared with that of pure vapor condensation, caused by an existed air layer which increases the thermal resistance during condensation process. As inlet volume fraction of air increases from 0.5% to 3%, the decline trend of heat transfer coefficient slows down.


Author(s):  
X. C. Li ◽  
J. Zhou ◽  
K. Aung

One of the most fundamental concepts in heat transfer is the convective heat transfer coefficient, which is closely related with the flow Reynolds number, flow geometry and the thermal conditions on the heat transfer surface. To define the heat transfer coefficient, a reference temperature is needed besides the surface temperature and heat flux. The reference temperature can be chosen differently, such as the fluid bulk mean temperature (for internal flows) and the temperature at the far field (for external flows). For complicated flows, the adiabatic wall temperature, defined as the wall temperature when the surface heat flux is zero, is commonly adopted as the reference temperature. Other options can also be applied to complicated flows. This paper analyzed some of the potential selections of the reference temperature for different flow settings, including film cooling, jet impingement with cross flows and a mixing flow in a straight duct with or without internal heat source. Both laminar and turbulent flows are considered with different boundary conditions. Dramatic changes of heat transfer coefficient are observed with different reference temperatures. In some special conditions the heat transfer coefficient becomes negative, which means the heat flux has a different direction with the driving temperature difference defined. An innovative method is proposed to calculate the heat transfer coefficient of complicated flows with constant surface temperature.


2019 ◽  
Vol 36 (9) ◽  
pp. 3180-3207
Author(s):  
Vishweshwara P.S. ◽  
Harsha Kumar M.K. ◽  
N. Gnanasekaran ◽  
Arun M.

Purpose Many a times, the information about the boundary heat flux is obtained only through inverse approach by locating the thermocouple or temperature sensor in accessible boundary. Most of the work reported in literature for the estimation of unknown parameters is based on heat conduction model. Inverse approach using conjugate heat transfer is found inadequate in literature. Therefore, the purpose of the paper is to develop a 3D conjugate heat transfer model without model reduction for the estimation of heat flux and heat transfer coefficient from the measured temperatures. Design/methodology/approach A 3 D conjugate fin heat transfer model is solved using commercial software for the known boundary conditions. Navier–Stokes equation is solved to obtain the necessary temperature distribution of the fin. Later, the complete model is replaced with neural network to expedite the computations of the forward problem. For the inverse approach, genetic algorithm (GA) and particle swarm optimization (PSO) are applied to estimate the unknown parameters. Eventually, a hybrid algorithm is proposed by combining PSO with Broyden–Fletcher–Goldfarb–Shanno (BFGS) method that outperforms GA and PSO. Findings The authors demonstrate that the evolutionary algorithms can be used to obtain accurate results from simulated measurements. Efficacy of the hybrid algorithm is established using real time measurements. The hybrid algorithm (PSO-BFGS) is more efficient in the estimation of unknown parameters for experimentally measured temperature data compared to GA and PSO algorithms. Originality/value Surrogate model using ANN based on computational fluid dynamics simulations and in-house steady state fin experiments to estimate the heat flux and heat transfer coefficient separately using GA, PSO and PSO-BFGS.


Author(s):  
Huasheng Wang ◽  
John W. Rose

The paper presents a theoretical model to predict film condensation heat transfer from a vapor flowing in a horizontal tube with equilateral triangular section minichannels or microchannels. The model is based on fundamental analysis which assumes laminar condensate flow on the channel walls and takes account of surface tension, vapor shear stress and gravity. The case considered here is where the channel wall temperature is uniform and the vapor is saturated at inlet. Sample numerical results are given for the channel size (side of triangle) of 1.0 mm and for refrigerant R134a. The general behaviour of the condensate flow pattern (spanwise and streamwise profiles of the condensate film), as well as streamwise variation in quality and local mean (over section perimeter) heat-transfer coefficient, are qualitatively in accord with expectations on physical grounds.


Author(s):  
E.-L. Pelletier ◽  
L. K. H. Leung ◽  
A. Teyssedou ◽  
R. Girard

Two methodologies to predict film boiling heat-transfer coefficient have been assessed against experimental wall-temperature measurements obtained under steady-state conditions with water flow inside vertical tubes. One of these methodologies employs heat flux as the independent parameter while the other uses wall temperature as the independent parameter. The film boiling heat transfer consists of developing and fully developed film-boiling regions. Film boiling heat-transfer coefficients are predicted using the film boiling look-up tables for fully developed flow. Developing film-boiling effect is accounted for using modification factors to the fully developed heat-transfer coefficient. Wall-temperature distributions along uniformly heated tubes were established using a semi-analytical scheme and compared against measurements. Both methodologies have provided good predictions. However, the overall prediction accuracy for the heat-flux-based correlation is slightly better than that for the wall-temperature-based correlation. Wall temperatures predicted with the heat-flux-based correlation follow closely measurements at the developing post-dryout region. The wall superheat correlation predicts a sharp temperature rise once the critical heat flux is exceeded, resulting in discrepancies between predictions and measurements of wall temperature and overpredictions of the maximum temperature. The wall-temperature-based modification factor for the developing film-boiling region has been revised using the tube heat-transfer database to improve the prediction accuracy of the wall temperature.


Sign in / Sign up

Export Citation Format

Share Document