Introduction & Context

For pseudoplastic (shear-thinning) fluids the viscosity measured in a Brookfield or cup test is meaningless at the high shear rates created by an impeller. The impeller instead “sees” an apparent viscosity that is a function of the local shear rate. Correcting the Reynolds number with this apparent viscosity lets the process engineer select the proper Power number and therefore size the motor, predict heat-of-agitation, and avoid under- or over-designing industrial fermenters, ketchup kettles, or slurry storage tanks.

Methodology & Formulas

All relations below are algebraic duplicates of the Python code; no empirical constants beyond those embedded in the code are introduced.

  1. Rotational speed conversion
    \(N=\dfrac{N_{\text{rpm}}}{60}\)
  2. Representative shear rate generated by a Rushton turbine
    \(\dot\gamma = k_s\,N\)
  3. Apparent viscosity for a power-law fluid (with consistency index \(K\) and flow index \(n\))
    \(\mu_{\text{app}}=K\bigl(\dot\gamma\bigr)^{\,n-1}\)
  4. Impeller Reynolds number corrected for shear-thinning behaviour
    \[\text{Re}' = \frac{\rho\,N\,D^{2}}{\mu_{\text{app}}}\]
  5. Mechanical power draw
    \[P=N_p\;\rho\;N^{3}\;D^{5}\]
Validation regime limits (from code)
Parameter Acceptable range Consequence if violated
Power-law index \(n\) \(0.2 < n < 1\) Code raises ValueError
‘Flat’ Power-number assumption (turbulent regime) \(\text{Re}' \geq 10,000\) Code raises ValueError; Power number may vary below this threshold

Units used in every expression are:

Consistent engineering units
Quantity Unit
\(K\)Pa·sn
\(n\)
\(N\)rev s-1
\(\rho\)kg m-3
\(D\)m
\(\mu_{\text{app}}\)Pa·s
\(P\)W

Constants embedded in the algorithm are:
\(k_s=11.5\) (Rushton turbine, fully-baffled standard tank)
\(N_p=5.5\) (turbulent regime flat value, valid only for \(\text{Re}' > 10,000\))