3.3 管道流动与泊肃叶定律

3.3 管道流动与泊肃叶定律

本节摘要：充分发展的圆管层流是 N-S 方程少数能精确求解的情形之一：速度剖面为抛物线，流量与单位长度压降、半径的四次方成正比（泊肃叶定律）。本节给出推导脉络与 Python 复算，讨论入口段长度，并给出从层流精确解过渡到湍流工程摩阻（达西公式、穆迪图）的路径。

医院输液、毛细血管灌注、色谱柱内的载液、精密仪器的润滑油膜——这些流动有一个共同点：尺度小、粘性大、速度低，处于层流区。而它们恰好落在 N-S 方程能够精确求解的那个窗口里。医生用泊肃叶定律换算输液滴速与滴注时间，并不是在"套经验公式"，用的是有严格解析解的定律。

一、从 N-S 方程到抛物线剖面

设定物理场景：水平等径圆管，流动充分发展（速度剖面沿流向不再变化，当地与沿程加速度均为零），不可压缩牛顿流体。此时 N-S 方程里几乎所有项清零，只剩压力梯度与粘性项的平衡——轴向压力推力被管内剪应力平衡。对这个简化方程在圆截面上一积分，得到：

速度剖面 u(r) = (Δp / (4 μ L)) (R² − r²)——一条抛物线：管中心最快，壁面为零（无滑移）。最大速度是平均速度的两倍。再对剖面积分得泊肃叶定律：

Q = π R⁴ Δp / (8 μ L)。

这个四次方是本节最重要的数字敏感性：管径加倍，同压差下流量涨 16 倍。血管造影术后压迫止血、喷墨打印的墨滴量控制，背后都是这条四次方律——管径毫米级的变化就是数倍的流量变化。

import numpy as np

mu, L, dp = 3.5e-3, 1.0, 5e3      # 血液近似粘度 3.5 mPa s
for R_mm in [1.0, 1.5, 2.0, 3.0]:
    R = R_mm * 1e-3
    Q = np.pi * R**4 * dp / (8 * mu * L)
    print(f"半径 {R_mm:.1f} mm  流量 {Q*1e6:8.1f} mL/s  "
          f"最大流速 {2*Q/(np.pi*R**2)*1e3:6.1f} mm/s")
​

跑一遍输出：半径 1 mm 时约 0.11 mL/s，3 mm 时 9 mL/s——半径三倍、流量 81 倍。微循环里小动脉半径的微小舒缩，就是血流量最有效的调节旋钮，比血压调节敏感得多。

二、壁面剪应力与摩擦系数

对结果再做一步加工：壁面剪应力 τw = Δp R / (2L)，把它无量纲化得到范宁摩擦系数，层流结果是 f = 16/Re（达西摩阻系数则是 λ = 64/Re，只差四倍口径，用谁的表就要配谁的公式）。层流摩擦只依赖雷诺数、与壁面粗糙度无关——粗糙凸起完全淹没在层流里，根本"探"不到主流。这一点与湍流形成鲜明对比。

三、入口段与"充分发展"的代价

泊肃叶解成立的前提是充分发展。流体从入口进入直管后，边界层从壁面逐渐长厚，要走过一段入口段长度 Le ≈ 0.06 Re D（层流）后才布满全管、剖面定格为抛物线。Re = 2000 时 Le 约 120 倍管径——很多实验台架的直管段不够长，测到的其实是还在发育的剖面，压降比泊肃叶公式偏高。做粘度测定或微流控标定时，要么加长直管，要么修正入口段损失。

四、层流之外的工程账：达西公式与穆迪图

流速升高进入湍流区后（临界雷诺数约 2300，详见第 6 章），泊肃叶定律失效，压降增长比线性快得多。工程标准算法是达西-魏斯巴赫公式：

hf = fD · (L/D) · U²/2g，摩擦系数 fD 由雷诺数与相对粗糙度从穆迪图（或 Colebrook 迭代）查出。层流区它在图中是一条直线 fD = 64/Re，与泊肃叶完全一致——精确解与工程图在这一点无缝衔接。

# 同一条管 层流与湍流压降对比（粗糙度 0.05 mm 钢管 Swamee-Jain 近似）
import numpy as np

rho, mu, D, eps, L, g = 998., 1.0e-3, 0.05, 5e-5, 50., 9.81
for U in [0.2, 1.0, 3.0]:
    Re = rho * U * D / mu
    if Re < 2300:
        fD = 64.0 / Re
        regime = "层流"
    else:
        fD = 0.25 / (np.log10(eps/(3.7*D) + 5.74/Re**0.9))**2
        regime = "湍流"
    hf = fD * (L/D) * U**2 / (2*g)
    print(f"U={U:4.1f} m/s  Re={Re:7.0f}  {regime}  摩阻水头 {hf:6.2f} m")
​

输出里能看到残酷的对比：流速从 0.2 提到 3 m/s（15 倍），层流段压降涨约 15 倍（线性），跨入湍流后同管压降涨到近 40 米水头——接近二次律再乘一个粗糙度敏感的摩擦系数。这就是为什么长距离输水靠大管径低流速，而不是小管径高流速：泵的电费按三次方涨（流量一次方乘压降二次方）。

⚠️ 常见坑：拿层流公式算湍流管路。Re 一过 2300，真实压降可能比泊肃叶结果高出几倍，选泵必然偏小。任何管路计算第一步是算雷诺数、判流态，再选公式。

五、局部损失：配件的账单

管路里弯头、三通、阀门、突缩突扩各自产生"局部损失"，工程上折算成损失系数：h_m = K·U²/2g。典型值量级：90° 弯头 K 约 0.9，全开闸阀约 0.2，突缩约 0.4，止回阀约 2.5。粗算一段抽水管线：DN100、流速 1.2 m/s、直管 80 m（fD 约 0.025）、弯头 6 个、闸阀 2 个：

U, D, L, fD = 1.2, 0.1, 80.0, 0.025
hf = fD * L/D * U**2/(2*9.81)
Ksum = 6*0.9 + 2*0.2 + 0.4 + 2.5
hm = Ksum * U**2/(2*9.81)
print(f"沿程损失 {hf:.2f} m  局部损失 {hm:.2f} m  合计 {hf+hm:.2f} m")
​

沿程约 1.8 m，局部约 1.3 m——配件吃掉了四成损失。长输管线沿程占绝对主导，短管线、多阀件的机房管路则局部为王。设计优化方向因此完全不同：前者降流速选大管径，后者减少弯头、优选低阻阀门、避免断面突变。

本节要点回顾

• 充分发展圆管层流是 N-S 方程的经典精确解：抛物线剖面，u_max = 2 ū；
• 泊肃叶定律 Q 与半径四次方成正比：管径的小变化被四次方放大成流量的大变化，微循环调控的核心机理；
• 入口段长度 Le ≈ 0.06 Re D，实验与标定必须确认"充分发展"成立；
• 层流摩擦 f = 16/Re 与粗糙度无关；湍流摩擦同时依赖 Re 与相对粗糙度；
• 达西公式加穆迪图是跨流态的工程通用算法，先判 Re 再选公式是铁律。

到此方程组有了第一个精确解。下一章换一条思路：干脆丢掉粘性，看理想化世界能解释什么、又在哪露馅。