Evaporation and flow dynamics of thin, shear-driven liquid films in microgap channels

Thin and ultra-thin shear-driven liquid films in a narrow channel are a promising candidate for the thermal management of advanced semiconductor devices in earth and space applications. Such flows experience complex, and as yet poorly understood, two-phase flow phenomena requiring significant advances in fundamental research before they could be broadly applied. This paper focuses on the results obtained in experiments with locally heated shear-driven liquid films in a flat mini-channel. A detailed map of the flow sub-regimes in a shear-driven liquid film flow of water and FC-72 have been obtained for a 2 mm channel operating at room temperature. While the water film can be smooth under certain liquid/gas flow rates, the surface of an intensively evaporating film of FC-72 is always distorted by a pattern of waves and structures. It was found, that when heated the shear-driven liquid films are less likely to rupture than gravity-driven liquid films. For shear-driven water films the critical heat flux was found of up to 10 times higher than that for a falling film, which makes shear-driven films (annular or stratified two-phase flows) more suitable for cooling applications than falling liquid films.     

Simulation and measurement of 3D shear-driven thin liquid film flow in a duct

Three-dimensional flow behavior of thin liquid film that is shear-driven by turbulent air flow in a duct is measured and simulated. Its film thickness and width are reported as a function of air velocity, liquid flow rate, surface tension coefficient, and wall contact angle. The numerical component of this study is aimed at exploring and assessing the suitability of utilizing the FLUENT-CFD code and its existing components, i.e. Volume of Fluid model (VOF) along with selected turbulence model, for simulating the behavior of 3D shear-driven liquid film flow, through a comparison with measured results. The thickness and width of the shear-driven liquid film are measured using an interferometric technique that makes use of the phase shift between the reflections of incident light from the top and bottom surfaces of the thin liquid film. Such measurements are quite challenging due to the dynamic interfacial instabilities that develop in this flow. The results reveal that higher air flow velocity decreases the liquid film thickness but increases its width, while higher liquid flow rate increases both its thickness and width. Simulated results provide good estimates of the measured values, and reveal the need for considering a dynamic rather than a static wall contact angle in the model for improving the comparison with measured values.     

Molecular dynamics simulation of liquid argon flow at platinum surfaces

The micro Poiseuille flow for liquid argon flowing in a nanoscale channel formed by two solid walls was studied in the present paper. The solid wall material was selected as platinum, which has well established interaction potential. We consider the intermolecular force not only among the liquid argon molecules, but also between the liquid argon atoms and the solid wall particles, therefore three regions, i.e. the liquid argon computation domain, the top and bottom solid wall regions are included for the force interaction. The present MD (Molecular Dynamics) simulation was performed without any assumptions at the wall surface. The objective of the study is to find how the flow and the slip boundaries at the wall surface are affected by the applied gravity force, or the shear rate. The MD simulations are performed in a nondimensional unit system, with the periodic boundary conditions applied except in the channel height direction. Once the steady state is reached, the macroscopic parameters are evaluated using the statistical mechanics approach. For all the cases tested numerically in the present paper, slip boundaries occur, and such slip velocity at the stationary wall surface increases with increasing the applied gravity force, or the shear rate. The slip length, which is defined as the distance that the liquid particles shall travel beyond the wall surfaces to reach the same velocity as the wall surface, sharply decreases at small shear rate, then slightly decreases with increasing the applied shear rate. We observe that the liquid viscosity remains nearly constant at small shear rates, and the Newtonian flow occurs. However, with increasing the shear rate, the viscosity increases and the non-Newtonian flow appears.     

Buoyancy Convection in Low Prandtl Number Liquids with Large Temperature Variation

Two-dimensional laminar convection in low Prandtl number liquids driven by the buoyancy force is studied. The liquid is contained in a closed square cavity with isothermal vertical walls kept at different temperatures. The top and bottom walls are assumed to be insulated. The thermal conductivity of the liquid is assumed to depend on temperature. ADI and SOR schemes are employed. The heat transfer is found to decrease appreciably across the cavity with a decrease in thermal conductivity.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Natural convection flow of water near its density maximum in a rectangular enclosure having isothermal walls with heat generation

We consider unsteady laminar natural convection flow of water subject to density inversion in a rectangular cavity formed by isothermal vertical walls with internal heat generation. The top and bottom horizontal walls are considered to be adiabatic, whereas the temperature of the left vertical wall is assumed to be greater than that of the right vertical wall. The equations are non-dimensionalized and are solved numerically by an upwind finite difference method together with a successive over-relaxation (SOR) technique. The effects of both heat generation and variations in the aspect ratio on the streamlines, isotherms and the rate of heat transfer from the walls of the enclosure are presented. Investigations are performed for water taking Prandtl number to be Pr=11.58 and the Rayleigh number to be Ra=10⁵.     

