The retraction of jetted slender viscoelastic liquid filaments

Long and slender liquid filaments are produced during inkjet printing, which can subsequently either retract to form a single droplet, or break up to form a primary droplet and one or more satellite droplets. These satellite droplets are undesirable since they degrade the quality and reproducibility of the print. Existing strategies for the suppression of satellite droplet formation include, among others, adding viscoelasticity to the ink. In the present work, we aim to improve the understanding of the role of viscoelasticity in suppressing satellite droplets in inkjet printing. We demonstrate that very dilute viscoelastic aqueous solutions (concentrations $\sim$ 0.003\% wt. polyethylene oxide (PEO), corresponding to nozzle Deborah number $De_{n}$ $\sim$ 3) can suppress satellite droplet formation. Furthermore, we show that, for a given driving condition, upper and lower bounds of polymer concentration exist, within which satellite droplets are suppressed. Satellite droplets are formed at concentrations below the lower bound, while jetting ceases for concentrations above the upper bound (for fixed driving conditions). Moreover, we observe that, with concentrations in between the two bounds, the filaments retract at velocities larger than the corresponding Taylor-Culick velocity for the Newtonian case. We show that this enhanced retraction velocity can be attributed to the elastic tension due to polymer stretching, which builds up during the initial jetting phase. These results shed some light on the complex interplay between inertia, capillarity, and viscoelasticity for retracting liquid filaments, which is important for the stability and quality of inkjet printing of polymer solutions.


