Fluid thread breakup

Fluid thread breakup is the process by which a single mass of fluid breaks into several smaller fluid masses. The process is characterized by the elongation of the fluid mass forming thin, thread-like regions between larger nodules of fluid. The thread-like regions continue to thin until they break, forming individual droplets of fluid.

Thread breakup occurs where two fluids or a fluid in a vacuum form a free surface with surface energy. If more surface area is present than the minimum required to contain the volume of fluid, the system has an excess of surface energy. A system not at the minimum energy state will attempt to rearrange so as to move toward the lower energy state, leading to the breakup of the fluid into smaller masses to minimize the system surface energy by reducing the surface area. The exact outcome of the thread breakup process is dependent on the surface tension, viscosity, density, and diameter of the thread undergoing breakup.

History
The examination of droplet formation has a long history, first traceable to the work of Leonardo da Vinci who wrote:

"'How water has tenacity in itself and cohesion between its particles. […] This is seen in the process of a drop becoming detached from the remainder, this remainder being stretched out as far as it can through the weight of the drop which is extending it; and after the drop has been severed from this mass the mass returns upwards with a movement contrary to the nature of heavy things.'"

He thus correctly attributed the fall of droplets to gravity and the mechanism which drives thread breakup to the cohesion of water molecules.

The first correct analysis of fluid thread breakup was determined qualitatively by Thomas Young and mathematically by Pierre-Simon Laplace between 1804 and 1805. They correctly attributed the driver of thread breakup to surface tension properties. Moreover, they also deduced the importance of mean curvature in the creation of excess pressure in the fluid thread. Through their analysis, they showed that surface tension can behave in two ways: an elastic mechanism that can support a hanging droplet and a pressure mechanism due to capillary pressure that promotes thread breakup.

In the 1820s, the Italian physicist and hydraulic engineer Giorgio Bidone studied the deformation of jets of water issuing from orifices of various shapes. Félix Savart followed in 1833 with experimental work, utilizing the stroboscopic technique to quantitatively measure thread breakup. He noted that breakup is a spontaneous process, occurring without an external stimuli. This work allowed him to determine that droplets are produced from a jet flowing from a tank at a distinct rate inversely proportional to the nozzle radius and proportional to pressure in the tank. These observations facilitated Joseph Plateau's work that established the relationship between jet breakup and surface energy. Plateau was able to determine the most unstable disturbance wavelength on the fluid thread, which was later revised by Lord Rayleigh to account for jet dynamics.

As the surface disturbance becomes large, non-linear theory must be applied. The behavior of jets with large disturbances was examined experimentally by Magnus and Lenard. Their experiments helped to characterize satellite droplets, droplets that are produced in addition to the large main droplet, through the introduction of high speed photography. High speed photography is now the standard method for experimentally analyzing thread breakup.

With the advent of greater computational power, numerical simulations have begun to replace experimental efforts as the chief means of understanding fluid breakup. However, difficulty remains in accurately tracking the free surface of many liquids due to its complex behavior. The most success has occurred with fluids of low and high viscosity where the boundary integral method can be employed as the Green's function for both cases is known. Dommermuth and Yue characterized irrotational, inviscid flow by this method as did Schulkes. Youngren and Acrivos considered the behavior of a bubble in a high viscosity liquid. Stone and Leal expanded this initial work to consider the dynamics of individual drops. For fluids of middling viscosity, full simulations using the Navier-Stokes equations are required with methods determining the free surface such as level-set and volume of fluid. The earliest work with full Navier-Stokes simulations was done by Fromm which focused on inkjet technology. Such simulations remain an active area of research.

Physical mechanism of thread breakup


The breakup process in a fluid thread or jet begins with the development of small perturbations on the free surface of the fluid. This is known as the linear theory of fluid thread breakup. These perturbations are always present and can be generated by numerous sources including vibrations of the fluid container or non-uniformity in the shear stress on the free surface. In general, these disturbances take an arbitrary form and are thus difficult to consider rigorously. It is therefore helpful to take a Fourier transform of the disturbances to decompose the arbitrary disturbances into perturbations of various single wavelengths on the surface of the thread. In doing so, this allows one to determine which wavelengths of the disturbance will grow and which will decay in time.

