Hostname: page-component-5d59c44645-mrcq8 Total loading time: 0 Render date: 2024-02-28T19:39:51.635Z Has data issue: false hasContentIssue false

Laser-induced forward transfer of viscoplastic fluids

Published online by Cambridge University Press:  09 October 2019

Maziyar Jalaal*
Affiliation:
Physics of Fluids Group, Max Planck Center for Complex Fluid Dynamics, MESA+ Institute and J.M. Burgers Centre for Fluid Dynamics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
Martin Klein Schaarsberg
Affiliation:
Physics of Fluids Group, Max Planck Center for Complex Fluid Dynamics, MESA+ Institute and J.M. Burgers Centre for Fluid Dynamics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
Claas-Willem Visser
Affiliation:
Thermal and Fluid Engineering Department, University of Twente, Enschede, P.O. Box 217, 7500 AE Enschede, The Netherlands
Detlef Lohse
Affiliation:
Physics of Fluids Group, Max Planck Center for Complex Fluid Dynamics, MESA+ Institute and J.M. Burgers Centre for Fluid Dynamics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands Max Planck Institute for Dynamics and Self-Organisation, 37077 Göttingen, Germany
*
Email addresses for correspondence: m.jalaal@utwente.nl, mazi@alumni.ubc.ca

Abstract

Laser-induced forward transfer (LIFT) is a nozzle-free printing technology that can be used for two- and three-dimensional printing. In LIFT, a laser pulse creates an impulse inside a thin film of material that results in the formation of a liquid jet. We experimentally study LIFT of viscoplastic materials by visualizing the process of jetting with high-speed imaging. The shape of the jet depends on the laser energy, focal height, surface tension and material rheology. We theoretically identify the characteristic jetting velocity and how it depends on the control parameters, and define non-dimensional groups to classify the regimes of jetting. Based on the results, we propose the optimal conditions for printing with LIFT technology.

Type
JFM Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s) 2019

1 Introduction

Laser-induced forward transfer (LIFT) is a direct-write method of printing. LIFT exploits the energy of a laser (continuous or pulsed) to locally transfer material from one surface to another. The various materials, such as metals (Pohl et al. Reference Pohl, Visser, Römer, Lohse, Sun and Huis2015; Visser et al. Reference Visser, Pohl, Sun, Römer, Huis in ’t Veld and Lohse2015; Zenou, Sa’ar & Kotler Reference Zenou, Sa’ar and Kotler2015), pastes (Wang et al. Reference Wang, Auyeung, Kim, Charipar and Piqué2010; Munoz-Martin et al. Reference Munoz-Martin, Brasz, Chen, Morales, Arnold and Molpeceres2016) or liquids (Colina et al. Reference Colina, Duocastella, Fernández-Pradas, Serra and Morenza2006), are initially deposited on a transparent surface (the donor). The laser passes through this substrate and is absorbed by the material. Given large enough laser energy, a small part of the material is ejected and jetted towards the receiver substrate. Repeating the process, and moving the donor and the receiver, leads to two-dimensional (2-D) or three-dimensional (3-D) printing of a structure. A LIFT configuration can feature a metallic or polymeric dynamic release layer (DRL) between the transparent substrate and the material (Turkoz, Fardel & Arnold Reference Turkoz, Fardel and Arnold2018a ). A DRL absorbs the majority of the laser energy and provides momentum for the jetted material. LIFT is non-lithographic, nozzle-free and does not require a cleanroom for fabrication. For an overview of the research on LIFT, we refer the reader to the recent comprehensive review of Serra & Piqué (Reference Serra and Piqué2018) and references therein.

The majority of LIFT studies up to now have focused on Newtonian liquids. However, a wide range of applications are based on complex materials, such as hydrogels in 3-D bio-printers (Xiong et al. Reference Xiong, Zhang, Chai, Huang and Chrisey2015) and metallic pastes in electronic printing devices (Wang et al. Reference Wang, Auyeung, Kim, Charipar and Piqué2010; Munoz-Martin et al. Reference Munoz-Martin, Brasz, Chen, Morales, Arnold and Molpeceres2016). Very recently, Turkoz et al. (Reference Turkoz, Perazzo, Kim, Stone and Arnold2018b ) have studied LIFT of shear-thinning and viscoelastic jets. They used xanthan gum solutions of different concentrations in a LIFT system with a DRL and categorized the jet morphologies, identifying (modified) Deborah and Ohnesorge numbers as control parameters. They show that these non-dimensional groups span a phase space, with a no-droplet regime, a jetting regime and a fragmentation regime.

In the present work, we study LIFT of viscoplastic (VP) fluids. The main feature of such a material is the yield stress. Below this critical stress, the material behaves like an elastic solid. If the applied stress exceeds the yield stress, the material flows like a viscous fluid, typically featuring a nonlinear shear-dependent viscosity. For more information on VP fluids, see the reviews of Balmforth, Frigaard & Ovarlez (Reference Balmforth, Frigaard and Ovarlez2014) and Bonn et al. (Reference Bonn, Denn, Berthier, Divoux and Manneville2017).

The paper is organized as follows. In § 2, we provide information on the experimental conditions, including the materials and the set-up. In § 3, we present the phenomenology and report on our high-speed observations of the regimes of jetting. In § 4, we will then theoretically derive the jetting velocity and present a dimensional analysis, identifying the non-dimensional groups governing the process. Finally, § 5 concludes the results and gives an outlook.

2 Experiments

2.1 Materials

For the experiments, we use mixtures of commercial hair gel from Europrofit and milli-Q water. The material is mainly a pH-neutralized aqueous solution of Carbopol (with triethanolamine). The material also includes a small amount of additional polymers (polyvinylpyrrolidone (PVP) and propyleneglycol) and features elastoviscoplastic properties (see Dinkgreve et al. Reference Dinkgreve, Paredes, Denn and Bonn2016). We use an MCR 502 Anton-Paar rheometer with a cone-and-plate configuration to characterize the material properties. To avoid slip (see Roberts & Barnes (Reference Roberts and Barnes2001) and also Meeker, Bonnecaze & Cloitre (Reference Meeker, Bonnecaze and Cloitre2004), Jalaal, Balmforth & Stoeber (Reference Jalaal, Balmforth and Stoeber2015)), we use sand-blasted surfaces with a roughness length scale of $4.2\pm 0.3~\unicode[STIX]{x03BC}\text{m}$ . The roughness size is close to the characteristic length of the microstructures of the solutions, being the diameter of the soft blobs of the polymers (see e.g. Kim et al. Reference Kim, Song, Lee and Park2003; Géraud et al. Reference Géraud, Jørgensen, Ybert, Delanoë-Ayari and Barentin2017). We first perform oscillatory tests for a stress range of $0.03~\text{Pa}\lesssim \unicode[STIX]{x1D70F}\lesssim 100~\text{Pa}$ . Figure 1(a) shows the values of the elastic storage modulus $G^{\prime }$ and of the loss modulus $G^{\prime \prime }$ both versus the stress $\unicode[STIX]{x1D70F}$ at an operating frequency of 1 Hz. At low stresses, the material is predominantly elastic ( $G^{\prime }\gtrsim 10G^{\prime \prime }$ ). The elastic moduli, however, decrease significantly at high stresses (when stress exceeds the yield stress), while the loss modulus becomes larger than the storage modulus.

We perform shear-rate-controlled tests to measure the viscous properties above the yield stress. The protocol was as follows: After placing the sample on the plate, we first pre-shear by ramping up the shear rate ( $\dot{\unicode[STIX]{x1D6FE}}$ ) from $0.01~\text{s}^{-1}$ to $500~\text{s}^{-1}$ . We then decreased the shear rate for the same range. The data were recorded during this decreasing shear-rate ramp. Note that a subsequent increasing shear-rate sweep test resulted in the same reported flow curve. Each data point is an average value over 10 s. The changes in the flow curves were negligible for waiting times longer than 10 s – in contrast to some other yield-stress fluids that require much longer waiting time (see e.g. Jalaal et al. Reference Jalaal, Cottrell, Balmforth and Stoeber2017; Hopkins & de Bruyn Reference Hopkins and de Bruyn2019). Figure 1(b) shows the flow curves of the used material. The flow curves approach a plateau when $\dot{\unicode[STIX]{x1D6FE}}\rightarrow 0$ , i.e. the stress approaches the yield stress $\unicode[STIX]{x1D70F}_{0}$ and the apparent viscosity blows up. We quantify the properties of the materials using the following constitutive model:

