Visualize Couette flow and viscosity curves for Newtonian, power-law, Bingham, and Herschel-Bulkley fluids in real time. Includes Arrhenius temperature correction and Reynolds number calculation.
What exactly is rheology? I've heard of viscosity, but this sounds more complex.
🎓
Basically, rheology is the study of how materials flow and deform. Viscosity is a key part of it, but rheology covers all kinds of fluid behavior—not just simple, constant-thickness liquids. In practice, think of ketchup that won't pour until you shake it, or paint that spreads easily but doesn't drip. Try selecting the "Bingham Plastic" model in the simulator above to see a fluid that behaves just like that ketchup.
🙋
Wait, really? So a fluid's "thickness" can change? How does that work in a simple model like Couette flow, where one wall moves?
🎓
Exactly! In Couette flow, the shear rate—how fast the fluid layers slide past each other—is constant. For a Newtonian fluid like water, viscosity stays the same. But for a non-Newtonian one, the viscosity you measure depends on that shear rate. For instance, select the "Power Law" model and watch the velocity profile. Changing the model parameters will show you fluids that thin (like paint) or thicken (like cornstarch paste) as the shear increases.
🙋
The FAQ mentioned numerical instability with Bingham fluids and "regularization." What's that about, and why is it important for simulation?
🎓
Great question! A perfect Bingham plastic has a yield stress—it won't flow at all until you push hard enough. That sharp "on/off" switch creates a mathematical discontinuity that crashes simulations. In practice, CAE tools like OpenFOAM use a smoothed model, like Papanastasiou's. It's a clever trick that adds an exponential term to make the transition smooth for the solver. Try the "Bingham (Regularized)" option in the simulator to see a stable, computable version of that ketchup-like behavior.
Physical Model & Key Equations
The core concept is the relationship between shear stress ($\tau$) and shear rate ($\dot{\gamma}$). For Newtonian fluids, this is a simple linear law defined by a constant viscosity ($\mu$).
$$ \tau = \mu \dot{\gamma}$$
Here, $\tau$ is the shear stress [Pa], $\mu$ is the dynamic viscosity [Pa·s], and $\dot{\gamma}$ is the shear rate [1/s]. This is the baseline behavior of fluids like water and air.
Non-Newtonian models introduce a variable "apparent viscosity" ($\eta(\dot{\gamma})$). A common and powerful model is the generalized Newtonian fluid formulation, where the stress is still proportional to shear rate, but the proportionality factor changes with flow conditions.
$$ \tau = \eta(\dot{\gamma}) \dot{\gamma}$$
The function $\eta(\dot{\gamma})$ defines the rheology. For example, the Power Law model is $\eta = K \dot{\gamma}^{\,n-1}$, where $K$ is consistency and $n$ is the flow index. If $n \lt 1$, the fluid thins with shear (shear-thinning); if $n \gt 1$, it thickens (shear-thickening).
Frequently Asked Questions
If you select 'Newtonian fluid' in the model selection on the left, the viscosity becomes a constant straight line regardless of the shear rate. On the other hand, if you select 'Power law' or 'Bingham plastic', the viscosity curve changes, and the velocity profile also changes to a parabola or plug flow, allowing you to visually compare the differences in real time.
In the 'Temperature compensation' tab, enter the reference temperature and the current temperature, and the viscosity will be automatically adjusted based on the Arrhenius-type equation (viscosity ∝ exp(Ea/RT)). This allows you to easily estimate the effect of temperature fluctuations in actual processes on flow characteristics.
It assumes a steady shear flow between parallel plates (fixed gap). The upper plate moves at a constant speed, while the lower plate is fixed. The velocity distribution is analytically obtained from the constitutive equations of each model and plotted, and you can also check the shear rate and stress distribution simultaneously.
Set it according to the shape of the flow channel. For example, enter the diameter for pipe flow, or the gap width for flow between parallel plates. The Reynolds number is automatically calculated from this value, along with the flow velocity, density, and viscosity, allowing you to check the laminar/turbulent flow criterion (Re < about 2000) in real time.
Real-World Applications
Food Processing: Designing pumps and pipes for products like yogurt, mayonnaise, or chocolate requires precise rheology models. A shear-thinning fluid like ketchup needs different handling than a Newtonian fluid like cooking oil to ensure consistent filling and packaging.
Polymer Processing & 3D Printing: Molten plastics are strongly shear-thinning. Accurate rheology data is critical for simulating injection molding or extrusion processes to predict flow into molds and ensure part quality without defects.
Drilling Mud in Oil & Gas: Drilling muds are often modeled as Bingham plastics. They must have a yield stress to suspend rock cuttings when circulation stops, but flow easily when pumped. CAE simulations of wellbore hydraulics rely on these models for safety and efficiency.
Biomedical Flows: Blood is a complex non-Newtonian fluid, exhibiting shear-thinning behavior. Understanding its rheology is vital for simulating blood flow in arteries, designing heart assist devices, and developing diagnostic equipment.
Common Misconceptions and Points to Note
When you start using this tool, there are a few common pitfalls to watch out for. The first is thinking that the power-law exponent 'n' alone tells you everything about the fluid's properties. It's true that n<1 indicates shear-thinning, but that only means there's "a tendency for viscosity to decrease as shear rate increases." Real materials, like polymer melts, often deviate from the power law at extremely low or high shear rates. Try comparing n=0.3 and n=0.8 in the tool. The trends are similar, but the way viscosity drops is completely different, right? In practice, it's rare for a single model to fit all measured data points; sometimes you need to switch models for different shear rate regimes.
The second point is the interpretation of the Reynolds number. The tool calculates it assuming pipe flow, but this value is only a "guideline." For instance, a Reynolds number exceeding 2300 doesn't guarantee turbulent flow. It varies greatly with channel geometry and inlet conditions. Use it as preliminary information for choosing CFD meshing or analysis methods, nothing more.
Finally, avoid simplistically equating Bingham fluids with ketchup. While ketchup does have a yield stress, it also has strong time-dependent (thixotropic) behavior. The tool's Bingham model is a highly idealized model that "behaves like a Newtonian fluid once it starts flowing." The real thing is more complex; it can solidify again after being stirred and left to stand. The tool teaches you the entry point to the concept of "yield stress." Keep in mind that actual product design often requires more complex models.
Select fluid type: Newtonian (constant viscosity) or non-Newtonian (power-law or Bingham)
Set viscosity (muVal) in Pa·s; for non-Newtonian fluids, adjust flow behavior index (nVal, typically 0.5–1.5) and yield stress (tau0Val in Pa) for Bingham models
Enter temperature (tempVal in °C); observe real-time viscosity curve updates and Couette flow velocity profile changes in the shear-rate vs. shear-stress plot
Monitor the strain-rate gradient and resulting shear thinning or thickening behavior on the rheology diagram
Worked Example
Simulate SAE 10W-30 engine oil at 40°C: set muVal=100 mPa·s (0.1 Pa·s), nVal=1.0 (Newtonian). At shear rate γ̇=1000 s⁻¹, shear stress τ=100 Pa. Now heat to tempVal=100°C; viscosity drops to ~35 mPa·s, yielding τ=35 Pa at same γ̇. For ketchup (Bingham fluid) at 25°C: tau0Val=50 Pa, muVal=5 Pa·s, nVal=0.8; no flow initiates until shear stress exceeds 50 Pa threshold, then velocity increases non-linearly with applied shear.
Practical Notes
Polymer melts exhibit strong shear thinning (nVal <1): increasing nIdxNum from 0.3 to 0.7 flattens viscosity curves, critical for injection molding at γ̇=100–10,000 s⁻¹
Temperature sensitivity: every 10°C rise reduces viscosity by ~20–50% for oils and polymers; use tempVal slider to validate Arrhenius predictions