2.1 Key concept of the experiment

In both set-ups, a drop falls down to the laser-impact position while it relaxes to a spherical shape with radius $R_{0}$ (see figure 2). On its route the drop intercepts a horizontal light sheet that generates a synchronisation signal. This signal is used to trigger the impact of the drop by the main laser, the acquisition of the laser-pulse energy $E_{L}$ by an energy meter, as well as a beam profiler and two cameras for the visualisation. The complete arrangement of the synchronisation laser, photodiode and equipment for the drop generation can be moved in the $yz$-plane to adjust the drop trajectory relative to the laser focus. The delay between the trigger and the laser pulse is tuned to align the drop with the pulse. The pulse enters from the left through a focusing lens $\text{f}_{1}$, hits the drop at $x=y=z=0$ and exits to the right through the imaging lens $\text{f}_{2}$, which allows us to characterise the pulse and the drop irradiation (see § 2.2).

Figure 2.

(a) Side-view sketch of the drop-impact experiment at the moment of laser impact ($t=0$). The laser pulse is focused with a lens of effective focal length $\text{f}_{1}$, hits the drop and is redirected with an imaging lens $\text{f}_{2}$ onto a charge-coupled device (CCD) for its characterisation. The drop centre at the impact location defines the origin of our coordinate system, which is sketched in (b) from a back view ($\boldsymbol{e}_{z}$-direction). The experiment is repeated each time a new drop reaches $x=0$. The technical equipments of the set-ups are described in detail in Klein et al. (Reference Klein, Lohse, Versluis and Gelderblom2017) and Kurilovich et al. (Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016), respectively.

The response of the drop to the laser impact is observed from two orthogonal views: the side view, aligned with $\boldsymbol{e}_{y}$, and the back view, aligned with the pulse and drop propagation ($\boldsymbol{e}_{z}$), see figure 2(b). Stroboscopic image sequences are obtained by performing a new impact experiment and incrementing the time delay between the laser impact and the pulsed light source that illuminates the scene for each image. Image analysis yields the drop centre-of-mass position in all three coordinate directions as a function of time, which is used to calculate the velocity $U$ along $\boldsymbol{e}_{z}$. For $t>\unicode[STIX]{x1D70F}_{e}$ this velocity is constant (Klein et al. Reference Klein, Bouwhuis, Visser, Lhuissier, Sun, Snoeijer, Villermaux, Lohse and Gelderblom2015). The equivalent sheet radius $R$ is determined as the radius of the circle with the same projected area as the sheet (in the $xy$-plane). Experiments that suffer considerably from a laser-to-drop misalignment or variations in the laser energy are excluded of our analysis. We typically filter out the worst 10 % of all experimental realisations.

The technical equipments used for the MEK and tin experiments differ and are described in detail in Klein et al. (Reference Klein, Lohse, Versluis and Gelderblom2017) and Kurilovich et al. (Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016), respectively. In the current work, the backlighting in the tin set-up has been improved: a pulsed dye laser pumped by the second harmonic wavelength of a Nd:YAG laser emitting an approximately 5 ns pulse of 560 nm light with a spectral width of ${\sim}4~\text{nm}$ is used. This lighting reduces the detrimental effects due to temporal coherence, such as speckle, which enables the visualisation of small features of the expanding tin sheets.

2.2 Laser–matter interaction

The nature of the laser–matter interaction is a key difference between the two systems. As this interaction will turn out to be important for understanding the late-time fragmentation of the sheet (see § 6), we summarise the difference here, while more details can be found in Klein et al. (Reference Klein, Bouwhuis, Visser, Lhuissier, Sun, Snoeijer, Villermaux, Lohse and Gelderblom2015, Reference Klein, Lohse, Versluis and Gelderblom2017) and Kurilovich et al. (Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016).

In the MEK system the driving mechanism for the drop acceleration and deformation is a local boiling that is induced by the absorption of laser energy in a superficial layer of the drop. The thickness $\unicode[STIX]{x1D6FF}$ of this layer is determined by the amount of dye dissolved in the liquid and the absorption coefficient of the dye at the laser wavelength (Klein et al. Reference Klein, Lohse, Versluis and Gelderblom2017). The laser–dye combination is chosen such that $\unicode[STIX]{x1D6FF}/R_{0}\sim 10^{-2}\ll 1$, which is also the case for the opaque tin drops (Cisneros, Helman & Wagner Reference Cisneros, Helman and Wagner1982). On a time scale $\unicode[STIX]{x1D70F}_{e}\sim 10~\unicode[STIX]{x03BC}\text{s}$ this layer vaporises and is ejected at the thermal velocity $u$. On the same time scale, the resulting recoil pressure $p_{e}$ accelerates the remainder of the drop to the centre-of-mass velocity (Klein et al. Reference Klein, Bouwhuis, Visser, Lhuissier, Sun, Snoeijer, Villermaux, Lohse and Gelderblom2015)

(2.2)$$\begin{eqnarray}U\sim \frac{E_{\mathit{abs}}-E_{\mathit{th}}}{\unicode[STIX]{x1D70C}R_{0}^{3}\unicode[STIX]{x0394}H}u,\end{eqnarray}$$

where $E_{\mathit{abs}}$ is the energy absorbed by the drop, $E_{\mathit{th}}$ is the threshold energy that is needed to heat the liquid layer to the boiling point and $\unicode[STIX]{x0394}H$ is the latent heat of vaporisation. The scaling law (2.2) motivates our choice to use the solvent MEK for the current study. The low value of $\unicode[STIX]{x0394}H$ results in large drop velocities for a given laser energy, which translates into a large range of accessible Weber numbers.

For the tin drops the local fluence of the laser exceeds the ionisation threshold. A plasma forms within a fraction of the laser-pulse duration $\unicode[STIX]{x1D70F}_{p}=10~\text{ns}$, after which inverse-bremsstrahlung absorption strongly decreases the initially high reflectivity of the metallic surface to negligible values (Kurilovich et al. Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016). Any further laser radiation is absorbed by the plasma cloud. The expanding plasma exerts a pressure $p_{e}$ on the drop surface that accelerates the drop. The time scale of this acceleration is set by the plasma dynamics, which is of the same order as the laser-pulse duration, i.e. $\unicode[STIX]{x1D70F}_{e}\sim \unicode[STIX]{x1D70F}_{p}=10~\text{ns}$. Hence, as for the vapour-driven MEK drops, the tin drops are propelled by a short recoil pressure $p_{e}$. Similarly, the centre-of-mass velocity $U$ for tin scales with the absorbed energy, that is $U\sim (E_{\mathit{abs}}-E_{\mathit{th}})^{0.59}$, where $E_{\mathit{abs}}$, $E_{\mathit{th}}$ and the exponent now have their origin in the plasma dynamics (Kurilovich et al. Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016).

To obtain the local laser fluence experienced by the drops, we characterise the laser beam in each system in the absence of the drop using the lens $\text{f}_{2}$ that images the incident fluence $F_{\mathit{inc}}$ in the impact plane (figure 2). First, the total radiative energy $E_{L}$ of the pulse is measured with an energy meter capturing the whole beam of light. Second, a CCD records the relative fluence $f(x,y,z=0)$, which is translated into absolute terms using

(2.3)$$\begin{eqnarray}\displaystyle F_{\mathit{inc}}=F(x,y,z=0)=\frac{f(x,y,z=0)}{\displaystyle \int f(x,y,z=0)\,\text{d}x\,\text{d}y}E_{L}. & & \displaystyle\end{eqnarray}$$

Using the position of the drop on impact obtained with the same CCD, we then compute the fluence $F_{\mathit{abs}}$ that is actually absorbed by the drop, as shown in figure 3(b). From the same arguments underlying (2.2), the local recoil pressure $p_{e}$ on the drop surface is expected to follow the spatial variations in $F_{\mathit{abs}}$ according to

(2.4)$$\begin{eqnarray}p_{e}(r,\unicode[STIX]{x1D719})\sim \frac{F_{\mathit{abs}}(r,\unicode[STIX]{x1D719})-F_{\mathit{th}}}{\unicode[STIX]{x0394}H}\frac{u}{\unicode[STIX]{x1D70F}_{e}}.\end{eqnarray}$$

Figure 3.

(a) Planar laser-beam profile for the MEK system as recorded without a drop ($y/R_{0}\leqslant 0$) and with a drop (for $y/R_{0}\geqslant 0$). The latter yields the drop radius $R_{0}$ and position in the beam profile as indicated by the red solid line. The quantity $F_{\mathit{inc}}$ is the average fluence incident on the drop as given by (2.3). (b) Fluence $F_{\mathit{abs}}$ absorbed by the drop considering the losses due to Fresnel reflection at the liquid–air interface (Hecht Reference Hecht2002). (c,d) Laser profile (red solid line) in radial (c, azimuthally averaged) and azimuthal directions (d, radially averaged) obtained from ${\sim}100$ recordings of the planar profile. The black solid line indicates a perfect flat-top beam profile (denoted as $F_{\mathit{abs,FT}}$ in d). (e) Planar laser-beam profile measured for the tin system. The red solid line indicates the drop location on impact. The colour bar is the same as in (a), which illustrates the smoother and more uniform irradiation of the drop compared to the MEK case. (f) Radial beam profile obtained from (e).

Given the spatial variation in fluence observed in figure 3(d) this suggests that the MEK drops are subject to a driving force that varies along the azimuthal direction $\unicode[STIX]{x1D719}$ by about $\pm 10\,\%$. Importantly, since $f$ is found to be independent of $E_{L}$, these spatial variations in the driving force are independent of $E_{L}$ and fixed in the laboratory frame.

Figure 4.

Sequence of events following the laser-pulse impact on a MEK drop for $\mathit{We}=330$. Images are recorded stroboscopically (i.e. on different drops) from side and back views. The former are shown in a frame co-moving with the propulsion speed $U$. At $t=0.27~\text{ms}$, the drop has deformed into a semi-transparent sheet with radius $R(t)$ and non-uniform thickness $h(r,\unicode[STIX]{x1D719},t)$ that is bordered by a rim. The pointers in the three subsequent pictures indicate the onset of fragmentation of the sheet. First, rim breakup occurs by the radial expulsion of ligaments (at $t=0.54~\text{ms}$) that subsequently destabilise. Second, corrugations of the sheet appear that finally pierce holes. This sheet breakup occurs close to the rim, leading to neck breakup at $t=1.1~\text{ms}$, and close to the centre of the sheet leading to centre breakup at $t=1.7~\text{ms}$. A final web of ligaments is shown for $t=2.5~\text{ms}$.

By contrast, the tin drops experience a smooth and highly symmetric driving force. The lens $\text{f}_{1}$ (with a focal length of 1 m) forms a Gaussian beam aligned with the drop with a diffraction-limited waist $\unicode[STIX]{x1D714}_{0}\sim \unicode[STIX]{x1D706}_{L}\,\text{f}_{1}/d_{0}$, where $\unicode[STIX]{x1D706}_{L}=1064~\text{nm}$ is the wavelength and $d_{0}$ the beam diameter before lens $\text{f}_{1}$ (Hecht Reference Hecht2002). In our optical arrangement $\unicode[STIX]{x1D714}_{0}\approx 100~\unicode[STIX]{x03BC}\text{m}$ is much larger than the drop size $R_{0}=24~\unicode[STIX]{x03BC}\text{m}$, which results in a homogeneous irradiation of each drop (see figure 3e,f). Moreover, the tin drops are shielded from direct laser illumination by their own plasma cloud, which smooths all spatial fluctuations in the laser fluence on scales smaller than $R_{0}$. As a consequence, the deforming tin drops obey a high degree of rotational symmetry, as we will see in § 3.

3.1 Sequence of events for MEK drops

The MEK experiment in figure 4 illustrates the response of a drop to the laser impact. First, the drop accelerates on the time scale $\unicode[STIX]{x1D70F}_{e}\sim 10~\unicode[STIX]{x03BC}\text{s}$ after which it moves in the $\boldsymbol{e}_{z}$-direction with a velocity $U$ while it expands radially. At $t=0.27~\text{ms}$, which is close to the inertial time $\unicode[STIX]{x1D70F}_{i}=R_{0}/U=0.28~\text{ms}$, the drop already resembles a thin sheet. The semi-transparent liquid reveals a thinner outer region of the sheet that is bordered by a thicker and hence darker rim. Likewise, the centre of the sheet is thick compared to the outer region. As the sheet further expands, its thickness decreases, as shown by the brightening of the sheet from $t=0.54$ to 1.7 ms. The spatial variations of the grey level indicate that the thickness also varies in space. However, in spite of these modulations, the sheet preserves a near-circular shape during the expansion.

While it expands, the sheet destabilises and fragments. Two types of breakup can be identified in figure 4. First, the breakup of the bordering rim: tiny ($\ll R$) corrugations are visible on the rim at $t=0.27~\text{ms}$ and grow over time to form ligaments (observed for the first time at $t=0.54~\text{ms}$, see pointer), which are expelled radially outward. These ligaments break into droplets that continue to move outward at a constant speed comparable to the rim velocity ${\dot{R}}$ at the moment of detachment. As a result of this rim breakup at $t=1.1~\text{ms}$, the sheet is surrounded by a cloud of tiny drops.

Second, sheet breakup occurs through the nucleation of holes. Corrugations on the sheet are visible at $t=1.1~\text{ms}$ (a pointer at the top highlights a patch with high spatial frequency components). We observe that such disturbances on the sheet precede any hole nucleation, including events with multiple holes piercing a single patch of corrugations. Figure 4 shows two cases where a single hole nucleates in a corrugated region. At 1.1 ms the lower pointer marks a hole shortly after it has pierced the sheet close to the outer rim ($r/R\sim 1$), which we term neck breakup. At 1.7 ms the same process is captured in the centre of the sheet ($r/R<0.5$, centre breakup). Once a hole nucleates on the sheet it continues to grow, thereby collecting the surrounding liquid mass into ligaments. The last frame at $t=2.5~\text{ms}$ in figure 4 shows the result of multiple holes growing and eventually merging over time. The liquid of the sheet is finally collected in a (quasi) two-dimensional structure of ligaments that breaks into droplets.

3.2 Comparison of MEK and tin drops

A comparison of the fragmentation in the MEK and tin systems is presented in figure 5. The first row (a,d) shows rim breakup for an unpierced sheet at low Weber number. In both systems ligaments are expelled and break into droplets. In the tin sheet, the rim itself cannot be observed directly because of the tin opacity at the chosen wavelength for visualisation (Cisneros et al. Reference Cisneros, Helman and Wagner1982).

Figure 5.

Fragmentation regimes for the vapour-driven MEK drops (a–c, $R_{0}=0.9~\text{mm}$) and plasma-driven tin drops (d–f, $R_{0}=24~\unicode[STIX]{x03BC}\text{m}$). In both systems drop fragmentation initiates at three distinct locations: the bordering rim (a,d), the neck (b,e) and the centre of the sheet (c,f). The neck and centre breakup are not consecutive processes, especially the neck breakup can occur multiple times. The apparent elliptical shape of the tin sheets is caused by the weak parallax angle of the camera relative to the propulsion direction ($\boldsymbol{e}_{z}$) and is corrected for in image analysis. The white glow in (e,f) to the left of the sheet centre is an artefact of the plasma that propels the tin drops.

While rim breakup is observed for MEK and tin at comparable Weber numbers, more than one order of magnitude in $\mathit{We}$ separates the sheet breakup for the two systems (figure 5b,c versus e,f). However, the qualitative features of the sheet breakup are similar. In both systems the sheet breaks by the nucleation of holes in two distinct regions: neck breakup (b,e) and centre breakup (c,f). Neck breakup occurs before centre breakup and may repeat several times during the sheet expansion.

The observation of the neck breakup requires a high spatial and temporal resolution. The process is strongly localised in space and difficult to separate from other breakup events. Indeed, once growing holes reach the outer rim of the sheet, the rim detaches and breaks up, leaving no other trace behind than a new corrugated rim and tiny droplets. These detached drops contribute to the cloud of droplets surrounding the sheet from the rim breakup. In figure 5(c) for instance, neck breakup has already taken place.

By contrast, the growth of holes during the centre breakup is much easier to observe experimentally. In both MEK and tin sheets holes nucleate in the centre of the sheet, merge and collect mass in a web of ligaments that breaks up into droplets. The opaque tin sheets prevent a further comparison of the two systems in terms of the corrugations that are visible for MEK in figure 5(b,c).

3.3 Some comments on jetting

In addition to the rim and sheet breakups, one observes the ejection of mass on the opposite side of the laser impact in the form of a liquid crown (see figure 4). This ejected mass moves at a speed larger than $U$, collapses on the $\boldsymbol{e}_{z}$-axis ($t=0.54~\text{ms}$) to form a jet that detaches from the sheet and finally breaks up ($t=1.1$–2.5 ms). A similar jetting is observed in the tin system, as shown in figure 6(a,b).

Figure 6.

(a,b) Side-view images showing the formation of a jet in the centre of the drop in the MEK (a) and tin (b) systems. (c) Sheet contour obtained from a boundary integral simulation illustrating the cross-section of the axisymmetric shape for $\mathit{We}=790$ (adapted from Gelderblom et al. (Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016)). (d) Sketch of the sheet showing the bordering rim and the tampered neck and centre regions.

This early jetting is not a direct consequence of the pressure pulse driving the drop expansion. Boundary integral (BI) simulations of the drop-shape evolution after pressure-pulse impact (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016), which are capable to reproduce jetting phenomena in principle (Peters et al. Reference Peters, Tagawa, Oudalov, Sun, Prosperetti, Lohse and van der Meer2013), do not show this feature (see figure 6c).

