Hostname: page-component-78c5997874-j824f Total loading time: 0 Render date: 2024-11-14T22:42:00.663Z Has data issue: false hasContentIssue false

Interfacial dynamics of a confined liquid–vapour bilayer undergoing evaporation

Published online by Cambridge University Press:  15 October 2018

Dipin S. Pillai*
Affiliation:
Department of Chemical Engineering, University of Florida, Gainesville, FL 32611, USA
R. Narayanan
Affiliation:
Department of Chemical Engineering, University of Florida, Gainesville, FL 32611, USA
*
Email address for correspondence: dipinsp@ufl.edu

Abstract

The dynamics of an interface between a thin liquid–vapour bilayer undergoing evaporation is studied. Both phases are considered to be hydrodynamically and thermally active, with momentum and thermal inertia taken into account. A reduced-order model based on the weighted-residual integral boundary layer method is used to investigate the dynamical behaviour for two cases, viz., phase change in the absence of gravity and then phase change in the presence of gravity. In the first case, it is shown that evaporative instability may cause rupture of either liquid or vapour layer depending on system parameters. Close to interfacial rupture, the disjoining pressure due to intermolecular forces results in the formation of drops (bubbles) separated by a thin film for low liquid (vapour) hold-up. Momentum inertia is shown to have a stabilizing effect, while thermal inertia has a destabilizing effect. In the second case, evaporative suppression of Rayleigh–Taylor (R–T) instability shows emergence of up to two neutral wavenumbers. Weak nonlinear analysis of these neutral wavenumbers suggests that the instability may be either supercritical or subcritical depending on the rate of evaporation. At high rates of evaporation, both neutral wavenumbers are supercritical and computations on the interface evolution lead to nonlinear saturated steady states. Momentum inertia slows down the rate of interface deformation and results in an oscillatory approach to saturation. Thermal inertia results in larger interface deformation and the saturated steady state is shifted closer to the wall. At very low evaporation rates, only one neutral wavenumber of subcritical nature exists. The nonlinear evolution of the interface in this case is then similar to pure R–T instability, exhibiting spontaneous lateral sliding as it approaches the wall.

