Visco-elastic jets

{{Citations needed|date=January 2025}}

File:Saliva Thread.jpg 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

{{multiple image

| width =

| image1 = Draining.jpg

| caption1 = Drop draining

| image2 = Merging.jpg

| caption2 = Drop merging

| image3 = Collision1.jpg

| caption3 = Drop collision

| total_width = 300

}}

{{multiple image

| width =

| image1 = Oscillation.jpg

| image2 = Oscillation1.jpg

| footer = Drop oscillation {{br}}(right image is a continuation of left image)

| total_width = 200

}}

A jet of a Newtonian fluid, such as honey poured from a bottle, thins continuously and coils regularly.{{Cite web |last=McKinley |first=Gareth |date=Nov 18, 2013 |title=Viscoelastic Jet |url=https://web.mit.edu/nnf/research/phenomena/viscoelastic_jet.html |access-date=2025-01-24 |website=Gareth McKinley's Non-Newtonian Fluid Dynamics Research Group |publisher=MIT}} 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 draining{{Snd}}a small bead between two beads shrinks as its fluid particles move towards the adjacent beads ("drains away");
  • drop merging{{Snd}}a smaller bead and a larger bead move close to each other and merge to form a single bead;
  • drop collision{{Snd}}a moving bead collides and combines with an adjacent bead;
  • drop oscillation{{Snd}}two 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.{{Needs update|date=January 2025|reason=Site claims "This research is in the process of being published" as of 2013. The resultant published paper, if any, would be a more reliable (and informative) source.}}

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 ({{Mvar|Oh}})

: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 ({{Mvar|De}}), 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:

{{NumBlk|:|\frac{\partial \ R}{\partial t} + \frac{\partial \ vR^2}{\partial z}=0|{{EquationRef|1}}}}

{{NumBlk|:|\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 (R^2\frac{\partial v}{\partial z})}{\partial z} + \frac{\frac{1}{R^2}\partial (R^2(\sigma_{zz}-\sigma_{rr}))}{\partial z}|{{EquationRef|2}}}}

{{NumBlk|:|\kappa=\frac{1}{R(1+R_z^2)^\frac{1}{2}}-\frac{R_{zz}}{(1+R_{zz}^2)^\frac{3}{2}}|{{EquationRef|3}}}}

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 ({{EquationNote|1}}) represents mass conservation, and Equation ({{EquationNote|2}}) 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.

References