Introduction
Drop-on-demand inkjet printing is known for its capability of highly-controlled, noncontact deposition of picoliters of liquid material (Basaran 2002;Wijshoff 2010;Derby 2010;Hoath 2016;Lohse 2022). Recent advances have enabled and enhanced the ability to deposit liquids over a wide range of surface tensions and viscosities (Castrejón-Pita et al. 2013). The high degree of reproducibility has led to inkjet printing applications in a diverse array of applications, including text or graphical printing on paper, fabrication of displays in electronics (Shimoda et al. 2003), electronics printing (Majee et al. 2016(Majee et al. , 2017, and in the life sciences (Villar et al. 2013;Daly et al. 2015;Simaite et al. 2016).
A typical inkjet printhead primarily consists of an ink reservoir, a piezo-acoustic transducer, and a dispensing nozzle (Wijshoff 2010), while the most simple driving waveform is a monopolar trapezoidal pulse (Castrejón-Pita et al. 2008), with a pulse width equal to half the period corresponding to the resonance frequency. Inkjet printheads are usually operated in the 'pull-push' mode, where the liquid is first pulled into the nozzle during the rise time of the trapezoidal pulse and then pushed out during the fall time of the pulse (Fraters et al. 2020). This results in the creation of a slender liquid jet of finite length and after pinch-off from the nozzle, a finite liquid ligament with a relatively large head droplet and a long tail. As the ligament is traveling towards the substrate, the tail retracts into the head droplet due to surface tension. However, during such motion, the tail may also break up due to the Rayleigh-Plateau instability (Fraters et al. 2020). This breakup leads to the formation of satellite droplets, which travel at a velocity lower than that of the head droplet. Thus the head droplet and the satellite droplet(s) reach the substrate at different times, potentially resulting in misalignment and substantially reduced print quality (Wijshoff 2010;Derby 2010), and contamination within the print device.
The detrimental effect of the formation of satellite droplets has resulted in an increased focus on the development of techniques to suppress such satellites. These techniques are either based on changing the driving waveform (Dong et al. 2006;Fraters et al. 2020) or modifying the properties of the ink. For the latter, since the satellite droplets primarily result from a Rayleigh-Plateau instability of the retracting tail, an obvious choice is to increase the viscosity of the ink in order to stabilize the tail filament (Notz & Basaran 2004;Castrejón-Pita et al. 2012;Driessen et al. 2013;Wang et al. 2019;Planchette et al. 2019;Anthony et al. 2019).
Another strategy (Christanti & Walker 2002) for the suppression of satellite droplet formation is to include polymer additives, which impart viscoelasticity, in the liquid being jetted. Early work has shown that viscoelasticity can stabilize a capillary jet against breakup (Goldin et al. 1969), and can also suppress satellite droplets when the liquid is jetted by a forced disturbance (such as in inkjet printing) (Christanti & Walker 2002). Shore & Harrison (2005) experimentally demonstrated that the addition of polymers to water-based inks suppresses satellite droplet formation in an inkjet printing configuration, while the same observation was also reported in the numerical investigation by Morrison & Harlen (2010). Yan et al. (2011) studied the effect of adding polymers, specifically polyethylene oxide (PEO), to water-based inks in inkjet printing, and concluded that the addition of low molecular weight polymers has no significant effect on the overall dynamics of the jetting behavior. Hoath et al. (2012) identified experimentally scaling laws relating the maximum jettable concentration to the molecular weight of the polymer additive for polystyrene (PS) in diethyl phthalate (DEP) solutions. For these polymer solutions, it was also reported (Hoath et al. 2014) that there is a delay in the pinch-off of the liquid filament from the nozzle, as compared to a purely Newtonian ink.
Despite the recent surge in the investigations of inkjet printing with polymer solutions, and the extensive literature (Bazilevskii et al. 1990;Anna & McKinley 2001;Amarouchene et al. 2001;Clasen et al. 2006a;Eggers et al. 2020;Zhou & Doi 2018 on the thinning and breakup of polymeric jets and pendant drops, there seems to be a dearth in the quantitative understanding of the fundamental physical mechanisms responsible for the suppression of satellite droplets in such conditions. Here, we quantitatively study the breakup (or retraction) of the long ligaments produced during (polymeric) ink jetting and the subsequent satellite droplet formation (or suppression). In particular, we identify the operating range where no satellite droplets are observed. In this range, the temporal retraction behavior of the jetted ligament length is characterized, and explained by a simple theoretical model. The model also quantifies the forces responsible for the suppression of satellite droplets, and demonstrates reasonable agreement with the experimental observations. The paper is organized as follows: § 2 describes the experimental procedure. In § 3 the experimental results for different values of the control parameters are shown, culminating in the phase diagram (Fig. 4). In § 4 we present detailed and quantitative experimental measurements of the filament retraction, which are theoretically explained in § 5. The paper ends with conclusions and an outlook in § 6.

Experimental procedure
Along the lines of previous work (Christanti & Walker 2002;Shore & Harrison 2005;Yan et al. 2011), PEO (average molecular weight 10 6 a.u., Sigma-Aldrich, henceforth referred to as PEO1M) was chosen to be the polymer additive in this work. Aqueous solutions of PEO1M, of concentrations (c, by mass) ranging from 0.001% to 0.009%, were prepared by adding the required amount of polymer powder to purified water (Milli-Q). For each polymer concentration, a 100 mL solution is made at first and all experiments were carried out using the same solution. The required polymer powder amount was measured using a precision laboratory balance (Secura 224-1S, Sartorius) with an accuracy of 0.1 mg, and then added to purified water (Milli-Q). Each solution was stirred with a magnetic stirrer for 24 h prior to use in order to ensure homogeneity of concentration. This results in a high degree of repeatability of the experiments (as evident later from the small error bars in Fig. 5a). In the present experiments, we wanted to have polymer solutions that can act as Boger fluids (James 2009), i.e., solutions whose shear viscosities are independent of the shear rate. This allows for the development of the simplified model, shown in § 5, where the shear rate dependence on viscosity is not taken into account. Dilute aqueous solutions of PEO1M act as Boger fluids (Cooper-White et al. 2002). Thus PEO1M was chosen as the polymer to be tested. This was also verified independently in the present work. The rheological characterization of the test liquids was performed on a rotational rheometer (MCR 502, Anton-Paar) with a cone-andplate configuration (1 • angle, 50 mm diameter, and mean gap of 0.1 mm). The measured viscosities are mentioned in Table 1, where c * is the critical overlap concentration (Clasen et al. 2006b). That being said, the results of the present work, in particular the simplified theoretical model ( § 5), can be extended to any Boger fluid composed of a different polymer. This can be easily facilitated by characterizing the rheological properties of the test fluid such as relaxation time (λ), solvent (η s ) and polymer (η p ) viscosities, along with knowing the density (ρ) and surface tension (γ). These material properties can then be incorporated into the theoretical approach ( § 5). The relaxation times were measured from the extensional thinning of the liquid filaments in a pendant droplet configuration (Deblais et al. 2018(Deblais et al. , 2020) (see also relaxation time measurement in a CaBER (Capillary Breakup Extensional Rheometer) device (Bazilevskii et al. 1990;Anna & McKinley 2001;Amarouchene et al. 2001;Clasen et al. 2006a)). The surface tension was measured by the pendant droplet method in an optical contact angle measurement and contour analysis instrument (OCA15, Dataphysics). It was noted that the surface tensions and the viscosities for the solutions of different polymer concentrations remained practically unchanged, while an appreciable change was observed in the relaxation time measurements (also observed by Yan et al. (2011)). It is known from literature (Noskov et al. 2000) that changing the PEO concentration can change the surface properties of the solutions. However, in the very narrow concentration range (0.001% -0.009% by wt.) used in the present work, the surface properties are practically independent of the polymer concentration (Noskov et al. 2000). Hence the surface properties of the solution were not taken into account. Nonetheless, it would be an interesting study to see how the surface properties of the solution affects the filament stability, and some recent numerical studies (Wee et al. 2020;Kamat et al. 2020;Wee et al. 2021) have already provided an indication of significant effects of surface properties. However, a detailed experimental study investigating the effects of the surface properties is beyond the scope of the present work. A recent study (Walls et al. 2015) also suggests that the viscosity of the outer gas can play a role in the breakup of a liquid filament. However, the focus of the present work is to study the effect of viscoelasticity on the retraction of a liquid filament. In that context, the effect of the outer gas viscosity is also not within the scope of the current study. Therefore, experiments with different surrounding gases were not performed. The schematic of the experimental setup is shown in Fig. 1a (a photograph of the setup has been provided as part of the supplementary information). A microdrop dispenser (AD-K-501, Microdrop Technologies GmbH), with a nozzle inner diameter of 50 µm was used to generate the liquid ligaments. The dispenser consists of a cylindrical piezoacoustic transducer glued around a glass capillary connected to a fluid reservoir. A detailed description of such a dispenser can be found in Dijksman (1984); Fraters et al. (2021). The piezoacoustic element is driven by an electrical pulse supplied from an arbitrary waveform generator (WW1072, Tabor Electronics) and amplified (50×) by a high-voltage amplifier (WMA-300, Falco Systems). A typical driving electrical pulse is trapezoidal in shape (Wijshoff 2010), as shown in Fig. 1b, where V A denotes voltage and t denotes time. The rise and fall times of the pulse are kept constant in the present experiments at 1 µs, while the amplitude (V A,max ) and the pulse width (t w ) were varied in the ranges 50 -75 V and 30 -50 µs, respectively. The lower limits of these ranges are set by the minimum driving required for jetting, while the upper limits are dictated by a bubble entrainment phenomenon associated with meniscus destabilization at strong driving conditions (Fraters et al. 2019(Fraters et al. , 2021. The shape of the supply waveform was verified by an oscilloscope (DPO 4034B, Tektronix). The liquid in the dispenser was supplied from a transparent plastic syringe (5 mL, Becton-Dickinson) connected via a flexible plastic PEEK tubing (Upchurch Scientific). The microdrop dispenser was driven continuously at a drop-on-demand frequency of 100 Hz to minimize selective evaporation from the nozzle, thereby ensuring a constant liquid composition. High-speed imaging of the jetting behavior was performed at 10 5 fps (frames-per-second), with a 600 ns exposure time, by a high-speed camera (Fastcam SA-Z, Photron) connected to a modular microscope (U-ECA and BXFM-F, Olympus) with multiple objectives (5×, 10×, and 20×, MPlanFL N, Olympus). This allowed for a spatial resolution of as low as 1 µm per pixel. The experiments were illuminated by a LED light source (70% intensity, KL 2500 LED, Schott), with the light beam being collimated onto the imaging plane by a collimation lens (Thorlabs). The light source used in the present work is a cold light source, and was only used at 70% of its maximum intensity. Furthermore, the light source would only be switched on during image capturing, and remained switched off when the camera was not capturing any frames (such as during changing or refilling the liquid, and saving the captured frames on the computer hard drive). Thus it is expected that the temperature of the liquid in the device does not change appreciably due to the presence of the light source. The dispenser and the camera were triggered simultaneously (with nanosecond precision) by a programmable pulse-delay generator (BNC 575, Berkley Nucleonics Corp.). The captured images were further analyzed using an OpenCV-based Python script developed in-house. The script utilized a Canny edge detection method to find the edges of the jetted ligament outside the nozzle and that of the retracting meniscus within the nozzle. Knowing the pixel location of the edges, one can find out the top-most and bottom-most point of the jetted ligament, and thus its length (defined as L(t) in Fig. 5a). Similarly, knowing the top-most point of the retracting meniscus within the nozzle, z m (Fig. 3d), is tracked.

Jetting liquids with polymer additives
The jetting behavior with water at V A,max = 50 V, t w = 40 µs is shown in Fig. 1c (and Movie 1 of the supplementary information). Here, and in all subsequent figures, the time t = 0 is defined as the moment when the dispenser is triggered. At that moment, the recording with the camera is also started. A detailed description of the inkjet droplet formation can be found in Hoath (2016); Wijshoff (2010). Briefly, a liquid ligament is formed following pinch-off from the nozzle. As the ligament propagates, the tail retracts due to capillarity, and merges with the head droplet to form a single droplet. When the pulse amplitude (V A,max ) is further increased to 60 V, a longer ligament is initially produced (as seen in Fig. 1d and Movie 1 of the supplementary information); however, in this case, the tail does not retract into the head droplet. Instead, the ligament breaks up into multiple smaller droplets through the Rayleigh-Plateau instability. These smaller droplets in Fig. 1d then merge to form a larger satellite droplet, which never coalesces with the primary droplet owing to the satellite's lower speed.
Next we add the polymer (PEO1M) to the fluid being jetted. The typical jetting behavior for three different concentrations is shown in Fig. 2. In all the three cases shown in the figure, the driving conditions are kept constant at V A,max = 60 V, t w = 40 µs; only the concentration of PEO1M is varied. The jetting behavior with water at this condition is shown in Fig. 2a, reproduced from Fig. 1d. For the 0.003% PEO1M solution ( Fig. 2b and Movie 2 of the supplementary information), the jetting behavior shows a stark difference as compared to that of water. Now the jetted ligament consists again of a spherical head droplet and a long slender tail, similar to the pure water case. Note that the tail here is significantly thinner than the one for the weaker driving case with water ( Fig. 1c). A small spherical tail droplet is also observed in this case, which grows in size as the tail retracts towards the head droplet. When the head and the tail droplet are sufficiently close, they merge to form a single droplet, without any satellite droplets being formed. Thus, the addition of a very small quantity of PEO1M (0.003% by mass) is sufficient to suppress the formation of satellite droplets. Addition of the long chain polymer imparts viscoelasticity to the aqueous solution. The slender tail acts like a stretched filament, being forced by elasticity (and capillarity) to retract its length. This results in the retraction of the whole liquid filament without any intermediate break up. By observing the time stamps, it can also be identified that the pinch-off from the nozzle happens at a later time for the 0.003% PEO1M solution (0.22 ms) as compared to that for water (0.15 ms). This has been reported in Hoath et al. (2014) also, albeit for a different polymer-solvent combination. We note that the suppression of satellite droplet formation by the addition of viscoelasticity comes at the cost of some jetting velocity. For example, the velocity of the jetted 0.003% PEO1M droplet (Fig. 2b) was measured to be 3.11 m/s, as compared to the 3.33 m/s measured for the primary droplet in the pure water case (Fig. 2a). A hypothesis for this decrease in velocity is as follows: In a piezo-acoustic fluid dispenser such as the one used in the present work, the electrical pulse input to the transducer induces a mechanical deformation of the piezo-actuator, resulting in a pressure pulse in the liquid within the dispenser. A detailed description of the associated mechanisms can be found in Dijksman (1984); Wijshoff (2010). From an energetics point of view, the pressure pulse imparts an additional energy in the fluid volume. For a Newtonian liquid, this energy is converted into the surface energy of the jetted droplet and its kinetic energy. However, for a viscoelastic liquid, there is an additional elastic energy associated with stretching the liquid ligament. Since the droplet sizes for the Newtonian and the viscoelastic cases are almost the same, resulting in similar surface energies (as surface tension remains unchanged at such low polymer concentration, as seen in Table 1), the resulting kinetic energy for the viscoelastic droplet is lower. This results in a lower jetting velocity. However, this is only a qualitative understanding of the physical phenomenon, and a detailed quantitative understanding is beyond the scope of the present work.
If the polymer concentration is further increased to 0.009% by mass (still with the same driving), the jetting behavior again changes dramatically, as seen in Fig. 2c (and Movie 2 of the supplementary information). In this case, a filament of the liquid, while still connected to the liquid inside the nozzle, appears downstream of the nozzle exit. However, now the additional polymer content has increased the elasticity to such an extent that the filament does not pinch-off from the nozzle. Instead, it retracts back into the nozzle, thus suppressing jetting altogether. Morrison & Harlen (2010) refer to such observations as the 'bungee jumper' in their numerical study. Obviously, for inkjet printing or similar droplet deposition applications, this phenomenon is undesirable, as no droplet is being produced.
The vastly different jetting behaviors observed in Fig. 2 suggest that, for a given driving condition, there exists a concentration range where jetting without the formation of satellites is observed, and on either side of that range the jetting deteriorates, namely by the occurrence of satellites for lower concentrations, or by suppressing the jetting altogether at high concentrations. Hence it is imperative to delineate these three regimes in an appropriate two-dimensional parameter space, with suitable variables representing the solution concentrations and the driving conditions.
The choice of the control parameter to represent the driving condition is not straightforward, as both pulse amplitude (V A,max ) and pulse width (t w ) are varied in the present experiments. Since the optically transparent microdrop dispenser is operated in the 'pull-push' mode, the motion of the liquid meniscus inside the nozzle, just prior to filament pinch-off, can be observed in order to gauge the effect of the driving condition. The meniscus motion for three different solution concentrations, at the same driving condition, are shown in Fig. 3 (and Movie 3 of the supplementary information). A detailed description of such motion, and how it is affected by changing the driving conditions, can be found in Fraters et al. (2021). The position of the meniscus, z m , is tracked as shown in Fig. 3a, and its temporal evolution is plotted in Fig. 3d. From Fig. 3d, it is observed that the maximum meniscus position, z m,max , is the same for the three different solution concentrations. In other words, z m,max , is only affected by the pulse amplitude and pulse width (Fraters et al. 2021), and not by the amount of polymer present in the solutions used in the present work. Thus, z m,max can be considered to be a suitable variable to represent the effect of the driving conditions. It is non-dimensionalized with the nozzle diameter, d n , asz m = z m,max d n . (3.1) The second dimensionless parameter for the two-dimensional parameter space must be related to the composition of the jetted liquid. The solution concentration is associated  resulting in the nozzle Deborah number, De n , defined as where ρ is the density of the fluid and γ is the surface tension. The regime map in thez m -De n parameter space is shown in Fig. 4. As expected from the observations made in Fig. 2, the desirable 'no satellite' (green markers) regime lies in the middle of the phase space, flanked by the 'satellite' (blue markers) and 'no jetting' (red markers) regimes. For a particular polymer concentration (constant De n ), a stronger driving (higherz m ) is required to overcome the elastic effects and traverse from the 'no  jetting' regime to the 'no satellite' regime. However, if the driving is too strong, the stabilizing effect of elasticity on the retracting tail is lost, resulting in the tail to break, forming satellites.

Filament retraction: experiments
The instantaneous retraction length, L(t), is measured between the extremities of the head droplet and the tail droplet, as shown in the inset of Fig. 5. The temporal evolution of L(t) is plotted in Fig. 5a for different solution concentrations (De n ) and driving conditions (z m ). The retraction length appears to decrease linearly with time even for the viscoelastic liquids (De n = 0) (dashed lines in Fig. 5a); a trend that is expected for Newtonian liquids (De n 0) (Planchette et al. 2019). The retraction velocity that is extracted from the slopes of these curves is approximately constant in time for each dataset shown in Fig. 5a, but different from that of a Newtonian filament (see Fig. 5a), where it is equal to the Taylor-Culick velocity (Keller 1983;Hoepffner & Paré 2013;Pierson et al. 2020): where R 0 is the radius of the retracting filament. In the present experiments, R 0 is measured at a location on the contracting filament close to the head droplet, and at a time instant t = t 0 after pinch-off from the nozzle. In certain experiments, post pinch-off, spurious effects such as out-of-plane oscillations were observed. These were adjudged to be initial transients, and removed from the datasets plotted in Fig. 5a. t 0 thus corresponds to the first datapoint of each dataset in Fig. 5a. The error bars in Fig. 5 represent ± one standard deviation across approximately 25 experiments for each dataset. The ratio of the retraction velocities (v ret ) of the viscoelastic filaments in the current experiments (measured from the linear fits in Fig. 5a) to the corresponding Newtonian Taylor-Culick velocity (v T C ) is plotted in Fig. 5b against the filament Deborah number, defined as In the present experiments, De 0 /De n ∼ O(10). Furthermore, R 0 was observed to vary over a very narrow range (4 -13 µm) in the present experiments, resulting in more than one De 0 value in some cases for a given De n . Figure 5b shows that the viscoelastic filaments have a higher retraction velocity than the Newtonian ones. Moreover, it can be observed from Fig. 5b that for the same De 0 , there are more than one v ret value. This suggests that the relaxation time, λ, is not the only parameter affecting the retraction of these filaments; the retraction is also affected by the driving conditions at which the solution is jetted. To characterize the effect of the driving conditions on v ret , we plot v ret /v T C as a function ofz m for different De n in Fig. 5c. However, no clear trend is observed in the plot. This is probably due to the complex flow profile associated with the motion of the meniscus, and is essentially a limitation of the experimental setup.

Filament retraction: theoretical model
In order to identify the role of viscoelasticity on the retraction velocity of liquid filaments, a simplified theoretical model is proposed. The retraction dynamics of viscoelastic liquid films have been studied in other geometries (Evers et al. 1997;Dalnoki-Veress et al. 1999;Villone et al. 2017;Tammaro et al. 2018;Villone et al. 2019), but not for a slender liquid filament. We follow the lines of Pierson et al. (2020) for Newtonian liquid filaments, but now account for the viscoelasticity due to the polymers. A careful examination of the retraction phenomenon (Fig. 6a) reveals that during the retraction, the size of the head droplet does not change noticeably (< 3%), while the slender tail is pulled towards the head droplet. During this retraction, the spherical tail droplet grows in size as the tail length decreases. This behavior is modeled by the geometry shown in Fig. 6b.
At time t = t 0 , the tail in Fig. 6a is modeled as a long cylindrical filament with a rounded end (shown by the black dashed lines in Fig. 6b), having radius R 0 and length L 0 . We consider an axisymmetric coordinate system that is co-moving with the head droplet (the head droplet moves with a constant velocity in the present experiments), which implies that the fluid within the tail is initially at rest in this coordinate system. One can then proceed via a momentum balance over the control volume that is indicated by the red dashed rectangle. The surface and elastic tensions pull the filament towards the head droplet. For t > t 0 , the rounded edge of the filament therefore starts to retract in the negative z-direction, with an instantaneous length L(t). As the tail retracts, liquid from the tail feeds the tail droplet (Savva & Bush 2009;Pierson et al. 2020), resulting in an increase of its size. Below we will use that most of the liquid momentum is localized inside the tail droplet, as it moves inwards with an instantaneous velocity dL/dt. Note that L(t) is defined slightly differently to the definition used in the experimental data. However, since the diameter of the head droplet does not change appreciably, the following analysis holds well for describing the retraction. We formalize these ideas using the slender jet approximation (R 0 L), within that the mass and momentum conservations can be written as (Eggers 1993;Shi et al. 1994;Eggers & Dupont 1994;Clasen et al. 2006a) Here v(z, t) and r(z, t), respectively are the axial velocity and the filament radius, prime denotes a derivative along z, while η s is the solvent viscosity, and κ is the curvature of the filament, given by The viscoelasticity is accounted for by σ zz and σ rr , the components of the polymer stress tensor σ, for which a separate constitutive equation needs to be specified (Bird et al. 1987). In the slender jet geometry, the predominant stretching and viscoelastic contribution is along the z-direction, so σ rr can be omitted in the remainder (Clasen et al. 2006a;Eggers & Fontelos 2015). In order to perform the control volume analysis, one can bring the slender jet equations to a conservative form by multiplying Eq. (5.1) by ρv and Eq. (5.2) by r 2 , and then adding them up to obtain (Clasen et al. 2006a;Eggers & Villermaux 2008): where K = 1 r(1 + r 2 ) 1/2 + r (1 + r 2 ) 3/2 . (5.5) The right hand side of Eq. (5.4) can readily be integrated from z = 0 to z = L(t), over the control volume in Fig. 6b. We integrate the equation using the same assumptions as in Pierson et al. (2020): (i) r vanishes at z = L(t), as that is the tip of the filament; (ii) the filament radius is uniform (r = R 0 and r = 0) at z = 0, which is at an arbitrary location close to the head droplet; (iii) the fluid is at rest (v = v = 0) at z = 0. Defining the total momentum P = πρ L(t) 0 r 2 vdz, we can then indeed integrate Eq. (5.4) from z = 0 to z = L(t) as Savva & Bush (2009);Pierson et al. (2020): We recover the anticipated momentum balance, with capillary and elastic forces pulling the liquid tail towards the head drop. Now, the retraction process is feeding the tail droplet (Savva & Bush 2009;Pierson et al. 2020), while the fluid between the head and the tail droplets remains at rest (Pierson et al. 2020). Therefore, P is essentially the momentum of the tail droplet P = M T dL/dt (assuming that the fluid velocity inside the droplet is constant (Pierson et al. 2020)). The mass M T (t) of the tail drop increases over time by the mass flow rate −πρR 2 0 dL/dt, such that In the analysis that follows, we omit the initial mass (2πρR 3 0 /3) of the edge of the filament (Pierson et al. 2020), which is negligible in the experimental data with which we compare our theoretical calculations. Integration of Eq. (5.6) in time then gives where we have introduced a change of variable with T = t − t 0 . For a Newtonian fluid (σ zz = 0), this equation can be integrated to (L 0 − L) 2 = γ ρR0 T 2 and one recovers a retraction with a constant (Taylor-Culick) velocity (Keller 1983;Hoepffner & Paré 2013;Pierson et al. 2020). It is clear that the presence of elastic stress will offer an extra contribution that speeds up the retraction, as observed in experiments. Importantly, however, the relaxation of σ zz (T ) will lead to a nonlinear evolution of L(T ), so that the retraction velocity is no longer constant.
To close the problem, we need a constitutive relation for the polymeric stress. Here we use the Oldroyd-B fluid that has been successfully used to describe the thinning of jets (Clasen et al. 2006a;Eggers et al. 2020). In terms of the conformation tensor A, the stress is then given by a constitutive relation (Bird et al. 1987 where η p is the polymer viscosity. In the Oldroyd-B fluid, the conformation tensor evolves by a linear relaxation dynamics, which in the slender jet approximation reads (Fontelos & Li 2004) (5.10) Using the same assumptions used for deriving Eq. (5.6) and using that A zz 1 (large stretching of polymer chains along the axis), Eq. (5.10) can be reduced to The polymer stress then follows as The initial condition A 0 = A zz (z = 0, T = 0) is not determined from the present analysis, but is an inherited condition from the jetting phase, where the polymers are deformed by the stretching flow. We now return to Eq. (5.8) with the polymer stress given by Eq. (5.12), so that (5.13) This can be integrated with the initial condition L(0) = L 0 , to yield the variation of the contracting length with time, given by This is the central result of the analysis. Although the resulting L(T ) is nonlinear, the variation of L with T appears nearly linear (see result plotted in Fig. 7a). To highlight the effect of viscoelasticity, the corresponding Newtonian Taylor-Culick prediction for each dataset is also shown in Fig. 7a by dashed lines. It can be clearly observed from Fig. 7a that viscoelastic retraction (discrete data points from experiments and solid lines from fitting Eq. (5.14)) is faster than the corresponding Newtonian Taylor-Culick retraction. Given the nearly linear appearance of L(T ), it is therefore instructive to expand Eq. (5.14) for early times T /λ 1, which gives (5.15) This illustrates the enhanced retraction velocity v 2 T C + ηpA0 λρ 1/2 during the initial stage. A result of this form can be obtained even more generally, beyond the assumptions underlying the Oldroyd-B fluid. Namely, evaluating the stress integral at short time in Eq. (5.8), one finds the initial retraction velocity incremented by elastic stress that is initially in the filament. Finally, one may use Eq. (5.14) to estimate A 0 in our experiments, which is otherwise difficult a priori. The fitting was performed by matching L 0 (at T = 0) and L from the experiments at the final T instant for each data set, as shown in Fig. 7a by the continuous lines. We rewrite Eq. (5.15) as When plotting the left hand side of Eq. (5.17), expressed as ψ, against T /λ, a straight line through (0,0) is expected at small values of T /λ, with the slope depicting A 1/2 0 . This is observed in Fig. 7b for the experimental data (discrete datapoints), with deviations from linear behavior (continuous lines) observed at larger T /λ values. Now, A 0 = A zz (z = 0, T = 0) is the axial component of the conformation tensor A, which itself is defined as A = X X /X 2 e (Bird et al. 1987), where X is the stretched length of each individual polymer molecule, and X e its equilibrium length. Here, the polymer molecules are pictured to be two spherical beads connected by a spring. Therefore, one can argue that A 1/2 0 is proportional to the stretched length of the polymer molecules at t = t 0 (or T = 0) and z = 0. The stretched length of the polymer is proportional to the local polymer stretching, l , while the initial aspect ratio, L 0 /R 0 , of the filament may be proportional to the stretching of the filament, g . It may be expected that l and g are correlated under strong axial tension for a slender liquid filament. Hence, A 1/2 0 can be assumed to be linearly correlated to L 0 /R 0 . The variation of A 1/2 0 with L 0 /R 0 is plotted in Fig. 7c; the dependence is not inconsistent with the assumed linear behavior. The values of A 1/2 0 obtained from fitting with the experiments are O(100), and one may wonder whether the finite extensibility of polymers (ignored in the Oldroyd-B model) may play a role. Determining the finite extensibility from rheological experiments is a challenge. Lindner et al. (2003) report, for higher molecular weight PEO (2 × 10 6 a.u. and 4 × 10 6 a.u.), maximum polymer stretched lengths in the range of O(10) to O(100). The deviations in A 1/2 0 from the assumed linear trend might therefore be attributed to the limitations of the Oldroyd-B model. We remark, however, that any analysis with nonlinear constitutive relations will come with additional (unknown) fitting parameters. In addition, there are other factors that may play a role in the retraction dynamics such as polydispersity and multiple relaxation time scales of the polymer molecular chains (Entov & Hinch 1997;Wagner et al. 2005), non-uniform radius of the filament along the axial direction, non-axisymmetric effects at the nozzle exit (van der Meulen et al. 2020), and wetting effects at the nozzle exit (Beulen et al. 2007;de Jong et al. 2007), which we have not considered in the present study. Figure 7c also suggests that A 1/2 0 depends on the polymer concentration, which is expected since the higher the polymer concentration, the longer is the relaxation time of the solution, resulting in the ability of the polymers to be stretched longer. Hence, intuitively, one can come to the realization that A 1/2 0 increases with the polymer concentration, thus the polymer relaxation time λ, and subsequently De 0 . This is also observed in Fig. 7d. However, the exact trend for the variation of A 1/2 0 with De 0 is rather difficult to predict in the present experiments. We note that this may be due to the complex flow profile associated with the motion of the meniscus (Fig. 5c), and a simplified approach to estimate the average elongation rate of the polymer macromolecules from the motion of the meniscus itself is unsatisfactory. This can be ascribed to being a limitation of the experimental setup, and more controlled experiments in a simplified geometry can shed light into the elongation rate of the macromolecules (thus A 0 ) prior to the commencement of retraction.
We further clarify the applicability of two of the assumptions used for developing the theoretical model: In the present experiments, the filament radius in Fig. 6a is not uniform but decreases along its distance from the head droplet. However, given the optical resolution of our imaging system, it is extremely difficult to characterize this variation. Hence in the theoretical model in § 5, this variation has not been taken into account. Furthermore, as highlighted by van der Bos et al. (2014), the fluid velocity inside the droplet produced by an inkjet printhead is not constant. However, the primary goal of the present work was to study the retraction of a viscoelastic filament. A simplified model § 5 for that purpose was developed, where it was assumed, for simplicity, that the fluid velocity inside the droplet is constant. Despite these aforementioned simplifications, the model was found to describe the experimentally-observed retraction speeds being faster than the Newtonian Taylor-Culick velocity reasonably well.
The regime map (Fig. 4) indicates a transition from the no satellite to satellite regime at low De n (0 -3) and highz m (4 -6). The exact location of this transition probably can be estimated by comparing the retraction timescale obtained from the theoretical model in the present work with the timescale associated with pinch-off. A starting point for such an analysis could be the recent work by Eggers et al. (2020). Furthermore, the theoretical model described in the present work is based on the Oldroyd-B model, which cannot capture pinch-off. One would need to consider nonlinear models such as FENE-P to suitably describe pinch-off. Such nonlinear models require additional (unknown) fitting parameters. However, a detailed calculation is beyond the scope of the present work, and has been left as a future exercise.

Conclusions and outlook
The present work demonstrates that adding small amounts of long chain polymers to water-based inks can result in the suppression of satellite droplets formation in inkjet printing. These polymers impart viscoelasticity to the liquid being jetted, resulting in stabilization of the slender finite-length filament against a Rayleigh-Plateau instability. Due to the action of both capillarity and viscoelasticity, the tail droplet is pulled towards the head droplet, while being connected by a thin liquid thread, resulting in the formation of a single droplet without any satellites. However, if the polymer concentration is too high, no droplet is jetted as the increased elasticity prevents pinch-off from the dispensing nozzle. A regime map in terms of the experimental operating parameters was identified, thus delineating the 'satellite', 'no satellite', and 'no jetting' regimes. Moreover, it was observed that, in the 'no satellite' regime, the filament essentially retracts linearly with time, with the retraction velocity being higher than the Newtonian Taylor-Culick velocity. A simple theoretical model was developed to derive the retraction velocity for a slender viscoelastic liquid filament, which works reasonably well in modeling the experimental observations. These results are helpful in revealing the complex interplay between inertia, capillarity, and viscoelasticity during the retraction of slender liquid filaments, and are highly useful to predict the behavior during inkjet printing of polymer solutions.
The novelty of the present work does not lie in demonstrating the suppression of satellite droplets by the addition of polymers. Such suppression behavior has been shown before (see Christanti & Walker (2002); Shore & Harrison (2005); Morrison & Harlen (2010); Yan et al. (2011);Hoath et al. (2012Hoath et al. ( , 2014). However, a quantitative understanding of the physical mechanism for the suppression of satellite droplets by the addition of viscoelasticity to the ink has been lacking in the scientific literature. We attempt to address these knowledge gaps through this present manuscript. Systematic experiments have been performed over a range of polymer concentrations and driving conditions that result in observations ranging from satellite droplets to no satellites and ultimately to the suppression of jetting altogether. To the best of our knowledge, a similar investigation has not been hitherto reported. The retraction behavior of viscoelastic filaments reported in this manuscript is also new; in particular the higher retraction speeds as compared to the Newtonian Taylor-Culick velocity. Furthermore, the manuscript presents a simplified theoretical model that explains the apparently linear retraction behavior and shows reasonable agreement with the experimentally-observed retraction speeds being faster than the Newtonian Taylor-Culick velocity. Until now there has not yet been such a comparison between experiment and theory, which we provide, though the complexity of the problem unfortunately does not allow for a quantitative one-to-one comparison.