Constitutive Relation for Newtonian Fluids

The intuitive feeling of resistance when stirring a liquid—its "thickness" or "stickiness"—is a phenomenon we all understand. This property, known as viscosity, is the defining feature of real-world fluids. However, moving from this intuitive concept to a precise, predictive scientific law requires a more rigorous language. How do we mathematically capture the internal friction that governs the flow of air over a wing, blood through an artery, or magma deep within the Earth? The answer lies in a foundational model of fluid dynamics first proposed by Isaac Newton.

This article addresses the gap between the qualitative experience of fluid friction and its quantitative description. It unpacks the constitutive relation for Newtonian fluids, the elegant rule that forms the bedrock of our understanding of a vast category of flows. By exploring this principle, you will gain insight into the fundamental physics that connects the forces within a fluid to its motion.

We will begin by examining the "Principles and Mechanisms," developing the language of stress and strain tensors to formally define the linear relationship that characterizes a Newtonian fluid. We will also explore its surprising consequences, from the generation of heat to the origins of the no-slip condition. Following this, the article will journey through "Applications and Interdisciplinary Connections," showcasing how this single, simple law provides critical insights into phenomena across engineering, geophysics, and biology, revealing the profound unity of physics in a complex world.

Imagine stirring a cup of tea versus a jar of honey. You feel a resistance in both, but it's vastly different. This resistance to flow, this internal friction, is what we call ​​viscosity​​. It’s the defining characteristic that separates the idealized, perfectly slippery fluids of introductory physics from the real fluids that shape our world, from the air flowing over a wing to the blood coursing through our veins. But how do we capture this sticky, messy, and beautiful property in the precise language of physics? How does the simple act of stirring lead to the complex dance of eddies and the inevitable warming of the fluid?

The journey to understand this requires us to build a new language, one that can describe how forces are transmitted and how fluid elements deform.

When you push on a solid block, the force is straightforward. But a fluid is a continuous medium. A force applied at one point creates a complex pattern of internal forces throughout. To describe this, we can’t just talk about a single force vector. We need a more powerful concept: the ​​stress tensor​​.

Imagine a tiny imaginary cube placed anywhere inside our moving fluid. Forces are acting on each of its six faces. These forces aren't necessarily perpendicular to the faces; they can be both pushing (or pulling) and scraping. To characterize the state of force at that point, we need to specify the force on each face.

Let's use a coordinate system $(x, y, z)$. Consider the face of our cube whose outward normal points in the $x$-direction. The fluid on the outside of this face exerts a force on the fluid inside. This force has three components: one in the $x$-direction (a normal push/pull), one in the $y$-direction (a tangential scrape), and one in the $z$-direction (another tangential scrape). We denote these forces-per-unit-area as $\sigma_{xx}$, $\sigma_{xy}$, and $\sigma_{xz}$. The first index tells you the orientation of the surface (a surface "facing" $x$), and the second index tells you the direction of the force component.

So, $\sigma_{xx}$ is a ​​normal stress​​, the kind of stress associated with pressure. But $\sigma_{xy}$ is a ​​shear stress​​: it's the force per unit area in the $y$-direction acting on a surface whose normal is in the $x$-direction. This is the mathematical embodiment of the scraping or rubbing friction between adjacent layers of fluid moving at different speeds. The collection of all nine such components forms the ​​Cauchy stress tensor​​, $\boldsymbol{\sigma}$, a machine that tells you the force vector on any surface you can imagine slicing through the fluid.

Now, what causes these stresses (apart from the ever-present static pressure)? The answer is motion—or more precisely, the relative motion of different parts of the fluid. A fluid that translates or rotates as a rigid block experiences no internal friction. Friction only arises when the fluid is being deformed. To describe this deformation, we look at the ​​velocity gradient​​, $\nabla \mathbf{u}$, a tensor that tells us how the velocity vector $\mathbf{u}$ changes from point to point.

Amazingly, any complex infinitesimal motion described by $\nabla \mathbf{u}$ can be cleanly split into two parts: a pure rotation and a pure strain (deformation). The rotational part is captured by the ​​spin tensor​​ $\boldsymbol{\Omega}$, which is related to the ​​vorticity​​ of the flow (the local spinning motion). The deformational part is captured by the ​​rate-of-strain tensor​​, $\boldsymbol{S}$. A crucial insight of fluid mechanics is that for most fluids, the viscous stress does not depend on the local rotation, only on the rate of strain. It doesn’t matter how fast a fluid element is spinning, only how much it is being stretched or sheared.

