Abstract
The flow in an evaporating glycerol-water binary submillimeter droplet with a Bond number Bo 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by microparticle image velocimetry 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.
| Original language | English |
|---|---|
| Article number | 114501 |
| Number of pages | 6 |
| Journal | Physical Review Letters |
| Volume | 122 |
| Issue number | 11 |
| DOIs | |
| Publication status | Published - 18 Mar 2019 |
Funding
Li Yaxing 1 Diddens Christian 1,2 Lv Pengyu 1 Wijshoff Herman 2,3 Versluis Michel 1 Lohse Detlef 1,4 ,* Physics of Fluids group, Department of Science and Technology, 1 Mesa+ Institute , Max Planck Center for Complex Fluid Dynamics and J. M. Burgers Centre for Fluid Dynamics, University of Twente , P.O. Box 217, 7500 AE Enschede, Netherlands Department of Mechanical Engineering, 2 Eindhoven University of Technology , P.O. Box 513, 5600 MB Eindhoven, Netherlands 3 Océ Technologies B.V. , P.O. Box 101, 5900 MA Venlo, Netherlands 4 Max Planck Institute for Dynamics and Self-Organization , 37077 Göttingen, Germany * [email protected] 18 March 2019 22 March 2019 122 11 114501 5 September 2018 © 2019 American Physical Society 2019 American Physical Society The flow in an evaporating glycerol-water binary submillimeter droplet with a Bond number Bo ≪ 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by microparticle image velocimetry 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. Nederlandse Organisatie voor Wetenschappelijk Onderzoek 10.13039/501100003246 Océ-Technologies Technische Universiteit Eindhoven 10.13039/501100003005 The evaporation of a microdroplet on a flat substrate has attracted a lot of attention because of its beautiful and phenomenologically rich fluid dynamics [1–14] and its relevance in various technological applications, such as medical diagnostics [15] and the fabrication of electronic devices [16] . For many of these applications, an understanding of the internal flow structure is crucial. One example is the so-called “coffee stain problem” [2] , i.e., an evaporating colloidal drop in which an outward capillary flow along the substrate carries the dispersed material from the interior towards the pinned contact line. This seminal study opened up a new line of research for surface coatings and patterning technologies, which is crucial for various applications in inkjet printing [17] , 3D printing technology [18] , and molecular biology [19] . However, in nearly all of these applications, the droplet liquid is not pure, but a binary or even ternary liquid. As is well known, then Marangoni flow, which is driven by surface tension gradients, is coming into play [6,20–23] , strongly affecting the evaporative behavior. The variation of the surface tension originates from two mechanisms or the combination of both, namely, a temperature gradient [20,21] or a solute concentration gradient [13,22–26] , due to the spatially varying local evaporation rates at the droplet surface. The conventional understanding is that the flows within submillimeter droplets can only be attributed to capillary and Marangoni convections, while natural convection is considered to be negligible as the surface tension force is dominant compared to gravity forces [21] . The Bond number of such a small-sized droplet system reads Bo = ρ g R 2 / γ ≪ 1 , which normally implies a gravity-independent system [27] . However, this only holds as a measure of the importance of the gravity force compared to the surface tension while relating to the droplet shape [28] . In recent years, studies of evaporating aqueous NaCl droplets revealed that natural convection can be driven by evaporation-induced density gradients inside a colloidal droplet [29–33] . However, to our best knowledge, an internal flow controlled by gravitational effects has never been observed, nor experimentally confirmed in a drying microdroplet consisting of two miscible liquid components. In this work, we investigate the flow field inside an evaporating glycerol-water binary miscible droplet. Glycerol is a very common liquid heavily used in industry [34] and laboratory experiments [35] , in particular due to its strong hygroscopic nature [36] . Recently, Shin et al. [37] found a spontaneous Bénard-Marangoni (BM) convection within the water-glycerol system: the hygroscopic absorption of water vapor at the interface generates thermal and solutal gradients, leading to a surface tension gradient, thus sustaining a BM instability. They confirmed that buoyancy is not the driving force by observing the same convection cells when orienting the system upside down. Here, however, we study the reverse process to absorption: the evaporation of a dilute aqueous glycerol droplet. In this system, glycerol can be considered nearly nonvolatile under room temperature as compared to water, implying a selective evaporation of the more volatile water, which leads to a concentration gradient of the components in the drying droplet. Such a concentration gradient generates surface tension and density gradients, which drives a convective flow in the droplet. Remarkably, the flow field is mainly controlled by the buoyancy force through the density gradient. We implemented microparticle image velocimetry ( μ PIV ) to measure the flow fields within the binary droplets with opposite configuration: a sessile and a pendant droplet. We deposited ( 0.18 ± 0.03 μ L ) binary droplets (glycerol:water, 50:50% by weight), seeded with fluorescent microspheres (520 nm diameter) at a concentration of 2 × 10 − 2 vol % on a OTS glass substrate (see Supplemental Material [38] ). Glycerol has a higher density ρ g = 1.261 × 10 3 kg / m 3 and lower surface tension γ g = 64 × 10 - 3 N m − 1 than water ( ρ w = 0.998 × 10 3 kg / m 3 , γ w = 75 × 10 - 3 N m − 1 ). During evaporation, the ambient temperature and relative humidity were stable, i.e., T = 21 ± 1 ° C and relative humidity = 50 ± 5 % . During the drying process the contact angles θ for both the sessile and the pendant droplet were between 90° to 104°. A typical droplet initially has a footprint radius R 0 = 400 μ m and a height h 0 = 490 μ m . We investigate the flow fields for both sessile and pendant droplets by adjusting the optical focus plane at different heights within the droplets: at 10 μ m (i) and at 200 μ m (ii) away from the substrate, as labeled in Fig. 1 . For sessile droplets, we observe an inward radial flow in the middle of the droplet (a1) and an outward radial flow close to the base of droplet (a2). Figure 1 presents the schematics of the flow velocity within the axis plane of the droplet from side view: the black arrows represent the measured flow velocities and the red arrows indicate the flow pattern reconstructed from the measurement. In contrast, Figs. 1 show a completely reversed flow fields for pendant droplets. The flow velocity near the substrate (b2) is radially outward, from center towards the contact line, which is opposite to the one of sessile droplet (b2). In the middle of the droplet (b1), it shows relatively weaker inward radial flow only in the outer region and asymmetric annulus flow near the edge of the droplet. Figure 1 (b3) illustrates the flow pattern in a pendant droplet. 1 10.1103/PhysRevLett.122.114501.f1 FIG. 1. The micro-PIV measurement of flow fields in both sessile (a) and pendant (b) droplets. The scale bars represent 200 μ m . (a1) and (a2) Flow fields in a sessile binary droplet measured at different heights: 200 and 10 μ m away from the substrate, respectively. The white arrows represent the flow direction. (a2) The measurement near the substrate shows an inward radial flow, (a2) the one at larger height reveals an outward radial flow. (b1) and (b2) Flow fields of a pendant binary droplet measured by the same method as the sessile droplets. (b2) The flow near the substrate follows outward radial direction, (b1) but the flow at midheight reveals an annular flow with deviations from axisymmetry near the edge and irregular flow within the inner part. The four PIV images were taken with four different droplets. (a3) and (b3) The schematics of the flow pattern in side views of both sessile and pendant droplets. The opposite radial flow directions near the substrate of differently orientated droplets clearly indicates that the gravitational effect is dominating to control the flow structure. In our system, the gravitational effect is driven by the density gradient in the bulk of the droplet, which in turn is generated from concentration gradients induced by the selective evaporation of the more volatile water. As the density of the glycerol-water mixture monotonically increases with increasing glycerol concentration, the local density decreases from the outer layer towards the inner bulk. For the sessile droplet, the denser glycerol-rich part is collecting on the top of the droplet, resulting in an unstable situation: the denser glycerol-rich part suspends atop the lighter water-rich part. As sketched in Fig. 1 , the lighter liquid part rises up due to the buoyancy, pushing the denser liquid to sink along the outer layer; hence an inward flow is passively driven from the contact line towards the center in the bottom layer. In contrast, for pendant droplets, the glycerol-rich part is at the bottom of the liquid bulk which cannot drive a buoyancy flow in the center. Instead, the denser liquid near the contact line flows down along the outer layer, and is replenished by the outward radial flow close to the substrate, implying a much weaker gravitational effect in the pendant droplet. Moreover, due to the intense coupling between composition and flow, instabilities can arise which lead to axial symmetry breaking near the interface [22,43] . We can estimate the flow velocity in the evaporating binary droplet by scaling arguments. The typical velocity from the μ PIV results is U ≈ 10 μ m s − 1 so that the Reynolds number Re = ρ U R / μ m ≈ 10 - 4 , where μ m ≈ 5 mPa s is the viscosity of the binary mixture. Furthermore, the Péclet number Pe = U R / D ≈ 10 , with the mutual diffusion coefficient D ≈ 0.4 × 10 - 9 m 2 / s , implying advection is dominating over diffusion. Assuming a quasisteady flow, the gravitational force due to the density difference is balanced by the viscous shear stress, which scales as g Δ ρ m ∼ μ ∇ 2 U . The density ρ m of the mixture liquid varies with the relative concentration ϕ following the Boussinesq approximation, ρ m = ρ 0 [ 1 + β ( ϕ - ϕ 0 ) - α ( T - T 0 ) ] [46] , where ρ 0 , ϕ 0 and T 0 denote, the density, concentration ratio and temperature of the reference state, respectively, and where α and β indicate the thermal and solutal expansion coefficient, respectively. Note that the evaporation of water cools down the liquid temperature near the surface, which enhances the density difference between the surface area and the bulk. For now we neglected the thermal effect and only considered the density changes due to concentration differences. This assumption is validated by numerical simulations (see below) with and without the thermal contribution with only minor changes to the flow field. Thus the density difference is given by Δ ρ m = ρ 0 β Δ ϕ , with the solutal expansion coefficient is β ≈ 0.2 . The evaporation rate is controlled by diffusion of water vapor molecules to the surrounding air. Therefore, the mass loss Δ m within a typical timescale Δ t ∼ R / U can be estimated from Δ m ∼ D w , air Δ C w , air R Δ t , with D w , air ≈ 10 - 5 m / s is the diffusion coefficient of water vapor at room temperature and Δ C w , air ≈ 10 - 2 kg / m 3 is the water vapor concentration difference between the air-liquid interface and the surrounding air. Hence, the concentration ratio difference Δ ϕ can be calculated from Δ ϕ ∼ Δ m / m 0 ∼ ( D w , air Δ C w , air R Δ t ) / ( 2 π R 3 ρ 0 / 3 ) . Therefore, we obtain the characteristic steady-state flow velocity U c as U c ∼ ( 3 D w , air Δ C w , air R 2 π μ m g β ) 1 / 2 ∼ 1 0 − 5 ms − 1 , (1) thus in very good agreement with the velocity estimate from experiment. Note that Eq. (1) is independent of the ratio of nonvolatile glycerol to water. We test the applicability of the scaling approach on various initial mass concentration of glycerol, ranging from 10% to 60%. The measured temporal evolutions of the averaged radial velocity near the substrate for sessile droplets are shown in Fig. 2(a) . We then rescale the experimental data by the characteristic velocity U c and timescale τ c . The latter is estimated by the diffusion lifetime of a sessile droplet [10] , τ c = ρ m R 2 D w , air Δ C w , air . (2) As shown in Fig. 2(b) , with the rescaling all data collapse, except the early stages of the droplets with 10% and 20% initial concentration. We conclude that the gravitational effect is enhanced by increasing the concentration of glycerol, before reaching the concentration at which glycerol starts to absorb the humidity from the air. 2 10.1103/PhysRevLett.122.114501.f2 FIG. 2. (a) The evolutions of the averaged radial velocity near the substrate for sessile droplets with various initial glycerol mass concentrations, 10% (blue), 20% (red), 30% (orange), 40% (purple), 50% (green), and 60% (cyan). (b) The dimensionless averaged radial velocity plotted against the dimensionless time for the same data as in (a), which are scaled by the characteristic velocity and lifetime according to Eqs. (1) and (2) , respectively. All the scaled data collapse into each other except the early stages of the ones with 10% and 20% initial concentration. We now introduce the Archimedes number Ar = g h 3 ρ 0 Δ ρ / μ m 2 [47] (which in this context can also be called Grashof number Gr) as the control parameter for this problem, where h is the height of droplets. For large Ar ≫ 1 gravity plays a prominent role, whereas for small Ar ≪ 1 , the gravity can be neglected. To test the applicability of Ar, we varied the height of a series of sessile glycerol-water droplets, with h g = 608 , 514, 320, and 154 μ m . For our system, with a density difference between glycerol and water Δ ρ g = 263 kg / m 3 , we obtain Ar = 23.1 , 14.0, 3.4, and 0.37. A prominent inward radial flow in the center of the droplets is observed for the cases when Ar is greater than 1, which indicates gravity-dominating flow. However, for the case when Ar < 1 , this flow disappears. More details are given in the Supplemental Material [38] . For comparison, we investigated the occurrence of gravitational effects in a 1,2-propanediol-water binary droplet in the same experimental setup. The density for 1,2-propanediol is 1.036 × 10 3 kg / m 3 at room temperature, which is slightly greater than that of water (density difference Δ ρ p = 36 kg / m 3 ) and μ p ≈ 7 mPa s for the viscosity [48] . As with glycerol, 1,2-propanediol is also nearly nonvolatile at room temperature conditions, resulting in a preferential evaporation of water. The contact angles θ p for sessile and pendant droplets are between 70° and 75° due to the low surface tension of 1,2-propanediol [49] , indicating a stronger Marangoni effect than that of the glycerol-water droplets. We measured the flow field in an optical focal plane 10 μ m above the substrate for oppositely configured droplets (sessile and pendant) with two different droplet heights, h p ≈ 800 and 410 μ m , with Ar ≈ 3.7 and 0.5, respectively, as shown in Fig. 3 . We observed the opposite radial outward flows in large droplets with different orientations, similar to glycerol-water droplets. But for small droplets, we find the same radially outwards flow direction for both droplets: there is no detectable flow within the center of the droplet, but only outward radial flow in the outer regime, which clearly indicates that for these droplets with low Ar, only the Marangoni effect dominates the flow near the edge. The two liquid systems clearly show that the Archimedes number is indeed the crucial control parameter for the occurrence of the gravitational effects. 3 10.1103/PhysRevLett.122.114501.f3 FIG. 3. Micro-PIV measurements of the flow field for both sessile and pendant 1,2-propanediol-water droplets with different heights: h g = 410 μ m (a1), 425 μ m (a2), 820 μ m (b1) and 800 μ m (b2). The flow fields were measured at same height: 10 μ m away from the substrate. (a1) and (a2) The small droplets have low Ar < 1 . The flows within two different orientated droplets show same radial direction, which indicates no gravitational effect. (b1) and (b2) The large droplets have Ar > 1 . The different radial flow directions clearly show gravitational effect. To validate the experimental results, corresponding finite element simulations were performed. To that end, an axisymmetric model [43] was employed, which had already been successfully validated with binary water-ethanol and ternary ouzo droplets [44] . The model solves the diffusion equation for the water vapor concentration C w , air in the gas phase, assuming vapor-liquid equilibrium according to Raoult’s law. It includes the activity coefficient [50] of the two components at the liquid-gas interface and considers the ambient vapor concentration far away from the droplet. The diffusive water vapor flux at the interface determines the volume evolution and the normal component of the velocity in the droplet and is furthermore used as an interfacial sink term in the convection-diffusion equation for the local composition in the droplet. The liquid properties, including the surface tension γ , mass density ρ , viscosity μ , diffusivity D , and thermodynamic activity of water a w , are not constant, but coupled to the composition field. The composition-dependent properties of the glycerol-water mixture have been extracted from experimental data and are plotted in Refs. [43,45] . While the gravitational body force has been neglected in the original model [43] , the discussed experimental results clearly indicate the relevance of this term. Hence, the generalized model used here solves the Navier-Stokes flow inside the droplet with the composition-dependent body force ρ g → in the bulk and the Marangoni shear stress at the liquid-gas interface. Figure 4 shows snapshots of the simulations for the sessile and the pendant droplet using the experimental parameters. It is apparent that, as in the micro-PIV results, close to the substrate the flow is directed inward (outward) for the sessile (pendant) droplet, whereas it is reversed in the bulk layer approximately at the center of the droplet. Along the liquid-gas interface the simulations reveal a counterrotating vortex as is apparent from the Stokes stream function depicted in the right half of Fig. 4 . For the pendant droplet, the flow at the interface is in the same direction as predicted by pure Marangoni flow in the absence of natural convection, i.e., from the apex towards the contact line. For the sessile droplet, however, the strong natural convection in the bulk transports water from the bulk to the apex. Thereby, despite of the enhanced evaporation rate at the top, there is a higher water concentration at the apex as compared to the contact line. This results in a Marangoni flow from the contact line towards the apex, i.e., in the opposite direction as predicted when considering the Marangoni effect without natural convection. Hence, the impact of the gravity is not only able to decisively control the flow direction in the bulk, but can also reverse the interfacial composition gradient and the corresponding Marangoni flow. 4 10.1103/PhysRevLett.122.114501.f4 FIG. 4. Snapshot of the numerical results for sessile (upper) and pendant (lower) droplets, respectively. In the gas phase, the water vapor concentration (color-coded) is shown along with the evaporation rate (arrows). Inside the droplet, the glycerol mass fraction (left) and the Stokes stream function (right) are depicted, where the arrows indicate the local flow direction. A quantitative comparison of the experimental and numerical results is shown in Fig. 5 . Here, the temporal evolution of the angular-averaged radial velocity of the micro-PIV measurement in the plane close to the substrate is compared with the numerically obtained radial velocity evaluated at the same height. It is apparent that the radial velocity profile is in good agreement for the entire drying time. For the case of the pendant droplet, however, the simulation predicts an inward flow region close to the contact line, which was not found in the corresponding micro-PIV measurement. This difference can presumably be attributed to the assumption of axisymmetry in the model, a condition that may easily be broken for the pendant droplet system [cf. Fig. 1 ]. 5 10.1103/PhysRevLett.122.114501.f5 FIG. 5. Experimental (a1),(b1) and numerical (a2),(b2) results for the evolution of the radial velocity near the substrate for both sessile (a1),(a2) and pendant (b1),(b2) droplets. Positive values indicates outward flow. For sessile droplets (a1),(a2), the numerical simulation shows a great agreement with experiment. However, for pendant droplet (b1),(b2), the numerical simulation shows a stronger flow near the contact line than in the experiment, which is presumably due to the axial symmetry breaking of the flow in experiment. Gravitational effects triggered by density gradients due to selective evaporation can play a dominating role in controlling the flow in evaporating multicomponent droplets, even at the submillimeter scale and small Bond numbers. Misled by the small Bo, hitherto, most studies on this subject until now have disregarded the influence of natural convection. Our results show conclusively that natural convection can readily dominate the flow for μ L droplets. Thus, our findings stimulate a careful treatment of the interplay of natural convection and Marangoni flow in multicomponent droplets in future studies. Furthermore, the possibility to reverse the bulk flow by overturning the system opens new application perspectives for surface coating and particle patterning. The authors thank A. Prosperetti, G. Lajoinie, T. Segers, H. Reinten, and M. van den Berg for the invaluable discussions, and the two anonymous referees for their very constructive comments. This work is part of the FIP Industrial Partnership Programme (IPP) of the Netherlands Organization for Scientific Research (NWO). This research programme is co-financed by Océ-Technologies B.V., University of Twente and Eindhoven University of Technology.
Fingerprint
Dive into the research topics of 'Gravitational effect in evaporating binary microdroplets'. Together they form a unique fingerprint.Cite this
- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver