-
Evaporating sessile droplets: solutal Marangoni effects overwhelm thermal Marangoni flow
Authors:
Duarte Rocha,
Philip L. Lederer,
Pim J. Dekker,
Alvaro Marin,
Detlef Lohse,
Christian Diddens
Abstract:
When an evaporating water droplet is deposited on a thermally conductive substrate, the minimum temperature will be at the apex due to evaporative cooling. Consequently, density and surface tension gradients emerge within the droplet and at the droplet-gas interface, giving rise to competing flows from, respectively, the apex towards the contact line (thermal-buoyancy-driven flow) and the other wa…
▽ More
When an evaporating water droplet is deposited on a thermally conductive substrate, the minimum temperature will be at the apex due to evaporative cooling. Consequently, density and surface tension gradients emerge within the droplet and at the droplet-gas interface, giving rise to competing flows from, respectively, the apex towards the contact line (thermal-buoyancy-driven flow) and the other way around (thermal Marangoni flow). In small droplets with a diameter below the capillary length, the thermal Marangoni effects are expected to dominate over thermal buoyancy ("thermal Rayleigh") effects. However, contrary to these theoretical predictions, our experiments mostly show a dominant circulation from the apex towards the contact line, indicating a prevailing of thermal Rayleigh convection. Furthermore, our experiments often show an unexpected asymmetric flow that persisted for several minutes. We hypothesise that a tiny amount of contaminants, commonly encountered in experiments with water/air interfaces, act as surfactants and counteract the thermal surface tension gradients at the interface and thereby promote the dominance of Rayleigh convection. Our finite element numerical simulations demonstrate that, under our specified experimental conditions, a mere 0.5% reduction in the static surface tension caused by surfactants leads to a reversal in the flow direction, compared to the theoretical prediction without contaminants. Additionally, we investigate the linear stability of the axisymmetric solutions, revealing that the presence of surfactants also affects the axial symmetry of the flow.
△ Less
Submitted 22 October, 2024;
originally announced October 2024.
-
On the Consistency of Dynamic Wetting Boundary Conditions for the Navier-Stokes-Cahn-Hilliard Equations
Authors:
T. H. B. Demont,
S. K. F. Stoter,
C. Diddens,
E. H. van Brummelen
Abstract:
We investigate the limiting behavior of the Navier-Stokes-Cahn-Hilliard model for binary-fluid flows as the diffuse-interface thickness passes to zero, in the presence of fluid-fluid-solid contact lines. Allowing for motion of such contact lines relative to the solid substrate is required to adequately model multi-phase and multi-species fluid transport past and through solid media. Even though di…
▽ More
We investigate the limiting behavior of the Navier-Stokes-Cahn-Hilliard model for binary-fluid flows as the diffuse-interface thickness passes to zero, in the presence of fluid-fluid-solid contact lines. Allowing for motion of such contact lines relative to the solid substrate is required to adequately model multi-phase and multi-species fluid transport past and through solid media. Even though diffuse-interface models provide an inherent slip mechanism through the mobility-induced diffusion, this slip vanishes as the interface thickness and mobility parameter tend to zero in the so-called sharp-interface limit. The objective of this work is to present dynamic wetting and generalized Navier boundary conditions for diffuse-interface models that are consistent in the sharp-interface limit. We concentrate our analysis on the prototypical binary-fluid Couette-flow problems. To verify the consistency of the diffuse-interface model in the limit of vanishing interface thickness, we provide reference limit solutions of a corresponding sharp-interface model. For parameter values both at and away from the critical viscosity ratio, we present and compare the results of both the diffuse- and sharp-interface models. The close match between both model results indicates that the considered test case lends itself well as a benchmark for further research.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Role of surfactants on droplet formation in piezoacoustic inkjet printing across microsecond-to-second timescales
Authors:
Maaike Rump,
Christian Diddens,
Uddalok Sen,
Michel Versluis,
Detlef Lohse,
Tim Segers
Abstract:
In piezo acoustic drop-on-demand inkjet printing a single droplet is produced for each piezo driving pulse. This droplet is typically multicomponent, including surfactants to control the spreading and drying of the droplet on the substrate. However, the role of these surfactants on the droplet formation process remains rather elusive. Surfactant concentration gradients may manifest across microsec…
▽ More
In piezo acoustic drop-on-demand inkjet printing a single droplet is produced for each piezo driving pulse. This droplet is typically multicomponent, including surfactants to control the spreading and drying of the droplet on the substrate. However, the role of these surfactants on the droplet formation process remains rather elusive. Surfactant concentration gradients may manifest across microsecond-to-second timescales, spanning both the rapid ejection of ink from the nozzle exit and the comparatively slower idling timescale governing the firing of successive droplets. In the present work, we study the influence of surfactants on droplet formation across 6 orders of magnitude in time. To this end, we visualize the microsecond droplet formation process using stroboscopic laser-induced fluorescence microscopy while we vary the nozzle idle time. Our results show that increasing the idle time up to O(1) s affects only the break-up dynamics of the inkjet but not its velocity. By contrast, for idle times $>$ O(1) s, both the break-up dynamics are altered and the velocity of the inkjet increases. We show that the increased velocity results from a decreased surface tension of the ejected droplet, which we extracted from the observed shape oscillations of the jetted droplets in flight. The measured decrease in surface tension is surprising as the $μ$s timescale of droplet formation is much faster than the typical ms-to-s timescale of surfactant adsorption. By varying the bulk surfactant concentration, we show that the fast decrease in surface tension results from a local surfactant concentration increase to more than 200 times the CMC. Our results suggest that a local high concentration of surfactant allows for surfactant adsorption to the interface of an inkjet at the $μ$s-to-ms timescale, which is much faster than the typical ms-to-s timescale of surfactant adsorption.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
Deforming ice with drops
Authors:
Duco van Buuren,
Pallav Kant,
Jochem G. Meijer,
Christian Diddens,
Detlef Lohse
Abstract:
A uniform solidification front undergoes non-trivial deformations when encountering an insoluble dispersed particle in a melt. For solid particles, the overall deformation characteristics are primarily dictated by heat transfer between the particle and the surroundings, remaining unaffected by the rate of approach of the solidification front. In this Letter, we show that, conversely, when interact…
▽ More
A uniform solidification front undergoes non-trivial deformations when encountering an insoluble dispersed particle in a melt. For solid particles, the overall deformation characteristics are primarily dictated by heat transfer between the particle and the surroundings, remaining unaffected by the rate of approach of the solidification front. In this Letter, we show that, conversely, when interacting with a droplet or a bubble, the deformation behaviour exhibits entirely different and unexpected behaviour. It arises from an interfacial dynamics which is specific to particles with free interfaces, namely thermal Marangoni forces. Our study employs a combination of experiments, theory, and numerical simulations to investigate the interaction between the droplet and the freezing front and unveils its surprising behaviour. In particular, we quantitatively understand the dependence of the front deformation $Δ$ on the front propagation velocity, which, for large front velocities, can even revert from attraction ($Δ< 0$) to repulsion ($Δ> 0$).
△ Less
Submitted 29 February, 2024;
originally announced February 2024.
-
Non-monotonic surface tension leads to spontaneous symmetry breaking in a binary evaporating drop
Authors:
Christian Diddens,
Pim J. Dekker,
Detlef Lohse
Abstract:
The evaporation of water/1,2-hexanediol binary drops shows remarkable segregation dynamics, with hexanediol-rich spots forming at the rim, thus breaking axisymmetry. While the segregation of hexanediol near the rim can be attributed to the preferential evaporation of water, the symmetry-breaking and spot formation could not yet be successfully explained. With three-dimensional simulations and azim…
▽ More
The evaporation of water/1,2-hexanediol binary drops shows remarkable segregation dynamics, with hexanediol-rich spots forming at the rim, thus breaking axisymmetry. While the segregation of hexanediol near the rim can be attributed to the preferential evaporation of water, the symmetry-breaking and spot formation could not yet be successfully explained. With three-dimensional simulations and azimuthal stability analysis of a minimal model, we investigate the flow and composition in the drop. We show that a slightly non-monotonic surface tension causes the emergence of a counter-rotating Marangoni vortex in the hexanediol-rich rim region, which subsequently becomes azimuthally unstable and forms the observed spots. Accurate measurements reveal that the surface tension is indeed non-monotonic (~0.3 mN/m). This work provides valuable insight for applications like inkjet printing or spray cooling.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Evaporation of acoustically levitated bicomponent droplets: mass and heat transfer characteristics
Authors:
Yuki Wakata,
Xing Chao,
Chao Sun,
Christian Diddens
Abstract:
Evaporation of multicomponent droplets is important in a wide range of applications, albeit complex, and requires a careful investigation. We experimentally and numerically investigate the evaporation characteristics of spherical, ethanol-water droplets with different initial concentration ratios in the acoustic levitation field. Imaging techniques and infrared thermometry are used for acquiring v…
▽ More
Evaporation of multicomponent droplets is important in a wide range of applications, albeit complex, and requires a careful investigation. We experimentally and numerically investigate the evaporation characteristics of spherical, ethanol-water droplets with different initial concentration ratios in the acoustic levitation field. Imaging techniques and infrared thermometry are used for acquiring volume and surface temperature variations of droplets, reflecting their mass and heat transfer characteristics. Numerical simulations are conducted using modified parameters based on a theoretical model to consider the effect of the acoustic field. The calculation results show good agreement with the experimental data. The concentration and temperature distribution within the droplet is further investigated based on the numerical results.
△ Less
Submitted 24 February, 2024;
originally announced February 2024.
-
Enhanced bubble growth near an advancing solidification front
Authors:
Jochem G. Meijer,
Duarte Rocha,
Annemarie M. Linnenbank,
Christian Diddens,
Detlef Lohse
Abstract:
Frozen water might appear opaque since gas bubbles can get trapped in the ice during the freezing process. They nucleate and then grow near the advancing solidification front, due to the formation of a gas supersaturation region in its vicinity. A delicate interplay between the rate of mass transfer and the rate of freezing dictates the final shapes and sizes of the entrapped gas bubbles. In this…
▽ More
Frozen water might appear opaque since gas bubbles can get trapped in the ice during the freezing process. They nucleate and then grow near the advancing solidification front, due to the formation of a gas supersaturation region in its vicinity. A delicate interplay between the rate of mass transfer and the rate of freezing dictates the final shapes and sizes of the entrapped gas bubbles. In this work, we experimentally and numerically investigate the initial growth of such gas bubbles that nucleate and grow near the advancing ice front. We show that the initial growth of these bubbles is governed by diffusion and is enhanced due to a combination of the presence of the background gas concentration gradient and the motion of the approaching front. Additionally, we recast the problem into that of mass transfer to a moving spherical object in a homogeneous concentration field, finding good agreement between our experimental data and the existing scaling relations for that latter problem. Lastly, we address how fluid flow around the bubble might further affect this growth and qualitatively explore this through numerical simulations.
△ Less
Submitted 12 September, 2024; v1 submitted 9 February, 2024;
originally announced February 2024.
-
Bifurcation tracking on moving meshes and with consideration of azimuthal symmetry breaking instabilities
Authors:
Christian Diddens,
Duarte Rocha
Abstract:
We present a black-box method to numerically investigate the linear stability of arbitrary multi-physics problems. While the user just has to enter the system's residual in weak formulation, i.e. by a finite element method, all required discretized matrices are automatically assembled based on just-in-time generated and compiled highly performant C code. Based on this method, entire phase diagrams…
▽ More
We present a black-box method to numerically investigate the linear stability of arbitrary multi-physics problems. While the user just has to enter the system's residual in weak formulation, i.e. by a finite element method, all required discretized matrices are automatically assembled based on just-in-time generated and compiled highly performant C code. Based on this method, entire phase diagrams in the parameter space can be obtained by bifurcation tracking and continuation within minutes. Particular focus is put on problems with moving domains, e.g. free surface problems in fluid dynamics, since a moving mesh introduces a plethora of complicated nonlinearities to the system. By symbolic differentiation before the code generation, however, these moving mesh problems are made accessible to bifurcation tracking methods. In a second step, our method is generalized to investigate symmetry-breaking instabilities of axisymmetric stationary solutions by effectively utilizing the symmetry of the base state. Each bifurcation type is validated on the basis of results reported in the literature on versatile fluid dynamics problems.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
On the rising and sinking motion of bouncing oil drops in strongly stratified liquids
Authors:
Jochem G. Meijer,
Yanshen Li,
Christian Diddens,
Detlef Lohse
Abstract:
When an immiscible oil drop is immersed in a stably stratified ethanol-water mixture, the Marangoni flow on the surface of the drop can experience an oscillatory instability, so that the drop undergoes a transition from levitating to bouncing. The onset of the instability and its mechanisms have been studied previously, yet the bouncing motion of the drop itself, which is a completely different pr…
▽ More
When an immiscible oil drop is immersed in a stably stratified ethanol-water mixture, the Marangoni flow on the surface of the drop can experience an oscillatory instability, so that the drop undergoes a transition from levitating to bouncing. The onset of the instability and its mechanisms have been studied previously, yet the bouncing motion of the drop itself, which is a completely different problem, has not yet been investigated. Here we study how the bouncing characteristics (jumping height, rising and sinking time) depend on the control parameters (drop radius, stratification strength, drop viscosity). We first record experimentally the bouncing trajectories of drops of different viscosities in different stratifications. Then a simplified dynamical analysis is performed to get the scaling relations of the jumping height and the rising and sinking times. The rising and sinking time scales are found to depend on the drag coefficient of the drop $C_D^S$ in the stratified liquid, which is determined empirically for the current parameter space. For low viscosity (5 cSt) oil drops the results on the drag coefficient match the ones from the literature. For high viscosity (100 cSt) oil drops the parameter space had not been explored and the drag coefficients are not readily available. Numerical simulations are therefore performed to provide external verification for the drag coefficients, which well match with the experimental results.
△ Less
Submitted 29 June, 2023;
originally announced June 2023.
-
Evaporation of binary liquids from a capillary tube
Authors:
Lijun Thayyil Raju,
Christian Diddens,
Javier Rodríguez-Rodríguez,
Marjolein N. van der Linden,
Xuehua Zhang,
Detlef Lohse,
Uddalok Sen
Abstract:
Evaporation of multi-component liquid mixtures in confined geometries, such as capillaries, is crucial in applications such as microfluidics, two-phase cooling devices, and inkjet printing. Predicting the behaviour of such systems becomes challenging because evaporation triggers complex spatio-temporal changes in the composition of the mixture. These changes in composition, in turn, affect evapora…
▽ More
Evaporation of multi-component liquid mixtures in confined geometries, such as capillaries, is crucial in applications such as microfluidics, two-phase cooling devices, and inkjet printing. Predicting the behaviour of such systems becomes challenging because evaporation triggers complex spatio-temporal changes in the composition of the mixture. These changes in composition, in turn, affect evaporation. In the present work, we study the evaporation of aqueous glycerol solutions contained as a liquid column in a capillary tube. Experiments and direct numerical simulations show three evaporation regimes characterised by different temporal evolutions of the normalised mass transfer rate (or Sherwood number, $Sh$), namely $Sh (\tilde{t}) = 1$, $Sh \sim 1/\sqrt{\tilde{t}}$, and $Sh \sim \exp\left(-\tilde{t}\right)$. Here $\tilde{t}$ is a normalised time. We present a simplistic analytical model which shows that the evaporation dynamics can be expressed by the classical relation $Sh = \exp \left( \tilde{t} \right) \mathrm{erfc} \left( \sqrt{\tilde{t}}\right)$. For small and medium $\tilde{t}$, this expression results in the first and second of the three observed scaling regimes, respectively. This analytical model is formulated in the limit of pure diffusion and when the penetration depth $δ(t)$ of the diffusion front is much smaller than the length $L(t)$ of the liquid column. When $δ\approx L$, finite length effects lead to $Sh \sim \exp\left(-\tilde{t}\right)$, i.e. the third regime. Finally, we extend our analytical model to incorporate the effect of advection and determine the conditions under which this effect is important. Our results provide fundamental insight into the physics of selective evaporation from a multi-component liquid column.
△ Less
Submitted 18 November, 2023; v1 submitted 11 November, 2022;
originally announced November 2022.
-
Interfacial aggregation of self-propelled Janus colloids in sessile droplets
Authors:
Maziyar Jalaal,
Borge ten Hagen,
Hai le The,
Christian Diddens,
Detlef Lohse,
Alvaro Marin
Abstract:
Living microorganisms in confined systems typically experience an affinity to populate boundaries. The reason for such affinity to interfaces can be a combination of their directed motion and hydrodynamic interactions at distances larger than their own size. Here we will show that self-propelled Janus particles (polystyrene particles partially coated with platinum) immersed in droplets of water an…
▽ More
Living microorganisms in confined systems typically experience an affinity to populate boundaries. The reason for such affinity to interfaces can be a combination of their directed motion and hydrodynamic interactions at distances larger than their own size. Here we will show that self-propelled Janus particles (polystyrene particles partially coated with platinum) immersed in droplets of water and hydrogen peroxide tend to accumulate in the vicinity of the liquid/gas interface. Interestingly, the interfacial accumulation occurs despite the presence of an evaporation-driven flow caused by a solutal Marangoni flow, which typically tends to redistribute the particles within the droplet's bulk. By performing additional experiments with passive colloids (flow tracers) and comparing with numerical simulations for both particle active motion and the fluid flow, we disentangle the dominating mechanisms behind the observed interfacial particle accumulation. These results allow us to make an analogy between active Janus particles and some biological microswimmers concerning how they interact with their environment.
△ Less
Submitted 1 December, 2022; v1 submitted 6 October, 2022;
originally announced October 2022.
-
Evaporation-driven liquid flow in sessile droplets
Authors:
Hanneke Gelderblom,
Christian Diddens,
Alvaro Marin
Abstract:
The evaporation of a sessile droplet spontaneously induces an internal capillary liquid flow. The surface-tension driven minimisation of surface area and/or surface-tension differences at the liquid-gas interface caused by evaporation-induced temperature or chemical gradients set the liquid into motion. This flow drags along suspended material and is one of the keys to control the material deposit…
▽ More
The evaporation of a sessile droplet spontaneously induces an internal capillary liquid flow. The surface-tension driven minimisation of surface area and/or surface-tension differences at the liquid-gas interface caused by evaporation-induced temperature or chemical gradients set the liquid into motion. This flow drags along suspended material and is one of the keys to control the material deposition in the stain that is left behind by a drying droplet. Applications of this principle range from the control of stain formation in the printing and coating industry, to the analysis of DNA, to forensic and medical research on blood stains, and to the use of evaporation-driven self-assembly for nanotechnology. Therefore, the evaporation of sessile droplets attracts an enormous interest from not only the fluid dynamics, but also the soft matter, chemistry, biology, engineering, nanotechnology and mathematics communities. As a consequence of this broad interest, knowledge on evaporation-driven flows in drying droplets has remained scattered among the different fields, leading to various misconceptions and misinterpretations. In this review we aim to unify these views, and reflect on the current understanding of evaporation-driven liquid flows in sessile droplets in the light of the most recent experimental and theoretical advances. In addition, we outline open questions and indicate promising directions for future research.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
Evaporation of a sessile colloidal water-glycerol droplet: Marangoni ring formation
Authors:
Lijun Thayyil Raju,
Christian Diddens,
Yaxing Li,
Alvaro Marin,
Marjolein N. van der Linden,
Xuehua Zhang,
Detlef Lohse
Abstract:
The transport and aggregation of particles in suspensions is an important process in many physicochemical and industrial processes. In this work, we study the transport of particles in an evaporating binary droplet. Surprisingly, the accumulation of particles occurs not only at the contact line (due to the coffee-stain effect) or at the solid substrate (due to sedimentation), but also at a particu…
▽ More
The transport and aggregation of particles in suspensions is an important process in many physicochemical and industrial processes. In this work, we study the transport of particles in an evaporating binary droplet. Surprisingly, the accumulation of particles occurs not only at the contact line (due to the coffee-stain effect) or at the solid substrate (due to sedimentation), but also at a particular radial position near the liquid-air interface, forming a "ring", which we term as the Marangoni ring. The formation of this ring is primarily attributed to the solutal Marangoni flow triggered by the evaporation dynamics of the water-glycerol droplet. Experiments and simulations show fair agreement in the volume evolution and the general structure of the solutal Marangoni flow, that is, the Marangoni vortex. Experiments show that the location of the Marangoni ring is strongly correlated with the Marangoni vortex. However, finite element numerical simulations fail to describe the particle distribution seen in the experiments. Interestingly, the particles not only accumulate to form the Marangoni ring but also assemble as colloidal crystals close to the liquid-air interface, yielding iridescence. The formation of the colloidal crystals in the experiments is strong evidence that non-hydrodynamic interactions, which are not represented in the simulations, also play a significant role in our system.
△ Less
Submitted 20 September, 2022; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Physiochemical hydrodynamics of the phase segregation in an evaporating binary microdroplet
Authors:
Yaxing Li,
Pengyu Lv,
Christian Diddens,
Detlef Lohse
Abstract:
Phase segregation triggered by selective evaporation can emerge in multicomponent systems, leading to complex physiochemical hydrodynamics. Recently, Li et al. (Phys. Rev. Lett., vol. 120, 2018, 224501) and Kim & Stone (J. Fluid Mech., vol. 850, 2018, pp. 769-783) reported a segregative behavior (i.e., demixing) in an evaporating binary droplet. In this work, by means of experiments and theoretica…
▽ More
Phase segregation triggered by selective evaporation can emerge in multicomponent systems, leading to complex physiochemical hydrodynamics. Recently, Li et al. (Phys. Rev. Lett., vol. 120, 2018, 224501) and Kim & Stone (J. Fluid Mech., vol. 850, 2018, pp. 769-783) reported a segregative behavior (i.e., demixing) in an evaporating binary droplet. In this work, by means of experiments and theoretical analysis, we investigate the flow dynamics after the occurrence of the phase segregation. As example, we take the 1,2-hexanediol-water binary droplet system. First, we experimentally reveal the overall physiochemical hydrodynamics of the evaporation process, including the segregative behavior and the resulting flow structure close to the substrate. By quantifying the evolution of the radial flow, we identify three successive life stages of the evaporation process. At Stage I, a radially outward flow is observed. It is driven by the Marangoni effect. At the transition to Stage II, the radial flow partially reverses, starting from the contact line. This flow breaks the axial symmetry and remarkably is driven by the segregation itself. Finally at Stage III, the flow decays as the evaporation gradually ceases. At this stage the segregation has grown to the entire droplet, and the flow is again controlled by the Marangoni effect. The resulting Marangoni flow homogenizes the distribution of the entrapped volatile water over the whole droplet.
△ Less
Submitted 16 August, 2022;
originally announced August 2022.
-
Sessile drop evaporation in a gap -- crossover between diffusion-limited and phase transition-limited regime
Authors:
S. Hartmann,
C. Diddens,
M. Jalaal,
U. Thiele
Abstract:
We consider the time evolution of a sessile drop of volatile partially wetting liquid on a rigid solid substrate. Thereby, the drop evaporates under strong confinement, namely, it sits on one of the two parallel plates that form a narrow gap. First, we develop an efficient mesoscopic thin-film description in gradient dynamics form. It couples the diffusive dynamics of the vertically averaged vapou…
▽ More
We consider the time evolution of a sessile drop of volatile partially wetting liquid on a rigid solid substrate. Thereby, the drop evaporates under strong confinement, namely, it sits on one of the two parallel plates that form a narrow gap. First, we develop an efficient mesoscopic thin-film description in gradient dynamics form. It couples the diffusive dynamics of the vertically averaged vapour density in the narrow gap to an evolution equation for the profile of the volatile drop. The underlying free energy functional incorporates wetting, interface and bulk energies of the liquid and gas entropy. The model allows us to investigate the transition between diffusion-limited and phase transition-limited evaporation for shallow droplets. Its gradient dynamics character also allows for a full-curvature formulation.
Second, we compare results obtained with the mesoscopic model to corresponding direct numerical simulations solving the Stokes equation for the drop coupled to the diffusion equation for the vapour as well as to selected experiments. In passing, we discuss the influence of contact line pinning.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Selective evaporation at the nozzle exit in piezoacoustic inkjet printing
Authors:
Maaike Rump,
Uddalok Sen,
Roger Jeurissen,
Hans Reinten,
Michel Versluis,
Detlef Lohse,
Christian Diddens,
Tim Segers
Abstract:
In practical applications of inkjet printing the nozzles in a printhead have intermittent idle periods, during which ink can evaporate from the nozzle exit. Inks are usually multicomponent where each component has its own characteristic evaporation rate resulting in concentration gradients within the ink. These gradients may directly and indirectly (via Marangoni flows) alter the jetting process a…
▽ More
In practical applications of inkjet printing the nozzles in a printhead have intermittent idle periods, during which ink can evaporate from the nozzle exit. Inks are usually multicomponent where each component has its own characteristic evaporation rate resulting in concentration gradients within the ink. These gradients may directly and indirectly (via Marangoni flows) alter the jetting process and thereby its reproducibility and the resulting print quality. In the present work, we study selective evaporation from an inkjet nozzle for water-glycerol mixtures. Through experiments, analytical modeling, and numerical simulations, we investigate changes in mixture composition with drying time. By monitoring the acoustics within the printhead, and subsequently modeling the system as a mass-spring-damper system, the composition of the mixture can be obtained as a function of drying time. The results from the analytical model are validated using numerical simulations of the full fluid mechanical equations governing the printhead flows and pressure fields. Furthermore, the numerical simulations reveal that the time independent concentration gradient we observe in the experiments is due to the steady state of water flux through the printhead. Finally, we measure the number of drop formation events required in this system before the mixture concentration within the nozzle attains the initial (pre-drying) value, and find a stronger than exponential trend in the number of drop formations required. These results shed light on the complex physiochemical hydrodynamics associated with the drying of ink at a printhead nozzle, and help in increasing the stability and reproducibility of inkjet printing.
△ Less
Submitted 27 February, 2023; v1 submitted 13 May, 2022;
originally announced May 2022.
-
A Robust and Accurate Adaptive Approximation Method for a Diffuse-Interface Model of Binary-Fluid Flows
Authors:
T. H. B. Demont,
G. J. van Zwieten,
C. Diddens,
E. H. van Brummelen
Abstract:
We present an adaptive simulation framework for binary-fluid flows, based on the Abels-Garcke-Grün Navier-Stokes-Cahn-Hilliard (AGG NSCH) diffuse-interface model. The adaptive-refinement procedure is guided by a two-level hierarchical a-posteriori error estimate, and it effectively resolves the spatial multiscale behavior of the diffuse-interface model. To improve the robustness of the solution pr…
▽ More
We present an adaptive simulation framework for binary-fluid flows, based on the Abels-Garcke-Grün Navier-Stokes-Cahn-Hilliard (AGG NSCH) diffuse-interface model. The adaptive-refinement procedure is guided by a two-level hierarchical a-posteriori error estimate, and it effectively resolves the spatial multiscale behavior of the diffuse-interface model. To improve the robustness of the solution procedure and avoid severe time-step restrictions for small-interface thicknesses, we introduce an $\varepsilon$-continuation procedure, in which the diffuse interface thickness ($\varepsilon$) are enlarged on coarse meshes, and the mobility is scaled accordingly. To further accelerate the computations and improve robustness, we apply a modified Backward Euler scheme in the initial stages of the adaptive-refinement procedure in each time step, and a Crank--Nicolson scheme in the final stages of the refinement procedure. To enhance the robustness of the nonlinear solution procedure, we introduce a partitioned solution procedure for the linear tangent problems in Newton's method, based on a decomposition of the NSCH system into its NS and CH subsystems. We conduct a systematic investigation of the conditioning of the monolithic NSCH tangent matrix and of its NS and CH subsystems for a representative 2D model problem. To illustrate the properties of the presented adaptive simulation framework, we present numerical results for a 2D oscillating water droplet suspended in air, and we validate the obtained results versus those of a corresponding sharp-interface model.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Droplet dissolution driven by emerging thermal gradients and Marangoni flow
Authors:
Binglin Zeng,
Yuliang Wang,
Christian Diddens,
Harold J. W. Zandvliet,
Detlef Lohse
Abstract:
The lifetime $τ$ of an isothermal and purely diffusively dissolving droplet in a host liquid scales as $τ\sim R_0^2$ with its initial radius $R_0$ [Langmuir, Phys. Rev. 12, 368 (1919)]. For a droplet dissolving due to natural convection driven by density differences, its lifetime scales as $τ\sim R_0^{5/4}$ [Dietrich et al., J. Fluid Mech. 794, 45 (2016)]. In this paper we experimentally find and…
▽ More
The lifetime $τ$ of an isothermal and purely diffusively dissolving droplet in a host liquid scales as $τ\sim R_0^2$ with its initial radius $R_0$ [Langmuir, Phys. Rev. 12, 368 (1919)]. For a droplet dissolving due to natural convection driven by density differences, its lifetime scales as $τ\sim R_0^{5/4}$ [Dietrich et al., J. Fluid Mech. 794, 45 (2016)]. In this paper we experimentally find and theoretically derive yet another droplet dissolution behavior, resulting in $τ\sim R_0^4$. It occurs when the dissolution dynamics is controlled by local heating of the liquid, leading to a modified solubility and a thermal Marangoni flow around the droplet. The thermal gradient is achieved by plasmonic heating of a gold nanoparticle decorated sample surface, on which a sessile water droplet immersed in water-saturated 1-butanol solution is sitting. The resulting off-wall thermal Marangoni flow and the temperature dependence of the solubility determine the droplet dissolution rate, resulting in a shrinkage $R(t) \sim (τ-t )^{1/4}$ of the droplet radius and thus in $τ\sim R_0^{4}$.
△ Less
Submitted 11 February, 2022;
originally announced February 2022.
-
Oscillatory droplet dissolution from competing Marangoni and gravitational flows
Authors:
Ricardo Arturo Lopez de la Cruz,
Christian Diddens,
Xuehua Zhang,
Detlef Lohse
Abstract:
The dissolution or growth of a droplet in a host liquid is an important part for processes like chemical extraction, chromatography or emulsification. In this work we look at the dissolution of a pair of vertically aligned droplets immersed in water, both experimentally and with numerical simulations. The liquids used for the droplets are long chain alcohols with a low but finite solubility in wat…
▽ More
The dissolution or growth of a droplet in a host liquid is an important part for processes like chemical extraction, chromatography or emulsification. In this work we look at the dissolution of a pair of vertically aligned droplets immersed in water, both experimentally and with numerical simulations. The liquids used for the droplets are long chain alcohols with a low but finite solubility in water and a significantly lower density than that of the host liquid. Therefore, a solutal plume is formed above of the bottom droplet and natural convection dominates the dissolution process. We monitor the volume of the droplets and the velocity field around them over time. When the liquids of the two droplets are the same, our previously found scaling laws for the Sherwood and Reynolds numbers as functions of the Rayleigh number (Dietrich et al., 2016, J. Fluid Mech.) can be applied to the lower droplet. However, remarkably, when the liquid of the top droplet is different than that of the bottom droplet the volume as function of time becomes non-monotonic, and an oscillatory Marangoni flow at the top droplet is observed. We identify the competition between solutal Marangoni flow and density driven convection as the origin of the oscillation, and numerically model the process.
△ Less
Submitted 22 January, 2022;
originally announced January 2022.
-
The asymmetric coalescence of two droplets with different surface tensions is caused by capillary waves
Authors:
Michiel A. Hack,
Patrick Vondeling,
Menno Cornelissen,
Detlef Lohse,
Jacco H. Snoeijer,
Christian Diddens,
Tim Segers
Abstract:
When two droplets with different surface tensions collide, the shape evolution of the merging droplets is asymmetric. Using experimental and numerical techniques, we reveal that this asymmetry is caused by asymmetric capillary waves, which are the result of the different surface tensions of the droplets. We show that the asymmetry is enhanced by increasing the surface tension difference, and suppr…
▽ More
When two droplets with different surface tensions collide, the shape evolution of the merging droplets is asymmetric. Using experimental and numerical techniques, we reveal that this asymmetry is caused by asymmetric capillary waves, which are the result of the different surface tensions of the droplets. We show that the asymmetry is enhanced by increasing the surface tension difference, and suppressed by increasing the inertia of the colliding droplets. Furthermore, we study capillary waves in the limit of no inertia. We reveal that the asymmetry is not directly caused by Marangoni forces. In fact, somehow counterintuitive, asymmetry is strongly reduced by the Marangoni effect. Rather, the different intrinsic capillary wave amplitudes and velocities associated with the different surface tensions of the droplets lie at the origin of the asymmetry during droplet coalescence.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
Time-resolved velocity and pressure field quantification in a flow-focusing device for ultrafast microbubble production
Authors:
Sarah Cleve,
Christian Diddens,
Tim Segers,
Guillaume Lajoinie,
Michel Versluis
Abstract:
Flow-focusing devices have gained great interest in the past decade, due to their capability to produce monodisperse microbubbles for diagnostic and therapeutic medical ultrasound applications. However, up-scaling production to industrial scale requires a paradigm shift from single chip operation to highly parallelized systems. Parallelization gives rise to fluidic interactions between nozzles tha…
▽ More
Flow-focusing devices have gained great interest in the past decade, due to their capability to produce monodisperse microbubbles for diagnostic and therapeutic medical ultrasound applications. However, up-scaling production to industrial scale requires a paradigm shift from single chip operation to highly parallelized systems. Parallelization gives rise to fluidic interactions between nozzles that, in turn, may lead to a decreased monodispersity. Here, we study the velocity and pressure field fluctuations in a single flow-focusing nozzle during bubble production. We experimentally quantify the velocity field inside the nozzle at 100 ns time resolution, and a numerical model provides insight into both the oscillatory velocity and pressure fields. Our results demonstrate that, at the length scale of the flow focusing channel, the velocity oscillations propagate at fluid dynamical time scale (order of microseconds) whereas the dominant pressure oscillations are linked to the bubble pinch-off and propagate at a much faster time scale (order of nanoseconds).
△ Less
Submitted 5 November, 2021;
originally announced November 2021.
-
Marangoni instability triggered by selective evaporation of a binary liquid inside a Hele-Shaw cell
Authors:
Ricardo Arturo Lopez de la Cruz,
Christian Diddens,
Xuehua Zhang,
Detlef Lohse
Abstract:
Interfacial stability is important for many processes involving heat and mass transfer across two immiscible phases. When this transfer takes place in the form of evaporation of a binary solution with one component being more volatile than the other, gradients in surface tension can arise. These gradients can ultimately destabilise the liquid-gas interface. In the present work, we study the evapor…
▽ More
Interfacial stability is important for many processes involving heat and mass transfer across two immiscible phases. When this transfer takes place in the form of evaporation of a binary solution with one component being more volatile than the other, gradients in surface tension can arise. These gradients can ultimately destabilise the liquid-gas interface. In the present work, we study the evaporation of an ethanol-water solution, for which ethanol has a larger volatility. The solution is contained in a horizontal Hele-Shaw cell which is open from one end to allow for evaporation into air. A Marangoni instability is then triggered at the liquid-air interface. We study the temporal evolution of this instability by observing the effects that it has on the bulk of the liquid. More specifically, the growth of convective cells is visualized with confocal microscopy and the velocity field close to the interface is measured with micro-particle-image-velocimetry. The results of numerical simulations based on quasi 2D equations satisfactorily compare with the experimental observations, even without consideration of evaporative cooling, although this cooling can play an extra role in experiments. Furthermore, a linear stability analysis applied to a simplified version of the quasi 2D equations showed reasonably good agreement with the results from simulations at early times, when the instability has just been triggered and no coarsening has taken place. In particular, we find a critical Marangoni number below which a regime of stability is predicted.
△ Less
Submitted 13 June, 2021;
originally announced June 2021.
-
Periodic bouncing of a plasmonic bubble in a binary liquid by competing solutal and thermal Marangoni forces
Authors:
Binglin Zeng,
Kai Leong Chong,
Yuliang Wang,
Christian Diddens,
Xiaolai Li,
Marvin Detert,
Harold J. W. Zandvliet,
Detlef Lohse
Abstract:
The physicochemical hydrodynamics of bubbles and droplets out of equilibrium, in particular with phase transitions, displays surprisingly rich and often counterintuitive phenomena. Here we experimentally and theoretically study the nucleation and early evolution of plasmonic bubbles in a binary liquid consisting of water and ethanol. Remarkably, the submillimeter plasmonic bubble is found to be pe…
▽ More
The physicochemical hydrodynamics of bubbles and droplets out of equilibrium, in particular with phase transitions, displays surprisingly rich and often counterintuitive phenomena. Here we experimentally and theoretically study the nucleation and early evolution of plasmonic bubbles in a binary liquid consisting of water and ethanol. Remarkably, the submillimeter plasmonic bubble is found to be periodically attracted to and repelled from the nanoparticle-decorated substrate, with frequencies of around a few kHz. We identify the competition between solutal and thermal Marangoni forces as origin of the periodic bouncing. The former arises due to the selective vaporization of ethanol at the substrate's side of the bubble, leading to a solutal Marangoni flow towards the hot substrate, which pushes the bubble away. The latter arises due to the temperature gradient across the bubble, leading to a thermal Marangoni flow away from the substrate which sucks the bubble towards it. We study the dependence of the frequency of the bouncing phenomenon from the control parameters of the system, namely the ethanol fraction and the laser power for the plasmonic heating. Our findings can be generalized to boiling and electrolytically or catalytically generated bubbles in multicomponent liquids.
△ Less
Submitted 20 April, 2021;
originally announced April 2021.
-
Marangoni instability of a drop in a stably stratified liquid
Authors:
Yanshen Li,
Christian Diddens,
Andrea Prosperetti,
Detlef Lohse
Abstract:
Marangoni instabilities can emerge when a liquid interface is subjected to a concentration or temperature gradient. It is generally believed that for these instabilities bulk effects like buoyancy are negligible as compared to interfacial forces, especially on small scales. Consequently, the effect of a stable stratification on the Marangoni instability has hitherto been ignored. Here we report, f…
▽ More
Marangoni instabilities can emerge when a liquid interface is subjected to a concentration or temperature gradient. It is generally believed that for these instabilities bulk effects like buoyancy are negligible as compared to interfacial forces, especially on small scales. Consequently, the effect of a stable stratification on the Marangoni instability has hitherto been ignored. Here we report, for an immiscible drop immersed in a stably stratified ethanol-water mixture, a new type of oscillatory solutal Marangoni instability which is triggered once the stratification has reached a critical value. We experimentally explore the parameter space spanned by the stratification strength and the drop size and theoretically explain the observed crossover from levitating to bouncing by balancing the advection and diffusion around the drop. Finally, the effect of the stable stratification on the Marangoni instability is surprisingly amplified in confined geometries, leading to an earlier onset.
△ Less
Submitted 10 March, 2021; v1 submitted 27 January, 2021;
originally announced January 2021.
-
Competing Marangoni and Rayleigh convection in evaporating binary droplets
Authors:
Christian Diddens,
Yaxing Li,
Detlef Lohse
Abstract:
For a small sessile or pendant droplet it is generally assumed that gravity does not play any role once the Bond number is small. This is even assumed for evaporating binary sessile or pendant droplets, in which convective flows can be driven due to selective evaporation of one component and the resulting concentration and thus surface tension differences at the air-liquid interface. However, rece…
▽ More
For a small sessile or pendant droplet it is generally assumed that gravity does not play any role once the Bond number is small. This is even assumed for evaporating binary sessile or pendant droplets, in which convective flows can be driven due to selective evaporation of one component and the resulting concentration and thus surface tension differences at the air-liquid interface. However, recent studies have shown that in such droplets gravity indeed can play a role and that natural convection can be the dominant driving mechanism for the flow inside evaporating binary droplets (Edwards et al., Phys. Rev. Lett. 121, 184501 (2018); Li et al., Phys. Rev. Lett. 122, 114501 (2019)). In this study, we derive and validate a quasi-stationary model for the flow inside evaporating binary sessile and pendant droplets, which successfully allows to predict the prevalence and the intriguing interaction of Rayleigh and/or Marangoni convection on the basis of a phase diagram for the flow field expressed in terms of the Rayleigh and Marangoni numbers.
△ Less
Submitted 28 May, 2020;
originally announced May 2020.
-
Rayleigh-Taylor instability by segregation in an evaporating multi-component microdroplet
Authors:
Yaxing Li,
Christian Diddens,
Tim Segers,
Herman Wijshoff,
Michel Versluis,
Detlef Lohse
Abstract:
The evaporation of multi-component droplets is relevant to various applications but challenging to study due to the complex physicochemical dynamics. Recently, Li (2018) reported evaporation-triggered segregation in 1,2-hexanediol-water binary droplets. In this present work, we added 0.5 wt% silicone oil into the 1,2-hexanediol-water binary solution. This minute silicone oil concentration dramatic…
▽ More
The evaporation of multi-component droplets is relevant to various applications but challenging to study due to the complex physicochemical dynamics. Recently, Li (2018) reported evaporation-triggered segregation in 1,2-hexanediol-water binary droplets. In this present work, we added 0.5 wt% silicone oil into the 1,2-hexanediol-water binary solution. This minute silicone oil concentration dramatically modifies the evaporation process as it triggers an early extraction of the 1,2-hexanediol from the mixture. Surprisingly, we observe that the segregation of 1,2-hexanediol forms plumes, rising up from the rim of the sessile droplet towards the apex during the droplet evaporation. By orientating the droplet upside down, i.e., by studying a pendant droplet, the absence of the plumes indicates that the flow structure is induced by buoyancy, which drives a Rayleigh-Taylor instability (i.e., driven by density differences & gravitational acceleration). From micro-PIV measurement, we further prove that the segregation of the non-volatile component (1,2-hexanediol) hinders the evaporation near the contact line, which leads to a suppression of the Marangoni flow in this region. Hence, on long time scales, gravitational effects play the dominant role in the flow structure, rather than Marangoni flows. We compare the measurement of the evaporation rate with the diffusion model of Popov (2005), coupled with Raoult's law and the activity coefficient. This comparison indeed confirms that the silicone-oil-triggered segregation of the non-volatile 1,2-hexanediol significantly delays the evaporation. With an extended diffusion model, in which the influence of the segregation has been implemented, the evaporation can be well described.
△ Less
Submitted 28 May, 2020;
originally announced May 2020.
-
Gravitational effect in evaporating binary microdroplets
Authors:
Yaxing Li,
Christian Diddens,
Pengyu Lv,
Herman Wijshoff,
Michel Versluis,
Detlef Lohse
Abstract:
The flow in an evaporating glycerol-water binary sub-millimeter droplet with Bond number Bo $\ll$ 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by micro-PIV for both sessile and pendant droplets during the evaporation process, which surprisingly show opposite radial flow directions -- inward and outward, respectively. This observation clearl…
▽ More
The flow in an evaporating glycerol-water binary sub-millimeter droplet with Bond number Bo $\ll$ 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by micro-PIV for both sessile and pendant droplets during the evaporation process, which surprisingly show opposite radial flow directions -- inward and outward, respectively. This observation clearly reveals that in spite of the small droplet size, gravitational effects play a crucial role in controlling the flow fields in the evaporating droplets. We theoretically analyze that this gravity-driven effect is triggered by the lower volatility of glycerol which leads to a preferential evaporation of water then the local concentration difference of the two components leads to a density gradient that drives the convective flow. We show that the Archimedes number Ar is the nondimensional control parameter for the occurrence of the gravitational effects. We confirm our hypothesis by experimentally comparing two evaporating microdroplet systems, namely a glycerol-water droplet and a 1,2-propanediol-water droplet. We obtain different Ar, larger or smaller than a unit by varying a series of droplet heights, which corresponds to cases with or without gravitational effects, respectively. Finally, we simulate the process numerically, finding good agreement with the experimental results and again confirming our interpretation.
△ Less
Submitted 5 February, 2019;
originally announced February 2019.
-
A bouncing oil droplet in a stratified liquid and its sudden death
Authors:
Yanshen Li,
Christian Diddens,
Andrea Prosperetti,
Kai Leong Chong,
Xuehua Zhang,
Detlef Lohse
Abstract:
Droplets can self-propel when immersed in another liquid in which a concentration gradient is present. Here we report the experimental and numerical study of a self-propelling oil droplet in a vertically stratified ethanol/water mixture: At first, the droplet sinks slowly due to gravity, but then, before having reached its density matched position, jumps up suddenly. More remarkably, the droplet b…
▽ More
Droplets can self-propel when immersed in another liquid in which a concentration gradient is present. Here we report the experimental and numerical study of a self-propelling oil droplet in a vertically stratified ethanol/water mixture: At first, the droplet sinks slowly due to gravity, but then, before having reached its density matched position, jumps up suddenly. More remarkably, the droplet bounces repeatedly with an ever increasing jumping distance, until all of a sudden it stops after about 30 min. We identify the Marangoni stress at the droplet/liquid interface as responsible for the jumping: its strength grows exponentially because it pulls down ethanol-rich liquid, which in turn increases its strength even more. The jumping process can repeat because gravity restores the system. Finally, the sudden death of the jumping droplet is also explained. Our findings have demonstrated a type of prominent droplet bouncing inside a continuous medium with no wall or sharp interface.
△ Less
Submitted 4 December, 2018;
originally announced December 2018.
-
Microdroplets nucleation by dissolution of a multicomponent drop in a host liquid
Authors:
Huanshu Tan,
Christian Diddens,
Ali Akash Mohammed,
Junyi Li,
Michel Versluis,
Xuehua Zhang,
Detlef Lohse
Abstract:
Multicomponent liquid drops in a host liquid are very relevant in various technological applications. Their dissolution or growth dynamics is complex. Differences in solubility between the drop components combined with the solutal Marangoni effect and natural convections contribute to this complexity, which can be even further increased in combination with the ouzo effect, i.e., the spontaneous nu…
▽ More
Multicomponent liquid drops in a host liquid are very relevant in various technological applications. Their dissolution or growth dynamics is complex. Differences in solubility between the drop components combined with the solutal Marangoni effect and natural convections contribute to this complexity, which can be even further increased in combination with the ouzo effect, i.e., the spontaneous nucleation of microdroplets due to composition-dependent miscibilities in a ternary system. The quantitive understanding of this combined process is important for applications in industry, particularly for modern liquid-liquid microextraction processes. In this work, as a model system, we experimentally and theoretically explore water/ethanol drops dissolving in anethole oil. During the dissolution, we observed two types of microdroplets nucleation, namely water microdroplet nucleation in the surrounding oil at drop midheight and oil microdroplet nucleation in the aqueous drop, again at midheight. The nucleated oil microdroplets are driven by Marangoni flows inside the aqueous drop and evolve into microdroplets rings. A one-dimensional multiphase and multicomponent diffusion model in combination with thermodynamical equilibrium theory is proposed to predict the behavior of spontaneous emulsification, i.e. the microdroplet nucleation, that is triggered by diffusion. A scale analysis together with experimental investigations of the fluid dynamics of the system reveals that both the solutal Marangoni flow inside the drop and the buoyancy-driven flow in the host liquid influence the diffusion-triggered emulsification process. Our work provides a physical understanding of the microdroplet nucleation by dissolution of a multicomponent drop in a host liquid.
△ Less
Submitted 8 November, 2018;
originally announced November 2018.
-
Evaporation-triggered segregation of sessile binary droplets
Authors:
Yaxing Li,
Pengyu Lv,
Christian Diddens,
Huanshu Tan,
Herman Wijshoff,
Michel Versluis,
Detlef Lohse
Abstract:
Droplet evaporation of multicomponent droplets is essential for various physiochemical applications, e.g. in inkjet printing, spray cooling and microfabrication. In this work, we observe and study phase segregation of an evaporating sessile binary droplet, consisting of a mixture of water and a surfactant-like liquid (1,2-hexanediol). The phase segregation (i.e., demixing) leads to a reduced water…
▽ More
Droplet evaporation of multicomponent droplets is essential for various physiochemical applications, e.g. in inkjet printing, spray cooling and microfabrication. In this work, we observe and study phase segregation of an evaporating sessile binary droplet, consisting of a mixture of water and a surfactant-like liquid (1,2-hexanediol). The phase segregation (i.e., demixing) leads to a reduced water evaporation rate of the droplet and eventually the evaporation process ceases due to shielding of the water by the non-volatile 1,2-hexanediol. Visualizations of the flow field by particle image velocimetry and numerical simulations reveal that the timescale of water evaporation at the droplet rim is faster than that of the Marangoni flow, which originates from the surface tension difference between water and 1,2-hexanediol, eventually leading to segregation.
△ Less
Submitted 30 January, 2018;
originally announced January 2018.
-
Evaporating pure, binary and ternary droplets: thermal effects and axial symmetry breaking
Authors:
Christian Diddens,
Huanshu Tan,
Pengyu Lv,
Michel Versluis,
J. G. M. Kuerten,
Xuehua Zhang,
Detlef Lohse
Abstract:
The Greek aperitif Ouzo is not only famous for its specific anise-flavored taste, but also for its ability to turn from a transparent miscible liquid to a milky-white colored emulsion when water is added. Recently, it has been shown that this so-called Ouzo effect, i.e. the spontaneous emulsification of oil microdroplets, can also be triggered by the preferential evaporation of ethanol in an evapo…
▽ More
The Greek aperitif Ouzo is not only famous for its specific anise-flavored taste, but also for its ability to turn from a transparent miscible liquid to a milky-white colored emulsion when water is added. Recently, it has been shown that this so-called Ouzo effect, i.e. the spontaneous emulsification of oil microdroplets, can also be triggered by the preferential evaporation of ethanol in an evaporating sessile Ouzo drop, leading to an amazingly rich drying process with multiple phase transitions [H. Tan et al., Proc. Natl. Acad. Sci. USA 113(31) (2016) 8642]. Due to the enhanced evaporation near the contact line, the nucleation of oil droplets starts at the rim which results in an oil ring encircling the drop. Furthermore, the oil droplets are advected through the Ouzo drop by a fast solutal Marangoni flow. In this article, we investigate the evaporation of mixture droplets in more detail, by successively increasing the mixture complexity from pure water over a binary water-ethanol mixture to the ternary Ouzo mixture (water, ethanol and anise oil). In particular, axisymmetric and full three-dimensional finite element method simulations have been performed on these droplets to discuss thermal effects and the complicated flow in the droplet driven by an interplay of preferential evaporation, evaporative cooling and solutal and thermal Marangoni flow. By using image analysis techniques and micro-PIV measurements, we are able to compare the numerically predicted volume evolutions and velocity fields with experimental data. The Ouzo droplet is furthermore investigated by confocal microscopy. It is shown that the oil ring predominantly emerges due to coalescence.
△ Less
Submitted 21 June, 2017;
originally announced June 2017.
-
Self-wrapping of an ouzo drop induced by evaporation on a superamphiphobic surface
Authors:
Huanshu Tan,
Christian Diddens,
Michel Versluis,
Hans-Jürgen Butt,
Detlef Lohse,
Xuehua Zhang
Abstract:
Evaporation of multi-component drops is crucial to various technologies and has numerous potential applications because of its ubiquity in nature. Superamphiphobic surfaces, which are both superhydrophobic and superoleophobic, can give a low wettability not only for water drops but also for oil drops. In this paper, we experimentally, numerically and theoretically investigate the evaporation proce…
▽ More
Evaporation of multi-component drops is crucial to various technologies and has numerous potential applications because of its ubiquity in nature. Superamphiphobic surfaces, which are both superhydrophobic and superoleophobic, can give a low wettability not only for water drops but also for oil drops. In this paper, we experimentally, numerically and theoretically investigate the evaporation process of millimetric sessile ouzo drops (a transparent mixture of water, ethanol, and trans-anethole) with low wettability on a superamphiphobic surface. The evaporation-triggered ouzo effect, i.e. the spontaneous emulsification of oil microdroplets below a specific ethanol concentration, preferentially occurs at the apex of the drop due to the evaporation flux distribution and volatility difference between water and ethanol. This observation is also reproduced by numerical simulations. The volume decrease of the ouzo drop is characterized by two distinct slopes. The initial steep slope is dominantly caused by the evaporation of ethanol, followed by the slower evaporation of water. At later stages, thanks to Marangoni forces the oil wraps around the drop and an oil shell forms. We propose an approximate diffusion model for the drying characteristics, which predicts the evaporation of the drops in agreement with experiment and numerical simulation results. This work provides an advanced understanding of the evaporation process of ouzo (multi-component) drops.
△ Less
Submitted 5 June, 2017;
originally announced June 2017.
-
Evaporation-triggered microdroplet nucleation and the four life phases of an evaporating Ouzo drop
Authors:
Huanshu Tan,
Christian Diddens,
Pengyu Lv,
J. G. M. Kuerten,
Xuehua Zhang,
Detlef Lohse
Abstract:
Evaporating liquid droplets are omnipresent in nature and technology, such as in inkjet printing, coating, deposition of materials, medical diagnostics, agriculture, food industry, cosmetics, or spills of liquids. While the evaporation of pure liquids, liquids with dispersed particles, or even liquid mixtures has intensively been studied over the last two decades, the evaporation of ternary mixtur…
▽ More
Evaporating liquid droplets are omnipresent in nature and technology, such as in inkjet printing, coating, deposition of materials, medical diagnostics, agriculture, food industry, cosmetics, or spills of liquids. While the evaporation of pure liquids, liquids with dispersed particles, or even liquid mixtures has intensively been studied over the last two decades, the evaporation of ternary mixtures of liquids with different volatilities and mutual solubilities has not yet been explored. Here we show that the evaporation of such ternary mixtures can trigger a phase transition and the nucleation of microdroplets of one of the components of the mixture. As model system we pick a sessile Ouzo droplet (as known from daily life - a transparent mixture of water, ethanol, and anise oil) and reveal and theoretically explain its four life phases: In phase I, the spherical cap-shaped droplet remains transparent, while the more volatile ethanol is evaporating, preferentially at the rim of the drop due to the singularity there. This leads to a local ethanol concentration reduction and correspondingly to oil droplet nucleation there. This is the beginning of phase II, in which oil microdroplets quickly nucleate in the whole drop, leading to its milky color which typifies the so-called 'Ouzo-effect'. Once all ethanol has evaporated, the drop, which now has a characteristic non-spherical-cap shape, has become clear again, with a water drop sitting on an oil-ring (phase III), finalizing the phase inversion. Finally, in phase IV, also all water has evaporated, leaving behind a tiny spherical cap-shaped oil drop.
△ Less
Submitted 5 June, 2017;
originally announced June 2017.