Type
JFM Papers
Copyright
© 2018 Cambridge University Press 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Adham-Khodaparast, K., Kawaji, M. & Antar, B. N. 1995 The Rayleigh–Taylor and Kelvin–Helmholtz stability of a viscous liquid–vapor. Phys. Fluids 7, 359364.Google Scholar
Ajaev, V. S. 2005 Evolution of dry patches in evaporating liquid films. Phys. Rev. E 72, 031605.Google Scholar
Alexeev, A. & Oron, A. 2007 Suppression of the Rayleigh–Taylor instability of thin liquid films by the Marangoni effect. Phys. Fluids 19, 082101.Google Scholar
Berg, J. C., Boudart, M. & Acrivos, A. 1966 Natural convection in pools of evaporating liquids. J. Fluid Mech. 24, 721735.Google Scholar
Bestehorn, M. & Merkt, D. 2006 Regular surface patterns on Rayleigh–Taylor unstable evaporating films heated from below. Phys. Rev. Lett. 97, 127802.Google Scholar
Boos, W. & Thess, A. 1999 Cascade of structures in long-wavelength Marangoni instability. Phys. Fluids 11, 14841494.Google Scholar
Burelbach, J. P., Bankoff, S. G. & Davis, S. H. 1988 Nonlinear stability of evaporating/condensing liquid films. J. Fluid Mech. 195, 463494.Google Scholar
Chatterjee, A., Plawsky, J. L. & Wayner, P. C. Jr 2013 Constrained vapor bubble heat pipe experiment aboard the international space station. J. Thermophys. Heat Transfer 27, 309319.Google Scholar
Dietze, G., Picardo, J. R. & Narayanan, R. 2018 Sliding instability of draining fluid films. J. Fluid Mech. (in press).Google Scholar
Dietze, G. F. & Ruyer-Quil, C. 2013 Wavy liquid films in interaction with a confined laminar gas flow. J. Fluid Mech. 722, 348393.Google Scholar
Dietze, G. F. & Ruyer-Quil, C. 2015 Films in narrow tubes. J. Fluid Mech. 762, 68109.Google Scholar
Guo, W., Labrosse, G. & Narayanan, R. 2012 The Application of the Chebyshev-Spectral Method in Transport Phenomena. Springer.Google Scholar
Guo, W. & Narayanan, R. 2010 Interfacial instability due to evaporation and convection: linear and nonlinear analyses. J. Fluid Mech. 650, 363389.Google Scholar
Ho, S. 1980 Linear Rayleigh–Taylor stability of viscous fluids with mass and heat transfer. J. Fluid Mech. 101, 111127.Google Scholar
Hsieh, D. Y. 1972 Effects of heat and mass transfer on Rayleigh–Taylor instability. J. Basic Engng 94, 156160.Google Scholar
Hsieh, D. Y. 1978 Interfacial stability with mass and heat transfer. Phys. Fluids 21, 745748.Google Scholar
Huang, A. & Joseph, D. D. 1992 Instability of the equilibrium of a liquid below its vapor between horizontal heated plates. J. Fluid Mech. 242, 235247.Google Scholar
Israelachvili, J. N. 2011 Intermolecular and Surface Forces. Academic Press.Google Scholar
Johns, L. E. & Narayanan, R. 2002 Interfacial Instabilities. Springer.Google Scholar
Kalliadasis, S., Demekhin, E. A., Ruyer-Quil, C. & Velarde, M. G. 2003a Thermocapillary instability and wave formation on a film falling down a uniformly heated plane. J. Fluid Mech. 492, 303338.Google Scholar
Kalliadasis, S., Kiyashko, A. & Demekhin, E. A. 2003b Marangoni instability of a thin liquid film heated from below by a local heat source. J. Fluid Mech. 475, 377408.Google Scholar
Kalliadasis, S., Ruyer-Quil, C., Scheid, B. & Velarde, M. 2011 Falling Liquid Films. vol. 176. Springer Science & Business Media.Google Scholar
Kanatani, K. 2010 Interfacial instability induced by lateral vapor pressure fluctuation in bounded thin liquid–vapor layers. Phys. Fluids 22, 012101.Google Scholar
Kanatani, K. & Oron, A. 2011 Nonlinear dynamics of confined thin liquid–vapor bilayer systems with phase change. Phys. Fluids 23, 032102.Google Scholar
Kharangate, C. R., Oneill, L. E. & Mudawar, I. 2016 Effects of two-phase inlet quality, mass velocity, flow orientation, and heating perimeter on flow boiling in a rectangular channel. Part 1. Two-phase flow and heat transfer results. Intl J. Heat Mass Transfer 103, 12611279.Google Scholar
Kim, B. J. & Kim, K. D. 2016 Rayleigh–Taylor instability of viscous fluids with phase change. Phys. Rev. E 93, 043123.Google Scholar
Konovalov, V. V., Lyubimov, D. V. & Lyubimova, T. P. 2016 The Rayleigh–Taylor instability of the externally cooled liquid lying over a thin vapor film coating the wall of a horizontal plane heater. Phys. Fluids 28, 064102.Google Scholar
Konovalov, V. V., Lyubimov, D. V. & Lyubimova, T. P. 2017 Influence of phase transition on the instability of a liquid–vapor interface in a gravitational field. Phys. Rev. Fluids 2, 063902.Google Scholar
Labrosse, G. 2011 Méthodes numriques: Méthodes spectrale: Méthodes locales globales, méthodes globales, problèmes d’Helmotz et de Stokes, équations de Navier–Stokes. Ellipses Marketing.Google Scholar
Lin, Z. 2012 Evaporative Self-Assembly of Ordered Complex Structures. World Scientific.Google Scholar
Lister, J. R., Morrison, N. F. & Rallison, J. M. 2006 Sedimentation of a two-dimensional drop towards a rigid horizontal plate. J. Fluid Mech. 552, 345351.Google Scholar
Lister, J. R., Rallison, J. M. & Rees, S. J. 2010 The nonlinear dynamics of pendent drops on a thin film coating the underside of a ceiling. J. Fluid Mech. 647, 239264.Google Scholar
Margerit, J., Colinet, P., Lebon, G., Iorio, C. S. & Legros, J. C. 2003 Interfacial nonequilibrium and Benard–Marangoni instability of a liquid–vapor system. Phys. Rev. E 68, 041601.Google Scholar
McFadden, G. B. & Coriell, S. R. 2009 Onset of oscillatory convection in two liquid layers with phase change. Phys. Fluids 21, 034101.Google Scholar
McFadden, G. B., Coriell, S. R., Gurski, K. F. & Cotrell, D. L. 2007 Onset of convection in two liquid layers with phase change. Phys. Fluids 19, 104109.Google Scholar
Miller, C. A. 1973 Stability of moving surfaces in fluid systems with heat and mass transport II. Combined effects of transport and density difference between phases. AIChE J. 19, 909915.Google Scholar
Narendranath, A. D., Hermanson, J. C., Kolkka, R. W., Struthers, A. A. & Allen, J. S. 2014 The effect of gravity on the stability of an evaporating liquid film. Microgravity Sci. Technol. 26, 189199.Google Scholar
Oron, A., Davis, S. H. & Bankoff, S. G. 1997 Long-scale evolution of thin liquid films. Rev. Mod. Phys. 69, 931980.Google Scholar
Ozen, O. & Narayanan, R. 2004 The physics of evaporative and convective instabilities in bilayer systems: linear theory. Phys. Fluids 16, 46444652.Google Scholar
Ozen, O. & Narayanan, R. 2006 A note on the Rayleigh–Taylor instability with phase change. Phys. Fluids 18, 042110.Google Scholar
Palmer, H. J. 1976 The hydrodynamic instability of rapidly evaporating liquids at reduced pressure. J. Fluid Mech. 75, 487511.Google Scholar
Persad, A. H. & Ward, C. 2016 Expressions for the evaporation and condensation coefficients in the Hertz–Knudsen relation. Chem. Rev. 116, 77277767.Google Scholar
Rajabi, A. A. A. & Winterton, R. H. S. 1987 Heat transfer across vapour film without ebullition. Intl J. Heat Mass Transfer 30, 17031708.Google Scholar
Ruyer-Quil, C. & Manneville, P. 2000 Improved modeling of flows down inclined planes. Eur. Phys. J. B 15, 357369.Google Scholar
Ruyer-Quil, C. & Manneville, P. 2002 Further accuracy and convergence results on the modeling of flows down inclined planes by weighted-residual approximations. Phys. Fluids 14, 170183.Google Scholar
Shankar, P. N. & Deshpande, M. D. 1990 On the temperature distribution in liquid–vapor phase change between plane liquid surfaces. Phys. Fluids 2, 10301038.Google Scholar
Smith, M. K. & Vrane, D. R. 1999 Deformation and rupture in confined, thin liquid films driven by thermocapillarity. In Fluid Dynamics at Interfaces (ed. Shyy, W. & Narayanan, R.), pp. 221233. Cambridge University Press.Google Scholar
Theofanous, T. G., Tu, J. P., Dinh, A. T. & Dinh, T. N. 2002 The boiling crisis phenomenon. Part I. Nucleation and nucleate boiling heat transfer. Exp. Therm. Fluid Sci. 26, 775792.Google Scholar
Trevelyan, P. M. J. & Kalliadasis, S. 2004 Wave dynamics on a thin-film falling down a heated wall. J. Engng Maths 50, 177208.Google Scholar
Ward, C. A. & Stanga, D. 2001 Interfacial conditions during evaporation or condensation of water. Phys. Rev. E 64, 051509.Google Scholar
Zhang, N. 2006 Surface tension-driven convection flow in evaporating liquid layers. In Surface Tension-Driven Flows and Applications (ed. Savino, R.), pp. 147170. Research Signpost.Google Scholar

