Bingham plastic

Figure 1 shows a graph of the behaviour of an ordinary viscous (or Newtonian) fluid in red, for example in a pipe. If the pressure at one end of a pipe is increased this produces a stress on the fluid tending to make it move (called the shear stress) and the volumetric flow rate increases proportionally. However, for a Bingham Plastic fluid (in blue), stress can be applied but it will not flow until a certain value, the yield stress, is reached. Beyond this point the flow rate increases steadily with increasing shear stress. This is roughly the way in which Bingham presented his observation, in an experimental study of paints.^[3] These properties allow a Bingham plastic to have a textured surface with peaks and ridges instead of a featureless surface like a Newtonian fluid.

Figure 2 shows the way in which it is normally presented currently.^[2] The graph shows shear stress on the vertical axis and shear rate on the horizontal one. (Volumetric flow rate depends on the size of the pipe, shear rate is a measure of how the velocity changes with distance. It is proportional to flow rate, but does not depend on pipe size.) As before, the Newtonian fluid flows and gives a shear rate for any finite value of shear stress. However, the Bingham plastic again does not exhibit any shear rate (no flow and thus no velocity) until a certain stress is achieved. For the Newtonian fluid the slope of this line is the viscosity, which is the only parameter needed to describe its flow. By contrast, the Bingham plastic requires two parameters, the yield stress and the slope of the line, known as the plastic viscosity.

The physical reason for this behaviour is that the liquid contains particles (such as clay) or large molecules (such as polymers) which have some kind of interaction, creating a weak solid structure, formerly known as a false body, and a certain amount of stress is required to break this structure. Once the structure has been broken, the particles move with the liquid under viscous forces. If the stress is removed, the particles associate again.

The material is an elastic solid for shear stress ${\displaystyle \tau }$, less than a critical value ${\displaystyle \tau _{0}}$. Once the critical shear stress (or "yield stress") is exceeded, the material flows in such a way that the shear rate, ∂u/∂y (as defined in the article on viscosity), is directly proportional to the amount by which the applied shear stress exceeds the yield stress:

${\displaystyle {\frac {\partial u}{\partial y}}={\begin{cases}0,&\tau <\tau _{0}\\{\frac {\tau -\tau _{0}}{\mu _{\infty }}},&\tau \geq \tau _{0}\end{cases}}}$

In fluid flow, it is a common problem to calculate the pressure drop in an established piping network.^[4] Once the friction factor, f, is known, it becomes easier to handle different pipe-flow problems, viz. calculating the pressure drop for evaluating pumping costs or to find the flow-rate in a piping network for a given pressure drop. It is usually extremely difficult to arrive at exact analytical solution to calculate the friction factor associated with flow of non-Newtonian fluids and therefore explicit approximations are used to calculate it. Once the friction factor has been calculated the pressure drop can be easily determined for a given flow by the Darcy–Weisbach equation:

${\displaystyle f={2h_{\text{f}}gD \over LV^{2}}}$

where:

• ${\displaystyle f}$ is the Darcy friction factor (SI units: dimensionless)
• ${\displaystyle h_{\text{f}}}$ is the frictional head loss (SI units: m)
• ${\displaystyle g}$ is the gravitational acceleration (SI units: m/s²)
• ${\displaystyle D}$ is the pipe diameter (SI units: m)
• ${\displaystyle L}$ is the pipe length (SI units: m)
• ${\displaystyle V}$ is the mean fluid velocity (SI units: m/s)

An exact description of friction loss for Bingham plastics in fully developed laminar pipe flow was first published by Buckingham.^[5] His expression, the Buckingham–Reiner equation, can be written in a dimensionless form as follows:

${\displaystyle f_{\text{L}}={64 \over \operatorname {Re} }\left[1+{\operatorname {He} \over 6\operatorname {Re} }-{64 \over 3}\left({\operatorname {He} ^{4} \over f^{3}\operatorname {Re} ^{7}}\right)\right]}$

where:

