Numerical Simulation of Saturated Film Boiling on a Horizontal Surface

1997 ◽  
Vol 119 (3) ◽  
pp. 525-533 ◽  
Author(s):  
G. Son ◽  
V. K. Dhir

The past efforts in applying linear Taylor instability theory to the prediction of heat transfer during film boiling on a horizontal surface have suffered from the fact that empirical correlations must be used to define the shape of vapor-liquid interfaces and to determine the transport of mass and heat across these interfaces. The objective of this study is to clarify the physics of film boiling and to predict heat transfer coefficients through complete numerical simulation of the evolving interface between superposed layers of immiscible fluids. A coordinate transformation technique supplemented by a numerical grid generation method and a second-order projection method are combined to solve for the flow and temperature fields associated with an evolving interface. From the numerical simulation, the film thickness and, in turn, the heat transfer coefficient are found to vary both spatially and temporally. Increased wall superheat not only thickens the vapor film in the valley but also enlarges the vapor bulge. The effect of increased system pressure is to slow down the growth of the interface.

1999 ◽  
Vol 121 (3) ◽  
pp. 623-631 ◽  
Author(s):  
G. Son ◽  
V. K. Dhir ◽  
N. Ramanujapu

In this study, a complete numerical simulation of a growing and departing bubble on a horizontal surface has been performed. A finite difference scheme is used to solve the equations governing conservation of mass, momentum, and energy in the vapor-liquid layers. The vapor-liquid interface is captured by a level set method which is modified to include the influence of phase change at the liquid-vapor interphase. The disjoining pressure effect is included in the numerical analysis to account for heat transfer through the liquid microlayer. From the numerical simulation, the location where the vapor-liquid interface contacts the wall is observed to expand and then retract as the bubble grows and departs. The effect of static contact angle and wall superheat on bubble dynamics has been quantified. The bubble growth predicted from numerical analysis has been found to compare well with the experimental data reported in the literature and that obtained in this work.


Author(s):  
Qiusheng Liu ◽  
Masahiro Shiotsu ◽  
Akira Sakurai ◽  
Katsuya Fukuda

Forced convection film boiling heat transfer from a horizontal cylinder in water and Freon-113 flowing upward perpendicular to the cylinder under subcooled conditions was measured for the flow velocities from 0 to 1 m/s at the system pressures ranging from 100 to 500 kPa: the platinum horizontal cylinders with diameters ranging from 0.7 to 5 mm were used as the test heaters. The film boiling heat transfer coefficients were obtained for the surface superheats from about 800 K for water and from about 400 K for Freon-113 down to minimum film boiling surface superheats. These heat transfer coefficients increase with the increase in flow velocity, liquid subcooling, system pressure, and with the decrease in cylinder diameter. A correlation for subcooled forced convection film boiling heat transfer was presented, which can describe the experimental data obtained within ±20% for the flow velocities below 0.7 m/s, and within −30% to +20% for the higher flow velocities. The correlation also predicted well the data by Shigechi (1983), Motte and Bromley (1957), and Sankaran and Witte (1990) obtained for the larger diameter cylinders and higher flow velocities in various liquids at the pressures of near atmospheric. The Shigechi’s data were in the range from about −20% to +15%, the data of Motte and Bromley were about ±30%, and the data of Sankaran and Witte were within +20% of the curves given by the corresponding predicted values.


2013 ◽  
Vol 34 (1) ◽  
pp. 5-16 ◽  
Author(s):  
Jozef Cernecky ◽  
Jan Koniar ◽  
Zuzana Brodnianska

Abstract The paper deals with a study of the effect of regulating elements on local values of heat transfer coefficients along shaped heat exchange surfaces with forced air convection. The use of combined methods of heat transfer intensification, i.e. a combination of regulating elements with appropriately shaped heat exchange areas seems to be highly effective. The study focused on the analysis of local values of heat transfer coefficients in indicated cuts, in distances expressed as a ratio x/s for 0; 0.33; 0.66 and 1. As can be seen from our findings, in given conditions the regulating elements can increase the values of local heat transfer coefficients along shaped heat exchange surfaces. An optical method of holographic interferometry was used for the experimental research into temperature fields in the vicinity of heat exchange surfaces. The obtained values correspond very well with those of local heat transfer coefficients αx, recorded in a CFD simulation.