The growth and decay of wavelengths can be determined by examining the change in pressure a perturbation wavelength imposes on the interior of the fluid thread. Changes to the internal pressure of the thread are induced by capillary pressure as the free surface of the thread deforms. Capillary pressure is a function of the mean curvature of the interface at a given location at the surface, meaning the pressure is dependent on the two radii of curvature that give the shape of the surface. Within the thinned area of a fluid thread undergoing breakup, the first radius of curvature is smaller than the radius of curvature in the thickened area, leading to a pressure gradient that would tend to force liquid from the thinned to thickened areas. However, the second radius of curvature remains important to the breakup process. For some perturbation wavelengths, the effect of the second radius of curvature can overcome the pressure effect of the first radius of curvature, inducing a larger pressure in the thickened regions than the thinned regions. This would push fluid back toward the thinned regions and tend to return the thread to its original, undisturbed shape. However, for other perturbation wavelengths, the capillary pressure induced by the second radius of curvature will reinforce that of the first radius of curvature. This will drive fluid from the thinned to the thickened regions and further promote thread breakup.



The wavelength of the perturbation is therefore the critical parameter in determining whether a given fluid thread will breakup into smaller masses of fluid. Rigorous mathematical examination of the perturbation wavelengths can lead to a relation showing which wavelengths are stable for a given thread as well as which perturbation wavelengths will grow most rapidly. The size of the fluid masses resulting from the breakup of a fluid thread can be approximated by the wavelengths of the perturbation that grow most rapidly.

Non-linear behavior
While linear theory is useful in considering the growth of small disturbances on the free surface, when the disturbances grow to have a significant amplitude, non-linear effects begin to dominate breakup behavior. The non-linear behavior of the thread governs its final breakup and ultimately determines the final shape and number of the resulting fluid masses.

Nonlinearity is captured through the use of self-similarity. Self-similarity assumes that the behavior of the fluid thread as the radius approaches zero is the same as the behavior of the fluid thread when it has some finite radius. Detailed understanding of non-linear thread behavior requires the use of asymptotic expansions to generate the appropriate scaling behavior. Numerous solutions have been found for the non-linear behavior of fluid threads based on the forces that are relevant in particular circumstances.

Important parameters
How a fluid thread or jet undergoes breakup is governed by several parameters among which are the Reynolds number, the Weber number, Ohnesorge number, and the disturbance wavelength. While these numbers are common in fluid mechanics, the parameters selected as scales must be appropriate to thread breakup. The length scale most often selected is the radius of the fluid thread, while the velocity is most often taken to be the velocity of the bulk fluid motion. However, these scales can change based on the characteristics of the considered problem.

The Reynolds number is the ratio between inertia and viscous effects within the thread. For large Reynolds numbers, the effects of motion of the thread are much greater than viscous dissipation. Viscosity only has a minimal damping effect on the thread. For small Reynolds numbers, viscous dissipation is large and any disturbances are rapidly damped from the thread.

The Weber number is the ratio between inertia and surface tension effects within the thread. When the Weber number is large, the inertia of the thread is large which resists the tendency of surface tension to flatten bent surfaces. For small Weber numbers, the changes in the capillary pressure due to the surface disturbances is large and surface tension dominates thread behavior.

The Ohnesorge number is the ratio between viscous and surface tension effects within the thread. As it eliminates the effects of inertia and the need for a velocity scale, it is oftentimes more convenient to express scaling relationships in terms of the Ohnesorge number rather than the Reynolds and Weber number individually.

The perturbation wavelength is the characteristic length of the disturbance on the surface of the jet, assuming that any arbitrary disturbance can be decomposed via a Fourier transform into its constitutive components. The wavelength of the perturbation is critical in determining if a particular disturbance will grow or decay in time.

Linear stability of inviscid liquids
The linear stability of low viscosity liquids was first derived by Plateau in 1873. However, his solution has become known as the Rayleigh-Plateau instability due to the extension of the theory by Lord Rayleigh to include fluids with viscosity. Rayleigh-Plateau instability is often used as an introductory case to hydrodynamic stability as well as perturbation analysis.

