-
Evaporation-driven coalescence of two droplets undergoing freezing
Authors:
Sivanandan Kavuri,
George Karapetsas,
Chander Shekhar Sharma,
Kirti Chandra Sahu
Abstract:
We examine the evaporation-induced coalescence of two droplets undergoing freezing by conducting numerical simulations employing the lubrication approximation. When two sessile drops undergo freezing in close vicinity over a substrate, they interact with each other through the gaseous phase and the simultaneous presence of evaporation/condensation. In an unsaturated environment, the evaporation fl…
▽ More
We examine the evaporation-induced coalescence of two droplets undergoing freezing by conducting numerical simulations employing the lubrication approximation. When two sessile drops undergo freezing in close vicinity over a substrate, they interact with each other through the gaseous phase and the simultaneous presence of evaporation/condensation. In an unsaturated environment, the evaporation flux over the two volatile sessile drops is asymmetric, with lower evaporation in the region between the two drops. This asymmetry in the evaporation flux generates an asymmetric curvature in each drop, which results in a capillary flow that drives the drops closer to each other, eventually leading to their coalescence. This capillary flow, driven by evaporation, competes with the upward movement of the freezing front, depending on the relative humidity in the surrounding environment. We found that higher relative humidity reduces the evaporative flux, delaying capillary flow and impeding coalescence by restricting contact line motion. For a constant relative humidity, the substrate temperature governs the coalescence phenomenon, and resulting condensation can accelerate this process. Interestingly, lower substrate temperatures are observed to facilitate faster propagation of the freezing front, which, in turn, restricts coalescence.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
Flow Structure near Three Phase Contact Line of Low-Contact-Angle Evaporating Droplets
Authors:
Zhenying Wang,
George Karapetsas,
Prashant Valluri,
Chihiro Inoue
Abstract:
Flow structure near three phase contact line (TPCL) of evaporating liquids plays a significant role in liquid wetting and dewetting, liquid film evaporation and boiling, etc. Despite the wide focus it receives, the interacting mechanisms therein remain elusive and in specific cases, controversial. Here, we reveal the profile of internal flow and elucidate the dominating mechanisms near TPCL of eva…
▽ More
Flow structure near three phase contact line (TPCL) of evaporating liquids plays a significant role in liquid wetting and dewetting, liquid film evaporation and boiling, etc. Despite the wide focus it receives, the interacting mechanisms therein remain elusive and in specific cases, controversial. Here, we reveal the profile of internal flow and elucidate the dominating mechanisms near TPCL of evaporating droplets, using mathematical modelling, microPIV, and infrared thermography. We indicate that for less volatile liquids such as butanol, the flow pattern is dominated by capillary flow. With increasing liquid volatility, e.g., alcohol, the effect of evaporation cooling, under conditions, induces interfacial temperature gradient with cold droplet apex and warm edge. The temperature gradient leads to Marangoni flow that competes with outwarding capillary flow, resulting in the reversal of interfacial flow and the formation of a stagnation point near TPCL. The spatiotemporal variations of capillary velocity and Marangoni velocity are further quantified by mathematically decomposing the tangential velocity of interfacial flow. The conclusions can serve as a theoretical base for explaining deposition patterns from colloidal suspensions, and can be utilized as a benchmark in analyzing more complex liquid systems.
△ Less
Submitted 18 January, 2024;
originally announced January 2024.
-
Freezing of sessile droplet and frost halo formation
Authors:
Sivanandan Kavuri,
George Karapetsas,
Chander Shekhar Sharma,
Kirti Chandra Sahu
Abstract:
The freezing of a sessile droplet unveils fascinating physics, characterised by the emergence of a frost halo on the underlying substrate, the progression of the liquid-ice interface, and the formation of a cusp-like morphology at the tip of the droplet. We investigate the freezing of a volatile sessile droplet, focusing on the frost halo formation, which has not been theoretically explored. The f…
▽ More
The freezing of a sessile droplet unveils fascinating physics, characterised by the emergence of a frost halo on the underlying substrate, the progression of the liquid-ice interface, and the formation of a cusp-like morphology at the tip of the droplet. We investigate the freezing of a volatile sessile droplet, focusing on the frost halo formation, which has not been theoretically explored. The formation of the frost halo is associated with the inherent evaporation process in the early freezing stages. We observe a negative evaporation flux enveloping the droplet in the initial stages, which indicates that vapour produced during freezing condenses on the substrate close to the contact line, forming a frost halo. The condensate accumulation triggers re-evaporation, resulting in a temporal shift of the frost halo region away from the contact line. Eventually, it disappears due to the diffusive nature of the water vapour far away from the droplet. We found that increasing the relative humidity increases the lifetime of the frost halo due to a substantial reduction in evaporation that prolonged the presence of net condensate on the substrate. Increasing liquid volatility increases the evaporation flux, and condensation occurs closer to the droplet, as a higher amount of vapour is in the periphery of the droplet. We also found that decreasing the thermal conductivity of the substrate increases the total freezing time. The slower freezing process is accompanied by increased vaporized liquid, resulting in condensation with its concentration reaching supersaturation.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
Intricate Role of Thermal Properties and Volatility in Droplet Spreading: A Generalization to Tanner's Law
Authors:
Zhenying Wang,
George Karapetsas,
Prashant Valluri,
Chihiro Inoue
Abstract:
Droplet spreading is ubiquitous and plays a significant role in liquid-based energy systems, thermal management devices, and microfluidics. While the spreading of non-volatile droplets is quantitatively understood, the spreading and flow transition in volatile droplets remains elusive due to the complexity added by interfacial phase change and non-equilibrium thermal transport. Here we show, using…
▽ More
Droplet spreading is ubiquitous and plays a significant role in liquid-based energy systems, thermal management devices, and microfluidics. While the spreading of non-volatile droplets is quantitatively understood, the spreading and flow transition in volatile droplets remains elusive due to the complexity added by interfacial phase change and non-equilibrium thermal transport. Here we show, using both mathematical modeling and experiments, that the wetting dynamics of volatile droplets can be scaled by the spatial-temporal interplay between capillary, evaporation, and thermal Marangoni effects. We elucidate and quantify these complex interactions using phase diagrams based on systematic theoretical and experimental investigations. A spreading law of evaporative droplets is derived by generalizing Tanner's law (valid for non-volatile liquids) to a full range of liquids with saturation vapor pressure spanning from 10^1 to 10^4 Pa and on substrates with thermal conductivity from 10^{-1} to 10^3 W/m/K. Besides its importance in fluid-based industries, the conclusions also enable a unifying explanation to a series of individual works including the criterion of flow reversal and the state of dynamic wetting, making it possible to control liquid transport in diverse application scenarios.
△ Less
Submitted 23 August, 2023; v1 submitted 8 April, 2023;
originally announced April 2023.
-
Dynamics of the interaction of multiple evaporating droplets on compliant substrates
Authors:
Anna Malachtari,
George Karapetsas
Abstract:
The dynamics of the interaction of a system of two thin volatile liquid droplets resting on a soft viscoelastic solid substrate are investigated theoretically. The developed model fully considers the effect of evaporative cooling and the generated Marangoni stresses due to the induced thermal gradients, while also accounting for the effect of the gas phase composition and the diffusion of vapour i…
▽ More
The dynamics of the interaction of a system of two thin volatile liquid droplets resting on a soft viscoelastic solid substrate are investigated theoretically. The developed model fully considers the effect of evaporative cooling and the generated Marangoni stresses due to the induced thermal gradients, while also accounting for the effect of the gas phase composition and the diffusion of vapour in the atmosphere of the droplets. Using the framework of lubrication theory, we derive evolution equations for both the droplet profile and the displacement of the elastic solid, which are solved in combination with Laplace's equation for the vapour concentration in the gas phase. A disjoining-pressure/precursor-film approach is used to describe contact-line motion. The evolution equations are solved numerically, using the finite-element method, and we present a thorough parametric analysis to investigate the physical properties and mechanisms that affect the dynamics of droplet interactions. The results show that the droplets interact through both the soft substrate and the gas phase. In the absence of thermocapillary phenomena, the combined effect of non-uniform evaporation due to the increased vapour concentration between the two droplets and elastocapillary phenonema determines whether the drop-drop interaction is attractive or repulsive. The Marangoni stresses suppress droplet attraction at the early stages of the drying process and lead to longer droplet lifetimes. For substrates with intermediate stiffness, the emergence of spontaneous symmetry breaking at late stages of evaporation is found. The rich dynamics of this complex system is explored by constructing a detailed map of the dynamic regimes.
△ Less
Submitted 12 October, 2023; v1 submitted 1 February, 2023;
originally announced February 2023.
-
Physics-agnostic and Physics-infused machine learning for thin films flows: modeling, and predictions from small data
Authors:
Cristina P. Martin-Linares,
Yorgos M. Psarellis,
Georgios Karapetsas,
Eleni D. Koronaki,
Ioannis G. Kevrekidis
Abstract:
Numerical simulations of multiphase flows are crucial in numerous engineering applications, but are often limited by the computationally demanding solution of the Navier-Stokes (NS) equations. Here, we present a data-driven workflow where a handful of detailed NS simulation data are leveraged into a reduced-order model for a prototypical vertically falling liquid film. We develop a physics-agnosti…
▽ More
Numerical simulations of multiphase flows are crucial in numerous engineering applications, but are often limited by the computationally demanding solution of the Navier-Stokes (NS) equations. Here, we present a data-driven workflow where a handful of detailed NS simulation data are leveraged into a reduced-order model for a prototypical vertically falling liquid film. We develop a physics-agnostic model for the film thickness, achieving a far better agreement with the NS solutions than the asymptotic Kuramoto-Sivashinsky (KS) equation. We also develop two variants of physics-infused models providing a form of calibration of a low-fidelity model (i.e. the KS) against a few high-fidelity NS data. Finally, predictive models for missing data are developed, for either the amplitude, or the full-field velocity and even the flow parameter from partial information. This is achieved with the so-called "Gappy Diffusion Maps", which we compare favorably to its linear counterpart, Gappy POD.
△ Less
Submitted 29 January, 2023;
originally announced January 2023.
-
Open-source finite volume solvers for multiphase (n-phase) flows involving either Newtonian or non-Newtonian complex fluids
Authors:
Chris Dritselis,
George Karapetsas
Abstract:
Two new control volume solvers multiFluidInterFoam and rheoMultiFluidInterFoam are presented for the simulation of Newtonian and non-Newtonian n-phase flows, respectively, fully accounting for interfacial tension and contact-angle effects for each phase. The multiFluidInterFoam solver modifies certain crucial aspects of the regular multiphaseInterFoam solver provided by OpenFOAM for Newtonian flow…
▽ More
Two new control volume solvers multiFluidInterFoam and rheoMultiFluidInterFoam are presented for the simulation of Newtonian and non-Newtonian n-phase flows, respectively, fully accounting for interfacial tension and contact-angle effects for each phase. The multiFluidInterFoam solver modifies certain crucial aspects of the regular multiphaseInterFoam solver provided by OpenFOAM for Newtonian flows, improving its efficiency and robustness, but most importantly improving considerably its accuracy for surface tension driven flows. The new solver uses the volume-of-fluid (VOF) method based on the MULES method and artificial interface compression for the interface capturing, in combination with a suitable continuum surface force model; i.e. the interfacial tension coefficient decomposition method is employed to treat pairings of interfacial tension between the phases. VOF smoothing is also implemented by applying a Laplacian filter on the distribution of volume fractions. The rheoMultiFluidInterFoam solver is capable of fully taking into account complex non-Newtonian effects (e.g. yield stress, viscoelasticity, etc.) of the involved liquids. To this end, the developed solver incorporates the RheoTool toolbox (Pimenta and Alves (2017) J. Non-Newtonian Fluid Mech. 239, 85-104), utilizing a wealth of constitutive equations suitable for modelling different types of fluids with complex rheology. Our solvers are tested for a wide range of different flow setups, considering typical benchmark two-phase and three-phase flows; flows involving Newtonian, viscoplastic and viscoelastic liquids have been considered. Comparisons against analytical solutions or previously published numerical data are performed clearly demonstrating the capability of the new solvers to efficiently and accurately simulate interfacial tension dominated three-phase flows for both simple and complex liquids.
△ Less
Submitted 18 March, 2022;
originally announced March 2022.
-
Spreading and retraction dynamics of sessileevaporating droplets comprising volatile binary mixtures
Authors:
A. G. L. Williams. G. Karapetsas,
D. Mamalis,
K. Sefiane,
O. K. Matar,
P. Valluri
Abstract:
The dynamics of thin volatile droplets comprising of binary mixtures deposited on a heated substrate are investigated. Using lubrication theory, we develop a novel one-sided model to predict the spreading and retraction of an evaporating sessile axisymmetric droplet formed of a volatile binary mixture on a substrate with high wettability. A thin droplet with a moving contact line is considered, ta…
▽ More
The dynamics of thin volatile droplets comprising of binary mixtures deposited on a heated substrate are investigated. Using lubrication theory, we develop a novel one-sided model to predict the spreading and retraction of an evaporating sessile axisymmetric droplet formed of a volatile binary mixture on a substrate with high wettability. A thin droplet with a moving contact line is considered, taking into account the variation of liquid properties with concentration as well as the effects of inertia. The parameter space is explored and the resultant effects on wetting and evaporation are evaluated. Increasing solutal Marangoni stress enhances spreading rates in all cases, approaching those of superspreading liquids. To validate our model, experiments are conducted with binary ethanol-water droplets spreading on hydrophilic glass slides heated from below. The spreading rate is quantified, revealing that preferential evaporation of the more volatile component (ethanol) at the contact line drives superspreading, leading in some cases to a contact line instability. Good qualitative agreement is found between our model and experiments, with quantitative agreement being achieved in terms of spreading rate.
△ Less
Submitted 22 September, 2020;
originally announced September 2020.
-
The stability of slowly evaporating thin liquid films of binary mixtures
Authors:
R. K. Nazareth,
G. Karapetsas,
K. Sefiane,
O. Matar,
P. Valluri
Abstract:
We consider the evaporation of a thin liquid layer which consists of a binary mixture of volatile liquids. The mixture is on top of a heated substrate and in contact with the gas phase that consists of the same vapour as the binary mixture. The effect of thermocapillarity, solutocapillarity and the van der Waals interactions are considered. We derive the long-wave evolution equations for the free…
▽ More
We consider the evaporation of a thin liquid layer which consists of a binary mixture of volatile liquids. The mixture is on top of a heated substrate and in contact with the gas phase that consists of the same vapour as the binary mixture. The effect of thermocapillarity, solutocapillarity and the van der Waals interactions are considered. We derive the long-wave evolution equations for the free interface and the volume fraction that govern the two-dimensional stability of the layer subject to the above coupled mechanisms and perform a linear stability analysis. Our results demonstrate two modes of instabilities, a monotonic instability mode and an oscillatory instability mode. We supplement our results from stability analysis with transient simulations to examine the dynamics in the nonlinear regime and analyse how these instabilities evolve with time. More precisely we discuss how the effect of relative volatility along with the competition between thermal and solutal Marangoni effect defines the mode of instability that develops during the evaporation of the liquid layer due to preferential evaporation of one of the components.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Effect of substrate topography, material wettability and dielectric thickness on reversible electrowetting
Authors:
Nikolaos T. Chamakos,
George Karapetsas,
Athanasios G. Papathanasiou
Abstract:
Recent experiments by Kavousanakis et al., Langmuir, 2018 [1], showed that reversible electrowetting on superhydrophobic surfaces can be achieved by using a thick solid dielectric layer (e.g. tens of micrometers). It has also been shown, through equilibrium (static) computations, that when the dielectric layer is thick enough the electrostatic pressure is smoothly distributed along the droplet sur…
▽ More
Recent experiments by Kavousanakis et al., Langmuir, 2018 [1], showed that reversible electrowetting on superhydrophobic surfaces can be achieved by using a thick solid dielectric layer (e.g. tens of micrometers). It has also been shown, through equilibrium (static) computations, that when the dielectric layer is thick enough the electrostatic pressure is smoothly distributed along the droplet surface, thus the irreversible Cassie to Wenzel wetting transitions can be prevented. In the present work we perform more realistic, dynamic simulations of the electrostatically-induced spreading on superhydrophobic surfaces. To this end, we employ an efficient numerical scheme which enables us to fully take into account the topography of the solid substrate. We investigate in detail the role of the various characteristics of the substrate (i.e. the dielectric thickness, geometry and material wettability) and present relevant flow maps for the resulting wetting states. Through our dynamic simulations, we identify the conditions under which it is possible to achieve reversible electrowetting. We have found that not only the collapse (Cassie-Baxter to Wenzel) transitions but also the contact angle hysteresis of the substrate significantly affects the reversibility.
△ Less
Submitted 27 July, 2018;
originally announced July 2018.