Fast jetting often results from the implosion of a cavitation bubble (Crum Reference Crum1979; Ohl et al. Reference Ohl, Arora, Dijkink, Janve and Lohse2006; Thoroddsen et al. Reference Thoroddsen, Takehara, Etoh and Ohl2009; Utsunomiya et al. Reference Utsunomiya, Kajiwara, Nishiyama, Nagayama, Kubota and Nakahara2010; Tagawa et al. Reference Tagawa, Oudalov, Visser, Peters, van der Meer, Sun, Prosperetti and Lohse2012; Gonzalez Avila & Ohl Reference Gonzalez Avila and Ohl2016). In the opaque tin and MEK drops ($\unicode[STIX]{x1D6FF}/R_{0}\ll 1$) direct laser-induced cavitation is unlikely. However, pressure transients resulting from the ablation and thermoelastic effects (Sigrist & Kneubühl Reference Sigrist and Kneubühl1978; Wang & Xu Reference Wang and Xu2001; Vogel & Venugopalan Reference Vogel and Venugopalan2003; Masnavi et al. Reference Masnavi, Nakajima, Horioka, Araghy and Endo2011) and shock waves accompanying plasma generation (Clauer, Holbrook & Fairand Reference Clauer, Holbrook and Fairand1981; Marpaung et al. Reference Marpaung, Kurniawan, Tjia and Kagawa2001) travel through or may even focus inside the drop and induce potential cavitation spots (Reijers et al. Reference Reijers, Snoeijer and Gelderblom2017).

As the jet carries little mass, it has only a small effect on the overall response of the drop, and in particular on the late-time sheet dynamics. Therefore, a more detailed description of the jetting phenomenon is beyond the scope of the present study.

4.1 Model derivation

The description of the rim and sheet breakup requires a model for the deformation of the drop into an expanding sheet of radius $R$ and thickness $h$. Previous models have considered a sheet with uniform thickness (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016). However, from the MEK data it is clear that the sheet thickness has a radial dependency (see e.g. figure 4). Therefore, we employ here a slightly more sophisticated model that has previously been used for the sheet formed by an impact on a pillar (Villermaux & Bossa Reference Villermaux and Bossa2011)

(4.1)$$\begin{eqnarray}\frac{R(t)-R_{0}}{R_{0}}=\sqrt{3\mathit{We }_{d}}\frac{t}{\unicode[STIX]{x1D70F}_{c}}\left(1-\frac{\sqrt{3}}{2}\frac{t}{\unicode[STIX]{x1D70F}_{c}}\right)^{2},\end{eqnarray}$$

with

(4.2)$$\begin{eqnarray}\mathit{We}_{d}=\frac{E_{k,d}}{E_{k,cm}}\mathit{We},\end{eqnarray}$$

where $E_{k,d}/E_{k,cm}$ is the ratio of the deformation to the propulsion kinetic energies, which depends on the laser-beam profile (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016). The rescaled Weber number $\mathit{We}_{d}$ is only based on the fraction of the kinetic energy that is actually used for deformation. Its relation to $\mathit{We}$ accounts for the difference in impact conditions between the laser case and the pillar case, as derived in appendix A.

In the model by Villermaux & Bossa (Reference Villermaux and Bossa2011) the sheet thickness away from its axis has been described by $h(r,t)\sim R_{0}^{2}\mathit{We}_{d}^{-1/2}\unicode[STIX]{x1D70F}_{c}/(rt)$, which has been validated experimentally by Vernay, Ramos & Ligoure (Reference Vernay, Ramos and Ligoure2015). For the evolution of the sheet thickness in the centre region, which is required for the discussion on the sheet breakup in § 6, we use here a mass-averaged description, simply reflecting the conservation of mass:

(4.3)$$\begin{eqnarray}\frac{h}{R_{0}}\sim \left(\frac{R}{R_{0}}\right)^{-2}.\end{eqnarray}$$

The energy partition $E_{k,d}/E_{k,cm}$ differs between the MEK and tin cases. In the MEK system, the relative fluence $f$ in the impact plane is kept constant for all experiments and is directly related to the recoil pressure $p_{e}$ as expressed by (2.4). For the flat fluence profile observed experimentally, the energy partition can be obtained analytically (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016), which yields $E_{k,d}/E_{k,cm}=1.8$, independently of $E_{L}$.

By contrast, in the tin experiments we find that $E_{k,d}/E_{k,cm}$ follows a power-law dependence on $E_{L}$ (see figure 7). This power law expresses the fact that the plasma dynamics and hence the corresponding recoil pressure is a function of the incident laser energy, even at constant focusing conditions. A theoretical prediction of the plasma dynamics goes beyond the scope of this study. However, the trend with the laser energy can be explained qualitatively: a comparison of figure 7(b–d) shows that at lower laser energy the plasma cloud covers a smaller area of the drop surface, which results in an effective focusing of the recoil pressure to a confined region. A focussed pressure pulse in turn results in a larger $E_{k,d}/E_{k,cm}$ (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016). As a result, we expect $E_{k,d}/E_{k,cm}$ to increase with decreasing laser energy $E_{L}$, which is in agreement with the experimental observations in figure 7.

Figure 7.

(a) Energy partition as a function of laser energy for MEK (blue solid line) and tin drops (red square markers). For MEK $E_{k,d}/E_{k,cm}=1.8$, independently of $E_{L}$ as calculated analytically (Gelderblom et al. Reference Gelderblom, Lhuissier, Klein, Bouwhuis, Lohse, Villermaux and Snoeijer2016). The value for tin is determined for each experiment by the best fit of expression (4.1) to the experimental curves shown in figure 8. The black solid line is the power law $E_{k,d}/E_{k,cm}=0.19\,(E_{L}/E_{0})^{-0.27}$ with $E_{0}=1.0~\text{J}$ that follows from a linear regression. The three insets (b–d) show the white plasma clouds inducing the deformation of the tin drops (the initial undeformed tin drop is indicated in each inset by a red circle).

4.2 Comparison between model and experiments

The comparisons of Villermaux & Bossa’s analytical model (4.1) to experiments with both MEK and tin are shown in figure 8(a) and (b), respectively. When the experimental data are rescaled by the deformation Weber number $\mathit{We}_{d}$ (figure 8c) they all collapse onto (4.1). The model accurately captures the expansion up to the maximum radius $R_{max}$, the moment when $R_{max}$ is reached at $t_{max}=2\unicode[STIX]{x1D70F}_{c}/\sqrt{27}\approx 0.38\unicode[STIX]{x1D70F}_{c}$, and the recoil of the sheet due to surface tension. Especially for tin the agreement between model and experiment holds over nearly four decades in Weber number (figure 8b). For MEK (figure 8a) the deviation between the model and the experimental data is larger, in particular at higher Weber numbers ($\mathit{We}>170$). As we will discuss below, the model deviates from the experimental results when the fragmentation severely affects the topology of the sheet.

Figure 8.

Sheet-radius evolution as a function of time for MEK (a, circle markers in c) and tin drops (b, square markers in c). The black solid lines represent the model (4.1). The experimental curves are shown with a reduced marker density in (a,b) for better visualisation and the curves are stopped when the sheet evolution becomes too much affected by the fragmentation (i.e. when ligaments detach or holes in the sheet reach the rim). (c) Rescaled experimental data comparing all experiments of (a,b) (grey markers) with the analytical prediction (4.1). The highlighted cases and insets illustrate the influence of rim breakup (tin drop at $\mathit{We}=130$) and sheet breakup (MEK drop at $\mathit{We}=330$, tin drop at $\mathit{We}=2600$) on the apparent sheet expansion. In the absence of fragmentation (tin drop at $\mathit{We}=5$) the agreement between the model and the experiments is excellent.

Figure 9.

Evolution of the rim breakup for $\mathit{We}=132$ with the dimensionless time $t/\unicode[STIX]{x1D70F}_{c}$ obtained from a tin experiment exhibiting a highly symmetric expansion. (a) Total number of ligaments $N_{r}$. Each marker (●) indicates a new realisation of the experiment (with a delayed measurement) and the black dashed line (– – –) is a running average. The inset (b) shows the sheet radius $R(t)$, the amplitude $\unicode[STIX]{x1D709}$, and the wavenumber $k_{r}$ of the corrugation as observed at $t/\unicode[STIX]{x1D70F}_{c}=0.2$. During the recoil of the sheet ($t>t_{max}$) two or more ligaments may merge as shown in insets (c) and (d). (e) Sheet-radius evolution. Measurements (●) and model (4.1) (——).

In the collapsed view of figure 8(c) a few cases are highlighted to illustrate how fragmentation affects the comparison between model and experiment. In the absence of fragmentation the experimental data follow the model closely (e.g. for tin at $\mathit{We}=5$). At $\mathit{We}=130$ (▪, blue) the sheet is subject to rim breakup. The ligaments, which are expelled outward, do not follow the recoil and lead to an apparent over-expansion of the sheet for $t>t_{max}$ (see inset in figure 8c) since our image analysis for $R$ excludes detached ligaments but not those connected to the sheet. The same behaviour is observed for MEK at $\mathit{We}=90$ (●, purple) (see figure 8a). Interestingly, the effect of the rim breakup on the sheet dynamics decreases with increasing Weber number. For $\mathit{We}=960$ (▪, green) the apparent over-expansion during the recoil phase is much smaller (figure 8b), although rim breakup is observed in the experiments. Indeed, the sheet model (4.1) predicts the rim diameter $b$ and hence the mass contained by the rim to decrease with Weber number as $b/R_{0}\sim \mathit{We}_{d}^{-1/4}$ (Villermaux & Bossa Reference Villermaux and Bossa2011).