Buoyancy-driven motion of a two-dimensional bubble or drop through a viscous liquid in the presence of a vertical electric field

The effect of an electric field on the buoyancy-driven motion of a two-dimensional gas bubble rising through a quiescent liquid is studied computationally. The dynamics of the bubble is simulated numerically by tracking the gas–liquid interface when an electrostatic field is generated in the vertical gap of the rectangular enclosure. The two phases of the system are assumed to be perfect dielectrics with constant but different permittivities, and in the absence of impressed charges, there is no free charge in the fluid bulk regions or at the interface. Electric stresses are supported at the bubble interface but absent in the bulk and one of the objectives of our computations is to quantify the effect of these Maxwell stresses on the overall bubble dynamics. The numerical algorithm to solve the free-boundary problem relies on the level-set technique coupled with a finite-volume discretization of the Navier–Stokes equations. The sharp interface is numerically approximated by a finite-thickness transition zone over which the material properties vary smoothly, and surface tension and electric field effects are accounted for by employing a continuous surface force approach. A multi-grid solver is applied to the Poisson equation describing the pressure field and the Laplace equation governing the electric field potential. Computational results are presented that address the combined effects of viscosity, surface tension, and electric fields on the dynamics of the bubble motion as a function of the Reynolds number, gravitational Bond number, electric Bond number, density ratio, and viscosity ratio. It is established through extensive computations that the presence of the electric field can have an important effect on the dynamics. We present results that show a substantial increase in the bubble’s rise velocity in the electrified system as compared with the corresponding non-electrified one. In addition, for the electrified system, the bubble shape deformations and oscillations are smaller, and there is a reduction in the propensity of the bubble to break up through increasingly larger oscillations.     

Effect of magnetic field on the buoyancy and thermocapillary driven convection of an electrically conducting fluid in an annular enclosure

The main objective of this article is to study the effect of magnetic field on the combined buoyancy and surface tension driven convection in a cylindrical annular enclosure. In this study, the top surface of the annulus is assumed to be free, and the bottom wall is insulated, whereas the inner and the outer cylindrical walls are kept at hot and cold temperatures respectively. The governing equations of the flow system are numerically solved using an implicit finite difference technique. The numerical results for various governing parameters of the problem are discussed in terms of the streamlines, isotherms, Nusselt number and velocity profiles in the annuli. Our results reveal that, in tall cavities, the axial magnetic field suppresses the surface tension flow more effectively than the radial magnetic field, whereas, the radial magnetic field is found to be better for suppressing the buoyancy driven flow compared to axial magnetic field. However, the axial magnetic field is found to be effective in suppressing both the flows in shallow cavities. From the results, we also found that the surface tension effect is predominant in shallow cavities compared to the square and tall annulus. Further, the heat transfer rate increases with radii ratio, but decreases with the Hartmann number.     

Numerical comparison of conjugate and non-conjugate natural convection for internally heated semi-circular pools

In this numerical investigation, a detailed comparison of the conjugate and non-conjugate natural convection within a semi-cylindrical cavity has been presented. The cavity is assumed to be filled with a fluid containing uniformly distributed internal heating sources. The bottom circular wall of the cavity is taken to be thick with finite conductive properties, while the top wall is considered to be isothermal. The Navier–Stokes and energy equations are solved numerically by using the SIMPLER algorithm. A Rayleigh number range from 3.2×10⁶ to 3.2×10¹¹ has been investigated and the effects of solid-to-fluid conductivity ratios of 1.0, 5.0 and 23.0 have been analysed. The present numerical results for a semi-circular cavity with entirely isothermal walls are compared with known results from the open literature. It was found that these results for the non-conjugate problem are in very good agreement. The present results for a conjugate cavity show a remarkable difference from the non-conjugate analysis. The average Nusselt number for the solid–fluid interface shows a decrease while the top wall average Nu number has increased. It has been concluded that these effects increase for a system with a low solid-to-fluid conductivity ratio. It is evident from the present conjugate results that the assumption of isothermal enclosing walls gives somewhat different results when the walls are thick and the solid-to-fluid conductivity ratio is small.     

Unsteady flow and heat transfer of porous media sandwiched between viscous fluids

The problem of unsteady oscillatory flow and heat transfer of porous medin sandwiched between viscous fluids has been considered through a horizontal channel with isothermal wall temperatures. The flow in the porous medium is modeled using the Brinkman equation. The governing partial differential equations are transformed to ordinary differential equations by collecting the non-periodic and periodic terms. Closed-form solutions for each region are found after applying the boundary and interface conditions. The influence of physical parameters, such as the porous parameter, the frequency parameter, the periodic frequency parameter, the viscosity ratios, the conductivity ratios, and the Prandtl number, on the velocity and temperature fields is computed numerically and presented graphically. In addition, the numerical values of the Nusselt number at the top and bottom walls are derived and tabulated.     

