## Viscosity

### Quick Start

We have three apps relating to viscosity.

- Calculating the viscosity of a mixture
- Calculating the change of viscosity with polymer concentration
- Calculating how viscosity depends on emulsion or particle concentration

Explore whichever ones interest you!

### Viscosity Mix

_{1}cP

_{1}g/cc

_{2}cP

_{2}g/cc

_{mix}R : η

_{mix}K-M

The first topic is the apparently simple question of what happens if you mix two solutions of different *dynamic* viscosities η. There are two methods, neither of which has much backing, based on the *kinematic* viscosities `ν=η/ρ`. The first is the Refutas equation. The procedure is to calculate Viscosity Blending Numbers, VBN, for each component, then create a VBN for the mix based on the mass fractions, x_{i}, then convert the mix VBN back to viscosity.

`VBN_i = 14.534.ln(ln(η+0.8))+10.975`

`VBN_(mix) = x_1.VBN_1 + x_2.VBN_2`

`ν_R = exp(exp((VBN_( mix)-10.975)/14.534))-0.8`

The second method is Kendall-Munroe.

`ν_R^(1/3) = x_1ν_1^(1/3) + x_2ν_2^(1/3)`

Given that the predicted values can differ greatly, it is hard to know what to do with these values. Just pick whichever seems to work for you.

Up to May 2022 the app used only *dynamic* viscosities till my error was pointed out Michał Kawa and the app updated.

### Viscosity dependence on polymer concentration

### Concentration Dependence

_{solv}cP

_{Entangled}g/l

_{Max}g/l

_{g}nm

Next we need to know what sorts of things happen when we change the concentration of a polymer in our coating/printing solution. The approach here is comparatively user friendly. A more sophisticated version is in Polymer Viscosity app in Practical Solubility.

The complicated behaviour is simplified in the app. You have the viscosity of the base solvent. You then set a concentration, C_{E} at which the polymer chains become fully entangled - at this point the viscosity increases rapidly. You select a concentration range of interest. Finally, you need to provide R_{g}, the radius of gyration of your polymer which you can estimate from a number, N, of pseudo bonds and a pseudo bond length (Kuhn length), b.

The viscosity increase depends on a pseudo viscosity which in turn depends on R_{g}³, then there is a Rousean-like 1.3 power dependency on concentration, C, along with a cubic dependency on C/C_{E}. With fudge factors A and B we have the inelegant formula:

η = η_{0}A(B.C.R_{g}³)^{1.3}(C/C_{E})³

The point of the app is to get a feel for what happens as you change your polymer, its concentration and the solvent used (good solvents give larger R_{g} values than bad ones).

If your polymer concentration is relatively small and you "only" increase the viscosity to a few 10's of cP this might seem harmless. But if you also have a dispersion then the effect of this "harmless" increase might be disasterous, as discussed below the Dispersion Viscosity app.

### Viscosity dependence on emulsion/particle concentration

### Dispersion Viscosity

_{0}cP

_{dispersed}cP

_{max}

_{max}

Finally, if we have an emulsion or fillers or pigments, how does the viscosity depend on the volume fraction? Actually, the question is about "low shear" viscosity as the high-shear viscosity decreases rapidly with shear as shown in the Particle Rheology app.

The overall viscosity depends on the initial viscosity of the bulk phase, η_{0} and the volume fraction of the dispersed component, φ. Many equations are used to describe this dependence and three are shown here. The Einstein equation obviously is the standard for solid particles.

η=η_{0}(1+2.5φ)

However, for this is the limiting factor when k, which is the ratio of the viscosity of the dispersed material to that of the bulk phase, η_{dispersed}/η_{0} reaches a very high value. For more modest k values, the Taylor equation is used:

η=η_{0}(1+φ(5k+2)/(2(k+1)))

The Dougherty-Krieger formula is popular but tends to overshoot for emulsions at high φ values. It incorporates a "close packing" fraction φ_{c} at which the viscosity becomes infinite. This is OK for solids but less suited for emulsions. In the app φ_{c} is assumed to be 0.74. The factor of 2.5 is the default value for the "intrinsic viscosity" term. Because the viscosity rises exceptionally above φ=0.61 (near the random close packing limit) the calculation is not performed above that value.

η=η_{0}(1-φ/φ_{c})^{-2.5φc}

Reading some early Krieger literature it seems that this curve is for viscosities at reasonably high shear stress. At low stress the curve behaves as if φ_{c} is ~0.58, so viscosities take off much faster as φ increases.

The most realistic, at least for emulsions, is also the most monstrous - the Yaron, Gal-Or. This uses the k from the Taylor model and instead of φ uses λ=φ^{0.333}.

η=η_{0}(1+5.5φ[4λ^{7}+10-(84/11)λ^{2}+(4/k)(1-λ^{7} )]/(10(1-λ^{10} )-25λ^{3}(1-λ^{4})+(10/k)(1-λ^{3} )(1-λ^{7})))

Another choice, based on an extensive evaluation of real emulsion viscosities (rather than dispersions) with very large ranges of k is from Pal^{1}:

η=η_{0}η_{r}

where η_{r} is given by the following with φ_{m}=0.637 (this is a general value which in principle can be fitted. For simplicity, it is fixed in the app and the plot stops at φ=0.6 before the curve gets too large):

η_{r}[(2η_{r}+5k)/(2+5k)]^{1.5}=(1-φ/φ_{m})^{-2.5φm}

None of the equations has any dependence on dispersed particle size. This is surprising but supported by the data.

All you have to do is choose a range over which to plot (limited by φ_{max}), enter η_{0} and, for Taylor and for Yaron, Gal-Or, η_{dispersed}. You can read out viscosity values with the mouse. For those with solid dispersed particles, set η_{dispersed} to its maximum value. As you can easily see, above a few 1000 cP the exact values doesn't matter.

As a quick check, the Wt% at your maximum φ is calculated from the density of your particle relative to the medium.

The most important point from all this is that apparently small changes in the base viscosity can make a huge difference to the dispersion's viscosity. Suppose the base viscosity is 1cP and at your chosen φ the viscosity is 100cP. Now increase the base viscosity by 4cP. This is trivial. But not for the dispersion which now has a viscosity of 500cP, not 104cP! My own experience out in the world of coating machines is that this effect of apparently trivial changes in the base viscosity is largely unknown.

^{1}Rajinder Pal,

*Novel viscosity equations for emulsions of two immiscible liquids*, J. Rheol., 4, 509-520, 2001