As the Weber number is further increased, sheet breakup in the neck region leads to a deviation between model and experiment, which is illustrated for MEK at $\mathit{We}=330$ (●, green) in figure 8(c). When holes nucleating in the neck region reach the outer rim, the latter partially detaches from the sheet and the measured radius $R$ decreases rapidly (see inset). This decrease in $R$ due to the neck breakup is also visible for tin sheets, e.g. for $\mathit{We}=2600$ (▪, purple) in figure 8(c). The onset of the sheet breakup occurs earlier for MEK than for tin, as we will show in § 6. Consequently, in figure 8(a) the MEK data deviate from relation (4.1) at earlier times than tin, especially for large Weber numbers where a severe neck breakup is observed.

5.1 Observations

A typical evolution of the rim breakup is illustrated in figure 9 for tin drops with $\mathit{We}=132$. Corrugations with an amplitude $\unicode[STIX]{x1D709}$ develop on the rim. Initially, these corrugations are visible in the experiments as mere noise. Later they form clear perturbations with a characteristic wavenumber $k_{r}$ from which ligaments evolve. We define the latter moment as the time $t_{r}$ of rim breakup, whereas the number $N_{r}$ of ligaments is obtained by counting.

Figure 9(a) shows that $N_{r}$ is initially constant but decreases for $t>t_{max}$ due to the compression of the rim during the recoil of the sheet. These ligaments, which are still attached to the sheet, get closer to each other and merge from their base, as shown in figure 9(c,d). The rim breakup time $t_{r}$ is plotted in figure 10 as a function of $\mathit{We}_{d}$. Ligaments form earlier for larger Weber numbers and always form before the sheet starts retracting ($t_{r}<t_{max}$). The maximum number of ligaments observed over $t\leqslant t_{max}$ is found to increase with increasing $\mathit{We}_{d}$ as illustrated in figure 11(a–d) with tin. This observation is confirmed by plotting $N_{r}$ versus $\mathit{We}_{d}$ in figure 11(e). For MEK drops neck breakup takes place much earlier than for tin and interacts with the formation of the rim ligaments. Therefore, neck breakup in MEK drops limits the range in $\mathit{We}_{d}$ for which reliable measurement of $N_{r}$ can be obtained. However, the two measurements we obtained are in quantitative agreement with the tin data at the same Weber number.

Figure 10.

Time $t_{r}$ when the rim corrugations become visible (see figure 9b) as a function of the Weber number $\mathit{We}_{d}$. The data are acquired manually from a subset of tin experiments that are recorded at identical camera and lighting settings to exclude any influence of the image resolution, for which we estimate the maximum error in $t_{r}/\unicode[STIX]{x1D70F}_{c}$ to be 0.02 as indicated by the error bars. The solid line is the scaling law (5.1) with a prefactor of 1.1.

Figure 11.

(a–d) Radial expulsion of ligaments during rim breakup for increasing Weber numbers (left to right). The back-view images are taken from experiments with tin drops exhibiting a highly symmetric expansion. When the depth of focus limits the detection of ligaments to a fraction $\unicode[STIX]{x0394}\unicode[STIX]{x1D719}/2\unicode[STIX]{x03C0}$ of the rim (see c) the total number of ligaments is estimated from $N_{r}=2\unicode[STIX]{x03C0}(\unicode[STIX]{x0394}N-1)/\unicode[STIX]{x0394}\unicode[STIX]{x1D719}$. (e) $N_{r}$ as a function of $\mathit{We}_{d}$ for tin (▪) and MEK drops (●). The data for MEK are limited to two experiments since the early hole nucleation in the neck region prevents an accurate measurement of the rim breakup for larger Weber numbers. The solid line is (5.2) with a prefactor of 4.4. The error bars indicate a relative error of 15 % based on figure 9.

5.2 Model derivation and comparison with experiments

Inspired by the similarity with the sheet dynamics following the impact on a pillar, we follow the approach of Villermaux & Bossa (Reference Villermaux and Bossa2011) to describe the rim breakup. We model the rim as a planar liquid cylinder of diameter $b\sim R_{0}\mathit{We}_{d}^{-1/4}$, which is justified since $k_{r}R\gg 1$ such that the curvature of the rim is negligible. The rim is subject to two destabilisation mechanisms. First, the Rayleigh–Plateau instability leads to a destabilisation of the rim on a time scale $(\unicode[STIX]{x1D70C}b^{3}/\unicode[STIX]{x1D6FE})^{1/2}\sim \unicode[STIX]{x1D70F}_{c}\mathit{We}_{d}^{-3/8}$ (Villermaux & Bossa Reference Villermaux and Bossa2011), which agrees with our experimental observation in figure 10. Second, the rim undergoes a time-dependent deceleration $-\ddot{R}(t)$, which induces a Rayleigh–Taylor instability with growth rate $\unicode[STIX]{x1D714}\sim (\unicode[STIX]{x1D70C}(-\ddot{R})^{3}/\unicode[STIX]{x1D6FE})^{1/4}$, because of the rim inertia.

For high Weber numbers and large rim decelerations reached in our experiments the instability is expected and found to develop at early times ($t_{r}\ll t_{max}$, see figure 10), in contrast to the experiments of Villermaux & Bossa (Reference Villermaux and Bossa2011). Using $t_{r}\ll \unicode[STIX]{x1D70F}_{c}$, the expansion of (4.1) into $\ddot{R}\sim -R_{0}\mathit{We}_{d}^{1/2}/\unicode[STIX]{x1D70F}_{c}^{2}$ gives the following time scale for the Rayleigh–Taylor instability

(5.1)$$\begin{eqnarray}t_{r}\sim \unicode[STIX]{x1D70F}_{c}\mathit{We}_{d}^{-3/8}.\end{eqnarray}$$

This time scale is identical to that of the Rayleigh–Plateau instability, as already observed for liquid sheet edges in a different context by Lhuissier & Villermaux (Reference Lhuissier and Villermaux2011). Figure 10 shows that (5.1) is in excellent agreement with the experimental data with a prefactor of 1.1. The scaling (5.1) differs from the breakup time ${\sim}\unicode[STIX]{x1D70F}_{c}$ proposed by Villermaux & Bossa (Reference Villermaux and Bossa2011) assuming that the stretching of the sheet delays the rim breakup.

The sheet radius at $t_{r}$ and the characteristic wavenumber $k_{r}$ at that time determine the number of ligaments according to $N_{r}\sim R(t_{r})k_{r}$. The fastest growing Rayleigh–Taylor mode is given by $k_{r}\sim (-\ddot{R}\unicode[STIX]{x1D70C}/\unicode[STIX]{x1D6FE})^{1/2}\sim \mathit{We}_{d}^{1/4}/R_{0}$, identical to the characteristic wavenumber of the Rayleigh–Plateau instability. Using again the early-time expansion of (4.1) we find $R/R_{0}\sim \mathit{We}_{d}^{1/2}t/\unicode[STIX]{x1D70F}_{c}$, which leads to

(5.2)$$\begin{eqnarray}N_{r}\sim R(t_{r})k_{r}\sim \mathit{We}_{d}^{3/8}.\end{eqnarray}$$

The derived scaling relations (5.1) and (5.2) are consistent with the work by Wang et al. (Reference Wang, Dandekar, Bustos, Poulain and Bourouiba2018), in which the same scaling for the rim diameter, breakup time and wavenumber of the fastest growing mode can be found. Figure 11(e) shows that (5.2), with a prefactor of 4.4, is in good agreement with the tin data. Although the $\mathit{We}_{d}$ dependence cannot be verified on the sole basis of the limited MEK data, the MEK data available are found to follow the scaling (5.2) with the same prefactor as the tin data. Hence, we conclude that the difference in rim breakup between MEK and tin is completely captured by the rescaled Weber number $\mathit{We}_{d}$ that accounts for the different driving mechanisms, in particular the effect of the plasma dynamics on the expansion of the tin sheets.

6.1 Observations

Figure 12.

Sheet breakup observed from a back view for MEK drops with increasing Weber numbers. (a) The sheet is smooth and starts to recoil from its maximum radius $R_{max}/R_{0}=6$ reached at $t/\unicode[STIX]{x1D70F}_{c}=2/\sqrt{27}\approx 0.38$, the moment the image is taken. Rim breakup leads to the formation of ligaments but breakup of the sheet itself is not observed. A slight increase in Weber number leads to a single piercing of the sheet (not shown). (b–e) The sheets are pierced near their neck and in the centre before $R_{max}$ is reached. The images are taken shortly after the first centre piercing event to allow for a characteristic hole density to develop. The resulting dimensionless time of each image ($t/\unicode[STIX]{x1D70F}_{c}=0.3$, 0.2, 0.15, 0.12) is decreasing with increasing Weber number. The shadowgraph visualisation with a small numerical aperture is sensitive to minute light refractions and reveals the sheet corrugations just before breakup. With increasing $\mathit{We}$ a larger hole density resulting in a finer web of ligaments is observed at the early moment of disintegration.