Plateau considered the stability of a thread of fluid when only inertial and surface tension effects were present. By decomposing an arbitrary disturbance on the free surface into its constitutive harmonics/wavelengths, he was able to derive the a condition for the stability of the jet in terms of the perturbation:


 * $$\omega^2 = \frac{\sigma k}{\rho a^2} \frac{I_1(ka)}{I_0(ka)} \left ( 1 - k^2 a^2 \right ),$$

where ω is the growth rate of the perturbation, σ is the surface tension of the fluids, k is the wavenumber of perturbation, ρ is the fluid density, a is the initial radius of the unperturbed fluid, and I is the modified Bessel function of the first kind. By computing the growth rate as a function of wavenumber, one can determine that the fastest growing disturbance wavelength occurs at:


 * $$ \lambda_\text{max} \approx 9.02a.$$

The wavelength of maximum instability increases as the radius of the fluid thread increases. As importantly, unstable modes are only possible when:


 * $$ ka < 1.$$

Linear stability of viscous liquids
Reynolds and later Tomotika extended Plateau's work to consider the linear stability of viscous threads. Rayleigh solved for the stability of a viscous thread of viscosity $$ \mu_A $$ without the presence of an external fluid. Tomokita solved for the stability of a fluid thread in the presence of an external fluid with its own viscosity $$ \mu_B $$. He considered three cases where the viscosity of the fluid thread was much greater than the external environment, the viscosity of the external environment was much greater than the fluid thread, and the general case where the liquids are of arbitrary viscosity.

Fluid thread highly viscous
For the limiting case where the fluid thread is much more viscous than the external environment, the viscosity of the external environment falls from the growth rate completely. The growth rate thus becomes only a function of the initial radius of thread, the perturbation wavelength, the surface tension of the thread, and the thread viscosity.


 * $$\omega = \frac{\sigma \left ( k^2 a^2-1 \right )}{2a \mu_A} \frac{1}{k^2 a^2 + 1 - k^2 a^2 I_0^2(ka) / I_1^2(ka)}$$

Plotting this, one finds that the longest wavelengths are the most unstable. As importantly, one can note that the viscosity of the fluid thread does not influence which wavelengths will be stable. Viscosity only acts to decrease how rapidly a given perturbation will grow or decay with time.

Examples of when this case would apply are when almost any liquid undergoes thread/jet breakup in an air environment.

External fluid highly viscous
For the limiting case where the external environment of the fluid thread is much more viscous than the thread itself, the viscosity of the fluid thread falls from the perturbation growth rate completely. The growth rate thus becomes only a function of the initial radius of the thread, the perturbation wavelength, the surface tension of the thread, the viscosity of the external environment, and the second order Bessel functions of the second kind.


 * $$ \omega = \frac{\sigma \left (1 - k^2 a^2 \right )}{2a \mu_B} \frac{1}{k^2 a^2 + 1 - k^2 a^2 K_0^2(ka) / K_1^2(ka)}$$

If one were to plot the growth rate as a function of the perturbation wavelength, one would find that the most unstable wavelengths again occur at the longest wavelengths and that the viscosity of the external environment would only act to decrease how rapidly a perturbation would grow or decay in time.

Examples of when this case would apply are when gas bubbles enter a liquid or when water falls into honey.

General case - arbitrary viscosity ratio
The general case for two viscous fluids is much more difficult to solve directly. Tomotika expressed his solution as:


 * $$\omega = \frac{\sigma \left ( 1 - k^2 a^2 \right ) }{2a \mu_B} \Phi(ka) $$

where $$ \Phi $$ was defined as:


 * $$ \begin{align}

\Phi(ka) ={} &\frac{N(ka)}{D(ka)} \\ N(ka)   ={} &I_1 (ka) \Delta_1 - \{ka I_0 (ka) - I_1 (ka)\} \Delta_2 \\ D (ka)  ={} &\frac{\mu_A}{\mu_B} \{ka I_0 (ka) - I_1 (ka)\} \Delta_1 - \frac{\mu_A}{\mu_B} \left \{\left (k^2a^2 +1 \right ) I_1(ka) - ka I_0(ka) \right \} \Delta_2 \\ &- \{kaK_0 (ka) + K_1(ka)\} \Delta_3 - \left \{ \left ( k^2 a^2 +1 \right ) K_1 (ka) + ka K_0 (ka) \right \} \Delta_4 \end{align}$$