1987 ◽  
Vol 109 (3) ◽  
pp. 599-605 ◽  
Author(s):  
An-Shik Yang ◽  
Ching-Chang Chieng

An anisotropic factor is carefully selected from eleven distributions and adopted to the k–ε two-equation model of turbulence to obtain detailed velocity and temperature fields for steady-state, fully developed turbulent flow through infinite triangular/square rod array. The present study covers the ranges of pitch-to-diameter ratio from 1.123 to 1.5, and Reynolds number from 2.4 × 104 to 106. Velocity and wall shear stress are calculated and compared to experimental data. Normalized fluid temperature, friction factor, and heat transfer coefficient are also computed. The correlations of friction factor and heat transfer coefficients for flow inside circular pipe and flow through finite rod arrays are compared with the results for flow through infinite rod arrays.


2014 ◽  
Vol 697 ◽  
pp. 235-238
Author(s):  
Gang Wu ◽  
Can Chao Huang ◽  
Hong Ling Qin ◽  
Chun Hua Zhao

Using the basic principle of heat transfer, tribology and numerical simulation, a two-dimensional heat transfer model of the three-layer composite brake pair materials were established. The temperature fields of brake pairs during the process of friction were analyzed. Applied given heat loads at different time node on the brake pair model, the temperatures of different bicycle brake pairs were compared and analyzed. Results show that the improved surface structures of brake pair have positive effect on decreasing the temperature of contact areas than that of ordinary surface structure.


2015 ◽  
Vol 764 ◽  
pp. 362-394 ◽  
Author(s):  
T. Dairay ◽  
V. Fortuné ◽  
E. Lamballais ◽  
L.-E. Brizzi

AbstractDirect numerical simulation (DNS) of an impinging jet flow with a nozzle-to-plate distance of two jet diameters and a Reynolds number of 10 000 is carried out at high spatial resolution using high-order numerical methods. The flow configuration is designed to enable the development of a fully turbulent regime with the appearance of a well-marked secondary maximum in the radial distribution of the mean heat transfer. The velocity and temperature statistics are validated with documented experiments. The DNS database is then analysed focusing on the role of unsteady processes to explain the spatial distribution of the heat transfer coefficient at the wall. A phenomenological scenario is proposed on the basis of instantaneous flow visualisations in order to explain the non-monotonic radial evolution of the Nusselt number in the stagnation region. This scenario is then assessed by analysing the wall temperature and the wall shear stress distributions and also through the use of conditional averaging of velocity and temperature fields. On one hand, the heat transfer is primarily driven by the large-scale toroidal primary and secondary vortices emitted periodically. On the other hand, these vortices are subjected to azimuthal distortions associated with the production of radially elongated structures at small scale. These distortions are responsible for the appearance of very high heat transfer zones organised as cold fluid spots on the heated wall. These cold spots are shaped by the radial structures through a filament propagation of the heat transfer. The analysis of probability density functions shows that these strong events are highly intermittent in time and space while contributing essentially to the secondary peak observed in the radial evolution of the Nusselt number.


1987 ◽  
Vol 109 (4) ◽  
pp. 912-918 ◽  
Author(s):  
J. R. Parsons ◽  
M. L. Arey

Experiments have been performed which describe the transient development of natural convective flow from both a single and two vertically aligned horizontal cylindrical heat sources. The temperature of the wire heat sources was monitored with a resistance bridge arrangement while the development of the flow field was observed optically with a Mach–Zehnder interferometer. Results for the single wire show that after an initial regime where the wire temperature follows pure conductive response to a motionless fluid, two types of fluid motion will begin. The first is characterized as a local buoyancy, wherein the heated fluid adjacent to the wire begins to rise. The second is the onset of global convective motion, this being governed by the thermal stability of the fluid layer immediately above the cylinder. The interaction of these two motions is dependent on the heating rate and relative heat capacities of the cylinder and fluid, and governs whether the temperature response will exceed the steady value during the transient (overshoot). The two heat source experiments show that the merging of the two developing temperature fields is hydrodynamically stabilizing and thermally insulating. For small spacing-to-diameter ratios, the development of convective motion is delayed and the heat transfer coefficients degraded by the proximity of another heat source. For larger spacings, the transient behavior approaches that of a single isolated cylinder.


Sign in / Sign up

Export Citation Format

Share Document