Herschel-Bulkley model

Many materials in our daily lives defy simple classification as either a solid or a liquid. Consider ketchup or toothpaste: they hold their shape under gravity but flow easily when squeezed. These materials are known as yield-stress fluids, and their behavior is captured by a powerful mathematical relationship called the Herschel-Bulkley model. This model addresses the gap in our understanding of substances that possess a solid-like stubbornness at rest but transition to a fluid-like state under sufficient force. This article will guide you through this fascinating topic. First, we will explore the "Principles and Mechanisms," dissecting the model's equation and its three key parameters, and examining the microscopic origins of yield stress. Following that, in "Applications and Interdisciplinary Connections," we will journey through its real-world relevance, from everyday engineering and large-scale geological events to the intricate mechanics of biological systems.

Imagine trying to get the last bit of ketchup out of a bottle. You turn it upside down. Nothing happens. You wait. Still nothing. The ketchup, under its own weight, refuses to flow. It possesses a stubbornness, a resistance to motion that a simple fluid like water doesn't have. You have to give it a good, sharp shake or a squeeze—you have to apply a force that overcomes its initial reluctance. Only then does it decide to move. This everyday struggle is a perfect introduction to the world of yield-stress fluids, and the key to understanding them is a beautifully versatile equation known as the ​​Herschel-Bulkley model​​.

At its heart, the Herschel-Bulkley model is a simple relationship between the force you apply to make a fluid flow (the ​​shear stress​​, $\tau$) and how fast it flows as a result (the ​​shear rate​​, $\dot{\gamma}$). For stresses greater than some threshold, it reads:

$\tau = \tau_y + K(\dot{\gamma})^n$

This equation may look a little dense, but it's really just a story with three characters, each controlling a different aspect of the fluid's behavior.

The first and most important character is $\tau_y$, the ​​yield stress​​. This is the ketchup's stubbornness, quantified. It's the minimum stress you must apply to get the fluid to "yield" and start moving. Below this value, the material acts like a solid—it might deform a little, but it won't flow. Toothpaste stays on your brush, and paint doesn't drip off the wall, precisely because of their yield stress.

Once you've paid the "entry fee" of $\tau_y$, the rest of the equation, $K(\dot{\gamma})^n$, kicks in to describe how the fluid behaves while it's flowing. Here we meet our other two characters.

The parameter $K$ is the ​​consistency index​​. You can think of it as the material's "heft" or "thickness" once it's in motion. A high $K$ means the fluid is very resistive, like trying to stir thick mud, while a low $K$ means it flows more easily, like a slurry. What's fascinating about $K$ is that its physical meaning, and even its units, are tied to the third parameter, $n$. It's not a simple viscosity like you'd find for water or honey; its role is more nuanced, which hints at the complex physics hiding beneath the surface.

The final character, $n$, is the ​​flow behavior index​​, and it's the most interesting of all. It describes the personality of the flow.

• If $\boldsymbol{n=1}$, the equation simplifies to $\tau = \tau_y + K\dot{\gamma}$. This describes a ​​Bingham plastic​​. Once it yields, the stress increases linearly with the shear rate, just like a normal Newtonian fluid, but with that initial offset $\tau_y$. Thick mud or clay slurries often behave this way.

• If $\boldsymbol{n 1}$, we have a ​​shear-thinning​​ fluid. This is the magic behind paint. As you shear it faster (by spreading it with a brush), it becomes less resistant—its apparent viscosity drops. This is why paint can be thick enough not to drip from the brush (low shear rate) but spread smoothly and easily on the wall (high shear rate).

• If $\boldsymbol{n>1}$, we get a ​​shear-thickening​​ fluid. This is the strange behavior of "oobleck" (a mix of cornstarch and water). If you stir it slowly, it's a liquid. But if you punch it or try to stir it quickly, it abruptly becomes almost solid.