(2.1) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D70F}=\unicode[STIX]{x1D6FE}G_{0}^{\prime }\quad \text{if}\quad \unicode[STIX]{x1D70F}<\unicode[STIX]{x1D70F}_{0},\quad \text{otherwise}\quad \unicode[STIX]{x1D70F}=\unicode[STIX]{x1D70F}_{0}+K\dot{\unicode[STIX]{x1D6FE}}^{n}, & & \displaystyle\end{eqnarray}$$

where we assume that the material is a simple linear elastic solid below the yield stress ( $G_{0}^{\prime }$ being the storage modulus at a small deformation limit) and a Herschel–Bulkley fluid above it. Using such a constitutive law, we ignore the viscous effects below the yield stress as well as elastic effects above the yield stress (viscoelasticity). We will later comment on the elastic effects in appendix B.

Figure 1. Rheology of the seven viscoplastic fluids analysed. Symbols and numbers refer to the different samples, whose material properties are given in table 1. (a) Storage modulus $G^{\prime }$ (symbols) and loss modulus $G^{\prime \prime }$ (lines) as functions of stress. (b) Flow curves: symbols are the experimental data and lines are the Herschel–Bulkley fits. (c) Viscous stress $\unicode[STIX]{x1D70F}-\unicode[STIX]{x1D70F}_{0}$ versus the shear rate (symbols) and their corresponding Herschel–Bulkley fits (lines).

In (2.1), $K$ and $n$ are the consistency and flow indices, respectively. We find these values by directly fitting the flow-curve data (solid lines in figure 1 b). To show the quality of the fits, we also plot the variation of the viscous stresses ( $\unicode[STIX]{x1D70F}-\unicode[STIX]{x1D70F}_{0}$ ) with $\dot{\unicode[STIX]{x1D6FE}}$ in figure 1(c) (see Katgert et al. Reference Katgert, Latka, Möbius and van Hecke2009; Lidon, Villa & Manneville Reference Lidon, Villa and Manneville2017). We note that our fits, although satisfactory for $\dot{\unicode[STIX]{x1D6FE}}>0.1~\text{s}^{-1}$ , include some error for small values of $\dot{\unicode[STIX]{x1D6FE}}$ . Alternatively, one can fit a linear curve to $\ln (\unicode[STIX]{x1D70F}-\unicode[STIX]{x1D70F}_{0})$ versus $\ln (\dot{\unicode[STIX]{x1D6FE}})$ with initially guessing and then searching for a yield stress that minimizes the error residuals (see appendix A in Katgert et al. (Reference Katgert, Latka, Möbius and van Hecke2009)). Such a method results in slightly larger values of the yield stress within the range of measurement errors ( ${\sim}5\,\%$ ). Table 1 lists the final coefficients for the samples used in our experiments. The rheological properties of our samples are similar to those reported for 0.1 %–1.5 % mass concentration of Ultrez 10 Carbopol in water (Dinkgreve et al. Reference Dinkgreve, Paredes, Denn and Bonn2016). For increasing polymer concentration, the elastic modulus and the yield stress increase, i.e. the gel becomes stiffer. Table 1 also lists the values of the yield stress from the stress sweep tests, where $G^{\prime }$ and $G^{\prime \prime }$ intersect (another way to measure the yield stress). The yield-stress values from the intersection of the elastic moduli curves in the stress sweep tests are larger than those found by the Herschel–Bulkley fits as also reported in Dinkgreve et al. (Reference Dinkgreve, Paredes, Denn and Bonn2016). In the text, we use the values obtained by the Herschel–Bulkley fit as that finds the most accurate values of yield stress (see the discussion in Dinkgreve et al. (Reference Dinkgreve, Paredes, Denn and Bonn2016)).

Table 1. Values of yield stress $\unicode[STIX]{x1D70F}_{0}$ from the Herschel–Bulkley (HB) fits and from the intersection of $G^{\prime }$ and $G^{\prime \prime }$ in stress sweep (SS) tests, consistency index $K$ , flow index $n$ and storage modulus ${\mathcal{G}}_{0}^{\prime }$ , from (2.1).

Measuring the surface tension of yield-stress materials is known to be challenging. The classic methods such as capillary rise or pendant drop methods fail when the yield stress is high enough to govern the interface shape (see Boujlel & Coussot Reference Boujlel and Coussot2013; Jørgensen et al. Reference Jørgensen, Le Merrer, Delanoë-Ayari and Barentin2015). Nonetheless, the values of surface tension reported for (different types of) Carbopol in the literature (Hu, Wang & Hartnett Reference Hu, Wang and Hartnett1991; Ishiguro & Hartnett Reference Ishiguro and Hartnett1992; Manglik, Wasekar & Zhang Reference Manglik, Wasekar and Zhang2001; Boujlel & Coussot Reference Boujlel and Coussot2013; Jørgensen et al. Reference Jørgensen, Le Merrer, Delanoë-Ayari and Barentin2015) show only small deviations from that of water. We hence use the standard value for water at $20\,^{\circ }\text{C}$ , $\unicode[STIX]{x1D70E}=0.072~\text{N}~\text{m}^{-1}$ . Note that we could only measure the surface tension at very low gel concentrations (using a pendant drop method), where $\unicode[STIX]{x1D70E}=0.069\pm 0.02~\text{N}~\text{m}^{-1}$ was obtained. The slightly smaller values could be due to the presence of the other polymers inside the hair gel, which might even result in smaller surface tension values at higher gel concentrations. Therefore, the value of $0.072~\text{N}~\text{m}^{-1}$ should be taken as an upper bound for the real surface tensions of our samples. The density of the gels is just slightly ( ${\sim}0.2\,\%$ ) larger than that of water (measured using a DMA-35 Anton-Paar density meter).

2.2 Experimental set-up and parameters

Figure 2 shows a schematic of our experimental set-up. We trigger a pulsed green Nd:YAG laser (Litron Nano S65-15PIV) with a delay generator (BNC Model 575). The emitted laser beam has a diameter of 4 mm, a pulse duration of $t_{p}=6~\text{ns}$ and a wavelength of $\unicode[STIX]{x1D706}=532~\text{nm}$ . We control the energy of the laser ( $0.5~\text{mJ}<E<6.5~\text{mJ}$ ) with a $\unicode[STIX]{x1D706}/2$ wave plate, a polarizing beam splitter (PBS) and a beam dump (BD). We measure the pulse energy for every single experiment using a Gentec QE8SP-B-BL energy meter (EM). Note that, using a second energy meter (Gentec QE12LP-H-MB), we found the calibration relationship between the energy of the pulse inside the film and those measured. The laser pulse passes through a 10 $\times$ objective (Thorlabs LMH-10X-532) and is focused inside the gel layer. Using a Gaussian beam approximation we find a spot size diameter of $3.4~\unicode[STIX]{x03BC}\text{m}$ . To form the gel layers, we use spacers of given height and flatten a previously deposited blob of the material, using a blade with a thickness of $H=1~\text{mm}$ (see figure 2 b). The substrate is a glass slide clipped into a motorized 2-D stage. The yield stress of the materials prevents further spreading of the liquids. We can adjust the vertical position of the objective lens with a motorized stage, to control the focal spot $z_{f}$ . The side views are imaged with a high-speed camera (Photron FASTCAM SA-X2) attached to a long-distance microscope (Navitar 12X). The imaging is performed at a frame rate of $100\,000$ frames per second and a resolution of $17.8~\unicode[STIX]{x03BC}\text{m}$ per pixel. The illumination uses a collimated white light source (Sumita LS-M352A), passing through a diffuser. The set-up also utilizes a low-power continuous red laser for alignment and calibration of the vertical position (not shown here). The control parameters are the focal height $z_{f}$ , the laser energy $E$ and the rheological properties, e.g. the yield stress  $\unicode[STIX]{x1D70F}_{0}$ .

Figure 2. (a) Experimental set-up for the investigation of LIFT of viscoplastic materials (dimensions are not to scale). A laser pulse passes through a $\unicode[STIX]{x1D706}/2$ plate and a polarizing beam splitter (PBS), which serve to control the beam energy. We measure the energy of the laser with an energy meter (EM). The laser is guided through a 10 $\times$ objective lens (OL) and focused on a specified focal position. The side view is imaged with a high-speed camera attached to a long-distance microscope (LDM). The test section is illuminated through backlighting. The camera is protected with a notch filter. In the schematic, BD is the beam dump. (b) Magnified view of the test section. The grey part with thickness ${\mathcal{L}}$ is the cavity at the focal height $z_{f}$ ; $D_{s}$ is the spot size; and the dashed line shows the deformed interface with a characteristic early velocity of  $U_{0}$ .

3 Phenomenology