Figure 12 illustrates the sheet breakup for MEK drops over one decade of Weber number. For Weber numbers up to 170 (panel a) the sheet remains smooth and intact at all times, and only fragments due to rim breakup. For slightly higher Weber numbers, single sheet breakup events are observed, which are preceded by corrugations on the sheet surface (see also § 3). For $\mathit{We}=330$ (panel b) and higher (panels c–e) the sheet is more and more corrugated and ruptures both in the neck and the centre regions before it reaches its maximum expansion. The images in panels (b–e) are taken just after the first piercing event. They show that, with increasing Weber number, the sheet breakup becomes more severe. The number of holes $N_{s}$ that pierce the sheet per unit area and the corresponding wavenumber $k_{s}\sim N_{s}^{1/2}$ increase with increasing $\mathit{We}$. In addition, the time scale of the breakup becomes shorter as $\mathit{We}$ is increased ($t/\unicode[STIX]{x1D70F}_{c}=0.30$ and 0.12 in (c) and (e), respectively).

Figure 13.

Corrugations and hole nucleation on MEK sheets at $\mathit{We}=440$. (a) Close-up view of the sheet in (b) illustrating the result of the algorithm used to detect the corrugations that precede sheet breakup (the local corrugations (shown for $\unicode[STIX]{x1D719}<\unicode[STIX]{x03C0}$ as ○) with a spatial frequency $1/k_{corr}$ are identified by cross-correlation of the image with a circular Gaussian image kernel having a standard deviation $\unicode[STIX]{x1D70E}\sim 1/k_{corr}$). (b) Preferred regions for sheet breakup as identified by the analysis shown in (c) (neck:

, yellow, centre:

, orange) on top of a typical sheet observed in the experiments. (c) Probability density function (PDF) for the radial location $r/R$ of the sheet corrugations obtained from approximately 100 realisations of the experiment. The PDF is approximated by $\mathit{PDF}=2r/R\,g(r)$, where $g(r)$ is a radial modulation that describes the deviation of the hole nucleation location from a spatially uniform distribution. The experimental data (—●—) are well described by a two-component Gaussian mixture model $g(r|\unicode[STIX]{x1D707}_{i},\unicode[STIX]{x1D70E}_{i})$ (——) with $\unicode[STIX]{x1D707}_{i}$ and $\unicode[STIX]{x1D70E}_{i}$ being the mean and standard deviations of the radial location of hole nucleation. The detection algorithm fails close to the rim due to neck breakup and the receding rim, which explains the apparent deviation for $r/R>0.9$. The highlighted areas, i.e. $\unicode[STIX]{x1D707}_{i}-\unicode[STIX]{x1D70E}_{i}\leqslant r/R\leqslant \unicode[STIX]{x1D707}_{i}+\unicode[STIX]{x1D70E}_{i}$, illustrate the preferred hole locations in the centre (

, orange, $\unicode[STIX]{x1D707}=0.37$, $\unicode[STIX]{x1D70E}=0.13$) and the neck region (

, yellow, $\unicode[STIX]{x1D707}=0.96$, $\unicode[STIX]{x1D70E}=0.18$) of the sheet. (d) PDF of the azimuthal position $\unicode[STIX]{x1D719}$ of preferred hole locations for the centre (——, orange) and neck regions (——, yellow).

Hole nucleation in MEK is always preceded by corrugations with a high $k_{corr}$ on the sheet surface. However, no direct relation between $k_{corr}$ and $k_{s}$ is found. Only a few holes pierce a corrugated area, such that $k_{corr}\gg k_{s}$. The corrugations can, however, be used as an indicator for the areas where holes are likely to nucleate. We verified this concept with an image-analysis algorithm that is sensitive to spatial frequencies much larger than the hole density, as shown in figure 13(a). From the data of approximately 100 experimental realisations we obtain the probability density function (PDF) of hole nucleation in the radial direction (figure 13c). Not surprisingly, the quantitative analysis recovers a bimodal Gaussian distribution with two preferred areas for hole nucleation, as already identified visually in § 3: the neck and centre region, which are marked in figure 13(b). For each region the PDF in the azimuthal direction is shown in figure 13(d). Again, there is a clear deterministic influence. Three preferred areas of hole nucleation are observed in the centre region and approximately six in the neck region. More strikingly, the final web of ligaments preserves these deterministic influences. As shown in figure 14, the web formed for a single sheet with $\mathit{We}=2000$ shows the same pattern as the overlay of 31 realisations of the same experiment.

Figure 14.

(a) Back view of a fragmented MEK sheet at $\mathit{We}=2000$ and $t/\unicode[STIX]{x1D70F}_{c}=0.15$. The nucleation, growth and merger of holes on the sheet lead to a web of ligaments. (b) Image overlay of 31 MEK sheets from 31 different drops under the same experimental condition as in (a) and at the same time $t/\unicode[STIX]{x1D70F}_{c}=0.15\pm 0.006$. The grey scale is proportional to the probability that a ligament is present at a given position and a black pixel means that a ligament is present at that particular position in 100 % of the 31 experiments. This superposition reveals the highly deterministic nature of the final web of ligaments.

For the opaque tin drops potential corrugations on the sheet cannot be visualised. However, as already mentioned, deterministic influences can be found in the radial location of hole nucleation by visual inspection (see figure 5). In figure 15 we analyse the centre breakup of a tin sheet at $\mathit{We}=30\,000$. At this Weber number the hole density and radial extent $L_{c}$ of the centre region are such that $k_{s}L_{c}\gg 1$. Hence, we can sample a large number of holes to obtain unbiased statistics, i.e. unaffected by large-scale radial variations in the sheet thickness. The distribution of holes follows a linearly increasing PDF in radial direction (figure 15b) and uniform PDF in azimuthal direction (figure 15c), which express a uniform surface density in the centre region.

Figure 15.

Hole nucleation in the centre of tin sheets at $\mathit{We}=30\,000$. (a) Tin sheet with individual holes (○) as detected by an image-analysis algorithm sensitive to grey scale variations. The circle with radius $L_{c}=2.7R_{0}$ encloses 90 % of the hole nucleation events observed over ${\sim}100$ realisations of the experiment. (b,c) Radial (b) and azimuthal (c) distributions of nucleation events over $r\leqslant L_{c}$. The experimental distribution (

) is close to uniform, i.e. $\mathit{PDF}_{r}=2rR_{0}/L_{c}^{2}$ and $\mathit{PDF}_{\unicode[STIX]{x1D719}}=1/(2\unicode[STIX]{x03C0})$ (——). The wavenumber of hole nucleation, $k_{s}R_{0}=(N_{s}R_{0}^{2}/(\unicode[STIX]{x03C0}nL_{c}^{2}))^{1/2}=0.86$, is obtained from the total number of holes $N_{s}$ observed over $n$ experimental realisations.

6.3 Model derivation

We now derive a model for the Rayleigh–Taylor instability-driven sheet breakup to obtain a prediction for the characteristic breakup time $t_{s}$ and wavenumber $k_{s}$. To this end we modify the model for sheet breakup by Bremond & Villermaux (Reference Bremond and Villermaux2005) to account for the formation of the sheet from the spherical drop (see figure 16). In our analysis the local thickness variations of the centre and the neck (marked by $L_{c}$ and $L_{n}$ in figure 16d) that lead to the preferred areas of hole nucleation discussed in § 6.2.2 are neglected. Instead, we focus on the underlying mechanism of destabilisation. Consistently, the global scaling (4.3) is used for the overall kinematics of the sheet, thereby ignoring thickness fluctuations present in the sheet.

Figure 16.

Sketch of the three-phase model for the evolution of the impulsive Rayleigh–Taylor instability of the deforming drop. (a) Phase 1: the drop is accelerated perpendicular to its surface by the ablation pressure $p_{e}$ on time scale $\unicode[STIX]{x1D70F}_{e}$. This acceleration amplifies the Fourier modes of initial amplitude $\unicode[STIX]{x1D702}_{0}$ and wavenumber $k$. (b) Phase 2: for $\unicode[STIX]{x1D70F}_{e}<t\lessapprox \unicode[STIX]{x1D70F}_{t}$ the drop deforms into a sheet in the absence of any external acceleration. (c) Phase 3: the sheet expands radially until it breaks at time $t_{s}$ when the perturbation amplitude is of the order of the sheet thickness $h_{s}$. (d) Detail of the sheet-thickness profile (black solid line) and the perturbation with characteristic wavenumber $k_{s}$ (red dashed line) that causes hole nucleation. The solid red line marks the average sheet thickness $h_{s}$ at the moment of breakup. In two regions where the sheet is thinnest, in the neck (marked as $L_{n}$) and in the centre (marked as $L_{c}$), the criterion for breakup is fulfilled first and holes nucleate.

We model the drop as a uniform sheet of initial thickness $h_{0}\sim R_{0}$ and density $\unicode[STIX]{x1D70C}$ that is surrounded by a gas phase of negligible density. The laser impact induces an axial acceleration of the sheet given by

