Hagen–Poiseuille law and vascular resistance. Since Q∝R⁴, halving the radius cuts the flow rate to 1/16. R: vessel radius [m], ΔP: pressure drop [Pa], μ: dynamic viscosity [Pa·s], L: length [m].
Wall shear stress τ_w [Pa]. A physiological range of 1–7 Pa is typical for atherosclerosis risk assessment.
What is Vascular Fluid Mechanics?
🙋
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.
Enter vessel radius (mm) and stenosis ratio (0–1, where 1 = complete occlusion) to define geometry
Set vessel length (mm) and pressure gradient (mmHg/mm) using sliders or direct input
Input blood dynamic viscosity (cP, typically 3.5–4.5 for normal blood) and click Calculate
Review flow rate Q in ml/min, wall shear stress τ_w in Pa, Reynolds and Womersley numbers, and stenosis risk flag
Worked Example
Femoral artery: radius 3 mm, stenosis ratio 0.6 (60% diameter reduction), length 80 mm, pressure gradient 0.05 mmHg/mm, viscosity 3.8 cP. Simulator outputs: Q = 180 ml/min, mean velocity 64 cm/s, τ_w = 2.8 Pa, Re = 320, Wo = 3.2. WSS Risk Assessment shows elevated risk because τ_w < 0.4 Pa indicates flow stasis; combined with high stenosis ratio triggers atherosclerosis warning for intervention planning.
Practical Notes
Womersley number (Wo) > 10 confirms pulsatile flow dominance in large arteries; Wo < 1 signals viscous-dominated capillary flow—use for vessel classification
Wall shear stress 1.0–3.0 Pa protects endothelium; τ_w < 0.4 Pa or > 4 Pa triggers atherosclerotic plaque formation—adjust stenosis inputs to evaluate intervention outcomes
Stenosis ratio 0.5+ combined with Re > 400 produces turbulent jets; model critical narrowings in carotid or coronary arteries for procedural planning
Blood viscosity rises in dehydration (5 cP) and sickle cell disease—rerun simulations to assess flow compromise