Numerical Simulation of Laminar Liquid Film Condensation in a Horizontal Circular Minichannel

2012 ◽  
Vol 134 (5) ◽  
Author(s):  
E. Da Riva ◽  
D. Del Col

A three-dimensional volume of fluid (VOF) simulation of condensation of R134a inside a 1 mm i.d. minichannel is presented. The minichannel is horizontally oriented and the effect of gravity is taken into account. Simulations have been run both with and without taking into account surface tension. A uniform interface temperature and a uniform wall temperature have been fixed as boundary conditions. The mass flux is G = 100 kg m−2 s−1 and it has been assumed that the flow is laminar inside the liquid phase while turbulence inside the vapor phase has been handled by a modified low Reynolds form of the k–ω model. The fluid is condensed till reaching 0.45 vapor quality. The flow is expected to be annular without the presence of waves, therefore the problem was treated as steady state. Computational results displaying the evolution of vapor–liquid interface and heat transfer coefficient are reported and validated against experimental data. The condensation process is found to be gravity dominated, while the global effect of surface tension is found to be negligible. At inlet, the liquid film is thin and evenly distributed all around the tube circumference. Moving downstream the channel, the film thickness remains almost constant in the upper half of the minichannel, while the film at the bottom of the pipe becomes thicker because the liquid condensed at the top is drained by gravity to the bottom.

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.


2007 ◽  
Vol 588 ◽  
pp. 75-129 ◽  
Author(s):  
ROBERT T. COLLINS ◽  
MICHAEL T. HARRIS ◽  
OSMAN A. BASARAN

Breakup of electrified jets is important in applications as diverse as electrospraying, electroseparations and electrospray mass spectrometry. Breakup of a perfectly conducting, incompressible Newtonian liquid jet surrounded by a passive insulating gas that is stressed by a radial electric field is studied by a temporal analysis. An initially quiescent jet is subjected to axially periodic shape perturbations and the ensuing dynamics are followed numerically until pinch-off by both a three-dimensional but axisymmetric (two-dimensional) and a one-dimensional slender-jet algorithm. Results computed with these algorithms are verified against predictions from linear theory for short times. Breakup times, ratios of the sizes of the primary to satellite drops formed at pinch-off, and the Coulombic stability of these drops are reported over a wide range of electrical Bond numbers, NE (ratio of electric to surface tension force), Ohnesorge numbers, NOh (ratio of viscous to surface tension force), and disturbance wavenumbers, k. Effect of surface charge on interface overturning is investigated. Furthermore, the influence of electrostatic stresses on the dynamics of pinch-off and the mechanisms of satellite drop formation is also addressed.


1974 ◽  
Vol 75 (3) ◽  
pp. 405-426 ◽  
Author(s):  
P. F. Rhodes-Robinson

AbstractIn this paper wave-maker theory including the effect of surface tension is determined for three-dimensional motion of water in a semi-infinite rectangular channel with outgoing surface wave modes allowed for at infinity; the motion is generated by a harmonically oscillating vertical plane wave-maker at the end of the channel and the cases of both infinite and finite constant depth are treated. The solution of the boundary-value problem for the velocity potential is more complicated in the presence of surface tension due mainly to the additional effect of the channel walls at which the normal free surface slopes are prescribed—as also is the slope at the wave-maker—to ensure uniqueness. The simpler three-dimensional solution for a semi-infinite region—obtained long ago by Sir Thomas Havelock in the absence of surface tension for the case of infinite depth—is also noted.


Author(s):  
E. Da Riva ◽  
D. Del Col ◽  
A. Cavallini

3D steady-state Volume of Fluid (VOF) simulations of condensation of R134a inside a 1 mm i.d. circular minichannel are proposed. The minichannel is horizontally oriented and both the effects of gravity and surface tension are taken into account. A uniform interface temperature, as well as a uniform wall temperature, were fixed as boundary conditions in order to model the phase change process. Simulations were performed at mass fluxes G = 100 kg m−2s−1 and G = 400 kg m−2s−1. It has been assumed that the flow was turbulent inside the vapour core, while for the condensate film two different computational approaches have been considered. The first approach (i.e. “laminar liquid film”) corresponds to the assumption that the flow is laminar inside the liquid phase and turbulent inside the vapour phase. For the second approach (i.e. “turbulent liquid film”), instead, a low Reynolds form of a turbulence model has been used through both phases. The aforementioned approaches are compared to each other and the computed heat transfer coefficients are compared against experimental data by Matkovic et al.: data at low mass flux is well predicted by the “laminar liquid film” approach and overpredicted by “turbulent liquid film” approach, while data at high mass flux is underpredicted by the former approach and well predicted by the latter. The evolution of the vapour-liquid interface along the minichannel, as well as the velocity field, are reported. Besides, the computed cross sectional void fraction is compared against empirical correlations available in the literature.


Sign in / Sign up

Export Citation Format

Share Document