The Herschel-Bulkley model is so powerful because, by just tuning these three parameters, it can describe this whole family of behaviors, from the simple Bingham plastic ($n=1$) to the quintessential power-law fluid (by setting $\tau_y=0$). It's a unified description of a vast range of materials we encounter every day.

This macroscopic model is wonderfully descriptive, but as physicists, we can't help but ask: why? Where does this yield stress come from? To find the answer, we have to zoom in and look at the microscopic structure of the material. It turns out that yield stress isn't one thing; it can arise from at least two different underlying physical mechanisms.

Imagine a dense suspension of particles, like a thick soup of microscopic beads.

​​Mechanism 1: The Sticky Network.​​ If the particles have some attraction to each other—a bit like tiny, weak magnets—they can form a connected, sprawling network that spans the entire container. This structure of interlocked particles is what gives the material its initial solidity. It can resist a small force by just stretching its bonds. The yield stress, $\tau_y$, is the force per unit area required to globally break this network. Once the network is broken, the particles can flow past each other, but they are still interacting, breaking and reforming bonds. The faster you shear the fluid, the more broken-up the network becomes, offering less resistance. This "structural weakening" is the physical origin of shear-thinning behavior ($n 1$). This mechanism—a cohesive, breaking network—is beautifully captured by the full Herschel-Bulkley model.

​​Mechanism 2: The Jammed Crowd.​​ Now, imagine the particles are not sticky, but are simply crowded together and held in place by an external pressure, like sand in a sandbag. They are jammed. To make them flow, you have to push hard enough to overcome the static friction between the particles so they can slide past one another. The yield stress, in this case, comes from this collective friction, and it's proportional to the confining pressure (just as it's harder to drag a heavy box across the floor than a light one). Once they start sliding, the resistance to flow comes mainly from the particles having to push the surrounding liquid out of the way (viscous drag). This leads to a behavior where the stress needed is the yield stress (from friction) plus a term proportional to the shear rate (from viscosity). This is the microscopic picture of a Bingham plastic ($n=1$).

These two pictures tell us that what we call a "yield-stress fluid" isn't a single category but a collection of behaviors with different origins. The Herschel-Bulkley model is the mathematical language that allows us to speak about them all.

This raises a subtle but profound question: Is a yield-stress fluid just an extremely thick liquid? Is toothpaste just very, very viscous honey? The answer is a definitive "no," and the mathematics of the Herschel-Bulkley model shows us why.

Let's define a fluid's ​​apparent viscosity​​ as $\eta_{app} = \tau/\dot{\gamma}$. For a simple fluid like water, this is a constant. For a shear-thinning fluid without a yield stress (like a polymer solution), this viscosity might be very high at low shear rates, but it's always a finite number. This means that if you apply any stress, no matter how small, it will flow, albeit very slowly.

But for a Herschel-Bulkley fluid, something dramatic happens. As the shear rate $\dot{\gamma}$ approaches zero, the apparent viscosity becomes:

$\eta_{app} = \frac{\tau_y + K\dot{\gamma}^n}{\dot{\gamma}} = \frac{\tau_y}{\dot{\gamma}} + K\dot{\gamma}^{n-1}$

Because $\tau_y$ is a positive constant, as $\dot{\gamma}$ gets infinitesimally small, the term $\tau_y/\dot{\gamma}$ explodes to ​​infinity​​. A true yield-stress fluid, in this idealized model, has infinite resistance to the very first inkling of flow. It's not just a very thick liquid; it's a different kind of beast altogether. It possesses a true solidity that must be broken. This mathematical divergence is the sharp, clear signature of a yield stress.

Armed with this model, we can start to explore how these fluids behave in real-world situations, and we quickly run into fascinating "paradoxes" that challenge our intuition and deepen our understanding.

The Stubborn Plug

Imagine our fluid being pumped down a long, straight pipe. The pressure pushes it forward. Our intuition for a shear-thinning fluid might suggest that its strange properties would affect the entire flow field. But here, physics neatly separates into two parts: force balance and material response.