Exchange Processes in a Channel with Two Vertical Emerged Obstructions

The flow around two vertical obstructions situated in a long open channel with lateral walls is simulated using Large Eddy Simulation (LES). The emerged obstructions (groynes) are oriented perpendicular to the channel side walls. The incoming channel flow is fully turbulent. The focus of the present paper is to examine the mass transfer between a contaminant situated initially in the (embayment) region between the obstructions and the main channel. The mass transfer is simulated using a passive (conserved) scalar transport equation. The scalar is introduced instantaneously inside the embayment. The eddy structures that populate the detached shear layer starting at the tip of the upstream obstruction and their interaction with the coherent structures inside the embayment area are shown to play an important role in the contaminant entrainment from the embayment area into the main channel. It is also found that the flow convected into the embayment from the region downstream of the second obstruction can substantially accelerate the removal of the contaminant from the embayment region. A detailed analysis of the contaminant flux variation in the top, middle and bottom layers inside the embayment region is carried out to better understand how the contaminant exits the embayment area, the role played by vertical motions within the embayment, and the effects induced by the presence of the bottom surface which delays the contaminant purging from the bottom layer. It is found that one half of the contaminant mass situated initially in the bottom third of the embayment volume is not leaving the embayment through the corresponding embayment–channel interface. The opposite is observed for the top third layer where the mass of contaminant leaving through the top interface area is 50% higher than the corresponding mass of contaminant initially situated in the top layer. This shows that the mass exchange is highly non-uniform over the depth and there is an overall contaminant flux within the embayment toward the free surface. The decay of contaminant mass within the embayment is calculated enabling estimation of a global 1D exchange coefficient based on dead-zone theory. It is found that though dead-zone models can relatively accurately describe the contaminant mass decay in time within the embayment, the mass exchange is not characterized by a unique value of the exchange coefficient over the whole length of the ejection process for the present geometry and flow conditions. Rather, two distinct phases of the decay process are identified. In the initial phase of decay over which about 68% of the total mass of contaminant leaves the embayment, the exchange coefficient is found to be about twice the value estimated for the final phase.     

Turbulent breaking of overturning gravity waves below a critical level

The interaction of an internal gravity wave with its evolving critical layer and the subsequent generation of turbulence by overturning waves are studied by three-dimensional numerical simulations. The simulation describes the flow of a stably stratified Boussinesq fluid between a bottom wavy surface and a top flat surface, both without friction and adiabatic. The amplitude of the surface wave amounts to about 0.03 of the layer depth. The horizontal flow velocity is negative near the lower surface, positive near the top surface with uniform shear and zero mean value. The bulk Richardson number is one. The flow over the wavy surface induces a standing gravity wave causing a critical layer at mid altitude. After a successful comparison of a two-dimensional version of the model with experimental observations (Thorpe [21]), results obtained with two different models of viscosity are discussed: a direct numerical simulation (DNS) with constant viscosity and a large-eddy simulation (LES) where the subgrid scales are modelled by a stability-dependent first-order closure. Both simulations are similar in the build-up of a primary overturning roll and show the expected early stage of the interaction between wave and critical level. Afterwards, the flows become nonlinear and evolve differently in both cases: the flow structure in the DNS consists of coherent smaller-scale secondary rolls with increasing vertical depth. On the other hand, in the LES the convectively unstable primary roll collapses into three-dimensional turbulence. The results show that convectively overturning regions are always formed but the details of breaking and the resulting structure of the mixed layer depend on the effective Reynolds number of the flow. With sufficient viscous damping, three-dimensional turbulent convective instabilities are more easily suppressed than two-dimensional laminar overturning.     

Time-dependent liquid metal flows with free convection and a deformable free surface

The finite element method is employed to investigate time-dependent liquid metal flows with free convection, free surfaces and Marangoni effects. The liquid circulates in a two-dimensional shallow trough with differentially heated vertical walls. The spatial formulation incorporates mixed Lagrangian approximations to the velocity, pressure, temperature and free surface position. The time integration is performed with the backward Euler and trapezoid rule methods with step size control. The Galerkin method is used to reduce the problem to a set of non-linear equations which are solved with the Newton–Raphson method. Calculations are performed for conditions relevant to the electron beam vaporization of refractory metals. The Prandtl number is 0·015 and Grashof number are in the transition range between laminar and turbulent flow. The results reveal the effects of flow intensity, surface tension gradients, mesh refinement and time integration strategy.     

Time-dependent behaviour of the liquid film in horizontal annular flow