In this section, we first report on the phenomenological behaviour of the jetting process for the variation of each control parameter and then classify the morphological changes with the help of the relevant non-dimensional groups. Figure 3(a) shows the process of jet formation for different laser energies for sample 5 and $z_{f}/H=0.5$ . Since we do not add any absorbing dye to our material, the gel layer is transparent to the laser. Therefore, we only see the jet formation when the laser fluence is above the optical breakdown threshold, $I=I_{p}\approx 1.3\times 10^{8}~\text{J}~\text{m}^{-2}$ . This experimentally measured value is close to the previously reported breakdown threshold for water (Vogel et al. Reference Vogel, Nahen, Theisen and Noack1996; Kennedy, Hammer & Rockwell Reference Kennedy, Hammer and Rockwell1997) and remains the same for all the solutions, regardless of the polymer concentration. Above the optical threshold, plasma forms and a bubble starts to grow (Vogel et al. Reference Vogel, Nahen, Theisen and Noack1996; Ando, Liu & Ohl Reference Ando, Liu and Ohl2012; Obreschkow et al. Reference Obreschkow, Tinguely, Dorsaz, Kobel, De Bosset and Farhat2013; Supponen et al. Reference Supponen, Obreschkow, Tinguely, Kobel, Dorsaz and Farhat2016). The growth of the bubble and its interaction with the free surface (Mézel et al. Reference Mézel, Hallo, Souquet, Breil, Hébert and Guillemot2009; Koukouvinis et al. Reference Koukouvinis, Gavaises, Supponen and Farhat2016) results in jet formation. Visualization of bubble dynamics when the film is thin is experimentally challenging and not the focus of this article. We provide more details on the laser-induced bubble in appendix A and focus on the jet formation in the main body of the article.