(6.1)$$\begin{eqnarray}a\approx \left\{\begin{array}{@{}ll@{}}{\displaystyle \frac{U}{\unicode[STIX]{x1D70F}_{e}}}=R_{0}/(\unicode[STIX]{x1D70F}_{c}\unicode[STIX]{x1D70F}_{e})\mathit{We}^{1/2}\quad & \text{for}~0\leqslant t\leqslant \unicode[STIX]{x1D70F}_{e},\\ 0\quad & \text{for}~t>\unicode[STIX]{x1D70F}_{e}.\end{array}\right.\end{eqnarray}$$

This acceleration amplifies any initial modulation of the surface, which can be represented by the Fourier modes (Bremond & Villermaux Reference Bremond and Villermaux2005)

(6.2)$$\begin{eqnarray}\unicode[STIX]{x1D702}(r,t)=\unicode[STIX]{x1D702}_{0}\,f(t)\text{e}^{\text{i}kr},\end{eqnarray}$$

with $k$ the wavenumber and $r$ a generalised coordinate system tangent to the sheet. The initial amplitude $\unicode[STIX]{x1D702}_{0}$, which can be as small as the thermal noise in the system (Eggers & Villermaux Reference Eggers and Villermaux2008), is assumed to be characteristic to each liquid system and independent of the wavenumber. The temporal evolution $f(t)$ follows from a potential flow analysis of the sheet and is given by $\ddot{f}(t)=-\unicode[STIX]{x1D714}^{2}f(t)$, with $\ddot{f}(t)=\text{d}^{2}f/\text{d}t^{2}$ and $\unicode[STIX]{x1D714}(k)$ the instantaneous growth rate (Keller & Kolodner Reference Keller and Kolodner1954; Bremond & Villermaux Reference Bremond and Villermaux2005).

We describe the evolution of the instability on the sheet in three consecutive phases, where we make use of the separation of time scales (2.1). In the first phase ($0\leqslant t\leqslant \unicode[STIX]{x1D70F}_{e}$, figure 16a) the drop is accelerated according to (6.1) and the modes are excited. In the second phase ($\unicode[STIX]{x1D70F}_{e}\leqslant t\lessapprox \unicode[STIX]{x1D70F}_{t}$, figure 16b) the acceleration is zero and the drop starts to deform. We define $\unicode[STIX]{x1D70F}_{t}$ as the time when the transition from a deforming drop to an expanding thin sheet takes place. In § 4 we observed that $\unicode[STIX]{x1D70F}_{t}\sim \mathit{We}_{d}^{-1/2}\unicode[STIX]{x1D70F}_{c}$, which implies that $\unicode[STIX]{x1D70F}_{t}\sim \unicode[STIX]{x1D70F}_{i}$. Even though during the second phase the drop no longer accelerates, the Fourier modes continue to evolve inertially, as they have acquired some velocity during the first phase. The third phase ($\unicode[STIX]{x1D70F}_{t}<t\leqslant t_{s}$, figure 16c) is characterised by a large radial expansion $R/R_{0}\gg 1$ of the sheet, which stretches the Fourier modes. During this phase the sheet gets pierced at a time $t_{s}$ when the amplitude of the evolving perturbations equals the sheet thickness (Bremond & Villermaux Reference Bremond and Villermaux2005). The sheet is not uniform in thickness but has thinner regions in the neck and centre, as illustrated in figure 16(d). Consequently, the perturbations can pierce the sheet in the neck and centre regions first, which thereby form preferred areas for hole nucleation.

In the following analysis, lengths and times are non-dimensionalised by the initial drop radius $R_{0}$ and capillary time $\unicode[STIX]{x1D70F}_{c}$,

(6.3a-c)$$\begin{eqnarray}\hat{\unicode[STIX]{x1D714}}=\unicode[STIX]{x1D714}\unicode[STIX]{x1D70F}_{c},\quad \hat{t}=\frac{t}{\unicode[STIX]{x1D70F}_{c}}\quad \text{and}\quad \hat{k}=kR_{0}.\end{eqnarray}$$

During phase 1 the capillary wavenumber $\hat{k}_{c}=\sqrt{\unicode[STIX]{x1D70C}aR_{0}^{2}/\unicode[STIX]{x1D6FE}}=\mathit{We}^{1/4}\hat{\unicode[STIX]{x1D70F}_{e}}^{-1/2}\gg 1$ and the sheet can be considered as thick with respect to the capillary length, which is set by the axial acceleration $a$ of the drop. The dispersion relation is then given by $\hat{\unicode[STIX]{x1D714}}_{1}^{2}=\hat{k}^{3}-\hat{k}_{c}^{2}\hat{k}$ (Bremond & Villermaux Reference Bremond and Villermaux2005). The modes of interest are the unstable ones that fit inside the sheet, i.e. $1\leqslant \hat{k}\leqslant \hat{k}_{c}$. As initial conditions for the shape function $f_{1}(t)$ in (6.2) we assume that all modes are initially excited at the same amplitude and zero initial velocity, such that

(6.4a,b)$$\begin{eqnarray}f_{1}(\hat{t}=0)=1\quad \text{and}\quad {\dot{f}}_{1}(\hat{t}=0)=0.\end{eqnarray}$$

Following Bremond & Villermaux (Reference Bremond and Villermaux2005), we treat the acceleration $a$ of the drop as impulsive, i.e. we assume $1/\hat{\unicode[STIX]{x1D714}}_{1}\gg \hat{\unicode[STIX]{x1D70F}_{e}}$. The sheet then behaves as a harmonic oscillator subject to an impulsive driving force, such that the shape function is given by

(6.5)$$\begin{eqnarray}\ddot{f}_{1}(t)=-\hat{k}^{3}f_{1}(t)+\frac{We^{1/2}}{\hat{\unicode[STIX]{x1D70F}}_{e}}\hat{k}.\end{eqnarray}$$

From (6.4), (6.5) we obtain

(6.6)$$\begin{eqnarray}f_{1}(\hat{t})=\cos (\hat{k}^{3/2}\hat{t})+\frac{\mathit{We}^{1/2}}{\hat{k}^{2}\hat{\unicode[STIX]{x1D70F}_{e}}}\{1-\cos (\hat{k}^{3/2}\hat{t})\}.\end{eqnarray}$$

To find the amplitude and growth rate of the modes at the end of phase 1, we again use the fact that the acceleration is impulsive and expand (6.6) for $\hat{\unicode[STIX]{x1D70F}}_{e}\rightarrow 0$ to obtain $f_{1}(\hat{\unicode[STIX]{x1D70F}_{e}})\approx 1$ and ${\dot{f}}_{1}(\hat{\unicode[STIX]{x1D70F}_{e}})\approx \hat{k}\mathit{We}^{1/2}$. Hence, by the end of phase 1 each mode has a specific growth rate that results from the impulsive acceleration while its amplitude is still equal to unity as it did not yet have time to grow.

In phase 2 the modes are no longer directly amplified by an acceleration ($\hat{k}_{c}=0$) but evolve inertially. The dispersion relation therefore simplifies to $\hat{\unicode[STIX]{x1D714}}_{2}^{2}=\hat{k}^{3}$. The initial conditions are obtained from a matching to phase 1 at $t=\unicode[STIX]{x1D70F}_{e}$. Again treating the acceleration as impulsive, i.e. letting $\hat{\unicode[STIX]{x1D70F}}_{e}\rightarrow 0$ (Bremond & Villermaux Reference Bremond and Villermaux2005), we find

(6.7a,b)$$\begin{eqnarray}f_{2}(\hat{t}=0)=1\quad \text{and}\quad {\dot{f}}_{2}(\hat{t}=0)=\hat{k}\mathit{We}^{1/2}.\end{eqnarray}$$

The shape function in phase 2 is then a free harmonic oscillator

(6.8)$$\begin{eqnarray}f_{2}(\hat{t})=\cos (\hat{k}^{3/2}\hat{t})+\frac{\mathit{We}^{1/2}}{\hat{k}^{1/2}}\sin (\hat{k}^{3/2}\hat{t}).\end{eqnarray}$$

As the drop expands into a thin sheet with $\hat{R}\gg 1$ we reach phase 3. The modes experience a stretch while at the same time the two interfaces of the sheet start to interact as their spacing ${\hat{h}}$ becomes of order $1/\hat{k}$, such that (6.8) is no longer valid. The mode development during phase 3 is described by the thin-sheet limit of the dispersion relation ($\hat{k}{\hat{h}}\ll 1$) in the absence of any acceleration, $\hat{\unicode[STIX]{x1D714}}_{3}={\hat{h}}\hat{k}^{4}/2$ (Bremond & Villermaux Reference Bremond and Villermaux2005). The expansion of the sheet causes a self-similar stretch of the modes. At the end of phase 2 (i.e. just before the stretch) each mode $k$ gets deformed according to $kR_{0}=k_{d}R$ (Villermaux & Bossa Reference Villermaux and Bossa2011), or $\hat{k}=\hat{k}_{d}\hat{R}$. Combining this self-similar stretch in wavenumber with the expression for the sheet thickness (4.3) we find for the instantaneous growth rate during phase 3: $\hat{\unicode[STIX]{x1D714}}_{3}\sim \hat{R}^{-6}\hat{k}^{4}$. Hence, while the sheet expands the growth (but also decay) rates of the modes decrease to zero much faster than the sheet thickness itself (recall (4.3)). As a consequence, the sheet expansion and simultaneous thinning freeze the exponential growth of the modes ($\ddot{f}_{3}(\hat{t})\approx 0$). The fastest growing mode $\hat{k}_{max}$ at the time $\hat{\unicode[STIX]{x1D70F}}_{t}$ of the transition from phase 2 to phase 3 therefore determines the shape function $f_{3}$ according to

(6.9)$$\begin{eqnarray}f_{3}(\hat{t})=f_{2,max}(\hat{\unicode[STIX]{x1D70F}}_{t})+{\dot{f}}_{2,max}(\hat{\unicode[STIX]{x1D70F}}_{t})\hat{t},\end{eqnarray}$$

where $f_{2,max}$ refers to the shape function (6.8) evaluated for $k=k_{max}$. As discussed above, the transition from phase 2 to phase 3 occurs when $\hat{\unicode[STIX]{x1D70F}}_{t}\sim \mathit{We}_{d}^{-1/2}$. The final result of our analysis is insensitive to the prefactor in this relation, which we take equal to unity.

To determine $\hat{k}_{max}(\hat{\unicode[STIX]{x1D70F}}_{t})$ and evaluate (6.9) we assume that the sheet expansion in phase 2 is fast in comparison with the oscillation period of $f_{2}$. Hence, at $\hat{\unicode[STIX]{x1D70F}}_{t}$ the sheet is thin while all Fourier modes are still in their first oscillation period. This condition requires $\hat{k}^{3/2}\mathit{We}_{d}^{-1/2}\ll 1$ (see (6.8)), which is justified for our experiments where $\mathit{We}_{d}\gg 1$. Therefore, we can expand (6.8) in the limit $\hat{k}^{3/2}\hat{t}\rightarrow 0$ to obtain

(6.10)$$\begin{eqnarray}\displaystyle & \displaystyle f_{2}(\hat{t})\approx 1+\mathit{We}^{1/2}\hat{k}\hat{t}-{\textstyle \frac{1}{2}}\hat{k}^{3}\hat{t}^{2}, & \displaystyle\end{eqnarray}$$

(6.11)$$\begin{eqnarray}\displaystyle & \displaystyle {\dot{f}}_{2}(\hat{t})\approx \mathit{We}^{1/2}\hat{k}-\hat{k}^{3}\hat{t}. & \displaystyle\end{eqnarray}$$

The fastest growing mode $\hat{k}_{max}$ at the end of phase 2 is then obtained from $\text{d}{\dot{f}}_{2}/\text{d}k=0$ and given by $\hat{k}_{max}=\mathit{We}^{1/4}/(3\hat{\unicode[STIX]{x1D70F}}_{t})^{1/2}$. The shape function in phase 3 for $\mathit{We}\gg 1$ then reads

(6.12)$$\begin{eqnarray}f_{3}(\hat{t})\sim \mathit{We}\left(\frac{E_{k,d}}{E_{k,cm}}\right)^{1/4}\hat{t}.\end{eqnarray}$$

The time $t_{s}$ of sheet breakup is reached when

(6.13)$$\begin{eqnarray}\hat{\unicode[STIX]{x1D702}}(t_{s})=\hat{\unicode[STIX]{x1D702}}_{0}\,f_{3}(\hat{t}_{s})={\hat{h}}(\hat{t}_{s}),\end{eqnarray}$$

i.e. when the corrugation amplitude equals the sheet thickness (Bremond & Villermaux Reference Bremond and Villermaux2005). Mass conservation dictates that the stretch in the wavelength due to sheet expansion is accompanied by a decrease in the corrugation amplitude $\hat{\unicode[STIX]{x1D702}}_{0}=\hat{\unicode[STIX]{x1D702}}_{d}/{\hat{h}}$, where $\hat{\unicode[STIX]{x1D702}}_{0}=\unicode[STIX]{x1D702}_{0}/R_{0}$. The breakup criterion (6.13) then reduces to $\hat{\unicode[STIX]{x1D702}}_{0}f_{3}(\hat{t}_{s})=1$ and the breakup time (expressed dimensionally for convenience) reads

(6.14)$$\begin{eqnarray}\frac{t_{s}}{\unicode[STIX]{x1D70F}_{c}}\sim \left(\frac{\unicode[STIX]{x1D702}_{0}}{R_{0}}\right)^{-1}\mathit{We}^{-1}\left(\frac{E_{k,d}}{E_{k,cm}}\right)^{-1/4}.\end{eqnarray}$$

Hence, our analysis predicts how the breakup time depends on the initial amplitude of the perturbation (large initial perturbation means early breakup), the Weber number, which measures the initial acceleration of the drop and the energy partitioning, which determines how fast the expanding sheet becomes thinner.

From the breakup time we can find the wavenumber $k_{s}$ at breakup, which sets the hole density. To this end, we use again the self-similar stretch $\hat{k}_{s}=\hat{k}_{\mathit{max}}/\hat{R}(\hat{t}_{s})$, with $\hat{R}(\hat{t}_{s})\sim \mathit{We}_{d}^{1/2}\hat{t}_{s}$ for $\hat{t}_{s}\ll 1$ and $\mathit{We}_{d}\gg 1$, to find (again dimensionally)

(6.15)$$\begin{eqnarray}k_{s}R_{0}\sim \frac{\unicode[STIX]{x1D702}_{0}}{R_{0}}\mathit{We}.\end{eqnarray}$$

Similar to the breakup time, we find that the characteristic wavenumber $k_{s}$ depends on the initial amplitude of the perturbation and the Weber number. However, the self-similar stretch of the wavenumber causes the dependence on the expansion dynamics, and hence the energy partitioning, to vanish.

6.4 Comparison between model and experiments

Figure 17 compares the scaling for the breakup time (6.14) with the experimental data of the centre breakup for the tin and MEK drops. As the kinetic-energy partition $E_{k,d}/E_{k,cm}$ differs between the tin and MEK experiments, the plot shows the breakup time $t_{s}$ rescaled by $E_{k,d}/E_{k,cm}$ to allow for a direct comparison.

Both the MEK and the tin data sets show good agreement with the predicted $\mathit{We}^{-1}$-scaling, over respectively one and two decades in Weber number. Strikingly, the absolute time at which the MEK and the tin sheets break differs by almost an order of magnitude, as was already observed in figure 5 in § 3.2. This difference can be explained by a difference in the initial noise level from which the instability grows. Assuming a prefactor of order unity in (6.14) we find for the noise level of the MEK sheets $\unicode[STIX]{x1D702}_{0}/R_{0}=1.2\times 10^{-2}$ and for tin $\unicode[STIX]{x1D702}_{0}/R_{0}=1.3\times 10^{-3}$. The noise level for MEK is hence much larger than the thermal noise, which is expected to be of nanometre scale (Eggers & Villermaux Reference Eggers and Villermaux2008). Unfortunately we are not able to determine the noise level from an independent experiment. However, we can qualitatively explain the difference between MEK and tin. As discussed in § 2.2, the MEK drops are subject to a much rougher beam profile and are furthermore propelled by vapour pockets bursting from their surface, while the tin drops interact with a smooth beam and plasma cloud. As a result, the initial noise in the MEK system is expected to be of macroscopic scale and much larger than for tin. These differences being accounted for by the parameter $\unicode[STIX]{x1D702}_{0}/R_{0}$, the two data sets obtained with different liquids at vastly different length scales both confirm scaling (6.14).

In figure 17 we also show the time of destabilisation in the neck region of MEK sheets. Again, the experimental data follow the scaling (6.14). Using the same noise level $\unicode[STIX]{x1D702}_{0}/R_{0}=1.2\times 10^{-2}$ we now obtain a prefactor of two. The different prefactor between the neck and centre region can be explained by the radial thickness profile $h(r,t)$. For a given Weber number the neck region reaches the critical thickness $h_{s}$ earlier than the centre region. However, the development of the corrugation amplitude is independent of the position on the sheet: a global Rayleigh–Taylor instability is responsible for the breakup in the neck and centre region alike. Consequently, the scaling exponents for the neck and centre breakup are identical and in agreement with our prediction.

Figure 17.

Onset time of centre breakup in MEK (●, blue) and tin drops (▪) as a function of $\mathit{We}$ ($t_{s}$ is defined as the time at which the first hole is observed over $r/R\leqslant 0.5$). The solid lines (——) are the prediction (6.14) with a prefactor unity, $\unicode[STIX]{x1D702}_{0}/R_{0}=1.2\times 10^{-2}$ for MEK and $\unicode[STIX]{x1D702}_{0}/R_{0}=1.3\times 10^{-3}$ for tin. For MEK the onset time for the neck breakup is also shown in grey (●). It follows the same scaling (6.14) with a prefactor of two instead of one (assuming the same noise $\unicode[STIX]{x1D702}_{0}/R_{0}=1.2\times 10^{-2}$).

Validation of the scaling for the characteristic wavenumber of breakup (6.15) requires experiments with sufficient holes in the centre region, i.e. $k_{s}L_{c}\gg 1$, to avoid influences of the global sheet-thickness fluctuations on the statistical analysis. These conditions are out of reach for the MEK drops, whereas for tin they can only be reached at very high Weber number (${\sim}10^{4}$). An example of such a case was shown in figure 15. As a consequence of these extreme conditions required for statistical analysis, we were unable to experimentally validate (6.15) for a broad range of Weber numbers. However, we find the order of magnitude of $k_{s}$ in figure 15 to be in agreement with (6.15), assuming a prefactor of order unity and using the same noise level as obtained from figure 17.

7.1 Phase diagram

After laser impact the drop goes through a series of stages, as described above. First, the drop expands radially according to (4.1). Then, at time $t_{r}$ given by (5.1), radial ligaments evolve from the sheet rim. Finally, holes nucleate on the sheet from time $t_{s}$ onward given by (6.14). The phase diagram in figure 18 summarises these different regimes as a function of the Weber number and the radial sheet expansion. The diagram is based on the scaling laws presented above with prefactors determined from the MEK experiments.

Figure 18.

Overview of the drop fragmentation regimes and radial expansion in terms of the Weber number. The parameter range of stable liquid sheets (

, blue) is separated from the unstable domain (

, orange) where rim breakup and sheet breakup occur. The radii $R(t_{r})$ and $R(t_{s})$ are determined from (4.1) and the scalings (5.1) and (6.14) for the breakup time in the rim (——, blue), neck (– – –) and centre region (——, orange). The prefactors are those obtained from the experimental MEK data in figures 10 and 17. The maximum radius $R_{max}$ (——, black) follows from the sheet kinematics (4.1) for $t_{max}/\unicode[STIX]{x1D70F}_{c}=2/\sqrt{27}$ with an energy partition for MEK of $E_{k,d}/E_{k,cm}=1.8$. The radius $R_{frag}$ (——, grey) is an estimate for the extent of the cloud of fragments (

) that originate from the rim breakup (see text).

The maximum sheet expansion $R_{max}/R_{0}$ that can be achieved follows from the sheet kinematics (4.1), as illustrated by the black solid line in figure 18. At low Weber number this sheet remains fully intact as the deceleration of the rim and the acceleration of the sheet are not strong enough to trigger breakup. As the Weber number increases the rim destabilises and radial ligaments form once the sheet has reached an expansion $R(t_{r})$ (blue solid line). The trajectory of fragments that originate from this rim breakup is set by the sheet expansion rate at the moment of detachment. The fragment position at the moment of maximum sheet expansion $t_{max}$ is therefore assumed to be given by $R_{frag}=R(t_{r})+{\dot{R}}(t_{r})(t_{max}-t_{r})$, as marked by the grey solid line. Sheet breakup in the neck and centre regions occurs from $R(t_{s})$ on, as marked by the red dashed and solid lines, respectively. Both breakup phenomena follow the same scaling law for the characteristic time of destabilisation $t_{s}$ but with a different prefactor, as discussed in § 6.4.

The shaded regions in figure 18 indicate the different fragmentation regimes. A fully intact sheet (blue zone,

$R/R_{0}\approx 3.4$

$\mathit{We}\approx 400$

$t_{max}$

$t>t_{max}$

The phase diagram presented in figure 18 is a practical tool to determine the radial mass distribution of the liquid after the first laser pulse. Such information is crucial for EUV lithography applications. To access the different regimes one can either vary the Weber number by tuning the laser-pulse energy or the radial expansion of the sheet by adjusting the timing of the main laser pulse.

7.2 Fragment sizes

Both the rim and sheet breakup give rise to a structure of elongated ligaments, which break up into droplets. The collection of all resulting fragments, which finally relax into a spherical shape, then leads to a size distribution of stable drops.

The elementary drop-size distribution coming from a single ligament breakup depends on the roughness and the mean diameter of the ligament and can be described by a gamma distribution (Villermaux Reference Villermaux2007). In our experiments at least five different sources of ligaments exist, as illustrated in figure 19. First, the rim gives rise to two types of ligaments: the radial ligaments that are expelled from the rim and the remnant of the rim itself that forms a thick circumferential ligament (Villermaux & Bossa Reference Villermaux and Bossa2011). Second, a web of ligaments results from holes opening on the sheet (Lhuissier & Villermaux Reference Lhuissier and Villermaux2013). As the thickness of the sheet formed after laser impact is far from uniform, one might expect at least two drop-size contributions originating from the neck and centre regions. Indeed, in figure 19 we observe that the mean ligament diameter varies considerably between the centre (panel b) and neck (panel c). Furthermore, the ligament diameter is also widely spread in each individual region (compare also figure 20a for an example of polydisperse ligament diameters). A final source of very small ligaments and drops originates from the collision of rims from neighbouring holes, as reported by Lhuissier & Villermaux (Reference Lhuissier and Villermaux2013). In our experiments this phenomenon is observed in particular in the neck region where the sheet is thinnest and therefore the hole-opening speed is largest. As illustrated in figure 19(c) rims may collide in an asymmetric fashion and form highly corrugated ligaments or even splash.

Figure 19.

Web of ligaments from a MEK drop at $\mathit{We}=750$ and $t/\unicode[STIX]{x1D70F}_{c}=0.27$ (a) and two close-up views from the centre (b) and neck region (c) of the sheet. The typical diameter of the ligaments $d_{\ell }$ varies considerably between the neck and centre as exemplified by the two pointers.

Because of the abovementioned effects, many simultaneous and sequential processes are at play in the general case, such that the analysis is de facto challenging. Therefore, we focus on a simpler case: the fragment sizes coming from a low Weber number impact where fragmentation occurs through the formation of radial ligaments. Figure 20 shows such an analysis for the MEK drops. To obtain a PDF of the fragments sizes, we analysed 200 separate MEK experiments. At times after the rim breakup had completed, we measured the fragment sizes $d$ using a large depth-of-focus setting for the imaging equipment in order to capture all rim fragments.

Figure 20(b) shows the PDF of $d/\langle d\rangle$ resulting from this analysis, with $\langle d\rangle =0.15~\text{mm}$ the mean fragment diameter. The PDF is cut at the combined resolution limit of our fragment detection algorithm and imaging system, which is then of the order of a few tens of $\unicode[STIX]{x03BC}\text{m}$. Clearly, the PDF is far from the expected bell-shaped, single gamma distribution, but exhibits, in particular, too many small drops. This broad, composite size distribution presumably results from two complicating factors. First, in the MEK drops azimuthal fluctuations in the sheet thickness resulting from the inhomogeneous laser-beam profile cause the radial ligaments to have a large spread in diameter; see figure 20(a). This distribution of ligaments sizes broadens the final drop-size distribution. Second, there is no clear separation between fragments from the rim and from the mist cloud that originates from the drop propulsion, as illustrated in figure 20(c). Hence, the many small fragments visible in the PDF do not come from the rim, but from the mist cloud. Unfortunately, we are not able to resolve the size distribution in the mist cloud separately, as we are limited by the resolution of our imaging system in combination with the high density of fragments in the mist cloud. Additionally, these fast tiny mist fragments may move out of focus, causing a bias in the analysis of the smallest fragments. The interference of the mist cloud with the rim fragments also affects the mean fragment diameter $\langle d\rangle /2R_{0}=0.088$, which is approximately a factor three smaller than one would expect for pure rim fragmentation, whose fragment size distribution does not present this small-size divergence (Villermaux & Bossa Reference Villermaux and Bossa2011).

Tin drops do not suffer from these complications. However, for these much smaller drops the analysis is severely limited by the optical resolution of the shadowgraph images. The small ligaments in figure 19(c) for the MEK system suggest a minimum ligament diameter of $d_{\ell }/R_{0}\sim 10^{-3}$. In absolute terms this size would translate to a few tens of nanometres for the tin system, i.e. far below the resolution of optical imaging in the visible spectrum.

The qualitative analysis presented here underlines the difficulty of establishing in a non-ambiguous manner a drop-size distribution arising from a non-trivial fragmentation process when the origin of each fragment cannot be traced back to a precisely identified intermediate mechanism. This is sometimes possible (Lhuissier & Villermaux Reference Lhuissier and Villermaux2013; Vledouts et al. Reference Vledouts, Quinard, Vandenberghe and Villermaux2016) and, when it is not, the analysis is often bound to invoke general principles to quantify an unspecified and therefore untraceable fragmentation mechanism (see e.g. He et al. (Reference He, Xin, Zhao, Chu, Xi, Shui, Lu and Gu2017) in the context of laser-pulse fragmentation, and the other examples discussed in § 7.2 of Vledouts et al. (Reference Vledouts, Quinard, Vandenberghe and Villermaux2016)), a pitfall we conscientiously avoid here.

Figure 20.

Fragmentation by rim breakup at $\mathit{We}=90$. (a) The polydispersity in ligament diameter is revealed by the two highlighted ligaments with vastly different diameters ($t/\unicode[STIX]{x1D70F}_{c}=0.6$). (b) Probability density function of the final fragment sizes normalised by the mean fragment diameter $\langle d\rangle$. The PDF is obtained on frames similar to (a) from 200 realisations of the experiment. The solid line is a gamma distribution of order 5 (Villermaux & Bossa Reference Villermaux and Bossa2011). (c) Initial mist cloud interfering with the remaining drop to generate a cloud of very small fragments ($t/\unicode[STIX]{x1D70F}_{c}=0.02$).