4.2 库埃特流动与泊肃叶流动

4.2 库埃特流动与泊肃叶流动

本节摘要：库埃特流与泊肃叶流是纳维-斯托克斯方程极少数的精确解：前者由平板拖动产生线性速度剖面，后者由管道两端压差产生抛物线剖面。两者合起来构成粘性流动的"字母表"——边界层剖面、润滑理论、管路压降计算，都是这两个解的变奏。本节完整走一遍求解过程，让你第一次看到"方程吐出一个函数"长什么样。

3.3 节说过，N-S 方程的精确解需要"几何简单到速度梯度结构固定"。现在兑现承诺，把最容易的两道解完整算出来。求解路径相同：先分析几何定出哪些速度分量为零，再用连续方程砍方程，最后剩下的常微分方程一积就完。

库埃特流：被拖动的流体

两块无限大平行板，间距 h，下板固定，上板以速度 U 拖动，流体充填其间。几何的对称性告诉我们：速度只有 x 方向分量 u，且它只随 y 变化（u = u(y)），压强沿 x 均匀。代回运动方程后，粘性项与惯性项自动平衡，方程退化为：

μ d²u/dy² = 0，边界条件 u(0) = 0，u(h) = U
解：u(y) = U · y / h —— 一条直线
​

没有压差、没有惯性，流体完全靠上板"搓"动，剪切率 d u/d y 处处等于 U/h，剪应力处处等于 μU/h——直线剖面意味着均匀剪应力，这是库埃特流的指纹。油膜轴承、旋转粘度计、两片玻璃间的液体流，全是库埃特流的现实版本。第 5 章的边界层内，靠近壁面的速度剖面也近似线性（那一段剪切率几乎恒定），等于把库埃特流当局部样板。

泊肃叶流：被压差推动的流体

固定不动的一根直圆管，半径 R，两端压强差 Δp 驱动流体。此时速度剖面由压差力与粘性力平衡决定，对圆管求解（把直角坐标换成柱坐标，积分两次，用管壁无滑移与轴线对称定常数）：

u(r) = (Δp / 4μL) · (R² − r²) —— 抛物线
平均速度 V_avg = R²Δp / 8μL = 最大速度的一半
​

管壁处速度为零（无滑移），轴线处最快，全剖面是旋转抛物面。最大速度恰是平均速度的两倍，这个系数值得记——用轴线流速估算流量时会差一倍。压降与流量的关系（哈根-泊肃叶定律）从平均速度直接导出：

压降 Δp = 8 μ L Q / (π R⁴)
​

最后那个四次方是管路设计的核心教训：管径减半，同样流量下压降增为十六倍。血管支架、微流控芯片、长距离输液管，全是这条四次方定律的受害者或受益者。

图：库埃特线性剖面与泊肃叶抛物线剖面对照

动手算：一支滴眼液的压降

滴眼液瓶口细管内径约 1 mm（R = 0.5 mm），管长 15 mm，药液粘度接近水（μ ≈ 0.001 Pa·s），希望一秒钟挤出 0.05 mL（Q = 5×10⁻⁸ m³/s）。需要的压差：

import math
mu, R, L, Q = 0.001, 0.0005, 0.015, 5e-8
dp = 8 * mu * L * Q / (math.pi * R**4)
print(dp)   # 约 1224 Pa
​

约 1.2 kPa——轻松超出。把 R 改成 0.4 mm 重算，压降跳到约 3 kPa：口径缩两成，阻力翻倍还多，四次方定律在毫米尺度上毫不留情。同时算一下雷诺数（ρ ≈ 1000 kg/m³，管径 0.001 m，平均流速约 0.064 m/s）：Re ≈ 64，远低于管流转捩的两千上下——泊肃叶解合法。这个"先判雷诺数、再套公式"的顺序，是使用一切层流精确解的前提。

两解合体与工程变奏

上下板同时拖动又有压差时，剖面是线性与抛物线的叠加，总流量可正可负（压差不足以克服拖板的倒吸时出现"回流区"）。这个"合体解"是润滑理论的种子：轴承间隙里楔形油膜的承重能力，正来自楔形缝隙中类似叠加剖面产生的压强场。

变奏：两个解的合体与血液动力学的账

上下板一个拖动、两端又加压差时，速度剖面是线性与抛物线的代数和，形状随压差正负在"钝三角"与"带拐点的 S 形"之间切换。这个合体解有个漂亮的工程应用——滑块轴承：楔形缝隙里的油膜相当于间距缓慢变化的库埃特流叠加泊肃叶流，缝隙收窄处产生远高于环境的压强场，把几千克的转轴稳稳托起。大型汽轮机轴承的油膜压强可达兆帕量级，托举的全部资金来自"两个精确解的代数和"。

生物流体里泊肃叶定律是主角。小动脉管径减半，同样的血流量需要十六倍的压差去推（四次方定律），心脏的后负荷随之飙升——高血压与血管狭窄的流体力学解释就在这一行公式里。注意血液并非牛顿流体（红细胞使其剪切变稀），大血管里泊肃叶解只是首阶近似，微循环里则要换用非牛顿模型——精确解的适用前提时刻要跟，这一条与雷诺数判层流的纪律同源。

问题：管流转捩的雷诺数为什么和翼型转捩差那么多？

管道流有"现成的壁面约束"：扰动在整个截面内被壁面挤压放大，失稳门槛低，雷诺数两千上下即可转捩。平板边界层只在一侧受限，扰动需要更长距离成长，转捩雷诺数数十万起步。两个数字相隔百倍，但判据逻辑相同——"扰动增长率与耗散速率的竞争"。比较任何两个转捩雷诺数之前，先确认几何可比。

问题：库埃特流在空气里有工程意义吗？

有，而且很日常。硬盘磁头与盘片之间的气膜（厚度几十纳米）就是典型的空气库埃特流，磁头靠这层气膜的压强浮力悬浮在盘面上方——写进去的"飞行高度"就是用修正的库埃特理论算的。气象上的风剖面（贴近地面风速从零爬升）也近似库埃特式的剪切层。别被"平板拖水"的课堂图景限制，剪切流动的舞台比想象中大得多。

问题：这两个精确解在 CFD 教学里为什么地位特殊？

因为它们是“代码考卷”。任何新写的粘性求解器，第一步就用泊肃叶流验证：算出的抛物线剖面、最大与平均速度之比、压降梯度，与解析解逐点对比；第二步用库埃特流验证剪切边界条件。解析解给数值代码提供了无需实验设备的标定基准——精确解的价值一半在物理，一半在验证体系，8.3 节的验证流程直接调用它们。

本节要点

• 库埃特流：拖板驱动、线性剖面、剪应力均匀，是边界层近壁段的局部样板。
• 泊肃叶流：压差驱动、抛物线剖面，最大速度是平均速度的两倍。
• 哈根-泊肃叶定律：压降反比于管径四次方，管路设计的头号约束。
• 适用前提：层流（管流雷诺数约两千以下）、定常、不可压——先判雷诺数再套解。
• 两解叠加给润滑理论铺路，楔形缝隙的承重来自剖面叠加产生的压强场。

精确解到此为止。下面的流动不再那么"乖"，但运动学约束（散度为零）给了我们一套新的组合拳——势流方法。