参数设置
流体假定为水（ρ = 1000 kg/m³）。公式在 Re < 2300 的层流区严格成立。
圆管纵断面与速度分布
左：入口（高压）/ 右：出口（低压）/ 管壁的两条线 = 管半径 R / 箭头 = 层流的抛物线速度分布（中心最大，管壁为零）
体积流量 Q 与压力差 ΔP
横轴 = ΔP [kPa] / 纵轴 = Q [mL/s] / 层流中 Q ∝ ΔP 的直线关系（黄点 = 当前值）
理论与主要公式
哈根-泊肃叶定律为圆管内完全发达层流的体积流量给出解析表达式。它从纳维-斯托克斯方程推导而得，在牛顿流体、层流、圆管、恒定截面的条件下严格成立。
体积流量与平均流速：
$$Q = \frac{\pi R^{4}\,\Delta P}{8\,\mu\,L},\qquad V_{\mathrm{avg}} = \frac{Q}{\pi R^{2}} = \frac{\Delta P\,R^{2}}{8\,\mu\,L}$$
抛物线速度分布（最大流速是平均的2倍）：
$$u(r) = V_{\max}\!\left(1-\frac{r^{2}}{R^{2}}\right),\qquad V_{\max} = 2\,V_{\mathrm{avg}}$$
层流判定的雷诺数：
$$Re = \frac{\rho\,V_{\mathrm{avg}}\,(2R)}{\mu}$$
$R$ 为管半径 [m]，$L$ 为管长 [m]，$\Delta P$ 为两端压力差 [Pa]，$\mu$ 为粘性系数 [Pa·s]，$\rho$ 为密度 [kg/m³]，$r$ 为自管中心向径向的距离。
哈根-泊肃叶流 模拟器介绍
🙋
听说输液管的滴速跟针头的粗细、药液的粘性和瓶子高度有关，有计算公式吗？
🎓
没错，那就是哈根-泊肃叶定律 $Q = \pi R^{4}\Delta P/(8\mu L)$。流量与针头半径 R 的4次方成正比，与压力差 ΔP 成正比，与针头长度 L 和液体粘性 μ 成反比。R⁴ 这一项很关键——针头直径减半，流量就会减少到 1/16。这就是为什么在医疗现场"选择粗针还是细针"对流量影响这么大。
🙋
流量与半径的4次方成正比！这样的话，管子稍微细一点流量就会大幅下降呢。
🎓
完全正确。动脉硬化导致血管稍微变窄，血流也会大幅减少，道理是一样的。用模拟器试试把 R 从 2 mm 改成 1 mm，你会看到流量精确地变成 1/16。反过来说，压力差 ΔP 只有1次方，翻倍压力也只能让流量翻倍——"加粗管子"比"增加压力"效率高得多。
🙋
速度分布图中，中心的箭头最长，靠近壁面的箭头变短，这表示什么呢？
🎓
这就是著名的"抛物线速度分布"。公式是 $u(r) = V_{\max}(1 - (r/R)^{2})$，在管中心 r=0 处最大，在管壁 r=R 处为0。速度在壁面为0叫做"无滑移条件"，是流体力学的基本假设。把这个分布在截面上求平均，你会得到平均流速是最大流速的一半，即 $V_{\max} = 2\,V_{\mathrm{avg}}$。默认参数（R=2 mm, L=10 m, ΔP=10 kPa, μ=0.001 Pa·s）下，V_avg = 0.5 m/s，V_max = 1.0 m/s。
🎓
很好的问题。哈根-泊肃叶公式只在"层流"的前提下才能精确成立。所以要时刻监测雷诺数，确认"现在真的还是层流吗？"。经验值是 Re < 2300 为层流，2300～4000 为过渡区，4000 以上为湍流。默认参数下 Re 约为 2000，接近层流的临界线。当你增加 ΔP 时，Re 也会成比例增加，最后会进入湍流区。进入湍流后，需要用另外的公式——达西-韦斯巴赫式和穆迪图。建议也试试"穆迪图"模拟器。
常见问题
在层流中，纳维-斯托克斯方程里的惯性项可以忽略，运动方程只需平衡粘性项和压力项，这样就能解析求解。解出来就是抛物线速度分布，对截面积分后得到 Q = πR⁴ΔP/(8μL)。一旦进入湍流，惯性项主导产生的湍流应力就不可忽略，方程失去解析解。通常当 Re 超过约2300时流动就会遷入湍流，此时 Q ∝ ΔP 的线性关系被破坏，变成 Q ∝ √ΔP 的关系。此时需用实验相关式或数值计算。
在大动脉等大血管中流动是脉动流，严格说不适用，但毛细血管、细动脉、细静脉等低 Re（数到数十）的细血管中，这个公式是很好的近似，医学和生理学中广泛应用。不过血液不是牛顿流体，有非牛顿特性（特别是在低剪切速率时黏度增大），需要加入见观粘度修正。现代血流模拟通常以哈根-泊肃叶式为出发点，逐步加入非牛顿项、脉动项、血管壁弹性项等来提高精度。
在管入口处，速度分布接近均匀（塞流）。随着流体向下游，壁面的粘性阻力作用范围逐渐扩展到管中心，速度分布逐渐演变为抛物线。经过一定的距离（入口段长度）后，速度分布稳定不变，这种状态就叫"完全发达"。层流中入口段长度约为 L_e ≈ 0.06·Re·D（D是管径），当 Re=2000、D=4 mm 时，L_e ≈ 0.48 m。本模拟器计算的是完全发达区的值，入口附近的值略有不同。实际设计时要么选择足够长的管，要么单独计算入口损失系数。
直接用不了，但层流中每种截面形状都有各自的解析解。比如边长为 a 的正方形风道，流量公式为 Q = K·a⁴·ΔP/(μ·L)，其中 K ≈ 0.0351。实际工作中常用水力直径 D_h = 4A/P（A为截面积，P为润湿周长）来近似替换圆管公式，但层流中误差通常在10～20%。精密计算还是要用针对具体形状的解析解或者数值计算（CFD）。
实际应用
医疗：输液、输血、导管设计：输液滴速由针头内径和瓶子高度（静压）控制，设计基础公式就是哈根-泊肃叶式。由于半径效应 R⁴，规格化的针号（18G、20G、22G等）选择是决定药液到达速度的关键因素。肾透析血路、麻醉药持续注射、CT/MRI造影剂注入等都需用这个公式估算流量和导管径的关系。
微流体和芯片实验室（LOC）：微流体芯片通常有几十到几百微米的流路，雷诺数基本 < 1，层流是前提。哈根-泊肃叶式（以及矩形截面的解析解）是流路设计的标准公式，应用于 PCR、细胞培养、蛋白质结晶等生物器件。流路越细，ΔP 的 R⁴ 项增长越快，泵压力和流路尺寸的权衡是设计中心。
地下水和多孔介质渗流：多孔介质中的流动用达西定律描述 q = -(k/μ)∇P，这个定律实际上可以从多孔介质中孔隙的哈根-泊肃叶流平均推导出来（科泽尼-卡门公式的出发点）。地下水流动分析、油气储层模拟、土壤污染物扩散、CO₂地层封存等都本质上是众多哈根-泊肃叶流的集合体。
润滑和摩擦学：轴承、齿轮间隙内的润滑油流动也是薄隙内的层流，用哈根-泊肃叶式和同样的框架处理（库埃特-泊肃叶复合流）。它是雷诺方程的出发点，用于从油膜厚度、粘度、转速估算摩擦损失和发热量。
常见误区与注意事项
最常见的误区是直观地认为"增加压力流量就无限增加"。看哈根-泊肃叶式，Q ∝ ΔP 似乎是线性的，但这只在"维持层流"的前提下成立。压力增大使 Re 增大，一旦超过2300进入湍流，流量与压力的关系就变成 Q ∝ √ΔP，即便压力翻倍，流量也只增加约1.4倍。模拟器中当 statRe 超过2300后，显示的 Q 值是"如果还是层流的话"的理论值，实际流量会更小。
第二个常见误会是把"半径"和"直径"搞混，把"平均速度"和"最大速度"搞混。哈根-泊肃叶式里的 R 是半径，不是直径 D。直径形式的公式 Q = πD⁴ΔP/(128μL) 和半径形式的系数相差16倍，用错了后果很大。同样，Re = ρVD/μ 用的是直径和平均速度，V_max 是 V_avg 的2倍，这些细节的定义必须确认清楚。本工具统一用半径 R 输入、平均速度 V_avg、直径制雷诺数。
最后一个误区是"水粘度小，不用管"的想法。水在20°C时粘度约1.0 mPa·s，很小，但温度和添加剂能大幅改变。甘油1500 mPa·s（水的1500倍），发动机油100～500 mPa·s。哈根-泊肃叶公式中流量与 1/μ 成正比，粘度差10倍流量就差10倍。用模拟器把 μ 从0.0010（水）改到0.01、0.1（油），你会看到流量剧减。CAE分析中粘度系数输入错误是导致结果偏离的典型原因。