Visco-elastic jets

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

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 Reynolds number or Deborah number. 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:Mvar)

$Oh=\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:Mvar), defined as

$De=\frac{\lambda }{\sqrt[]{\rho R_0^3/\gamma }}$

where $t_r=\sqrt[]{\rho R_0^3/\gamma }$ is the "Rayleigh time scale" for inertio-capillary breakup of an inviscid jet. In these expressions, $\rho$ is the fluid density, ${\eta }_0$ is the fluid zero shear viscosity, $\gamma$ is the surface tension, $R_0$ is the initial radius of the jet, and $\lambda$ 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: Template:NumBlk Template:NumBlk Template:NumBlk where $(z,t)$ is the axial velocity; ${\eta }_s$ and ${\eta }_p$ are the solvent and polymer contribution to the total viscosity, respectively (total viscosity ${\eta }_0={\eta }_s+{\eta }_p$); $R_z$ indicates the partial derivative $\frac{\partial R}{\partial z}$; and ${\sigma }_{zz}$ and ${\sigma }_{rr}$ are the diagonal terms of the extra-stress tensor. Equation (Template:EquationNote) represents mass conservation, and Equation (Template:EquationNote) represents the momentum equation in one dimension. The extra-stress tensors ${\sigma }_{zz}$ and ${\sigma }_{rr}$ can be calculated as follows:

${\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}+\frac{\partial v}{\partial z}{\sigma }_{rr})+\frac{\alpha \lambda }{{\eta }_p}{\sigma }_{rr}^2=-{\eta }_p\frac{\partial v}{\partial z}$

where $\lambda$ is the relaxation time of the liquid, and $\alpha$ is the mobility factor, a positive dimensionless parameter corresponding to the anisotropy of the hydrodynamic drag on the polymer molecules.