Pillai and Narayanan supplementary movie 1

Spatio-temporal evolution of the interface profile for evaporative instability when heated from the liquid side showing vapour rupture

Download Pillai and Narayanan supplementary movie 1(Video)
Video 822.1 KB
Supplementary material: File

Pillai and Narayanan supplementary material

Supplementary material

Download Pillai and Narayanan supplementary material(File)
File 15.8 KB

Pillai and Narayanan supplementary movie 2

Spatio-temporal evolution of the interface profile for evaporative instability when heated from the liquid side showing liquid rupture

Download Pillai and Narayanan supplementary movie 2(Video)
Video 909.9 KB

Pillai and Narayanan supplementary movie 3

Spatio-temporal evolution of the interface profile for evaporative instability when heated from the liquid side showing liquid rupture

Download Pillai and Narayanan supplementary movie 3(Video)
Video 1.1 MB

Pillai and Narayanan supplementary movie 4

Nonlinear saturation of the interface profile for a Rayleigh-Taylor unstable system heated from the vapour side; H=0.3, E=6.1×10-5, k=0.2

Download Pillai and Narayanan supplementary movie 4(Video)
Video 770.1 KB

Pillai and Narayanan supplementary movie 5

Nonlinear saturation of the interface profile for a Rayleigh-Taylor unstable system heated from the vapour side; H=0.3, E=1.22×10-5, k=0.2.

Download Pillai and Narayanan supplementary movie 5(Video)
Video 748.7 KB

Pillai and Narayanan supplementary movie 6

Velocity profile in each phase for the R-T unstable configuration, exhibiting a flow reversal in the liquid phase as the interface attains its steady state

Download Pillai and Narayanan supplementary movie 6(Video)
Video 1.3 MB

Pillai and Narayanan supplementary movie 7

Steady interface profile close to the left neutral wavenumber (k = 1.01 kcL); kcL = 0.081

Download Pillai and Narayanan supplementary movie 7(Video)
Video 1.8 MB

Pillai and Narayanan supplementary movie 8

Oscillatory approach to saturation of the interface in the presence of momentum inertia in evaporative suppresion of R-T instability

Download Pillai and Narayanan supplementary movie 8(Video)
Video 689.3 KB
Supplementary material: File

Pillai and Narayanan supplementary material

Supplementary data

Download Pillai and Narayanan supplementary material(File)
File 15.8 KB