So, we have the stress $\boldsymbol{\tau}$ (the viscous part of $\boldsymbol{\sigma}$) representing the internal friction, and we have the rate of strain $\boldsymbol{S}$ representing the deformation. What is the relationship between them? In the late 17th century, Isaac Newton proposed the simplest possible relationship: they are directly proportional.

This is the essence of a ​​Newtonian fluid​​: the stress that arises from viscosity is a linear function of the rate of strain.

For an ​​incompressible fluid​​—one whose density doesn't change, like water or oil under typical conditions—this "Newtonian bargain" takes on a beautifully simple form:

Here, $\boldsymbol{\tau}$ is the viscous stress tensor, and $\boldsymbol{S}$ is the rate-of-strain tensor. The constant of proportionality, $\mu$, is the ​​dynamic viscosity​​, the number that tells us whether we have water or honey. The factor of 2 is a convention that makes things work out nicely in simple cases. For instance, in a simple shear flow between two plates, this grand tensor equation elegantly reduces to the familiar form relating shear stress $\tau$ to the shear rate $\dot{\gamma}$: $\tau = \mu \dot{\gamma}$.

The full stress tensor for an incompressible Newtonian fluid is then the sum of the isotropic pressure and the viscous stress:

where $\boldsymbol{I}$ is the identity tensor. A subtle but profound point lies hidden here. For an incompressible fluid, density is constant, so pressure can no longer be determined from an equation of state like the ideal gas law. What, then, determines the pressure $p$? The pressure becomes a kind of enforcer. It adjusts itself instantaneously throughout the fluid to whatever value is needed to ensure the fluid remains incompressible ($\nabla \cdot \mathbf{u} = 0$). In the language of advanced mechanics, the pressure field $p$ acts as a ​​Lagrange multiplier​​ for the incompressibility constraint. It is a mechanical pressure, not a thermodynamic one.

This simple linear law is not a universal truth, but a model. However, it is an astonishingly successful one, accurately describing a vast range of fluids, including air, water, oil, and many gases. For ​​compressible fluids​​, where density can change, the law is slightly more complex, involving a second coefficient of viscosity, $\lambda$, related to resistance to volume change. The full relation becomes $\boldsymbol{\sigma} = -p\boldsymbol{I} + 2\mu\boldsymbol{S} + \lambda(\nabla \cdot \mathbf{u})\boldsymbol{I}$. Notice that if the flow is incompressible, $\nabla \cdot \mathbf{u}=0$, and we recover the simpler form, revealing the beautiful unity of the underlying physics.

The linear relationship $\boldsymbol{\tau} = 2\mu\boldsymbol{S}$ seems humble, but its consequences are rich and often counter-intuitive.

First, viscosity is not just about shear. Consider a flow that is purely stretching, like a fluid being pulled apart in the $x$-direction and squeezed in the $y$-direction. In such a flow, there might be no shearing motion at all. Yet, because the fluid is being deformed (stretched), the rate-of-strain tensor $\boldsymbol{S}$ is non-zero. The constitutive law then predicts that the normal stresses will be different! For a planar extensional flow with velocity $\mathbf{u} = (\dot{\epsilon}x, -\dot{\epsilon}y, 0)$, the law tells us that $\sigma_{xx} - \sigma_{yy} = 4\mu\dot{\epsilon}$. An ideal, inviscid fluid could not support such a difference in normal pressures; it would just squeeze out. A viscous fluid, however, resists being pulled apart, generating extra tension in the direction of stretching. This is a crucial effect in processes like polymer extrusion and glass blowing.

Second, where does the energy from stirring your coffee go? It doesn't speed up the coffee indefinitely. The work you do against the viscous forces is converted into thermal energy, warming the fluid. This process is called ​​viscous dissipation​​. The rate of this energy conversion per unit volume is given by the dissipation function, $\Phi = \boldsymbol{\tau} : \boldsymbol{S}$. For a Newtonian fluid, this becomes $\Phi = 2\mu(\boldsymbol{S}:\boldsymbol{S})$. Since viscosity $\mu$ is positive and $\boldsymbol{S}:\boldsymbol{S}$ (a sum of squared terms) can never be negative, dissipation is always a one-way street: mechanical energy is irreversibly lost to heat. This is the second law of thermodynamics manifesting itself through fluid friction, the relentless march of entropy.

Finally, consider a fluid flowing past a solid boundary, like water in a pipe. We almost always assume the fluid right at the wall is stationary—the ​​no-slip boundary condition​​. Is this a fundamental law of nature? Not quite. It's an empirical observation that arises from the intense molecular forces between the fluid and the solid. We can model this by considering a finite interfacial friction. This reveals that the degree of slip depends on a "slip length," which is typically nanometers in scale for common fluid-solid pairs. Because this length is so much smaller than the size of pipes or airplanes, the friction at the interface is effectively infinite, and the no-slip condition becomes an extraordinarily accurate approximation of reality.

