Marvin
Detert
ab,
Yibo
Chen
a,
Harold J. W.
Zandvliet
b and
Detlef
Lohse
*ac
aPhysics of Fluids Group, Department of Science and Technology, Max Planck Center for Complex Fluid Dynamics and J. M. Burgers Centre for Fluid Dynamics, MESA+ Institute for Nanotechnology, University of Twente, 7500AE Enschede, The Netherlands. E-mail: d.lohse@utwente.nl
bPhysics of Interfaces and Nanomaterials, MESA+ Institute for Nanotechnology, University of Twente, 7500 AE Enschede, The Netherlands
cLaboratory for Fluid Physics, Pattern Formation and Biocomplexity, Max Planck Institute for Dynamics and Self-Organization, 37077 Göttingen, Germany
First published on 12th May 2022
Multi-component fluids with phase transitions show a plethora of fascinating phenomena with rich physics. Here we report on a transition in the growth mode of plasmonic bubbles in binary liquids. By employing high-speed imaging we reveal that the transition is from slow evaporative to fast convective growth and accompanied by a sudden increase in radius. The transition occurs as the three-phase contact line reaches the spinodal temperature of the more volatile component leading to massive, selective evaporation. This creates a strong solutal Marangoni flow along the bubble which marks the beginning of convective growth. We support this interpretation by simulations. After the transition the bubble starts to oscillate in position and in shape. Though different in magnitude the frequencies of both oscillations follow the same power law , which is characteristic of bubble shape oscillations, with the surface tension σ as the restoring force and the bubble's added mass as inertia. The transitions and the oscillations both induce a strong motion in the surrounding liquid, opening doors for various applications where local mixing is beneficial.
One example of such a system is a plasmonic bubble, which is created at immersed noble metal nanoparticles irradiated by a laser.7–13 In such a system light absorption via surface plasmon resonance leads to abundant heat causing the nucleation and growth of bubbles.1,6 These bubbles are initially almost solely composed of vapor, but eventually grow to be dominated by diffusion of dissolved gas from the surrounding liquid.14–17 Multiple recent studies elucidated a variety of fascinating phenomena occurring with plasmonic bubbles, such as a competition between solutal and thermal Marangoni forces leading to a bouncing bubble in the early growth stages,18 the formation of droplet plumes due to a Marangoni instability in ternary liquids19 or the deposition of gold nanoparticles from a suspension.20 Moreover, a multitude of potential applications of these bubbles has been demonstrated, e.g. in catalysis,21–23 medicine24–28 or microfluidics.29–32
This paper adds another surprising phenomenon to the rich phenomenology of physiochemical hydrodymanics involving plasmonic bubbles in binary liquids. Again, it builds on the preferential evaporation of one component of the binary liquid, which can lead to a bouncing bubble in the early growth stage via competing thermal and solutal Marangoni forces.18 The key idea here is to carry this on even further and to study what happens if we now take a drastic difference in the volatility of the two components of the binary liquids, much larger than we did before in Zeng et al. (2021).18 As we will see, this will lead to even stronger preferential evaporation and thereby larger concentration gradients, which gives rise to two new phenomena: one is a sudden, instantaneous transition in growth mode that is found for a variety of binary liquids. We study this experimentally via highspeed imaging and compare our results with simulations in order to identify the driving mechanism. The other phenomenon is a bubble oscillation right after the transition with a frequency that differs by an order of magnitude from that of the bouncing bubble. To rationalize our findings we propose a mechanism based on the competition between inertia and surface tension.
The paper is organised as follows: first, the experimental setup and the preparation of binary liquids are introduced in Section 2. Then, the transition in the growth mode is described in Section 3.1, followed by an detailed examination in Section 3.2. From this a mechanism is derived in Section 3.3 and supported by simulations. In Section 4 the bubble oscillations are discussed and compared to established theoretical predictions. Finally, the paper ends with conclusions and an outlook on further experiments and applications in Section 5.
Property | Benzyl alcohol | Cyclohexanone | 1-Heptanol | 1-Butanol | Acetone |
---|---|---|---|---|---|
Refractive index | 1.54 | 1.45 | 1.42 | 1.40 | 1.36 |
Boiling temperature (K) | 478.5 | 428.5 | 449.7 | 390.8 | 329.2 |
Density (kg m−3) | 1046 | 947 | 822 | 810 | 790 |
Thermal cond. (W (m K)−1) | 0.156 | 0.149 | 0.159 | 0.135 | 0.157 |
Heat capacity (J (mol K)−1) | 124.8 | 124.8 | 173.5 | 173.7 | 123.7 |
Latent heat of vap. (kJ mol−1) | 46.9 | 37.3 | 46.0 | 43.3 | 29.1 |
Surface tension (mN m−1) | 39.1 | 34.9 | 26.8 | 24.5 | 23.7 |
Viscosity (mPa s) | 5.93 | 2.26 | 6.92 | 2.94 | 0.32 |
Vapor pressure (kPa) | 0.008 | 0.40 | 0.02 | 0.64 | 24.50 |
Transition | Yes | Yes | Yes | No | — |
Initially, the bubble grows slowly (a1-a2), as expected, and a shadow forms close to the substrate (in Fig. 1(a2)). Similarly, in the bottom-view (in Fig. 1(b2)) the surrounding of the bubble shows a slight change in contrast. Moreover, a brighter plume above the bubble can be identified which corresponds to the bright structure in the bottom left of Fig. 1(b2). After 42 s, as shown in (a3), a sudden transition occurs and plume-like structures with rolls inside appear next to the bubble, while the shadowed region almost vanishes. The bottom-view in (b3) shows that this happens in every direction with three dominant plumes forming. Additionally, the transition is accompanied by a distinct increase in bubble size. At 43 s two plumes remain (see (b4)) and the rolls have developed further, while the depletion zone is completely gone (see (a4)). Then, the rolls and the bubble keep growing (see supplementary movie 1, ESI†) until the bubble reaches its maximum volume in (a5) and detaches. Afterwards, the whole process repeats itself. Additionally, after the transition, some bubble oscillations can be seen in supplementary movie 1 (ESI†). These oscillations will be discussed in detail in Section 4.
The contrast changes in the rolls and shadowed region only occur for mixtures and not in pure fluids. They are most pronounced if the difference in volatility and refractive index of the components of the binary liquid is large. Therefore, we interpret this as a change in the local refractive index due to a change in the local acetone concentration. Hence, the shadowed region in the beginning is considered to be a region of acetone depletion, whereas the plume above the bubble appears due to acetone enrichment caused by acetone vapor condensation. In addition, that means the rolls at the transition reflect the changes in the local acetone concentration and thereby indicate the onset of liquid motion along the rolls.
The evolution of the radius during the bubble growth is shown in Fig. 1(c) in yellow. The point at which the bubble undergoes a transition in growth mode (42 s) is clearly visible by a distinct increase in radius (red box). In movie Supplementary: Movie 2, ESI† and Fig. 1(d) the transition at τ = 41.23 s is shown in detail. Clearly, the jump in radius is almost instantaneous and the growth behavior changes drastically. Before the transition the growth of the radius is slow and dominated by the influx of dissolved air contained in the binary liquid as observed before for n-alkanes.14 However, after the transition the bubble grows rapidly until it detaches and the original growth mode is recovered for the new bubble. This, in conjunction with the formation of rolls around the bubble indicates that the bubble undergoes a transition from slow evaporative to convective growth. Prior studies from Zeng et al. (2021) support this as they have observed a strong flow near the bubble via particle image velocimetry measurements in the same system with an ethanol–water binary liquid.18
By systematically varying the mixing ratio of the binary liquid, the acetone fraction range in which the bubble undergoes a transition can be determined. Some exemplary cases are shown in Fig. 1(c). While the bubble undergoes a transition for intermediate acetone ratios (fA = 30–70%v) there is none in pure acetone (black), nor in pure benzyl alcohol (blue), nor in binary mixtures close to these pure liquids (fA = 90%v and fA = 10%v). Moreover, with increasing fA the bubble grows faster, irrespective of whether a transition takes place or not. This is readily explained by the fact that the bubble is mostly composed of air11,15,17 and the dissolved air entering the bubble during the evaporation of acetone is its main contributor. In comparison to water, where the air influx is purely diffusion driven,11 for the organic liquids the air solved in the fluid that enters when the fluid evaporates becomes more important. This is the case, because the solubility of air decrease in water with temperature, but increases in most organic liquids. Hence, the surrounding, heated fluid becomes slightly under-saturated in our case and hinders the growth via diffusion. Consequently, more available acetone implies that more acetone evaporates and more air enters the bubble, thereby accelerating the bubble growth.
To investigate whether this phenomenon is limited to a specific binary liquid, the experiments have been repeated for acetone mixed with either CHN or 1-heptanol or 1-butanol. While no transition could be observed for the case of 1-butanol, both CHN and 1-heptanol exhibit a transition and a similar behavior to benzyl alcohol. The results for these binary liquids are summarized in a phase diagram in Fig. 2. For all combinations shown, three different regimes can be identified. First, for low acetone ratios the bubble grows slowly and undergoes no transition (regime I, red). Then, for intermediate ratios the bubble undergoes a transition in growth mode (regime II, green). Finally, for high acetone ratios no transition occurs and the bubble grows fast and detaches (regime III, yellow). Regime II ranges from fA = 26–72%v for BnOH (a), fA = 30–67%v for 1-heptanol (b), fA = 30–50%v for CHN (c) and does not exist for 1-butanol. Clearly, the size of regime II depends on the fluid as it decreases from BnOH via 1-heptanol, and CHN to 1-butanol. Comparing the liquids properties, listed in Table 1, the boiling temperature (TB, right axis in Fig. 2) matches the decrease in regime II and hence we identify TB as the most probable determining parameter. In order to understand why TB should be relevant and what mechanism causes the bubble to undergo a transition, a more detailed examination of this phenomenon is required.
Another question which arises is whether the drastic change in bubble volume at the transition is caused by an influx of gas, evaporation or a combination of both? To find out, the bubble composition has been determined by studying the shrinkage of the bubble after exposure to a laser pulse, as done by Zaytsev et al.16 The method employed in that paper is based on the fact that the condensation of vapor is orders of magnitude faster than the dissolution of gas. Hence, the initial drop in bubble volume during its shrinkage reflects the vapor fraction. The remainder of the bubble content is gas (i.e. air). The result of this method is shown in Fig. 4(b) for CHN with fA = 43%v, where red, black and blue indicate a laser pulse length of 3 s, 4 s and 5 s, respectively. Without a transition the bubble consists mostly of gas with a vapor fraction fv,b ≈ 33%v, which is slightly lower for longer exposure (blue curve) as the bubble has more time to grow via diffusion. However, right after the transition the vapor fraction instantly jumps to fv,b = 50%v and then decreases within 2 s back to the pre-transition level of fv,b ≈ 33%v. So, the transition is accompanied by a drastic change in composition and the mechanism thus explains the initial burst of vapor production and bubble growth. With this we have collected enough hints to explain the mechanism.
Initially, the bubble grows fast and exceeds the area that is heated by the laser. Then, the bubble growth slows down and the footprint diameter D (t) remains constant. As D (t) stagnates, acetone replenishment only occurs via diffusion and hence the local acetone concentration decreases, thus reducing the evaporation rate. Since the evaporative cooling decreases, while the heat influx due to the continuous laser irradiation remains the same, the temperature at the contact line slowly increases. As long as the boiling temperature of the less volatile fluid is comparable to the spinodal temperature of acetone Ts (A) ≈ 0.9Tc = 457 K,35 it should not evaporate significantly. Thereby, it neither limits the temperature increase due to latent heat, nor contributes to the vapor fraction fv,b in the bubble. Once the temperature of the contact line reaches Ts (A), the latent heat of vaporization of acetone vanishes and almost all the remaining acetone in the mixture near the contact line will instantly evaporate without energy cost. This leads to a jump in bubble radius R (t) and also increases the vapor fraction fv,b inside the bubble as the additional volume is solely composed of vapor. At the same time, the rapid evaporation causes a sudden decrease in the local acetone concentration fA,l at the contact line. Hence, a strong gradient in surface tension to the ambient fA is created along the bubble, as fA,l ≪ fA, which leads to Marangoni flow. This flow ambient liquid along the bubble interface and marks the onset of streaming, as sketched in Fig. 5. The thermal Marangoni flow can be neglected, as the Marangoni numbers for the solutal and thermal Marangoni flow are and , i.e. Mas ≫ Mat. Here μ is the dynamic viscosity of the binary liquid, Δσc and ΔσT are the surface tension gradients caused by concentration and temperature and α and κ are the mass and thermal diffusivity. Once the flow sets in, it will supply fresh fluid with the ambient acetone concentration fA to the contact line, where acetone will preferentially evaporate. Consequently, a concentration gradient is maintained, allowing the streaming to persist. Besides that, the incoming fresh fluid will reduce the size of the depletion zone. As acetone contains more dissolved air than the less volatile fluids and evaporates easier, the growth of the bubble will accelerate once convection sets in and be similar to pure acetone (as observed in Fig. 1). As the growth is again dominated by the influx of dissolved air, the bubble composition shifts back towards the lower vapor fraction shortly after the transition.
Because this mechanism predicts a transition only when the boiling temperature Tb of the less volatile fluid is comparable to the spinodal temperature Ts (A) of the volatile one, it explains why no transition occurs for 1-butanol and why regime II increases with Tb. Besides this, it also explains why no transition occurs close to the pure fluids. For high fA the bubble grows too fast to reach Ts (A) at the contact line. On the other hand, for low fA the concentration gradient to the local fA,l is small and therefore the Marangoni flow becomes weak and not enough acetone can be resupplied to enhance the growth. Besides this, it provides a potential explanation to the surprising behavior of τ that seems to depend on the surface tension. It might be caused by the contact angle of the fluid, as it increases on the hydrophilic substrate from BnOH via CHN to HeOH, presumably due to their surface tensions. As the contact angle increases and approaches 90°, the heated liquid layer trapped between the bubble and the substrate shrinks and thereby reduces the distance between fA,l and fA. In turn, this could facilitate the onset of the solutal Marangoni flow and thereby decrease τ.
In order to verify this suggested mechanism, simulations for a Marangoni instability at a catalytic plane are performed (details of this simulation can be found in Li et al. (2021)19). The result of the simulations for Mas = 5 × 105 is shown in Fig. 6(a) and compared to the experimental data. Both cases show the bubble at t–τ = 0.5 s with Rτ = 300 μm. However, the simulations do not take into account bubble growth, so that the bubble is 5% larger in the experiment compared to the simulation with Re (t) = 315 μm and Rs (t) = 300 μm, respectively. For the simulations the local acetone concentration fA,l is shown via a color map, while the resulting flow is indicated by white arrows. The concentration gradients in the experimental data are indicated by different shades of gray. There is good qualitative agreement between the simulations and experiments as both show the same rolls and general structure. This further supports that the mechanism originates from a Marangoni instability.
The mechanism also predicts the same transition for other binary liquids as long as Tb of the less volatile fluid is comparable to Ts of the volatile one. One such combination is BnOH (Tb = 478.5 K) with ethanol (Ts (E) ≈ 464 K).35 Indeed, a transition could also be observed for this binary fluid as shown in supplementary movie 3 (ESI†) though no systematic analysis has been performed.
Summarizing this section, we conclude that the proposed mechanism is consistent with all our observations and explains why and when a bubble undergoes a transition in growth behavior.
In order to understand this phenomenon, we investigated the behavior on the much shorter timescale of milliseconds. The motion of the bubble is recorded in side-view and parameters such as the motion of its center-coordinates xc (t), yc (t), radius R (t) and the footprint diameter D (t) are extracted. Fig. 7 shows the bubble motion as deviation of the mean for an example of the oscillation (a) and the orbiting (b) from t0 = 7 s. The directions x and y are parallel and perpendicular to the substrate, respectively. In the case of an oscillating bubble, the bubble remains pinned to the laser spot (as xc (t) remains constant) and oscillates in shape as yc (t) shows a clear periodic signal. Interestingly, D (t) oscillates with the same frequency as yc (t), albeit with a phase shift of about 180°. In contrast, the orbiting motion is not accompanied by any other motion as only xc (t) oscillates. Additionally, the frequency of the oscillation and orbiting are clearly different. However, a bubble can switch from one kind of motion to another during its growth.
An example of this is shown in Fig. 8 for fA = 68%v, where the frequency during the whole life of the bubble is shown as a function of time. Closed symbols represent the main motion of the bubble, whereas open symbols also consider small amplitude oscillations. In both cases the bubble starts with shape oscillations yc (t) (black) which are accompanied by oscillations of the footprint diameter D (t) (red). Then, the motion switches to orbiting xc (t) (green) with a brief moment at t = 5.5 s where both motions are present and overlap. If only large amplitude motions are considered, there is no other motion present until shape oscillations appear again at t = 24 s and become dominant at t = 25 s. However, if small amplitude oscillations are not neglected oscillations in yc (t) and D (t) become visible during the orbiting. Surprisingly, those oscillations have about twice the frequency of xc (t) with . Additionally, they have a different frequency compared to the same motion with a large amplitude (closed symbols), as clearly visible at t = 24 s. Therefore, in the following analysis only the large amplitude oscillations will be considered. The reason why the bubble switches motion is not clear, however, one possible explanation could be that it is caused by an asymmetry. E.g., the local pinning could vary and lead to a small displacement of the bubble center from the laser spot. Then, the bubble would first move along a line and then undergo a transition into an orbiting motion. We have indeed observed such a behavior, though no detailed analysis has been performed.
Interestingly, all data for the oscillation and orbiting modes at different laser powers and concentration collapse on one curve, as shown in Fig. 9(a). Stars show the motion in D (t) if an oscillation (Dy (t)) or orbiting (Dx (t)) is detected. Similar to Fig. 7, the frequency of D (t) matches the oscillation, but not the orbiting. However, clearly both kinds of motion follow a power law with .
In order to understand this behavior we have another look at Fig. 5 where the situation is similar. In this case the bubble has grown rapidly and hence the temperature at the contact line is lower. Nevertheless, both thermal and solutal Marangoni flows are still present, since gradients in temperature and concentration remain (as sketched in Fig. 5). Similar to Section 3.3 it is reasonable to neglect the effect of the thermal Marangoni flow at the interface because Mas ≫ Mat. Then, the solutal Marangoni flow pushes the contact line inwards and in response the bubble elongates perpendicular to the substrate (along the y-axis). However, surface tension acts against this deformation and drives the bubble back to its spherical shape. Consequently, the contact line extends again and the whole process repeats itself, as for a spring-mass system. Basically, we have a small perturbation driven by solutal Marangoni flow and surface tension acts against it. Lamb determined the response frequency ω of a gas bubble in another inviscid fluid to a small perturbation of the interface using surface waves as36,37
Surprisingly, the orbiting motion also obeys the power law derived for shape oscillations. In order to understand this we have to take into account the previously neglected small amplitude oscillations (see Fig. 8). As their frequency is the second harmonic of the orbiting frequency and varies from the large amplitude oscillations, we assume they couple with the orbiting motion. We assume that they couple to the second harmonic of the dominant orbital motion, because it is the closest harmonic to their original frequency when shape oscillations are dominant. The small shape oscillations can be explained by the same mechanism used for the large amplitude oscillations, though the prefactor may change due to the coupling with the orbiting. Because of the coupling, the scaling of the small shape oscillation should transfer to the orbiting motion, too. Hence, it follows the same behavior as the large amplitude shape oscillations.
Both the transition and the oscillations cause motion in a volume more than an order of magnitude larger than the bubble and hence can be very beneficial for mixing, e.g. in catalysis or microfluidics. This is particularly relevant, because a wide range of fluids can be used as long as they satisfy that one is significantly more volatile (i.e. Tb ≲ Ts).
Further experiments that study the effect of other parameters, e.g. viscosity via silicon oils, are suggested. Another interesting aspect of the mechanism is that it does not require plasmonic heating, but only a localized heat source. Thus, we envision that our results could be reproduced with other heat sources, such as small resistive heaters, which removes the need of optical access to benefit from the enhanced mixing.
Footnote |
† Electronic supplementary information (ESI) available: Experimental videos of the bubble growth and transition as well as bubble oscillations in bottom- and side-view. See DOI: https://doi.org/10.1039/d2sm00315e |
This journal is © The Royal Society of Chemistry 2022 |