Bingham-Papanastasiou model explained

An important class of non-Newtonian fluids presents a yield stress limit which must be exceeded before significant deformation can occur – the so-called viscoplastic fluids or Bingham plastics. In order to model the stress-strain relation in these fluids, some fitting have been proposed such as the linear Bingham equation and the non-linear Herschel-Bulkley and Casson models.^[1]

Analytical solutions exist for such models in simple flows. For general flow fields, it is necessary to develop numerical techniques to track down yielded/unyielded regions. This can be avoided by introducing into the models a continuation parameter, which facilitates the solution process and produces virtually the same results as the ideal models by the right choice of its value.^[2]

Viscoplastic materials like slurries, pastes, and suspension materials have a yield stress, i.e. a critical value of stress below which they do not flow are also called Bingham plastics, after Bingham.^[3]

Viscoplastic materials can be well approximated uniformly at all levels of stress as liquids that exhibit infinitely high viscosity in the limit of low shear rates followed by a continuous transition to a viscous liquid. This approximation could be made more and more accurate at even vanishingly small shear rates by means of a material parameter that controls the exponential growth of stress. Thus, a new impetus was given in 1987 with the publication by Papanastasiou^[4] of such a modification of the Bingham model with an exponential stress-growth term. The new model basically rendered the original discontinuous Bingham viscoplastic model as a purely viscous one, which was easy to implement and solve and was valid for all rates of deformation. The early efforts by Papanastasiou and his co-workers were taken up by the author and his coworkers,^[5] who in a series of papers solved many benchmark problems and presented useful solutions always providing the yielded/unyielded regions in flow fields of interest. Since the early 1990s, other workers in the field also used the Papanastasiou model for many different problems.

Papanastasiou

Papanastasiou in 1987, who took into account earlier works in the early 1960s^[6] as well as a well-accepted practice in the modelling of soft solids^[7] and the sigmoidal modelling behaviour of density changes across interfaces.^[8] He introduced a continuous regularization for the viscosity function which has been largely used in numerical simulations of viscoplastic fluid flows, thanks to its easy computational implementation. As a weakness, its dependence on a non-rheological (numerical) parameter, which controls the exponential growth of the yield-stress term of the classical Bingham model in regions subjected to very small strain-rates, may be pointed. Thus, he proposed an exponential regularization of eq., by introducing a parameter m, which controls the exponential growth of stress, and which has dimensions of time. The proposed model (usually called Bingham-Papanastasiou model) has the form:

$$\backslash vec\; =\; \backslash left(\backslash mu\; +\; [1-\backslash exp(-m\backslash mid\backslash dot\{\backslash gamma\}\backslash mid)]\; \backslash right)\; (\backslash vec)$$

and is valid for all regions, both yielded and unyielded. Thus it avoids solving explicitly for the location of the yield surface, as was done by Beris et al.^[9]

Papanastasiou's modification, when applied to the Bingham model, becomes in simple shear flow (1-D flow):

Bingham-Papanastasiou model:

• $$\backslash tau\; =\; \backslash tau\_y\; [1-\backslash exp(-m\backslash dot\{\backslash gamma\})]\; +\; \backslash mu\; \backslash dot$$
  • $$\backslash eta\; =\; \backslash mu\; +\; [1-\backslash exp(-m\; \backslash mid\backslash dot\{\backslash gamma\}\backslash mid)]$$

where η is the apparent viscosity.

Notes and References

1. Soto . Hilda Pari . Martins-Costa . Maria Laura . Fonseca . Cleiton . Frey . Sérgio . A numerical investigation of inertia flows of Bingham-Papanastasiou fluids by an extra stress-pressure-velocity galerkin least-squares method . Journal of the Brazilian Society of Mechanical Sciences and Engineering . December 2010 . 32 . 5 . 450–460 . 10.1590/S1678-58782010000500004 . 10183/75845 . free . free .
2. Mitsoulis . Evan . Tsamopoulos . John . Numerical simulations of complex yield-stress fluid flows . Rheologica Acta . March 2017 . 56 . 3 . 231–258 . . 10.1007/s00397-016-0981-0 . 99126659 .
3. Book: Bingham . Eugene Cook . Fluidity and plasticity . 1922 . McGraw-Hill . 1118524319 .
4. Papanastasiou . Tasos C. . Flows of Materials with Yield . Journal of Rheology . July 1987 . 31 . 5 . 385–404 . 10.1122/1.549926 . 1987JRheo..31..385P .
5. Ellwood . K. R. J. . Georgiou . G. C. . Papanastasiou . T. C. . Wilkes . J. O. . Laminar jets of Bingham‐plastic liquids . Journal of Rheology . August 1990 . 34 . 6 . 787–812 . 10.1122/1.550144 . 1990JRheo..34..787E . 59521944 . free .
6. Shangraw . Ralph . Grim . Wayne . Mattocks . Albert M. . An Equation for Non‐Newtonian Flow . Transactions of the Society of Rheology . March 1961 . 5 . 1 . 247–260 . 10.1122/1.548898 . 1961JRheo...5..247S .
7. Gavrus . A. . E. . Ragneau . P. . Caestecker . Rheological behaviour formulation of solid metallic materials for dynamic forming processes simulation . Proceedings of the 4th International ESAFORM Conference on Material Forming . 2001 . 403–406 . 51843097 .
8. Hirt . C.W . Nichols . B.D . Volume of fluid (VOF) method for the dynamics of free boundaries . Journal of Computational Physics . January 1981 . 39 . 1 . 201–225 . 10.1016/0021-9991(81)90145-5 . 1981JCoPh..39..201H .
9. Beris . A. N. . Tsamopoulos . J. A. . Armstrong . R. C. . Brown . R. A. . Creeping motion of a sphere through a Bingham plastic . Journal of Fluid Mechanics . September 1985 . 158 . 219–244 . 10.1017/S0022112085002622 . 1985JFM...158..219B . 122781643 .