PDF version
TL;DR
The Herschel–Bulkley model unifies Newtonian, Bingham, and power-law fluids via a yield stress and a strain-rate-dependent viscosity. An 
-regularization ensures stable computations and recovers simpler models (Newtonian, Bingham) by tuning model parameters. Dimensionless groups (e.g., the plasto-capillary number  and the effective Ohnesorge) capture the interplay of fluid rheology, capillarity, and flow scales. Implementation details are provided, along with references, open-source code, and demonstrations of bubble-burst simulations in viscoplastic media.
Features:
- Yield stress
- Power law dependance on the strain rate
- Shear thinning for
. - Shear thickening for
.
- Shear thinning for
- Bingham model for
. - Newtonian fluid for
and .
-formulation
Normalizing stresses with
Here, the effective Ohnesorge is
The plasto-capillary number
One can easily see that putting Newtonian response.
More details on the implementation
Calculate the norm of the deformation tensor :
Note:
We use the formulation as given in Balmforth et al. (2013) [1], who use the strain rate tensor
Calculate the equivalent viscosity
Factorizing with
In this formulation,
Note: The fluid flows always, it is not a solid, but a very viscous fluid.
Reproduced from: P.-Y. Lagrée’s Sandbox. Here, we use a face implementation of the regularisation method, described here.
Further exploration:
Video showcasing a typical simulation of bubble bursting in a Herschel–Bulkley fluid medium
More resources
References
[1] N. J. Balmforth, I. A. Frigaard, and G. Ovarlez, “Yielding to Stress: Recent Developments in Viscoplastic Fluid Mechanics,” Annu. Rev. Fluid Mech., vol. 46, pp. 121–146, Jan. 2014, doi: 10.1146/annurev-fluid-010313-141424.
[2] V. Sanjay, D. Lohse, and M. Jalaal, “Bursting bubble in a viscoplastic medium,” J. Fluid Mech., vol. 922, p. A2, 2021.
Metadata
Author:: Vatsal Sanjay Date published:: Dec 31, 2024
Date modified:: Jan 26, 2025 at 11:50 CET
Back to main website