A crucial point still to be established in the prediction of the film thickness distribution in horizontal annular two-phase flow is the mechanism(s) for transporting liquid from the bottom to the top part of the tube. To resolve this issue, the time-dependent behaviour of the liquid film is studied. Wave characteristics such as velocity and frequency are measured around the circumference. It is inferred from the autospectral density functions of film thickness variation that disturbance waves play an important, but as yet unclear, role in the formation of a liquid film in the top part of the tube. A new mechanism, based on the shape of disturbance waves is proposed.     

Measurement of liquid film thickness in micro square channel

In micro channels, slug flow becomes one of the main flow regimes due to strong surface tension. In micro channel slug flow, elongated bubble flows with the thin liquid film confined between the bubble and the channel wall. Liquid film thickness is an important parameter in many applications, e.g., micro heat exchanger, micro reactor, coating process etc. In the present study, liquid film thickness in micro square channels is measured locally and instantaneously with the confocal method. Square channels with hydraulic diameter of D_h = 0.3, 0.5 and 1.0 mm are used. In order to investigate the effect of inertial force on the liquid film thickness, three working fluids, ethanol, water and FC-40 are used. At small capillary numbers, liquid film at the channel center becomes very thin and the bubble interface is not axisymmetric. However, as capillary number increases, bubble interface becomes axisymmetric. Transition from non-axisymmetric to axisymmetric flow pattern starts from lower capillary number as Reynolds number increases. An empirical correlation for predicting axisymmetric bubble radius based on capillary number and Weber number is proposed from the present experimental data.     

A direct boundary-layer stability analysis of steady-state cavity convection flow

Natural convection flow in cavities with insulated top and bottom and heated and cooled walls is known to exhibit travelling wave instabilities in the thermal boundary layers that form on the walls. In water (Pr = 7.5) at Rayleigh number Ra = 6 × 10⁸, these waves have been observed at start-up. However no such waves have been observed for the fully developed flow, although it may be assumed that the stability character of the boundary layers is at least approximately the same. The start-up waves are generated by perturbations to the system. In the present paper, an artificial perturbation is applied to the system to determine the stability character of the boundary layers in fully developed flow. It is shown that the thermal boundary layers in the fully developed flow have approximately the same stability character as the start-up flow.     

Effects of temperature-dependent viscosity variation on entropy generation,heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford [12], is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Peclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case by Ratts and Raut [14].     

Low-Reynolds-number motion of a deformable drop between two parallel plane walls

The motion of a three-dimensional deformable drop between two parallel plane walls in a low-Reynolds-number Poiseuille flow is examined using a boundary-integral algorithm that employs the Green’s function for the domain between two infinite plane walls, which incorporates the wall effects without discretization of the walls. We have developed an economical calculation scheme that allows long-time dynamical simulations, so that both transient and steady-state shapes and velocities are obtained. Results are presented for neutrally buoyant drops having various viscosity, size, deformability, and channel position. For nearly spherical drops, the decrease in translational velocity relative to the undisturbed fluid velocity at the drop center increases with drop size, proximity of the drop to one or both walls, and drop-to-medium viscosity ratio. When deformable drops are initially placed off the centerline of flow, lateral migration towards the channel center is observed, where the drops obtain steady shapes and translational velocities for subcritical capillary numbers. With increasing capillary number, the drops become more deformed and have larger steady velocities due to larger drop-to-wall clearances. Non-monotonic behavior for the lateral migration velocities with increasing viscosity ratio is observed. Simulation results for large drops with non-deformed spherical diameters exceeding the channel height are also presented.     

Cusped ripples at the plane surface of a viscous liquid

Summary We study the time-evolution of periodical ripples of a viscous liquid at the plane free surface under the action of a distant pure straining flow. We neglect inertial forces (Stokes flow) and include surface tension effects. The solutions for a contracting surface and constant strain rate show that the ripples may develop near-cusps during a stage of the evolution, though later the free surface inevitably asymptotically tends to a smooth plane with vanishing ripples due to the action of capillarity. We obtain the condition for cusp formation in this intermediate stage in terms of the initial capillary number and aspect ratio. If the capillary number is kept constant, the surface tends to shrink through a succession of self-similar trochoidal shapes, whose aspect ratio is given by the capillary number. Received 23 March 1998, accepted for publication 23 July 1998     

Effects of the side walls on unsteady flow of a second grade fluid over a plane wall

The effects of the side walls on unsteady flow of a second grade fluid over a plan wall are considered. The solution of the governing equation for velocity is obtained by the sine transform method. This gives a correct result for the shear stress at the bottom wall. The shear stress at the bottom wall is minimum at the middle of the plate and it increases near the side walls. It is shown that the mean thickness of the layer of the liquid over the plate increases with time and the ratio of the mean thickness to the distance between the side walls becomes ultimately 0.2714.