When increasing the laser energy slightly above the plasma threshold, the interface only deforms to a small degree and forms a bump. Further increasing the laser energy results in jet formation. Depending on the experimental conditions, the jet might feature a crown, which forms due to the bubble dynamics inside the thin film (see appendix A). We measure the position of the jets as a function of time by high-speed imaging and subsequent image processing. Figure 3(b) shows the variation of the height of the jets shown in panel (a). When the laser energy is small, the bump relaxes, going through an oscillating behaviour. This is due to the propagation and interaction of the surface waves (see the supplementary movies, which are available at https://doi.org/10.1017/jfm.2019.731). However, the exact mechanisms involved in the propagation of these waves (surface tension, viscosity or elasticity) are not clear to us. When higher laser energies are applied, we only see the growth of the jets, in the time frame of our recordings.

Increasing the laser energy even further results in a faster jet. We find the jet velocity at small times by fitting the jet peak position by a second-order polynomial (for $t\leqslant 0.1~\text{ms}$ ) and extracting the slope at $t=0$ (see the inset in figure 3 b). We later compare these velocities for different conditions and use them to validate our model.

Figure 3. (a) Snapshots showing the evolution of the jet for sample 5 at $z_{f}=0.5~\text{mm}$ for various laser energies. The snapshots show the evolution of the jet for $t=0.0$ , 0.01, 0.02, 0.05, 0.1, 0.18, 0.33, 0.61, 1.1 and 2.0 ms. The scale bar corresponds to 2 mm. (b) The vertical extent of the jets versus time. The inset shows a zoom-in of the region for the early stage of ejection. The red dashed line demarcates the region $t\leqslant 0.1~\text{ms}$ in which a quadratic fit is applied to extrapolate the jet velocity at $t=0$ . See the supplementary movies for the examples shown here.

Figure 4. (a) Snapshots for increasing focal height at $t=0$ , 0.01, 0.02, 0.04, 0.1, 0.2, 0.4, 1 and 2 ms for $E=2.1~\text{mJ}$ and sample 5. The scale bar represents 2 mm. See the supplementary movies for the examples shown here.

Figure 5. Snapshots for four different samples at $t=0$ , 0.01, 0.02, 0.04, 0.1, 0.2, 0.4, 1 and 2 ms where energy and focal height are fixed at $E=5.1~\text{mJ}$ , and $z_{f}=0.5~\text{mm}$ , respectively. Gels become stiffer from top to bottom (see figure 1). The scale bar represents 2 mm. See the supplementary movies for the examples shown here.

Figure 4 shows the jet evolution from sample 5 for different focal height $z_{f}$ , with fixed laser energy at 2.1 mJ. Focusing on the free surface ( $z_{f}=0$ ) results in spray formation. This regime is different from the others since no bubble formation takes place inside the film. In fact, it is similar to what occurs in laser ablation of soft materials and liquids (Vogel & Venugopalan Reference Vogel and Venugopalan2003; Thoroddsen et al. Reference Thoroddsen, Takehara, Etoh and Ohl2009). The high pressure of plasma leads to the formation of a sheet with a highly unstable front, which fragments into small droplets. The cylindrical sheet closes later and forms a thin jet that is unstable and disintegrates further in the process. By increasing the focal length $z_{f}$ , we suppress the fragmentation and transit from sprays to jets of smaller height and velocity and eventually bumps.

Besides the laser energy, we also change the rheology of the samples. Figure 5 demonstrates the effect of rheology (yield stress $\unicode[STIX]{x1D70F}_{0}$ and nonlinear viscosity) on the process of jet formation. In the experiments shown, the laser energy ( $E=5.1~\text{mJ}$ ) and focal height ( $z_{f}=0.5$ ) are fixed, and the gel stiffness varies (i.e. higher yield stress and apparent viscosity from top to bottom). It is shown that the ejection regimes can be significantly modified by increasing the effective viscosity. When the gel is very soft (sample 7) fragmentation occurs, producing small droplets. Additionally, the crown features a rim instability (fingers). Such instabilities, however, are suppressed by increasing the yield stress. Finally, the larger the yield stress, the shorter and slower the produced jet is.

As we have now shown, the jet morphologies vary from bumps to jets and sprays for different laser energies, focal heights and rheological parameters. Next, we will theoretically derive the jetting velocity and provide quantitative analysis to classify these regimes.

4 Theoretical analysis of the experimental results

4.1 Characteristic jet velocity

One of our main goals is to categorize the jetting regimes with non-dimensional groups to find optimal conditions for printing (a clean straight jet). To do so, a characteristic velocity is required that must be known prior to the specific experiments. Here, we provide a simple analysis to find such a velocity valid for small times, i.e. early on in the jetting process. We assume the absorbed energy contributes to optical breakdown (plasma formation) and subsequent bubble growth. Based on this assumption, the balance between the energies before and after the absorption reads

(4.1) $$\begin{eqnarray}\displaystyle (I(t)-I_{p})\,\text{d}t=p(t)\,\text{d}{\mathcal{L}}, & & \displaystyle\end{eqnarray}$$

where $I$ (measured in units of $\text{J}~\text{m}^{-2}$ ) is the energy absorbed by the film per unit surface, $I_{p}$ is the plasma threshold fluence, $p$ is the pressure and ${\mathcal{L}}$ is the length of the cavity (see figure 2 b). We assume that the cavity grows at a constant speed of $\text{d}{\mathcal{L}}/\text{d}t=p/(\unicode[STIX]{x1D70C}c)$ , where $\unicode[STIX]{x1D70C}$ and $c$ are the density and the speed of sound, respectively (Fabbro et al. Reference Fabbro, Fournier, Ballard, Devaux and Virmont1990; Asay & Shahinpoor Reference Asay and Shahinpoor2012). If we also assume that the laser fluence remains constant during the pulse, we can find a characteristic pressure, namely

(4.2) $$\begin{eqnarray}\displaystyle p_{c}=(I-I_{p})^{1/2}(\unicode[STIX]{x1D70C}c)^{1/2}. & & \displaystyle\end{eqnarray}$$

This equation indeed gives the same order of magnitude ( ${\sim}10~\text{MPa}$ ) of the maximum pressures as reported in experiments with nanosecond lasers (Lauterborn & Vogel Reference Lauterborn and Vogel2013). Here we assume $c=1500~\text{m}~\text{s}^{-1}$ , as for water, since the polymer concentration is low and therefore its effect on the speed of sound is negligible (Povey Reference Povey1997; Parker & Povey Reference Parker and Povey2012).

Note that $p_{c}$ is a constant characteristic pressure that gives a simple scaling with respect to the change of the laser energy. We will later elaborate on the simplifications made here. We now assume that $p_{c}$ , over the characteristic time $t_{p}$ , provides the momentum to a material cylinder above it. Hence, $p_{c}At_{p}=Az_{f}\unicode[STIX]{x1D70C}U_{0}$ , where $A$ is the area of the spot size of diameter $D_{s}$ . We eventually find the characteristic velocity $U_{0}$ to be

(4.3) $$\begin{eqnarray}\displaystyle & \displaystyle U_{0}=\unicode[STIX]{x1D6FD}t_{p}\left(\frac{c}{\unicode[STIX]{x1D70C}}\right)^{1/2}\frac{(I-I_{p})^{1/2}}{z_{f}}. & \displaystyle\end{eqnarray}$$

In (4.3), we have included an unknown parameter $\unicode[STIX]{x1D6FD}$ that will be found later by fitting the experimental results.

We should emphasize that, by using a momentum balance like the one above, we make simplifications in five main ways. First, in reality, the pressure signal is time-dependent and strongly non-monotonic: it rapidly rises, reaches a maximum and then exponentially decays (see e.g. Lauterborn & Vogel Reference Lauterborn and Vogel2013). Also, the lifetime of this pressure signal is much larger than the laser pulse duration ( $t_{p}$ ) (Lauterborn & Vogel Reference Lauterborn and Vogel2013; Wang et al. Reference Wang, Zaytsev, Lajoinie, Eijkel, van den Berg, Versluis, Weckhuysen, Zhang, Zandvliet and Lohse2018). Therefore, a larger total momentum is expected. The rather large fitting parameter that we obtain later is mostly due to this assumption. Second, we assume that the reflection, scattering and transmission of the laser energy are negligible (as suggested by Vogel et al. (Reference Vogel, Noack, Nahen, Theisen, Busch, Parlitz, Hammer, Noojin, Rockwell and Birngruber1999)). Third, upon the optical breakdown, an (almost spherical) bubble forms that pushes the fluid and the free surface above (see appendix A). This will create a more complicated flow field inside the films, in comparison to what is assumed here (Mézel et al. Reference Mézel, Hallo, Souquet, Breil, Hébert and Guillemot2009; Brasz et al. Reference Brasz, Arnold, Stone and Lister2015; Koukouvinis et al. Reference Koukouvinis, Gavaises, Supponen and Farhat2016; Jalaal et al. Reference Jalaal, Li, Klein Schaarsberg, Qin and Lohse2019b ). Fourth, we ignore the effects of temperature change inside the liquid, as it can influence the density and consequently the speed of sound. The latter is justified, as one expects that the temperature only changes within a few micrometres near the plasma point and also reaches the ambient temperature in a few microseconds (Vogel et al. Reference Vogel, Noack, Nahen, Theisen, Busch, Parlitz, Hammer, Noojin, Rockwell and Birngruber1999). Lastly, the rheological effects, such as elasticity, plasticity and viscosity, as well as the surface tension effects are assumed to be insignificant to our problem at small time. One can examine such an assertion through the corresponding time scales. The viscous time scale $t_{vis}=(z_{f}^{2}\unicode[STIX]{x1D70C}/\unicode[STIX]{x1D707})^{1/2}\approx 10^{-2}~\text{s}$ , the gravitational time scale $t_{g}=(z_{f}/g)^{1/2}\approx 10^{-3}~\text{s}$ and the capillary time scale $t_{cap}=(\unicode[STIX]{x1D70C}z_{f}^{3}/\unicode[STIX]{x1D70E})\approx 10^{-3}~\text{s}$ are all much larger than the considered early stage time scale (which is the advection time scale, $t_{ad}=z_{f}/U_{0}\approx 10^{-5}~\text{s}$ ). Similarly, the elastic relaxation time scale $t_{el}\approx (K/G_{0}^{\prime })^{1/n}\approx 10^{-3}~\text{s}$ is also much larger than the characteristic time scale of the early stage deformation. These time scales justify our inviscid assumption to find the early stage characteristic velocity (see e.g. Brasz et al. (Reference Brasz, Arnold, Stone and Lister2015) for a similar argument in the context of jet formation).

Despite the simplifications listed above, equation (4.3) successfully predicts the scaling laws previously seen in the laser–liquid interaction tests. The predicted power-law exponent $1/2$ is indeed close to the previously reported values in the context of droplet propulsion with a laser pulse (Basko, Novikov & Grushin Reference Basko, Novikov and Grushin2015; Kurilovich et al. Reference Kurilovich, Klein, Torretti, Lassise, Hoekstra, Ubachs, Gelderblom and Versolato2016). Also, the predicted $U_{0}\sim z_{f}^{-1}$ dependence is the same as experimentally observed for jet formation in a capillary tube (Tagawa et al. Reference Tagawa, Oudalov, Visser, Peters, van der Meer, Sun, Prosperetti and Lohse2012).

The application of (4.3) is restricted to the domain that excludes the boundaries. Focusing the laser close to the glass slide, i.e. $z_{f}/H\rightarrow 1$ , results in solid ablation. This became evident when we found cracks in the glass slide. Moreover, equation (4.3) becomes singular when $z_{f}/H\rightarrow 0$ (i.e. focusing close to the free surface) and obviously can then no longer hold. In our experiments, the divergence of the velocity appears as an explosive splash of the film in the sheet opening and splash regime (see figure 4). Figure 6(a) shows the variation of the early jet velocities $U_{0}$ for different samples at different laser energies. As expected, we do not observe any significant effect of the rheological properties on $U_{0}$ , meaning that inertia is dominating any possible elastoviscoplastic effects. Moreover, the value of the plasma threshold $I_{p}$ also remains insensitive to the concentration of polymers. The predicted power law (4.3) is in good agreement with the trends observed in the experiments $U_{0}\propto (I-I_{p})^{1/2}$ . Fitting through the mean values of the data (see figure 6 a), we find $\unicode[STIX]{x1D6FD}\approx 174$ .

We then compare the model with experiments of samples at different focal heights. An example is shown in figure 6(b), where $\unicode[STIX]{x1D6FD}=174$ is fixed. The results of the model are in good agreement, implying that $U_{0}\sim z_{f}^{-1}$ is satisfactory. Note that we cannot compare the model for $z_{f}=0$ , as then $U_{0}\rightarrow \infty$ . Furthermore, we use an average value for the sound speed $c$ for the case of $z_{f}=H$ , as the plasma forms at the interface of the glass and the gel. Nonetheless, the difference between the experiments and theory remains large, which might be because of additional complexities due to the glass ablation, as explained above.

Figure 6. Jet velocity at different laser fluences for (a) different samples, but fixed $z_{f}=0.5$ , and (b) sample 5, at different focal heights, $z_{f}$ (see legend). The yellow vertical dashed lines denote $I=I_{p}$ . In both panels (a) and (b), the solid lines are from (4.3) with $\unicode[STIX]{x1D6FD}=174$ . Error bars denote the standard deviation from multiple experiments.

4.2 Phase diagram and non-dimensional groups

We identify six different regimes of jets from our experiments, namely a bump regime (the jet is shallow with a maximum aspect ratio of less than unity), a jet regime (when a long straight jet with no significant crown forms), a crown jet regime (where a crown follows the straight jet), an unstable crown jet regime (where the edge of the crown becomes unstable and forms fingers), a fragmented jet regime (where the jet disintegrates into smaller droplets in the early stage of formation), and finally a spray regime (where the sheet opens, forms many small droplets and the sheet closure results in even more fragmentation). In figure 7, we show examples of these regimes and the dimensional phase space in these control parameters: laser energy, yield stress (gel stiffness) and focal height. When $z_{f}=0$ , we only observe sprays. As a trend, increasing the gel stiffness and/or the focal height, as well as decreasing the laser energy, results in a less developed morphology.

Figure 7. (a) Examples of the jetting regimes observed in the experiments. The experimental conditions from top to bottom are: bump, $E=2.1~\text{mJ}$ , $z_{f}/H=0.75$ , sample no. 2; jet, $E=2.7~\text{mJ}$ , $z_{f}/H=0.5$ , sample no. 1; jet with a crown, $E=2.6~\text{mJ}$ , $z_{f}/H=0.25$ , sample no. 5; jet with an unstable crown, $E=6.3~\text{mJ}$ , $z_{f}/H=0.75$ , sample no. 6; fragmented jet, $E=6.4~\text{mJ}$ , $z_{f}/H=0.25$ , sample no. 6; spray, $E=6.3~\text{mJ}$ , $z_{f}/H=0$ , sample no. 5. (b) Three-dimensional phase space of the regimes for different laser energy $E$ , yield stress $\unicode[STIX]{x1D70F}_{0}$ and focal height $z_{f}$ . Crosses denote the experiments in which no deformation was detected. See the supplementary material for movies for the examples shown here and the 3-D animation of the phase space.

In general, inertial forces due to the laser pulse compete with capillary, viscous and plastic forces to define the dynamics and morphology of the jets. With the characteristic velocity found above (4.3), we define the following two non-dimensional groups:

(4.4a,b ) $$\begin{eqnarray}\displaystyle & \displaystyle We=\frac{\unicode[STIX]{x1D70C}U_{0}^{2}z_{f}}{\unicode[STIX]{x1D70E}}\quad \text{and}\quad Re=\frac{\unicode[STIX]{x1D70C}U_{0}^{2}}{K(U_{0}/z_{f})^{n}+\unicode[STIX]{x1D70F}_{0}}. & \displaystyle\end{eqnarray}$$

The Weber number $We$ compares the inertial stresses and capillary pressure. The Reynolds number $Re$ compares the inertial stress and the total internal stress. The latter contains two parts due to the shear thinning (with a characteristic shear rate of $U_{0}/z_{f}$ ) and the plastic viscosities. A Reynolds number of this type has previously been used in the context of viscoplastic droplets (Blackwell et al. Reference Blackwell, Deetjen, Gaudio and Ewoldt2015; Jalaal, Kemper & Lohse Reference Jalaal, Kemper and Lohse2019a ) and other configurations (Thompson & Soares Reference Thompson and Soares2016). Note that by using such a non-dimensional group, we assume that the yield stress only contributes via the plastic viscosity. We also ignore the effect of the film thickness $H$ and assume that the height of the laser focus $z_{f}$ is the only important length scale.

Figure 8 shows the non-dimensional phase diagram. Note that, in figure 8, we do not include data for $z_{f}/H=1$ , due to the complexities from the glass ablation, nor those for $z_{f}/H=0$ , as they only result in spray formation and our theoretical velocity (therefore, both $Re$ and $We$ ) diverges in that asymptotic limit. At relatively low Weber and Reynolds numbers, the inertial effects are overwhelmed by the viscous and capillary forces and only a bump forms. Increasing the Reynolds and Weber numbers results in jets, crown formation, unstable crowns and eventually fragmentation.

Figure 8. Phase space in terms of Reynolds and Weber numbers. The dashed grey lines correspond to $Ca=0.1$ and 1. The horizontal black line shows $Re=200$ . The shaded region highlights the optimal printing condition with LIFT.

We inspect the non-dimensional map using two other non-dimensional numbers common in printing technologies: the capillary number $Ca=We/Re$ and the Ohnesorge number $Oh=We^{1/2}/Re$ . Our experimental data are bounded within the two bounds of $Ca=1$ and $Ca=0.1$ . This is a consequence of the range of experimental parameters employed here. Using (4.3), for a given gel, one finds $Ca\sim (I-I_{0})^{n/2}z_{f}^{1-2n}$ . For our gels, $2n\approx 1$ ; therefore, $Ca$ is a function of only the laser energy. For the range of laser energy we used here, $(I-I_{0})^{n/2}$ remains close to unity. Therefore, $Ca$ is almost constant for a given gel. Therefore the straight line envelopes show up. The dashed grey lines in figure 8 denote these bounds.

The Ohnesorge number (in its simple form for Newtonian liquids) is commonly used in inkjet printing application to find the regimes in which one single droplet forms (typically $0.1<Oh<1$ ). If $Oh$ is large, then the material is too viscous to be printed; and if $Oh\leqslant 1$ , then multiple droplets might form (Derby Reference Derby2010; McKinley & Renardy Reference McKinley and Renardy2011). We use the same analogy here to find the optimal regime of printing, where a straight jet forms. All the experiments shown in figure 8 feature $Oh<1$ (a common criterion in inkjet printing noted by the red dashed line in figure 8). For the range of $0.075<Oh<0.2$ , we obtain jet morphologies that can be used for LIFT printing. Within this range, the proper balance between total viscosity and inertia does not allow for fragmentation or significant crown formation. We also introduce another criterion of $Re>200$ to exclude bumps. When $Re$ is smaller than this critical value, the viscous dissipation is strong enough to resist jet formation. The shaded region in figure 8 shows the region we suggest for optimal printing with LIFT.

5 Conclusions and outlook

In this paper, we experimentally studied the LIFT of viscoplastic materials, in which the ejection is driven by direct absorption of the laser. Several parameters, such as the laser energy, focal height and the rheology of the material, determine the shape of the jet, which can vary from bumps, jets with stable or unstable crowns, and fragmented jets, to sprays. In most applications, a non-fragmenting straight jet is desirable. To look for suitable conditions, we reduce the parameter space to two non-dimensional groups, namely the Reynolds and Weber numbers, which include rheological properties such as the nonlinear shear dependence and the yield stress.

Our dimensional analysis for the current LIFT set-up suggests operating conditions that satisfy the two conditions of $0.075<Oh<0.2$ and $Re>200$ . Outside these ranges, the viscous dissipation (due to the shear-thinning and plastic viscosities) is either too large to avoid jetting or too small to make a stable straight jet. We note that the conditions above are for jet formation only. Depending on the technical details of printing, additional conditions can be introduced for optimal printing. For instance, one should avoid splashing when droplets are deposited on the surface (Derby Reference Derby2010).

Our study has direct applications in printing 2-D or 3-D structures with soft materials, featuring viscoplastic properties. In particular, it provides information for printing flexible electronics with metallic pastes and hydrogels in bio-printing. Moreover, the non-dimensional numbers used here can be obtained prior to printing, knowing the properties of the laser pulse and the material (ink). The only other requirement is to find a correct characteristic velocity. Here, we use an energy balance argument to obtain the latter. The type of analysis provided here can, in its general form, be used for any LIFT configuration or any other printing system that utilizes jets of inelastic generalized Newtonian fluids. For instance, for a blister-actuated LIFT (e.g. Turkoz et al. Reference Turkoz, Perazzo, Kim, Stone and Arnold2018b ), one only needs to use a different characteristic velocity related to the deformation of the blister.

We finally note that inks, such as silver pastes, can feature even more complex rheological properties. Moreover, the detailed mechanisms of optical absorption and breakdown can vary. A first attempt to study such a more complicated system could be to test LIFT for a granular suspension with Newtonian or non-Newtonian suspending fluids (Guazzelli & Pouliquen Reference Guazzelli and Pouliquen2018).

Acknowledgements

We would like to thank B. Benschop, M. Bos and G.-W. Bruggert for their technical assistance. The authors thank C. Seyfert, A. Klein, G. Arutinov, D. van der Meer and A. Prosperetti for useful discussions. This work was funded by the Netherlands Organisation for Scientific Research (NWO).

Supplementary movies

Supplementary movies are available at https://doi.org/10.1017/jfm.2019.731.

Appendix A. Bubble dynamics

To elaborate more on the underlying fluid mechanics of a LIFT process, we performed model experiments with thick films. In contrast to the main experiments, where $H=1~\text{mm}$ , here we are able to visualize the dynamics of the bubble as well as the jet formation. Instead of a thin film, we fill a container and focus the laser close to the interface. Such an experiment is not a realistic mimic of a common LIFT system with a thin ink layer (it is instead close to the so-called film-free LIFT, see e.g. Piqué et al. (Reference Piqué, Serra, Surdo, Diaspro, Duocastella, Piqué and Serra2018)). Nonetheless, it provides useful information on the underlying fluid mechanics of the problem. Figure 9 shows an example of jet formation for a water film with $H=4~\text{mm}$ and $z_{f}=1~\text{mm}$ .

Figure 9. Process of bubble, jet and crown formation in a thick film. Snapshots correspond to (a) $t=0~\unicode[STIX]{x03BC}\text{s}$ , (b) $t=5~\unicode[STIX]{x03BC}\text{s}$ , (c) $t=125~\unicode[STIX]{x03BC}\text{s}$ , (d $t=250~\unicode[STIX]{x03BC}\text{s}$ , (e) $t=350~\unicode[STIX]{x03BC}\text{s}$ , (f) $t=450~\unicode[STIX]{x03BC}\text{s}$ , (g) $t=575~\unicode[STIX]{x03BC}\text{s}$ and (h) $t=625~\unicode[STIX]{x03BC}\text{s}$ .

The process is as follows. Optical breakdown results in plasma formation (shown by the black arrow in figure 9 a) and as a consequence a bubble forms inside the liquid (figure 9 b) that contains vapour (see Jalaal et al. (Reference Jalaal, Li, Klein Schaarsberg, Qin and Lohse2019b ) for more information on the early stage of LIFT). The bubble grows, initiating the liquid jet (figure 9 c) and then collapses under hydrostatic pressure (shown by blue arrows in figure 9 d). Meanwhile, the liquid jet is extending vertically (shown by the orange arrow in figure 9 d). During the bubble collapse a downward micro-jet forms (due to gravity, the presence of the wall and the free surface). This micro-jet impacts on the bottom of the bubble (see Supponen et al. (Reference Supponen, Obreschkow, Tinguely, Kobel, Dorsaz and Farhat2016) for details), resulting in a complex toroidal shape of secondary bubbles (figure 9 f). The formation of the crown seems to be highly connected to the second growth of the bubble, when it pushes the perimeter around the jet and forms a crown. A complete explanation of the crown dynamics, however, requires more experiments and analysis.

The process explained above is expected to be more complicated when the film thickness becomes thin, as the interaction between the bubble, the free surface and the bottom wall is expected to be more pronounced. Exploring such a system requires more sophisticated visualization tools.

Appendix B. Elasticity

Our analysis does not include the effects of elasticity. Here, we briefly justify our assumption and note the possible effects of elasticity. Instead of (2.1), one might use a more complicated constitutive model that features viscoelasticity after yielding. An example of an elastoviscoplastic constitutive model is as follows (Luu & Forterre Reference Luu and Forterre2009; Saramito Reference Saramito2009; Fraggedakis, Dimakopoulos & Tsamopoulos Reference Fraggedakis, Dimakopoulos and Tsamopoulos2016):

(B 1) $$\begin{eqnarray}\displaystyle & & \displaystyle \displaystyle \unicode[STIX]{x1D70F}=G^{\prime }\unicode[STIX]{x1D6FE}_{el}:\nonumber\\ \displaystyle & & \displaystyle \displaystyle \frac{\text{d}\unicode[STIX]{x1D6FE}_{el}}{\text{d}t}=\dot{\unicode[STIX]{x1D6FE}}\quad \text{if}\quad |\unicode[STIX]{x1D6FE}_{el}|<\unicode[STIX]{x1D6FE}_{c},\end{eqnarray}$$
(B 2) $$\begin{eqnarray}\displaystyle & & \displaystyle \displaystyle \frac{\text{d}\unicode[STIX]{x1D6FE}_{el}}{\text{d}t}=\dot{\unicode[STIX]{x1D6FE}}-\frac{1}{\unicode[STIX]{x1D706}}\frac{|\unicode[STIX]{x1D6FE}_{el}|-\unicode[STIX]{x1D6FE}_{c}}{|\unicode[STIX]{x1D6FE}_{el}|}\unicode[STIX]{x1D6FE}_{el}\quad \text{if}\quad |\unicode[STIX]{x1D6FE}_{el}|>\unicode[STIX]{x1D6FE}_{c},\end{eqnarray}$$

where $\unicode[STIX]{x1D6FE}_{el}$ is the reversible elastic deformation and $\unicode[STIX]{x1D6FE}_{c}$ is the yield strain. The model above, similar to (2.1), assumes that the material is a linear elastic solid before yielding. After yielding, however, the rate of deformation includes an additional creep term, featuring a relaxation time $\unicode[STIX]{x1D706}\approx (K/G^{\prime })^{1/n}$ . For a fast transient experiment, looking at equation (B 2), one might expect that, if the deformation time scale is much shorter than the relaxation time scale, the response is always dominated by elasticity, even after yielding. An argument of this type has been made by Luu & Forterre (Reference Luu and Forterre2009) in the context of the impact of yield-stress droplets. In making such an argument, one should, however, consider the dependence of elastic properties on deformation history. Our rheological measurements (see figure 1) show clear non-constant elastic properties above the yield stress (or yield strain of $\unicode[STIX]{x1D6FE}_{c}\approx 25\,\%$ ), as the magnitude of the elastic modulus $G^{\prime }$ drops orders of magnitude when deformation is large. Therefore, for the time window of our study, we expect that the viscoplastic response is dominating. One should note that, much later in the process, when the jet, droplets or the crater are relaxing, the elastic effects can be important. Therefore, to investigate these processes, we should include the effects of elasticity, using dimensionless numbers such as the Deborah number (Luu & Forterre Reference Luu and Forterre2009; Turkoz et al. Reference Turkoz, Perazzo, Kim, Stone and Arnold2018b ).

References

Ando, K., Liu, A.-Q. & Ohl, C.-D. 2012 Homogeneous nucleation in water in microfluidic channels. Phys. Rev. Lett. 109 (4), 044501.Google Scholar
Asay, J. R. & Shahinpoor, M. 2012 High-Pressure Shock Compression of Solids. Springer Science & Business Media.Google Scholar
Balmforth, N. J., Frigaard, I. A. & Ovarlez, G. 2014 Yielding to stress: recent developments in viscoplastic fluid mechanics. Annu. Rev. Fluid Mech. 46, 121146.Google Scholar
Basko, M. M., Novikov, V. G. & Grushin, A. S. 2015 On the structure of quasi-stationary laser ablation fronts in strongly radiating plasmas. Phys. Plasmas 22 (5), 053111.Google Scholar
Blackwell, B. C., Deetjen, M. E., Gaudio, J. E. & Ewoldt, R. H. 2015 Sticking and splashing in yield-stress fluid drop impacts on coated surfaces. Phys. Fluids 27 (4), 043101.Google Scholar
Bonn, D., Denn, M. M., Berthier, L., Divoux, T. & Manneville, S. 2017 Yield stress materials in soft condensed matter. Rev. Mod. Phys. 89 (3), 035005.Google Scholar
Boujlel, J. & Coussot, P. 2013 Measuring the surface tension of yield stress fluids. Soft Matt. 9 (25), 58985908.Google Scholar
Brasz, C. F., Arnold, C. B., Stone, H. A. & Lister, J. R. 2015 Early-time free-surface flow driven by a deforming boundary. J. Fluid Mech. 767, 811841.Google Scholar
Colina, M., Duocastella, M., Fernández-Pradas, J. M., Serra, P. & Morenza, J. L. 2006 Laser-induced forward transfer of liquids: study of the droplet ejection process. J. Appl. Phys. 99 (8), 084909.Google Scholar
Derby, B. 2010 Inkjet printing of functional and structural materials: fluid property requirements, feature stability, and resolution. Annu. Rev. Mater. Res. 40, 395414.Google Scholar
Dinkgreve, M., Paredes, J., Denn, M. M. & Bonn, D. 2016 On different ways of measuring the yield stress. J. Non-Newtonian Fluid Mech. 238, 233241.Google Scholar
Fabbro, R., Fournier, J., Ballard, P., Devaux, D. & Virmont, J. 1990 Physical study of laser-produced plasma in confined geometry. J. Appl. Phys. 68 (2), 775784.Google Scholar
Fraggedakis, D., Dimakopoulos, Y. & Tsamopoulos, J. 2016 Yielding the yield stress analysis: a thorough comparison of recently proposed elasto-visco-plastic (evp) fluid models. J. Non-Newtonian Fluid Mech. 236, 104122.Google Scholar
Géraud, B., Jørgensen, L., Ybert, C., Delanoë-Ayari, H. & Barentin, C. 2017 Structural and cooperative length scales in polymer gels. Eur. Phys. J. E 40 (1), 5.Google Scholar
Guazzelli, É. & Pouliquen, O. 2018 Rheology of dense granular suspensions. J. Fluid Mech. 852, P1.Google Scholar
Hopkins, C. C. & de Bruyn, J. R. 2019 Gelation and long-time relaxation of aqueous solutions of pluronic f127. J. Rheol. 63 (1), 191201.Google Scholar
Hu, R. Y. Z., Wang, A. T. A. & Hartnett, J. P. 1991 Surface tension measurement of aqueous polymer solutions. Exp. Therm. Fluid Sci. 4 (6), 723729.Google Scholar
Ishiguro, S. & Hartnett, J. P. 1992 Surface tension of aqueous polymer solutions. Intl Commun. Heat Mass Transfer 19 (2), 285295.Google Scholar
Jalaal, M., Balmforth, N. J. & Stoeber, B. 2015 Slip of spreading viscoplastic droplets. Langmuir 31 (44), 1207112075.Google Scholar
Jalaal, M., Cottrell, G., Balmforth, N. & Stoeber, B. 2017 On the rheology of pluronic f127 aqueous solutions. J. Rheol. 61 (1), 139146.Google Scholar
Jalaal, M., Kemper, D. & Lohse, D. 2019a Viscoplastic water entry. J. Fluid Mech.; (to appear).Google Scholar
Jalaal, M., Li, S., Klein Schaarsberg, M., Qin, Y. & Lohse, D. 2019b Destructive mechanisms in laser induced forward transfer. Appl. Phys. Lett. 114 (21), 213703.Google Scholar
Jørgensen, L., Le Merrer, M., Delanoë-Ayari, H. & Barentin, C. 2015 Yield stress and elasticity influence on surface tension measurements. Soft Matt. 11 (25), 51115121.Google Scholar
Katgert, G., Latka, A., Möbius, M. E. & van Hecke, M. 2009 Flow in linearly sheared two-dimensional foams: from bubble to bulk scale. Phys. Rev. E 79 (6), 066318.Google Scholar
Kennedy, P. K., Hammer, D. X. & Rockwell, B. A. 1997 Laser-induced breakdown in aqueous media. Prog. Quantum Electron. 21 (3), 155248.Google Scholar
Kim, J.-Y., Song, J.-Y., Lee, E.-J. & Park, S.-K. 2003 Rheological properties and microstructures of carbopol gel network system. Colloid Polym. Sci. 281 (7), 614623.Google Scholar
Koukouvinis, P., Gavaises, M., Supponen, O. & Farhat, M. 2016 Simulation of bubble expansion and collapse in the vicinity of a free surface. Phys. Fluids 28 (5), 052103.Google Scholar
Kurilovich, D., Klein, A. L., Torretti, F., Lassise, A., Hoekstra, R., Ubachs, W., Gelderblom, H. & Versolato, O. O. 2016 Plasma propulsion of a metallic microdroplet and its deformation upon laser impact. Phys. Rev. A 6 (1), 014018.Google Scholar
Lauterborn, W. & Vogel, A. 2013 Shock wave emission by laser generated bubbles. In Bubble Dynamics and Shock Waves, pp. 67103. Springer.Google Scholar
Lidon, P., Villa, L. & Manneville, S. 2017 Power-law creep and residual stresses in a carbopol gel. Rheol. Acta 56 (3), 307323.Google Scholar
Luu, L.-H. & Forterre, Y. 2009 Drop impact of yield-stress fluids. J. Fluid Mech. 632, 301327.Google Scholar
Manglik, R. M., Wasekar, V. M. & Zhang, J. 2001 Dynamic and equilibrium surface tension of aqueous surfactant and polymeric solutions. Exp. Therm. Fluid Sci. 25 (1–2), 5564.Google Scholar
McKinley, G. H. & Renardy, M. 2011 Wolfgang von ohnesorge. Phys. Fluids 23 (12), 127101.Google Scholar
Meeker, S. P., Bonnecaze, R. T. & Cloitre, M. 2004 Slip and flow in soft particle pastes. Phys. Rev. Lett. 92 (19), 198302.Google Scholar
Mézel, C., Hallo, L., Souquet, A., Breil, J., Hébert, D. & Guillemot, F. 2009 Self-consistent modeling of jet formation process in the nanosecond laser pulse regime. Phys. Plasmas 16 (12), 123112.Google Scholar
Munoz-Martin, D., Brasz, C. F., Chen, Y., Morales, M., Arnold, C. B. & Molpeceres, C. 2016 Laser-induced forward transfer of high-viscosity silver pastes. Appl. Surf. Sci. 366, 389396.Google Scholar
Obreschkow, D., Tinguely, M., Dorsaz, N., Kobel, P., De Bosset, A. & Farhat, M. 2013 The quest for the most spherical bubble: experimental setup and data overview. Exp. Fluids 54 (4), 1503.Google Scholar
Parker, N. G. & Povey, M. J. W. 2012 Ultrasonic study of the gelation of gelatin: phase diagram, hysteresis and kinetics. Food Hydrocolloids 26 (1), 99107.Google Scholar
Piqué, A., Serra, P., Surdo, S., Diaspro, A. & Duocastella, M. 2018 Film-Free LIFT (FF-LIFT). In Laser Printing of Functional Materials (ed. Piqué, A. & Serra, P.). doi:10.1002/9783527805105.ch6.Google Scholar
Pohl, R., Visser, C. W., Römer, G.-W., Lohse, D., Sun, C., Huis, B. & others 2015 Ejection regimes in picosecond laser-induced forward transfer of metals. Phys. Rev. A 3 (2), 024001.Google Scholar
Povey, M. J. W. 1997 Ultrasonic Techniques for Fluids Characterization. Elsevier.Google Scholar
Roberts, G. P. & Barnes, H. A. 2001 New measurements of the flow-curves for carbopol dispersions without slip artefacts. Rheol. Acta 40 (5), 499503.Google Scholar
Saramito, P. 2009 A new elastoviscoplastic model based on the Herschel–Bulkley viscoplastic model. J. Non-Newtonian Fluid Mech. 158 (1–3), 154161.Google Scholar
Serra, P. & Piqué, A. 2018 Laser-induced forward transfer: fundamentals and applications. In Advanced Materials Technologies, p. 1800099.Google Scholar
Supponen, O., Obreschkow, D., Tinguely, M., Kobel, P., Dorsaz, N. & Farhat, M. 2016 Scaling laws for jets of single cavitation bubbles. J. Fluid Mech. 802, 263293.Google Scholar
Tagawa, Y., Oudalov, N., Visser, C. W., Peters, I. R., van der Meer, D., Sun, C., Prosperetti, A. & Lohse, D. 2012 Highly focused supersonic microjets. Phys. Rev. X 2 (3), 031002.Google Scholar
Thompson, R. L. & Soares, E. J. 2016 Viscoplastic dimensionless numbers. J. Non-Newtonian Fluid Mech. 238, 5764.Google Scholar
Thoroddsen, S. T., Takehara, K., Etoh, T. G. & Ohl, C.-D. 2009 Spray and microjets produced by focusing a laser pulse into a hemispherical drop. Phys. Fluids 21 (11), 112101.Google Scholar
Turkoz, E., Fardel, R. & Arnold, C. B. 2018a Advances in blister-actuated laser-induced forward transfer (ba-lift). In Laser Printing of Functional Materials: 3D Microfabrication, Electronics and Biomedicine, pp. 91121.Google Scholar
Turkoz, E., Perazzo, A., Kim, H., Stone, H. A. & Arnold, C. B. 2018b Impulsively induced jets from viscoelastic films for high-resolution printing. Phys. Rev. Lett. 120 (7), 074501.Google Scholar
Visser, C. W., Pohl, R., Sun, C., Römer, G.-W., Huis in ’t Veld, B. & Lohse, D. 2015 Toward 3D printing of pure metals by laser-induced forward transfer. Adv. Mater. 27 (27), 40874092.Google Scholar
Vogel, A., Nahen, K., Theisen, D. & Noack, J. 1996 Plasma formation in water by picosecond and nanosecond nd: Yag laser pulses. I. Optical breakdown at threshold and superthreshold irradiance. IEEE J. Selected Topics Quantum Electron. 2 (4), 847860.Google Scholar
Vogel, A., Noack, J., Nahen, K., Theisen, D., Busch, S., Parlitz, U., Hammer, D. X., Noojin, G. D., Rockwell, B. A. & Birngruber, R. 1999 Energy balance of optical breakdown in water at nanosecond to femtosecond time scales. Appl. Phys. B 68 (2), 271280.Google Scholar
Vogel, A. & Venugopalan, V. 2003 Mechanisms of pulsed laser ablation of biological tissues. Chem. Rev. 103 (2), 577644.Google Scholar
Wang, J., Auyeung, R. C. Y., Kim, H., Charipar, N. A. & Piqué, A. 2010 Three-dimensional printing of interconnects by laser direct-write of silver nanopastes. Adv. Mater. 22 (40), 44624466.Google Scholar
Wang, Y., Zaytsev, M. E., Lajoinie, G., Eijkel, J. C. T., van den Berg, A., Versluis, M., Weckhuysen, B. M., Zhang, X., Zandvliet, H. J. W. & Lohse, D. 2018 Giant and explosive plasmonic bubbles by delayed nucleation. Proc. Natl Acad. Sci. 115 (30), 76767681.Google Scholar
Xiong, R., Zhang, Z., Chai, W., Huang, Y. & Chrisey, D. B. 2015 Freeform drop-on-demand laser printing of 3D alginate and cellular constructs. Biofabrication 7 (4), 045011.Google Scholar
Zenou, M., Sa’ar, A. & Kotler, Z. 2015 Laser transfer of metals and metal alloys for digital microfabrication of 3D objects. Small 11 (33), 40824089.Google Scholar
Figure 0

Figure 1. Rheology of the seven viscoplastic fluids analysed. Symbols and numbers refer to the different samples, whose material properties are given in table 1. (a) Storage modulus $G^{\prime }$ (symbols) and loss modulus $G^{\prime \prime }$ (lines) as functions of stress. (b) Flow curves: symbols are the experimental data and lines are the Herschel–Bulkley fits. (c) Viscous stress $\unicode[STIX]{x1D70F}-\unicode[STIX]{x1D70F}_{0}$ versus the shear rate (symbols) and their corresponding Herschel–Bulkley fits (lines).

Figure 1

Table 1. Values of yield stress $\unicode[STIX]{x1D70F}_{0}$ from the Herschel–Bulkley (HB) fits and from the intersection of $G^{\prime }$ and $G^{\prime \prime }$ in stress sweep (SS) tests, consistency index $K$, flow index $n$ and storage modulus ${\mathcal{G}}_{0}^{\prime }$, from (2.1).

Figure 2

Figure 2. (a) Experimental set-up for the investigation of LIFT of viscoplastic materials (dimensions are not to scale). A laser pulse passes through a $\unicode[STIX]{x1D706}/2$ plate and a polarizing beam splitter (PBS), which serve to control the beam energy. We measure the energy of the laser with an energy meter (EM). The laser is guided through a 10$\times$ objective lens (OL) and focused on a specified focal position. The side view is imaged with a high-speed camera attached to a long-distance microscope (LDM). The test section is illuminated through backlighting. The camera is protected with a notch filter. In the schematic, BD is the beam dump. (b) Magnified view of the test section. The grey part with thickness ${\mathcal{L}}$ is the cavity at the focal height $z_{f}$; $D_{s}$ is the spot size; and the dashed line shows the deformed interface with a characteristic early velocity of $U_{0}$.

Figure 3

Figure 3. (a) Snapshots showing the evolution of the jet for sample 5 at $z_{f}=0.5~\text{mm}$ for various laser energies. The snapshots show the evolution of the jet for $t=0.0$, 0.01, 0.02, 0.05, 0.1, 0.18, 0.33, 0.61, 1.1 and 2.0 ms. The scale bar corresponds to 2 mm. (b) The vertical extent of the jets versus time. The inset shows a zoom-in of the region for the early stage of ejection. The red dashed line demarcates the region$t\leqslant 0.1~\text{ms}$ in which a quadratic fit is applied to extrapolate the jet velocity at $t=0$. See the supplementary movies for the examples shown here.

Figure 4

Figure 4. (a) Snapshots for increasing focal height at $t=0$, 0.01, 0.02, 0.04, 0.1, 0.2, 0.4, 1 and 2 ms for $E=2.1~\text{mJ}$ and sample 5. The scale bar represents 2 mm. See the supplementary movies for the examples shown here.

Figure 5

Figure 5. Snapshots for four different samples at $t=0$, 0.01, 0.02, 0.04, 0.1, 0.2, 0.4, 1 and 2 ms where energy and focal height are fixed at $E=5.1~\text{mJ}$, and $z_{f}=0.5~\text{mm}$, respectively. Gels become stiffer from top to bottom (see figure 1). The scale bar represents 2 mm. See the supplementary movies for the examples shown here.

Figure 6

Figure 6. Jet velocity at different laser fluences for (a) different samples, but fixed $z_{f}=0.5$, and (b) sample 5, at different focal heights, $z_{f}$ (see legend). The yellow vertical dashed lines denote $I=I_{p}$. In both panels (a) and (b), the solid lines are from (4.3) with $\unicode[STIX]{x1D6FD}=174$. Error bars denote the standard deviation from multiple experiments.

Figure 7

Figure 7. (a) Examples of the jetting regimes observed in the experiments. The experimental conditions from top to bottom are: bump, $E=2.1~\text{mJ}$, $z_{f}/H=0.75$, sample no. 2; jet, $E=2.7~\text{mJ}$, $z_{f}/H=0.5$, sample no. 1; jet with a crown, $E=2.6~\text{mJ}$, $z_{f}/H=0.25$, sample no. 5; jet with an unstable crown, $E=6.3~\text{mJ}$, $z_{f}/H=0.75$, sample no. 6; fragmented jet, $E=6.4~\text{mJ}$, $z_{f}/H=0.25$, sample no. 6; spray, $E=6.3~\text{mJ}$, $z_{f}/H=0$, sample no. 5. (b) Three-dimensional phase space of the regimes for different laser energy $E$, yield stress $\unicode[STIX]{x1D70F}_{0}$ and focal height $z_{f}$. Crosses denote the experiments in which no deformation was detected. See the supplementary material for movies for the examples shown here and the 3-D animation of the phase space.

Figure 8

Figure 8. Phase space in terms of Reynolds and Weber numbers. The dashed grey lines correspond to $Ca=0.1$ and 1. The horizontal black line shows $Re=200$. The shaded region highlights the optimal printing condition with LIFT.

Figure 9

Figure 9. Process of bubble, jet and crown formation in a thick film. Snapshots correspond to (a) $t=0~\unicode[STIX]{x03BC}\text{s}$, (b) $t=5~\unicode[STIX]{x03BC}\text{s}$, (c) $t=125~\unicode[STIX]{x03BC}\text{s}$, (d$t=250~\unicode[STIX]{x03BC}\text{s}$, (e) $t=350~\unicode[STIX]{x03BC}\text{s}$, (f) $t=450~\unicode[STIX]{x03BC}\text{s}$, (g) $t=575~\unicode[STIX]{x03BC}\text{s}$ and (h) $t=625~\unicode[STIX]{x03BC}\text{s}$.

Jalaal et al. supplementary movie 1

Jetting for sample 5 at $z_f=0.5\,$mm and $E=1.4\,$mJ. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 1(Video)
Video 637 KB

Jalaal et al. supplementary movie 2

Jetting for sample 5 at $z_f=0.5\,$mm and $E=1.7\,$mJ. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 2(Video)
Video 631 KB

Jalaal et al. supplementary movie 3

Jetting for sample 5 at $z_f=0.5\,$mm and $E=2.1\,$mJ. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 3(Video)
Video 654 KB

Jalaal et al. supplementary movie 4

Jetting for sample 5 at $z_f=0.5\,$mm and $E=3.3\,$mJ. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 4(Video)
Video 654 KB

Jalaal et al. supplementary movie 5

Jetting for sample 5 at $z_f=0.5\,$mm and $E=6.4\,$mJ. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 5(Video)
Video 657 KB

Jalaal et al. supplementary movie 6

Jetting for sample 5 at $E=2.1\,$mJ and $z_f=0\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 6(Video)
Video 804 KB

Jalaal et al. supplementary movie 7

Jetting for sample 5 at $E=2.1\,$mJ and $z_f=0.5\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 7(Video)
Video 654 KB

Jalaal et al. supplementary movie 8

Jetting for sample 5 at $E=2.1\,$mJ and $z_f=0.75\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 8(Video)
Video 581 KB

Jalaal et al. supplementary movie 9

Jetting for sample 5 at $E=2.1\,$mJ and $z_f=1\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 9(Video)
Video 657 KB

Jalaal et al. supplementary movie 10

Jetting for sample 7 at $E=5.1\,$mJ and $z_f=0.5\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 10(Video)
Video 913 KB

Jalaal et al. supplementary movie 11

Jetting for sample 6 at $E=5.1\,$mJ and $z_f=0.5\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 11(Video)
Video 630 KB

Jalaal et al. supplementary movie 12

Jetting for sample 3 at $E=5.1\,$mJ and $z_f=0.5\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 12(Video)
Video 645 KB

Jalaal et al. supplementary movie 13

Jetting for sample 1 at $E=5.1\,$mJ and $z_f=0.5\,$mm. The scale-bar corresponds to 2 mm.

Download Jalaal et al. supplementary movie 13(Video)
Video 612 KB

Jalaal et al. supplementary movie 14

Examples of the jetting regimes observed in the experiments. The experimental conditions from top to bottom: \textit{bump}: $E=2.1\,$mJ, $z_f/H=0.75$, Sample $\#=2$; \textit{jet}: $E=2.7\,$mJ, $z_f/H=0.5$, Sample $\#=1$; \textit{jet with a crown}:$E=2.6\,$mJ, $z_f/H=0.25$, Sample $\#=5$; \textit{jet with an unstable crown}: $E=6.3\,$mJ, $z_f/H=0.75$, Sample $\#=6$; \textit{fragmented jet}: $E=6.4\,$mJ, $z_f/H=0.25$, Sample $\#=6$; \textit{spray}: $E=6.3\,$mJ, $z_f/H=0$, Sample $\#=5$.

Download Jalaal et al. supplementary movie 14(Video)
Video 8 MB

Jalaal et al. supplementary movie 15

3D animation of the three-dimensional phase space of the regimes for different laser energy, $E$, yield stress, $\tau_0$, and the focal height, $z_f$. Crosses denote the experiments in which no deformation was detected.

Download Jalaal et al. supplementary movie 15(Video)
Video 647 KB
Supplementary material: PDF

Jalaal et al. supplementary captions list

Movie Captions

Download Jalaal et al. supplementary captions list(PDF)
PDF 66 KB