Visualization of the Evaporating Liquid-Vapor Interface in Micropillar Arrays

2016 ◽  
Vol 138 (2) ◽  
Author(s):  
Dion S. Antao ◽  
Solomon Adera ◽  
Edgardo Farias ◽  
Rishi Raj ◽  
Evelyn N. Wang

We captured interesting static and dynamic behavior of the liquid-vapor interface in well-defined silicon micropillar arrays during thermally driven evaporation of water from the microstructured surface. The 3-D shape of the meniscus was characterized via laser interferometry where bright and dark fringes result from the interference of incident and reflected monochromatic light due to a variable thickness thin liquid film (FIG. 1). During steady state evaporation experiments, water was supplied to the sample with a syringe pump at 10 μL/min. FIG. 2a and 2b show a SEM image of a typical fabricated micropillar array and a schematic of the experimental setup, respectively. When water wicks through the micropillar array, the meniscus in a unit cell (four pillars in FIG. 1) assumes an equilibrium shape depending on the location from the liquid source/reservoir and the ambient conditions (ambient evaporation at Qin = 0 W). At this point, the meniscus is pinned at the top of the pillars. As the evaporation rate increases due the applied heat flux, the meniscus increases in curvature, thus increasing the capillary pressure to sustain the higher evaporation rate. This is evidenced by the increasing number of fringes in the unit cell when Qin is increased (0 W, 0.11 W, 0.44 W, and 0.99 W, FIG. 1a-1d respectively). Beyond a maximum curvature, the meniscus de-pins from the pillar top surface and recedes within the unit cell. This occurs when the capillary pressure generated at this curvature, cannot balance the viscous loss resulting from flow through the micropillar array. We observed that this receding shape was independent of the applied heat, and only depended on the micropillar array geometry and the intrinsic wettability of the material. Representative meniscus profiles along the diagonal direction of the unit cell obtained from image analysis of FIG. 1 at various Qin are shown in FIG. 2c.

2020 ◽  
Vol 142 (12) ◽  
Author(s):  
Brian E. Fehring ◽  
Roman W. Morse ◽  
Jason Chan ◽  
Kristofer M. Dressler ◽  
Evan T. Hurlburt ◽  
...  

Abstract Instantaneous temperature measurements at the interface between a solid wall and a thin, unsteady liquid film are performed using thermoreflectance, a nonintrusive optical technique with high temporal resolution. A laser beam is directed at a wall–liquid interface, and the intensity of the light reflected at that interface is measured by a photodiode. The intensity of the reflected light varies with the index of refraction of the liquid at the wall. The index of refraction is a function of temperature, which enables the instantaneous measurement of the wall temperature. In the presence of thin liquid films, reflections from the liquid–vapor interface at the free surface of the film generate noise in the measurements. We demonstrate that orienting the laser beam at a large incident angle, close to total internal reflection, minimizes noise from the liquid–vapor interface while increasing the sensitivity of the measurement. The thermoreflectance technique is validated in an unsteady two-phase annular flow. Measurements of temperature fluctuations less than 1 K in amplitude are achieved, with an uncertainty of 0.1 K.


2014 ◽  
Vol 136 (8) ◽  
Author(s):  
Tingting Hao ◽  
Xuehu Ma ◽  
Zhong Lan ◽  
Nan Li ◽  
Yuzhe Zhao

The effects of superhydrophobic surface and superhydrophobic and superhydrophilic hybrid surface on the fluid flow and heat transfer of oscillating heat pipes (OHPs) were investigated in the paper. The inner surfaces of the OHPs were hydrophilic surface (copper), hybrid surface (superhydrophilic evaporation and superhydrophobic condensation section), and uniform superhydrophobic surface, respectively. Deionized water was used as the working fluid. Experimental results showed that superhydrophobic surface influenced the slug motion and thermal performance of OHPs. Visualization results showed that the liquid-vapor interface was concave in the OHP with copper surface. A thin liquid film existed between the vapor plug and the wall of the OHP. On the contrary, the liquid-vapor interface took a convex profile in the OHP with superhydrophobic surface and the liquid-vapor interface contact line length in the hybrid surface OHP was longer than that in the uniform superhydrophobic surface OHP. The liquid slug movements became stronger in the hybrid surface OHPs as opposed to the copper OHP, while the global heat transfer performance of the hybrid surface OHPs increased by 5–20%. Comparing with the copper OHPs, the maximum amplitude and velocity of the liquid slug movements in the hybrid surface OHPs increased by 0–127% and 0–185%, respectively. However, the maximum amplitude and velocity of the liquid slug movements in the uniform superhydrophobic OHPs was reduced by 0–100% and 0–100%, respectively. The partial dryout phenomenon took place in OHPs with uniform superhydrophobic surface. The liquid slug movements became weaker and the thermal resistance was increased by 10–35% in the superhydrophobic surface OHPs.


