What exactly is the Hagen-Poiseuille equation calculating in this simulator? I see it gives a flow rate Q.
🎓
Basically, it's the fundamental law for steady, laminar flow in a straight pipe. It tells you how much blood volume flows per second, given the vessel's geometry and the pressure pushing it. In practice, it shows that flow is incredibly sensitive to the radius—double the radius, and the flow increases by 16 times! Try moving the "Vessel Radius R" slider above to see this dramatic effect.
🙋
Wait, really? So a tiny plaque buildup that narrows the radius would massively reduce flow? What's the "Wall Shear Stress" number telling me then?
🎓
Exactly! That's the core of atherosclerosis risk. Wall Shear Stress (WSS), $\tau_w$, is the frictional force blood exerts on the vessel wall. Low or oscillating WSS can signal spots where plaque is likely to form. In the simulator, when you add a "Stenosis," watch how the WSS spikes in the narrow section—that's a dangerous mechanical signal for endothelial cell damage.
🙋
Okay, so this is for steady flow. But the heart beats! What does the "Heart Rate" parameter do? Does the Poiseuille equation still work?
🎓
Great question! For pulsatile flow, we need the Womersley number. It compares the pulse frequency to viscous effects. A high Womersley number (from high Heart Rate or large radius) means flow inertia dominates and the velocity profile is blunted, not parabolic. The simulator uses it to check if the steady-flow assumption is valid. Try cranking up the HR and see how the flow profile changes.
Physical Model & Key Equations
The primary equation governing steady, laminar flow in a long, straight, rigid cylindrical tube is the Hagen-Poiseuille Law. It derives from the Navier-Stokes equations and assumes a parabolic velocity profile.
$$Q = \frac{\pi R^4 \Delta P}{8 \mu L}$$
Q = Volumetric Flow Rate [m³/s] R = Inner Radius of the vessel [m] ΔP = Pressure drop along the vessel length [Pa] μ = Dynamic Viscosity of blood [Pa·s] L = Length of the vessel segment [m]
Wall Shear Stress (WSS) is the key hemodynamic force acting on the endothelial cells lining the vessel. It can be calculated directly from the flow rate or the pressure drop.
τ_w = Wall Shear Stress [Pa]. This is a critical biomarker. Physiological WSS (≈1-2 Pa) is protective. Low WSS (<0.5 Pa) promotes atherosclerosis, while very high WSS (>4 Pa) can cause mechanical damage.
Frequently Asked Questions
According to the Hagen-Poiseuille equation, the flow rate is proportional to the fourth power of the radius. Therefore, a 10% reduction in radius results in approximately a 34% decrease in flow rate. This dramatic change with just a slight adjustment of the stenosis slider is due to this physical law, which is an important point in assessing arteriosclerosis risk.
The Womersley number (α) is a dimensionless parameter that indicates the ratio of inertial forces to viscous forces in blood flow. When α is small, the flow approximates steady flow (Poiseuille flow), and when it is large, the influence of pulsatile flow becomes stronger. In this tool, α changes according to the vessel type and degree of stenosis, allowing evaluation of the flow characteristics.
When a vessel type (e.g., aorta, coronary artery, cerebral artery) is selected, standard values for vessel radius, length, and pressure difference are automatically set. This allows easy comparison of blood flow and wall shear stress across different regions, enabling simulation of sites prone to arteriosclerosis.
This tool is a simulator for educational and research purposes and cannot be directly used for actual diagnosis or treatment decisions. Since it assumes ideal steady laminar flow, it does not reflect individual patient vessel geometry or pulsatile effects. Please use it solely for qualitative understanding of arteriosclerosis risk.
Real-World Applications
Vascular Stent Design: Engineers use these exact calculations to optimize stent geometry. A poorly designed stent can create regions of low WSS downstream, inadvertently encouraging restenosis (re-narrowing). CFD simulations are validated against these hand calculations.
LVAD & Artificial Heart Development: When designing Left Ventricular Assist Devices, engineers must ensure the pump and connected grafts produce physiological WSS levels to prevent blood clot formation. This simulator's principles are used for first-pass sizing and analysis.
Aneurysm Risk Assessment: In cerebral or aortic aneurysms, the wall is weakened. High WSS at the aneurysm's "neck" can promote growth and rupture, while low, swirling flow inside the sac can lead to clot formation. Simple models help triage complex CFD cases.
Microfluidic Lab-on-a-Chip Design: Biomedical devices that manipulate blood in tiny channels for diagnostics rely on precise control of flow and shear. The Poiseuille equation is fundamental for designing these chips to ensure cells experience the correct mechanical forces without damage.
Common Misconceptions and Points to Note
When you start using this simulator, there are several points beginners often stumble on. First and foremost, understand that the radius-to-the-fourth-power law is not absolute. While Poiseuille's law $Q \propto R^4$ is powerful, it strictly applies under ideal conditions: "straight, rigid circular tubes," "Newtonian fluid," and "steady flow." Actual blood vessels are curved, elastic, and carry pulsatile flow. For instance, doubling the radius would theoretically increase flow 16-fold, but in living systems, autoregulatory mechanisms adjust to maintain constant pressure, preventing such a simple increase. Directly comparing calculation results with clinical data is risky.
Next, beware of the pitfall of "fixing the inlet-outlet pressure difference ΔP." The simulator treats ΔP as a fixed parameter, but in the actual circulatory system, cardiac output (≈ flow rate Q) is maintained nearly constant. When stenosis occurs, ΔP increases to maintain this flow (increasing cardiac workload). The result that "increasing stenosis reduces flow Q" assumes ΔP remains unchanged. In practice, you must always be conscious of "which parameter you consider the independent variable."
Finally, consider the treatment of blood viscosity μ. Here, viscosity is approximated simply using "hematocrit," but real blood is a non-Newtonian fluid, meaning its viscosity changes with shear rate (flow speed). Viscosity may differ near the vessel center versus the wall. This approximation is sufficient for a learning tool, but if used for research or advanced design, you need to keep the potential error from this assumption in mind.