First, a simple force balance (Newton's laws, really) tells us that the shear stress inside the pipe must be highest at the walls and zero at the very center. This is true for any fluid, whether it's water, paint, or pancake batter. The stress distribution is dictated by mechanics, not by the material.

Now, we bring in the material's personality. In the regions near the wall, the stress is high—higher than $\tau_y$—so the material yields and flows. But as we move toward the center of the pipe, the stress drops. At some point, the stress becomes equal to $\tau_y$. Inside this radius, the stress is below the yield stress. And what does a yield-stress fluid do when the stress is below $\tau_y$? It doesn't deform. It behaves like a solid.

The result is a "plug" of solid-like material moving down the center of the pipe, with layers of yielded fluid sliding along beside it. The size of this plug is determined only by where the stress equals the yield stress. It depends on the pressure gradient and $\tau_y$, but—and this is the surprise—it is completely independent of the flow parameters $K$ and $n$. The shear-thinning nature of the fluid affects how fast the yielded parts are flowing, but it doesn't change the boundary of the solid plug. This is a beautiful example of how separating principles—force balance and constitutive law—can lead to crystal-clear, if unexpected, conclusions.

The Sliding Solid

Let's consider another puzzle. What if our yield-stress fluid is in a channel where the walls are "slippery"? Instead of a no-slip condition, let's say the stress at the wall is proportional to the fluid's velocity at the wall. This is called a ​​Navier slip condition​​.

Now, let's apply only a very small pressure gradient, so small that the stress everywhere, even at the wall, is below the yield stress $\tau_y$. According to our model, the fluid cannot yield. It must behave as a solid. Our intuition screams that if it can't yield, it can't move.

But our intuition is wrong. The constitutive law, $\dot{\gamma} = 0$, doesn't say the velocity is zero; it says the velocity gradient is zero. This means the material can move, but it must do so as a single, rigid block, with a uniform velocity across the channel. This rigid translation is perfectly consistent with the "solid" state. But how fast does it move? The speed is determined by the boundary condition. The pressure force on the plug is balanced by the friction from the slippery walls. A non-zero pressure gradient results in a non-zero wall stress, which, through the slip law, requires a non-zero velocity.

So we have the remarkable situation of a solid block of unyielded material sliding down a channel, its motion enabled entirely by slip at the boundaries. It beautifully illustrates that the constitutive model governs internal deformation, while the boundary conditions govern how the object as a whole interacts with its surroundings.

The Herschel-Bulkley model is more than a theoretical curiosity; it's a workhorse in engineering and science. But using it effectively requires judgment and an appreciation for the art of modeling.

​​The Modeler's Dilemma.​​ Suppose you have experimental data that clearly shows a fluid is shear-thinning ($n \approx 0.7$). You need to run a computer simulation. Should you always use the more accurate Herschel-Bulkley model? Not necessarily. If your application only involves a very narrow range of shear rates, a simpler Bingham model ($n=1$) might be a perfectly adequate approximation within that range. A simpler model can be computationally cheaper and more robust, making it the smarter choice for that specific job. However, if your simulation involves a wide range of shear rates, using the wrong physics (like assuming $n=1$) can lead to dramatically wrong predictions. The lesson is that there is no single "best" model; the right choice is a trade-off between fidelity, cost, and the question you are trying to answer.

​​The Experimentalist's Trap.​​ How do we find the values of $\tau_y$, $K$, and $n$ from real data? A common trick for power-laws is to take the logarithm of the data to turn it into a straight line. One might try to guess a value for $\tau_y$, subtract it from the stress data, and then plot $\ln(\tau - \tau_y)$ versus $\ln(\dot{\gamma})$. If you've guessed $\tau_y$ correctly, you'll get a straight line with a slope of $n$. But what if your guess is slightly off? The math shows that your "straight line" will actually be a subtle curve. Fitting a line to this curve will give you a biased, incorrect value for $n$. The scientifically honest approach is to fit the full, nonlinear Herschel-Bulkley equation directly to the raw data. It's computationally harder, but it respects the integrity of the model and the data, a principle of utmost importance in any scientific endeavor.

From the stubbornness of ketchup to the design of advanced materials, the Herschel-Bulkley model provides a language to describe, understand, and predict the behavior of a fascinating class of materials. It reminds us that sometimes, the most profound insights are hidden in the most familiar phenomena, just waiting for a framework that can give them a voice.

Having grasped the fundamental principles of the Herschel-Bulkley model, we are now ready to embark on a journey. We will see that these ideas are not confined to the abstract world of equations but are, in fact, all around us, operating in our kitchens, inside our bodies, in the awesome and terrifying theater of nature, and at the very frontier of modern technology. The model's real beauty lies in its power to unify a vast range of seemingly disconnected phenomena, revealing a common physical principle at work: the fascinating duality of being solid-like at rest and liquid-like in motion.

Our first stop is the world of materials we use every day. Have you ever wondered why a good quality paint doesn't drip from the brush, yet spreads smoothly and effortlessly onto a wall? The secret is yield stress. At rest, under the low stress of gravity, the paint's yield stress $\tau_y$ is not exceeded, and it behaves like a soft solid, clinging to the bristles. But the moment you apply the high shear of a brushstroke, you overcome this yield stress. The paint begins to flow, and because it is also shear-thinning (with a flow index $n \lt 1$), its apparent viscosity plummets, allowing it to be spread into a thin, even layer.

This same principle is crucial in high-precision manufacturing. Consider the challenge of screen printing electronic circuits using a conductive ink. The ink must be liquid enough to be forced through a fine mesh by a squeegee—a high-shear process—but once it is on the circuit board, it must instantly "solidify" and hold its shape with perfect fidelity, without bleeding or smudging under the low stresses of gravity and surface tension. A material with a high yield stress and strong shear-thinning behavior exhibits a dramatically lower apparent viscosity during printing than when it is at rest, making it the perfect candidate for the job.

The journey of these materials doesn't end there. How do we transport thick substances like toothpaste, concrete, or food slurries through miles of pipeline? Pumping a simple, highly viscous fluid like honey requires enormous energy. Here, the Herschel-Bulkley model reveals a wonderfully counter-intuitive phenomenon: plug flow. When a yield-stress fluid is pumped through a pipe, the shear stress is highest at the pipe wall and zero at the center. This creates a fascinating situation where a layer of fluid near the wall yields and flows, while a large central core, where the stress is below $\tau_y$, moves as a single, solid "plug." This plug effectively slides along, lubricated by the flowing layer at the edge. The size of this solid core is directly related to the fluid's yield stress and inversely related to the pressure gradient driving the flow. This is not just a curiosity; it's a fundamental principle for efficiently engineering the transport of countless materials that define our modern world.

The same physics that keeps paint on your brush also governs some of the most powerful forces on our planet. Let's scale up our thinking from pipelines to entire landscapes. A lahar, or volcanic mudflow, is a terrifying slurry of water, ash, and rock. Why is it so destructive? A simple flood of water might move cars, but a lahar can transport enormous boulders, houses, and bridges. The reason is its immense yield stress.

A boulder sitting in a lahar is pulled down by gravity, but it is also pushed up by a supportive force from the fluid's internal strength. This supportive force is a direct consequence of the yield stress. For the boulder to be suspended and carried along, this supportive force must be greater than the boulder's submerged weight. This means there is a maximum size of rock that a lahar of a given yield stress $\tau_y$ can carry. Anything larger will be deposited. Geologists use this very principle, sometimes in reverse by examining the size of boulders left behind by ancient flows, to estimate the yield stress and destructive power of past and future volcanic hazards.

Descending from the surface deep into the Earth's crust, we find another dramatic application in the flow of magma. Magma is not the simple, uniform red liquid of cartoons; it is a complex suspension of solid crystals within a liquid silicate melt. This mixture behaves as a Herschel-Bulkley fluid. The rheological parameters—yield stress $\tau_y$, consistency $K$, and flow index $n$—are not constant but depend critically on the fraction of crystals present. As magma cools and crystallizes, its yield stress increases dramatically. This behavior governs everything from the speed of lava flows to the way magma intrudes into rock to form dikes. By modeling this crystal-dependent rheology, geophysicists can predict the width of the solid "plug" of magma in a dike and better understand the mechanics of volcanic eruptions and the very plumbing of our planet.

The Herschel-Bulkley model not only describes the world around us; it describes the world within us. Our bodies are master rheologists, constantly manipulating complex fluids. The simple act of swallowing is a profound biomechanical challenge. To move a bolus of pureed food from the mouth to the esophagus, the tongue must act like a piston, generating enough pressure to overcome the food's yield stress. If the pressure generated by the tongue is insufficient to make the food yield and flow, swallowing fails. This principle is not just academic; it is of vital clinical importance in designing foods for individuals with dysphagia (swallowing difficulties), ensuring the food is safe and easy to swallow while still being palatable.

Perhaps the most famous non-Newtonian fluid of all is our own blood. At a first glance, the Herschel-Bulkley model seems to fit blood quite well. At very low flow rates, red blood cells tend to clump together into stacks called "rouleaux," forming a network that gives blood an apparent yield stress. A certain amount of force is needed to break this network and initiate flow. Once flowing, these rouleaux break apart and the individual red blood cells deform and align with the flow, causing the apparent viscosity to drop—a classic shear-thinning behavior.

However, here we must learn a classic Feynman-esque lesson about the nature of physical models. While the Herschel-Bulkley model provides a valuable first approximation, it is not the whole story. Unlike our simple model, blood's properties are time-dependent (thixotropic)—it takes time for the rouleaux network to form or break. Furthermore, the model predicts that viscosity should continue to decrease indefinitely at high shear rates, whereas blood's viscosity approaches a constant plateau. This doesn't mean the model is "wrong"; it means it is an idealization that captures some, but not all, of the fluid's rich behavior. It reminds us that science progresses by creating models, testing their limits, and then building more refined ones.

The model finds a more direct home in other biological contexts. The "slime" produced by bacteria, known as a biofilm, is a city of microbes encased in a protective matrix of extracellular polymeric substances (EPS). This EPS matrix is a yield-stress fluid. This gives the biofilm a solid-like integrity, making it incredibly resistant to being washed away, whether from a ship's hull, a river stone, or your own teeth (as dental plaque). Conversely, the mucus lining our respiratory tracts and lubricating our joints uses its Herschel-Bulkley properties for a different purpose. Its yield stress allows it to form a stable, protective layer. Yet, when subjected to shear—by the sliding of a joint or the cilia in our lungs—it flows easily, acting as an exceptional lubricant. The yield stress, in this case, provides a form of static friction, or "stiction," that must be overcome to initiate sliding, a key principle in the field of biotribology.

We end our journey at the frontier, where our understanding of these materials allows us not just to describe the world, but to build it in new ways. The field of 3D bioprinting aims to construct living tissues, scaffolds, and even organs, layer by layer. The "ink" in this process is often a cell-laden hydrogel.

The success of bioprinting hinges on a delicate rheological balancing act. The bio-ink must have a low enough apparent viscosity under the high shear of the printing nozzle to flow smoothly without requiring pressures so high that they damage or kill the living cells within. Yet, the moment it is extruded onto the substrate, it must have a high enough yield stress to solidify and hold its shape, allowing the construction of complex, overhanging 3D structures. The goal is to find the perfect set of Herschel-Bulkley parameters $(\tau_y, K, n)$ that minimizes stress on the cells during printing while maximizing the shape fidelity of the final construct. Here, the model becomes more than a descriptive tool; it is a design equation for the fabrication of life itself.

From the mundane to the magnificent, from the geological to the biological, the Herschel-Bulkley model provides a unifying language. It shows us how a single, elegant physical concept can explain the behavior of a staggering diversity of materials, revealing the deep and beautiful connections that underlie the complexity of our world.