Author(s):  
Jingzhi Zhang ◽  
Wei Li ◽  
Tom I.-P. Shih ◽  
Yonghai Zhang ◽  
Yanping Shi ◽  
...  

Heat transfer and pressure drop characteristics of condensation for R410A inside horizontal tubes (dh = 3.78 mm) under normal and micro gravity are investigated numerically. The Volume of Fluid method is used to acquire liquid-vapor interface, while the low-Reynolds form of the Shear Stress Transport k∼ω (SST k∼ω) model is adopted to taking turbulent effect into account. The results indicate that the heat transfer coefficients decrease with increasing gravity accelerations, while the frictional pressure gradients increase with increases in gravity accelerations. The liquid film accumulates at the bottom of the tube, leading to a very thin liquid film attached to the upper part of inner tube wall. This accumulation effect decreases with decreases in gravitational accelerations. A more symmetrical liquid-vapor interface is obtained at lower gravity. The average liquid film thickness is nearly the same for different gravity accelerations at the same vapor quality (δave≈56 μm at x = 0.9 and δave≈230 μm at x = 0.5). The local heat transfer coefficients increase with increasing gravity at the top of the tube and decrease with increases in gravity at the bottom, while the bottom part of the tube has a limited contribution to the global heat transfer coefficient for stratified flow regime. The numerical data obtained under normal gravity agree well with well-known empirical correlations.


Author(s):  
Arnoldo A. Badillo

Although many years have past from the pioneer work of Lord Rayleigh [1] on bubble growth in the inertia controlled regime and later from Scriven [2], Plesset and Zwick [3] for the diffusion controlled regime, we are still missing mathematical model able to predict accurately both situations. Advances in computational power open the possibility of exploring up-close the transport phenomena in the vicinity of the liquid-vapor interface at an unprecedented resolution. Nonetheless, a high numerical resolution is not enough to fully solve the general problem of bubble growth. New models based on a sharp-interface interpretation of the liquid-vapor interface, have proven to provide accurate results in the diffusion controlled regime, however, these models must assume the interface temperature at the saturation value, restricting their application to physical situations where the evaporation rate satisfies the Stefan condition and bubbles are big enough as to neglect the curvature effects in the interface temperature. In an attempt to provide a more general framework to study bubble growth, a new phase-field model has recently been derived, where no assumption is made on the interface temperature. In this new model, the evaporation rate depends on the local interface temperature and not directly on the heat balance at the liquid-vapor interface. In principle, this particular feature of the model should allow us to simulate both, the inertia and diffusion controlled regimes, but the model has only been validated for the latter. The next step in the validation process is the simulation of bubble growth under convective conditions. Experiments of single bubbles growing and rising up under normal gravity conditions have shown that the growth exponent is about 0.8, in contrast to the value of 1.0 for the inertial controlled regime and 0.5 for the diffusion controlled regime. In this work, fully three dimensional phase-field simulations of bubble growth under convective conditions are presented, where the predicted bubble size and growth exponent compare very well to experimental observations.


Langmuir ◽  
2016 ◽  
Vol 32 (2) ◽  
pp. 519-526 ◽  
Author(s):  
Dion S. Antao ◽  
Solomon Adera ◽  
Yangying Zhu ◽  
Edgardo Farias ◽  
Rishi Raj ◽  
...  

2001 ◽  
Author(s):  
Qiang Bai ◽  
V. K. Dhir

Abstract Deposition of boron on the fuel rod cladding during boiling of water containing boron can depress the neutron flux and lead to a decrease in nuclear reactor power output. There is practically little precise information on the temperature field, the gradients of chemical concentration and deposition of boron on the cladding surface. The objective of the present work is to simulate the nucleate boiling process along with velocity, temperature and concentration fields of aqueous boron in the vicinity of the cladding of a fuel rod. As a first step in solving the complete problem, two-dimensional numerical simulation of a bubble growth on a horizontal surface is considered. A finite difference scheme is used to solve the equations governing conservation of mass, momentum, energy and species concentration. The calculation domain is divided into macro and micro regions. In macro-region, the governing equations are used to calculate the distributions of velocity, temperature, and concentration. The Level Set method is used to capture the evolving liquid-vapor interface. For micro-region, lubrication theory is used, which includes the disjoining pressure in the thin liquid film. The solutions for micro-region and macro-region are matched at the outer edge of the micro-layer. A dilute aqueous Boron solution is considered in the simulation. From numerical simulations, the dynamic change in concentration distribution of boron during the bubble growth shows that the precipitation of boron can occur near the advancing and receding liquid-vapor interface when the ambient boron concentration level is 0.003.


Author(s):  
Didier Jamet ◽  
Olivier Lebaigue ◽  
Jean-Marc Delhaye ◽  
N. Coutris

Sign in / Sign up

Export Citation Format

Share Document