The $$ \Delta $$ coefficients are most easily expressed as the determinants of the following matrices:


 * $$\begin{align}

\Delta_1 &= \begin{vmatrix} kaI_0 (ka) - I_1(ka) & K_1(ka) & -kaK_0(ka) - K_1(ka) \\ I_0(ka) + kaI_1(ka) & -K_0(ka) & -K_0(ka) + kaK_1(ka) \\ \frac{\mu_A}{\mu_B}kaI_0 (ka) & K_1(ka) & -kaK_0(ka) \end{vmatrix} \\[3pt] \Delta_2 &= \begin{vmatrix} I_1(ka) & K_1(ka) & -kaK_0 (ka) - K_1(ka) \\ I_0(ka) & - K_0(ka) & -K_0(ka) + ka K_1(ka) \\ \frac{\mu_A}{\mu_B}I_1(ka) & K_1(ka) & -kaK_0(ka) \end{vmatrix} \\[3pt] \Delta_3 &= \begin{vmatrix} I_1(ka) & ka I_0(ka) - I_1(ka) & -ka K_0(ka) - K_1(ka) \\ I_0(ka) & I_0(ka) + kaI_1(ka) & -K_0(ka) + ka K_1(ka) \\ \frac{\mu_A}{\mu_B} I_1(ka) & \frac{\mu_A}{\mu_B}kaI_0(ka) & -ka K_0(ka) \end{vmatrix} \\[3pt] \Delta_4 &= \begin{vmatrix} I_1(ka) & ka I_0(ka) - I_1(ka) & K_1(ka) \\ I_0(ka) & I_0(ka) + kaI_1(ka) & -K_0(ka) \\ \frac{\mu_A}{\mu_B}I_1 (ka) & \frac{\mu_a}{\mu_B} ka I_0(ka) & K_1(ka) \end{vmatrix} \end{align}$$

The resulting solution remains a function of both the thread and external environment viscosities as well as the perturbation wavelength. The most unstable combination of viscosities and perturbation occurs when $$ \mu_A/\mu_B \approx 0.28 $$ with $$ \lambda \approx 10.66a $$.

For most applications, use of the general case is unnecessary as the two fluids in question have significantly different viscosities which permits the use of one of the limiting cases. However, some instances such as the mixing of oils or oils and water may require the use of the general case.

Satellite drop formation


Satellite drops, also known as secondary droplets, are the drops produced during the thread breakup process in addition to the large main droplet. The drops result when the filament by which the main droplet in hanging from the larger fluid mass itself breaks off from the fluid mass. The fluid contained in the filament can stay as a single mass or breakup due to the recoil disturbances imposed on it by the separation of the main droplet. While the production of satellite droplets can be predicted based on fluid properties, their precise location and volume cannot be predicted.

In general, secondary droplets are an unwanted phenomenon, especially in applications where precise deposition of droplets is important. The production of satellite droplets is governed by the non-linear dynamics of the problem near the final stages of thread breakup.

Examples


Numerous examples of the breakup of fluid threads exist in daily life. It is one of the most common fluid mechanics phenomena one experiences and as such most give the process little thought.

Flow from a faucet
Dripping water is an everyday occurrence. As water leaves the faucet, the filament attached to the faucet begins to neck down, eventually to the point that the main droplet detaches from the surface. The filament cannot retract sufficiently rapidly to the faucet to prevent breakup and thus disintegrates into several small satellite drops.

Air bubbles
Air bubbles are another common breakup phenomenon. As air enters a tank of liquid, like a fish tank, the thread again necks at the base to produce a bubble. Blowing bubbles from a straw into a glass behaves in much the same manner.

Pitch drop experiment
The pitch drop experiment is a famous fluid breakup experiment using high viscous tar pitch. The rate of breakup is slowed to such a degree that only 11 drops have fallen since 1927.

Drops of honey
Honey is sufficiently viscous that the surface perturbations that lead to breakup are almost fully damped from honey threads. This results in the production of long filaments of honey rather than individual droplets.