The Newtonian model, in its elegant linearity, provides the foundation for much of fluid dynamics. But nature is full of rebels. Think of paint, which flows easily when brushed quickly but is thick otherwise, or ketchup, which refuses to leave the bottle until you shake it hard enough. These are ​​non-Newtonian fluids​​.

For these materials, the "Newtonian bargain" is broken. Viscosity is no longer a constant material property. Instead, they have an ​​apparent viscosity​​, $\eta$, that changes with the rate of strain. The constitutive law is generalized to $\boldsymbol{\tau} = 2\eta(\dot{\gamma})\boldsymbol{S}$, where $\dot{\gamma}$ is a measure of the overall magnitude of the strain rate, typically defined as $\dot{\gamma} = \sqrt{2\boldsymbol{S}:\boldsymbol{S}}$. For a shear-thinning fluid like paint, $\eta$ decreases as $\dot{\gamma}$ increases. For a shear-thickening fluid like a cornstarch-and-water mixture, $\eta$ increases.

By studying these more complex fluids, we gain an even deeper appreciation for the Newtonian model. It represents the simplest, most fundamental behavior of viscous fluids, a bedrock of linearity upon which our understanding of a vast and turbulent world is built. It is a testament to how a simple, well-posed physical hypothesis can unlock a universe of complex and beautiful phenomena.

In our exploration of physics, we often encounter laws of striking simplicity and breathtaking power. The linear relationship between stress and strain rate for a Newtonian fluid, $\tau = \mu \dot{\gamma}$, is a perfect example. At first glance, it might seem like a dry, technical definition confined to textbooks. But to think that is to miss the adventure. This simple rule is, in fact, a universal language spoken by a vast range of flowing matter, and understanding it allows us to read the secrets of phenomena from the microscopic to the planetary scale. It governs the blood in our veins, the magma beneath our feet, and the advanced materials in the devices we use every day. Let us embark on a journey across disciplines to witness this single, elegant principle in action.

Our journey begins in the world of human invention. Consider the burgeoning technology of 3D printing. An engineer designing a print head to extrude a molten polymer must know exactly how much material will flow for a given pressure. The answer is hidden in the constitutive relation. Inside the narrow nozzle, the fluid sticks to the walls and flows fastest at the center. This difference in speed creates a shear rate. By applying $\tau = \mu \frac{du}{dr}$, one can determine the shear stress at any point. This stress must balance the force from the pressure pushing the fluid through. By working backward from this fundamental balance, one can integrate across the nozzle's radius to build up the entire velocity profile and then integrate once more to find the total volumetric flow rate. The final result, a famous relationship known as the Hagen-Poiseuille equation, directly links the macroscopic flow rate to the fluid's viscosity and the nozzle's geometry—a powerful design tool born from a microscopic law of friction.

This principle is not just for plumbers and printers; it is at the heart of our most advanced technologies. In the sterile cleanrooms where microchips are made, processes like Plasma-Enhanced Chemical Vapor Deposition (PECVD) are used to lay down unimaginably thin films of material. To ensure that millions of transistors on a silicon wafer receive a perfectly uniform coating, engineers use complex computer simulations. These simulations must solve for the flow of precursor gases throughout the reactor chamber. The governing equations for that flow—the Navier-Stokes equations—are built upon the very foundation of the Newtonian constitutive relation, which defines how the gas resists being sheared. Getting this part right is the first and most crucial step in predicting how heat and chemical species are transported to the wafer surface, ultimately determining the quality and performance of the final microchip. From the humble pipe to the advanced reactor, engineering relies on understanding the internal friction of fluids.

Let's now turn our gaze from the engineered world to the natural one, expanding our scale from millimeters to kilometers. Deep within the Earth, molten rock, or magma, flows through conduits on its journey to the surface. Geophysicists seeking to understand and perhaps one day predict volcanic eruptions model this flow. As a first approximation, many types of magma behave as Newtonian fluids, albeit incredibly viscous ones. As the magma is forced upward, it drags against the stationary rock walls of the conduit. This drag is nothing more than a shear stress. By knowing the velocity profile of the flow—which is often nearly parabolic, like the flow in a pipe—scientists can use $\tau_w = \mu \left| \frac{du}{dr} \right|_{\text{wall}}$ to calculate the immense force exerted on the surrounding rock. This stress can contribute to seismic tremors and influence the dynamics of the eruption, reminding us that the same law that governs honey dripping from a spoon also describes the awesome power of a volcano.