• ${\displaystyle f_{\text{L}}}$ is the laminar flow Darcy friction factor (SI units: dimensionless)
• ${\displaystyle \operatorname {Re} }$ is the Reynolds number (SI units: dimensionless)
• ${\displaystyle \operatorname {He} }$ is the Hedstrom number (SI units: dimensionless)

The Reynolds number and the Hedstrom number are respectively defined as:

${\displaystyle \operatorname {Re} ={\rho VD \over \mu },}$ and

${\displaystyle \operatorname {He} ={\rho D^{2}\tau _{o} \over \mu ^{2}}}$

where:

• ${\displaystyle \rho }$ is the mass density of fluid (SI units: kg/m³)
• ${\displaystyle \mu }$ is the dynamic viscosity of fluid (SI units: kg/m s)
• ${\displaystyle \tau _{o}}$ is the yield point (yield strength) of fluid (SI units: Pa)

Darby and Melson developed an empirical expression^[6] that was then refined, and is given by:^[7]

${\displaystyle f_{\text{T}}=4\times 10^{a}\operatorname {Re} ^{-0.193}}$

where:

• ${\displaystyle f_{\text{T}}}$ is the turbulent flow friction factor (SI units: dimensionless)
• ${\displaystyle a=-1.47\left[1+0.146e^{-2.9\times {10^{-5}}\operatorname {He} }\right]}$

Note: Darby and Melson's expression is for a Fanning friction factor, and needs to be multiplied by 4 to be used in the friction loss equations located elsewhere on this page.

Although an exact analytical solution of the Buckingham–Reiner equation can be obtained because it is a fourth order polynomial equation in f, due to complexity of the solution it is rarely employed. Therefore, researchers have tried to develop explicit approximations for the Buckingham–Reiner equation.

The Swamee–Aggarwal equation is used to solve directly for the Darcy–Weisbach friction factor f for laminar flow of Bingham plastic fluids.^[8] It is an approximation of the implicit Buckingham–Reiner equation, but the discrepancy from experimental data is well within the accuracy of the data. The Swamee–Aggarwal equation is given by:

${\displaystyle f_{L}={64 \over \mathrm {Re} }+{64 \over \mathrm {Re} }\left({\mathrm {He} \over 6.2218\mathrm {Re} }\right)^{0.958}}$

Danish et al. have provided an explicit procedure to calculate the friction factor f by using the Adomian decomposition method.^[9] The friction factor containing two terms through this method is given as:

${\displaystyle f_{L}={\frac {K_{1}+{\dfrac {4K_{2}}{\left(K_{1}+{\frac {K_{1}K_{2}}{K_{1}^{4}+3K_{2}}}\right)^{3}}}}{1+{\dfrac {3K_{2}}{\left(K_{1}+{\frac {K_{1}K_{2}}{K_{1}^{4}+3K_{2}}}\right)^{4}}}}}}$

where

${\displaystyle K_{1}={16 \over \mathrm {Re} }+{16\mathrm {He} \over 6\mathrm {Re} ^{2}},}$

and

${\displaystyle K_{2}=-{16\mathrm {He} ^{4} \over 3\mathrm {Re} ^{8}}.}$

In 1981, Darby and Melson, using the approach of Churchill^[10] and of Churchill and Usagi,^[11] developed an expression to get a single friction factor equation valid for all flow regimes:^[6]

${\displaystyle f=\left[{f_{\text{L}}}^{m}+{f_{\text{T}}}^{m}\right]^{\frac {1}{m}}}$

where:

${\displaystyle m=1.7+{40000 \over \operatorname {Re} }}$

Both Swamee–Aggarwal equation and the Darby–Melson equation can be combined to give an explicit equation for determining the friction factor of Bingham plastic fluids in any regime. Relative roughness is not a parameter in any of the equations because the friction factor of Bingham plastic fluids is not sensitive to pipe roughness.

• Bagnold number
• Bernoulli's principle
• Bingham-Papanastasiou model
• Rheology
• Shear thinning