Visco-elastic jets

Saliva exhibits viscoelastic "beads-on-a-string" structure.

In fluid mechanics, a viscoelastic jet is a projected stream (jet) of a viscoelastic fluid (a fluid that disobeys Newton's law of viscosity). A viscoelastic fluid returns to its original shape after the applied stress is released.

Free surface continuous jets of viscoelastic fluids are relevant in engineering applications involving blood, paints, adhesives, and foodstuff as well as in industrial processes like fiber spinning, bottle filling, and oil drilling. In process engineering, it is essential to understand the instabilities a jet undergoes due to changes in fluid parameters like the Reynolds number (Template:Math) or Deborah number (Template:Math). With the advent of microfluidics, an understanding of the jetting properties of non-Newtonian fluids becomes essential from micro- to macro-length scales, and from low to high Reynolds numbers.

Description

A jet of a Newtonian fluid, such as honey poured from a bottle, thins continuously and coils regularly. In contrast, a viscoelastic jet breaks up much more slowly. Typically, it evolves into a "beads-on-a-string" structure, where large drops are connected by thin threads. The slow breakup process provides the viscoelastic jet sufficient time to exhibit other phenomena, including:

• drop drainingTemplate:Snda small bead between two beads shrinks as its fluid particles move towards the adjacent beads ("drains away");
• drop mergingTemplate:Snda smaller bead and a larger bead move close to each other and merge to form a single bead;
• drop collisionTemplate:Snda moving bead collides and combines with an adjacent bead;
• drop oscillationTemplate:Sndtwo adjacent beads start oscillating, their separation gradually decreases, and they eventually merge to form a single bead.

The behaviors of non-Newtonian fluids result from the interplay of non-Newtonian properties (e.g. viscoelasticity, shear-thinning) with gravitational, viscous, and inertial effects.Template:Needs update

The evolution of a viscoelastic fluid thread over time depends on the relative magnitude of the viscous, inertial, and elastic stresses and the capillary pressure. To study the inertio-elasto-capillary balance for a jet, two dimensionless parameters are defined: the Ohnesorge number (Template:Math) $$\mathrm{O}\mathrm{h}=\frac{{\eta }_0}{\sqrt[]{\rho \gamma R_0}}$$ which is the inverse of the Reynolds number based on a characteristic capillary velocity $\frac{\gamma }{{\eta }_0};$ and the intrinsic Deborah number (Template:Math), defined as $$\mathrm{D}\mathrm{e}=\frac{\lambda }{t_r}=\frac{\lambda }{\sqrt{\rho R_0^3/\gamma }}$$ where

Template:Mvar is the "Rayleigh time scale" for inertio-capillary breakup of an inviscid jet;

Template:Mvar is the fluid density;

Template:Math is the fluid zero shear viscosity;

Template:Mvar is the surface tension;

Template:Math is the initial radius of the jet;

Template:Mvar is the relaxation time associated with the polymer solution.

Equations

Like other fluids, when considering viscoelastic flows, the velocity, pressure, and stress must satisfy equations of mass and momentum, supplemented with a constitutive equation involving the velocity and stress.

The behaviors of weakly viscoelastic jets can be described by the following set of mathematical equations, with the first representing mass conservation, and the second representing the momentum equation in one dimension:

$$\begin{array}{cc}& \frac{\partial R}{\partial t}+\frac{\partial v{R}^2}{\partial z}=0\\ & \rho (\frac{\partial v}{\partial t}+\frac{v\partial }{\partial z})=-\gamma \frac{\partial \kappa }{\partial t}+\frac{3{\eta }_s}{R^2}\frac{\partial \left(R^2\frac{\partial v}{\partial z}\right)}{\partial z}+\frac{1}{R^2}\frac{\partial (R^2({\sigma }_{zz}-{\sigma }_{rr}))}{\partial z}\\ & \kappa =\frac{1}{R(1+R_z^2{)}^{\frac{1}{2}}}-\frac{R_{zz}}{(1+R_{zz}^2{)}^{\frac{3}{2}}}\end{array}$$

where

Template:Math is the axial velocity;

Template:Mvar and Template:Mvar are the solvent and polymer contribution to the total viscosity, respectively (total viscosity Template:Math);

Template:Mvar indicates the partial derivative $\frac{\partial R}{\partial z}$;

Template:Mvar and Template:Mvar are the diagonal terms of the extra-stress tensor, calculated as follows:$$\begin{array}{cc}{\sigma }_{zz}+\lambda (\frac{\partial {\sigma }_{zz}}{\partial t}+v\frac{\partial {\sigma }_{zz}}{\partial z}-2\frac{\partial v}{\partial z}{\sigma }_{zz})+\frac{\alpha \lambda }{{\eta }_p}{\sigma }_{zz}^2& =2{\eta }_p\frac{\partial v}{\partial z}\\ {\sigma }_{rr}+\lambda (\frac{\partial {\sigma }_{rr}}{\partial t}+v\frac{\partial {\sigma }_{rr}}{\partial z}+\phantom{2}\frac{\partial v}{\partial z}{\sigma }_{rr})+\frac{\alpha \lambda }{{\eta }_p}{\sigma }_{rr}^2& =-{\eta }_p\frac{\partial v}{\partial z}\end{array}$$

where Template:Mvar is the relaxation time of the liquid, and

Template:Mvar is the mobility factor, a positive dimensionless parameter corresponding to the anisotropy of the hydrodynamic drag on the polymer molecules.