Nowhere is the Newtonian constitutive relation more intimately and surprisingly at play than within the domain of living things. Life is bathed in, and built from, fluids. The principles of fluid dynamics are not external constraints; they are an integral part of biology's playbook.

Let's begin with the river of life: blood. The inner lining of our blood vessels is a delicate, single-cell-thick layer called the endothelium. These cells live in a world of perpetual flow, and their health is critically dependent on the mechanical forces exerted by the blood. In a healthy, straight artery, blood flows in smooth, orderly layers—a state known as laminar flow. The shear stress felt by the endothelial cells, typically around $1-7$ Pascals, is a vital, health-promoting signal. This gentle, steady rubbing stimulates the cells to produce molecules that keep the vessel relaxed and prevent clotting. It's a "good" stress that protects against atherosclerosis.

But what happens when this flow is disturbed? In a patient with a leaky heart valve, a high-pressure jet of blood can surge backward into a low-pressure chamber. This turbulent jet creates regions of ferociously high shear rates, sometimes thousands of times greater than normal. Our simple law, $\tau_w = \mu \dot{\gamma}_w$, tells us what happens next. A near-wall shear rate of, say, $9,000 \text{ s}^{-1}$ in blood with a viscosity of about $3.2 \times 10^{-3} \text{ Pa}\cdot\text{s}$ generates a shear stress of nearly $30 \text{ Pa}$. This is not a gentle, stimulating rub; it is a physical sandblasting. This pathologically high stress can physically damage and denude the endothelial lining, creating a raw spot. This injury site becomes a nidus for platelet aggregation and fibrin deposition, forming a sterile vegetation that, during a transient infection, can become colonized by bacteria, leading to the dangerous condition of infective endocarditis. The same physical law explains both the protective nature of healthy flow and the destructive mechanism of a disease.

The story gets even more subtle. Cells don't just survive stress; they read it and respond with astonishing precision. This process, mechanotransduction, is how physics gets translated into biochemistry. In a developing zebrafish embryo, the shear stress in the first primitive blood vessels, calculated directly from the blood's viscosity and the velocity gradient, provides the essential cue that tells arteries and veins how to form correctly. In a mature artery, the magnitude of the shear stress controls the expression of key protective genes like KLF2. In the laboratory, scientists can build microfluidic channels where they have complete control over the flow. By perfusing a culture medium over a layer of endothelial cells, they can dial in a precise wall shear stress. They observe that as the stress increases, the activation of genes like KLF2 follows a predictable, switch-like curve. This allows them to build quantitative models, such as Hill functions, that directly link the physical input, $\tau_w$, to the biological output, gene activation level. Physics becomes an input variable in the equations of life.

The influence of fluid shear permeates all of biology. During embryonic development, the very formation of our synovial joints—the knees, the elbows—depends on it. Fetal movements shear the primordial fluid between the nascent bones, and the resulting stress is a necessary signal to clear out the intervening cells and sculpt the joint cavity. A genetic defect that increases the fluid's viscosity would, according to our law, proportionally increase this shear stress, potentially disrupting this delicate developmental process. Even in the dental office, this physics is at work. During a root canal procedure, a key goal is to remove the bacterial biofilm clinging to the canal walls. The irrigation fluid is not just a chemical disinfectant; it is a mechanical tool. The effectiveness of the cleaning depends on whether the flow generates enough wall shear stress to overcome the biofilm's adhesion and scour it away.

And finally, consider a routine medical procedure: preparing a cytology slide. A technician takes a drop of sample and uses a second slide to smear it into a thin film. This simple action creates a planar Couette flow. The speed of the smear, $U$, and the height of the gap, $h$, define a shear rate, $\dot{\gamma} \approx U/h$. The cells in the sample experience a shear stress of $\tau = \mu (U/h)$. If the technician smears too quickly, this stress can exceed the critical threshold for membrane rupture, destroying the very cells they wish to examine and rendering the sample useless for diagnosis.

Our journey has taken us from the heart of a microchip to the heart of a volcano, and into the very heart of a living cell. In every case, the simple, linear relationship between stress and shear rate for a Newtonian fluid was not just an incidental feature; it was a central character in the story. It is a design principle in engineering, a diagnostic clue in geophysics, a signal for life and a mechanism of disease in medicine. The discovery that such a diverse and complex set of phenomena can all be understood through such an elegant and simple physical law is a profound testament to the unity and beauty of science. It reminds us that if we learn to speak the language of nature's fundamental rules, we can begin to understand worlds both strange and familiar.