学习目标 / Learning objectives
- 从速度梯度定义剪切速率，使用 Newtonian 本构关系并区分动力黏度 \(\mu\) 与运动黏度 \(\nu\)。
- 由圆管微元的力平衡推导 Hagen–Poiseuille 速度剖面、流量与流阻。
- 由 Navier–Stokes 尺度分析得到 Reynolds 数，并说明几何相似加无量纲参数匹配才构成动力相似。
- 解释无滑移、层流、湍流、边界层增长与逆压梯分离之间的关系。
- 在 Stokes 阻力和二次阻力模型间作条件判断，并用表面应力理解升力与阻力。
先修知识 / Prerequisites
- 连续性与 Bernoulli 方程。
- 受力平衡和一阶/二阶微分方程。
范围、约定与建模契约 / Scope, conventions, and modeling contract
- 固定惯性系，SI。壁面静止时采用 no-slip：流体切向速度在壁面等于零；若研究微纳尺度稀薄气体，滑移边界需另建模。
- \(\mu\) 是动力黏度，单位 \(\mathrm{Pa\,s}\)；\(\nu=\mu/\rho\) 是运动黏度，单位 \(\mathrm{m^2\,s^{-1}}\)。两者不可互换。
- Newtonian 流体的偏应力与应变率线性；水和空气在普通条件下常近似 Newtonian。血液、聚合物熔体等可能非 Newtonian，本课只作边界说明。
- 管流速度 \(u(r)\) 与截面平均速度 \(\bar u=Q/A\) 要显式区分。
- 升力 \(L\) 与阻力 \(D\) 相对于来流定义；表面压力和黏性剪应力共同贡献。
符号表 / Notation
| 符号 | 定义 | SI 单位 |
|---|---|---|
| \(\mu\) | 动力黏度 | \(\mathrm{Pa\,s}\) |
| \(\nu=\mu/\rho\) | 运动黏度 | \(\mathrm{m^2\,s^{-1}}\) |
| \(\dot\gamma\) | 剪切速率 | \(\mathrm{s^{-1}}\) |
| \(Q\) | 体积流率 | \(\mathrm{m^3\,s^{-1}}\) |
| \(R_h=\Delta p/Q\) | 液压阻力 | \(\mathrm{Pa\,s\,m^{-3}}\) |
| \(\mathrm{Re}=\rho VL/\mu\) | Reynolds 数 | 1 |
| \(C_D,C_L\) | 阻力、升力系数 | 1 |
Knowledge 1 — Viscosity and shear-rate relations / 黏度与剪切速率关系
Knowledge ID: P01-U08-L03-K01
1.1 动量扩散与 Newtonian 本构
两平行板间的简单剪切流 \(\mathbf v=u(y)\mathbf e_x\) 中，剪切速率为
\[\dot\gamma=\frac{du}{dy}.\]
Newtonian 流体满足
\[\boxed{\tau_{xy}=\mu\frac{du}{dy}}.\]
若上板以 \(U\) 匀速运动、板距 \(h\)、无压强梯度，则 \(u(y)=Uy/h\)，壁面对流体的剪应力大小为 \(\mu U/h\)。应力方向取决于所选截面的法向；上式的正负必须与坐标和面力约定一起解释。
图 1：平面 Couette 流的线性速度剖面。箭长按 \(u\propto y\) 绘制。
更一般的不可压缩 Newtonian 流体写成
\[\boldsymbol\sigma=-p\mathbf I+2\mu\mathbf D, \qquad \mathbf D=\frac12\left(\nabla\mathbf v+(\nabla\mathbf v)^{\mathsf T}\right).\]
\(\mathbf D\) 是对称应变率张量；刚体转动部分不应被误计为黏性形变。
1.2 非 Newtonian 边界
若 \(\tau\) 与 \(\dot\gamma\) 非线性或有屈服应力、记忆效应，就不能用单一常数 \(\mu\) 覆盖全部剪切速率。可定义表观黏度，但必须给出测量条件。
图 2：不同本构曲线具有不同教学含义；曲线为示意，未代表某一具体材料数据。
Worked Example 1 — Couette 剪切力与耗散功率
已知： 面积 \(A=0.200\ \mathrm{m^2}\) 的上板以 \(U=0.600\ \mathrm{m\,s^{-1}}\) 滑动，板距 \(h=2.00\ \mathrm{mm}\)，油的 \(\mu=0.450\ \mathrm{Pa\,s}\)。
研究对象： 上板；\(+x\) 为其运动方向。
模型： 稳恒、Newtonian、无滑移、无压强梯度、边缘效应忽略。
\[\tau=\mu\frac{U}{h} =(0.450)\frac{0.600}{2.00\times10^{-3}} =135\ \mathrm{Pa},\]
\[F=\tau A=27.0\ \mathrm N, \qquad P=FU=16.2\ \mathrm W.\]
流体对板的力方向为 \(-x\)，外力需以 \(+27.0\ \mathrm N\) 维持匀速。输入功率最终转化为内能；\(h\) 减半时剪切力加倍，符合梯度尺度检查。
Knowledge 2 — Poiseuille flow and hydraulic resistance / 泊肃叶流与流阻
Knowledge ID: P01-U08-L03-K02
2.1 圆管剖面的推导
考虑半径 \(R\)、长度 \(L\) 的直圆管，轴向 \(z\)，入口压强 \(p_1\)、出口 \(p_2\)，\(\Delta p=p_1-p_2>0\)。对半径 \(r\)、长度 \(L\) 的流体圆柱作轴向力平衡：
\[\Delta p\,\pi r^2+\tau_{rz}(r)\,2\pi rL=0,\]
因此
\[\tau_{rz}(r)=-\frac{\Delta p}{2L}r.\]
Newtonian 关系 \(\tau_{rz}=\mu\,du/dr\) 与壁面无滑移 \(u(R)=0\) 给出
\[u(r)=\frac{\Delta p}{4\mu L}(R^2-r^2).\]
图 3：对同轴流体圆柱作力平衡；剪切力沿侧面分布，方向反抗压强驱动。
积分流量为
\[Q=\int_0^R u(r)\,2\pi r\,dr =\frac{\pi R^4}{8\mu L}\Delta p,\]
\[\bar u=\frac{Q}{\pi R^2}=\frac{\Delta pR^2}{8\mu L}, \qquad u_{\max}=u(0)=2\bar u.\]
图 4：充分发展层流的抛物线剖面，壁面为零、中心最大；箭长按公式绘制。
2.2 液压电阻类比
定义
\[R_h=\frac{\Delta p}{Q}=\frac{8\mu L}{\pi R^4}.\]
同一流量依次通过串联管段时压降相加，\(R_{h,\rm eq}=\sum_iR_{h,i}\)；同一压差加在并联支路时流量相加，\(1/R_{h,\rm eq}=\sum_i1/R_{h,i}\)。这只是在线性层流范围内的结构类比，不代表流体变量与电路微观机制相同。
图 5：Poiseuille 线性关系允许串并联流阻核算；若进入湍流或非 Newtonian 区域，线性电阻模型失效。
Worked Example 2 — 细管流量与壁面剪应力
已知： \(R=0.500\ \mathrm{mm}\)、\(L=0.200\ \mathrm m\)，水黏度 \(\mu=1.00\times10^{-3}\ \mathrm{Pa\,s}\)，压降 \(\Delta p=12.0\ \mathrm{kPa}\)。
研究对象： 充分发展圆管流。
模型： 稳恒、不可压缩、Newtonian、层流、无滑移，入口段忽略。
\[Q=\frac{\pi(5.00\times10^{-4})^4}{8(1.00\times10^{-3})(0.200)}(1.20\times10^4) =1.47\times10^{-6}\ \mathrm{m^3\,s^{-1}}.\]
\[\bar u=\frac{Q}{\pi R^2}=1.88\ \mathrm{m\,s^{-1}}, \qquad \lvert\tau_w\rvert=\frac{\Delta pR}{2L}=15.0\ \mathrm{Pa}.\]
以水的 \(\rho\approx998\ \mathrm{kg\,m^{-3}}\) 检查，\(\mathrm{Re}_D=\rho\bar u(2R)/\mu\approx1870\)，仍可能为层流，但已接近对扰动敏感的范围；因此 Poiseuille 结果需结合入口和扰动条件谨慎使用。
Worked Example 3 — 并联毛细管
已知： 两根同长同液体毛细管并联，半径分别 \(R\) 与 \(2R\)，共同压差 \(\Delta p\)。
模型： 两支路均满足 Poiseuille。
因为 \(Q\propto R^4\)，
\[\frac{Q_2}{Q_1}=\left(\frac{2R}{R}\right)^4=16, \qquad Q_{\rm total}=17Q_1.\]
等效阻力为 \(R_{h,1}/17\)。半径只有两倍，流量却差十六倍，说明测量管径误差会被四次方放大；但在大 Reynolds 数下不能沿用这一比例。
Knowledge 3 — Reynolds number and dynamic similarity / 雷诺数与动力相似
Knowledge ID: P01-U08-L03-K03
以特征速度 \(V\)、长度 \(L\) 对惯性项和黏性项作尺度估计：
\[\rho(\mathbf v\cdot\nabla)\mathbf v\sim\frac{\rho V^2}{L}, \qquad \mu\nabla^2\mathbf v\sim\frac{\mu V}{L^2}.\]
比值为
\[\boxed{\mathrm{Re}=\frac{\rho VL}{\mu}=\frac{VL}{\nu}}.\]
\(\mathrm{Re}\ll1\) 表示黏性扩散相对强，惯性弱；\(\mathrm{Re}\gg1\) 表示惯性在主体流中强，但壁面附近黏性仍通过薄边界层决定无滑移和阻力。
图 6：Reynolds 数组织机制竞争，但不存在适用于所有几何与扰动的单一“临界 Re”。
3.1 动力相似
缩尺实验需先几何相似，再匹配控制方程中相关无量纲数。若重力、自由表面、压缩性或表面张力也重要，还可能需同时匹配 Froude、Mach、Weber 数；通常无法只靠调一个速度全部满足。
图 7：几何缩小后，若流体性质不变，匹配 Re 通常要求速度按长度反比增加；这可能与 Froude 或 Mach 相似冲突。
Worked Example 4 — 风洞 Reynolds 相似
已知： 原型长度 \(L_p=3.00\ \mathrm m\)，空气速度 \(V_p=25.0\ \mathrm{m\,s^{-1}}\)；\(1{:}5\) 模型在相同温度、压强的空气中测试。
目标： 只匹配 Reynolds 数。
模型： 两者 \(\rho,\mu\) 相同，几何相似。
\[\rho V_mL_m/\mu=\rho V_pL_p/\mu \quad\Rightarrow\quad V_m=V_p\frac{L_p}{L_m}=125\ \mathrm{m\,s^{-1}}.\]
这会把模型 Mach 数提高五倍；若接近可压缩范围，同时匹配 Re 与 Mach 可能需要改变气体或压力。答案不是“永远用 125 m/s”，而是暴露了相似性约束之间的冲突。
Knowledge 4 — Laminar flow, turbulence, and boundary layers / 层流、湍流与边界层
Knowledge ID: P01-U08-L03-K04
4.1 层流与湍流
层流并非“完全没有扰动”，而是速度场具有有序层状结构且扰动衰减或受控；湍流含多尺度、非规则但可统计的脉动。圆直管中常以 \(\mathrm{Re}_D\lesssim2300\) 作为层流经验范围、\(\gtrsim4000\) 作为通常湍流范围，但入口、粗糙度和扰动会改变转捩，不能把数字推广到所有流动。
4.2 边界层与分离
高 Re 外流中，主体可近似无黏，但无滑移要求壁面速度为零。速度从零过渡到外流 \(U_e(x)\) 的薄区为边界层。其厚度随下游增长，壁面剪应力与动量亏损在此产生。
图 8：平板边界层从前缘向下游增长；剖面箭长定性满足壁面零速与外缘接近 \(U_e\)。
若沿流向压强升高（逆压梯 \(dp/dx>0\)），近壁低动量流体可能先停滞并反向，流线脱离壁面形成分离和尾迹。此时理想附着流的压强预测与阻力估计可能严重失效。
图 9：逆压梯可使低动量边界层分离；尾迹中的回流箭头是局部平均流向示意。
Knowledge 5 — Drag and lift on bodies / 物体所受阻力与升力
Knowledge ID: P01-U08-L03-K05
物体所受流体力是表面总应力积分：
\[\mathbf F=\oint_S\boldsymbol\sigma\mathbf n\,dA =\underbrace{-\oint_Sp\mathbf n\,dA}_{\text{pressure}} +\underbrace{\oint_S\boldsymbol\tau\mathbf n\,dA}_{\text{viscous shear}}.\]
相对于来流分解为阻力和升力：
\[D=\frac12\rho V^2C_DA_{\rm ref}, \qquad L=\frac12\rho V^2C_LA_{\rm ref}.\]
\(C_D,C_L\) 不是普适常数；它们依赖 Reynolds 数、Mach 数、几何、攻角、粗糙度和非稳恒效应。\(A_{\rm ref}\) 的定义也必须说明。
极低 Reynolds 数、无滑移球的 Stokes 阻力为
\[D=6\pi\mu aV, \qquad \mathrm{Re}_d=\frac{\rho V(2a)}{\mu}\ll1.\]
高 Re 外流常用二次阻力形式；不能把 Stokes 线性式与固定 \(C_D\) 二次式在同一条件下同时使用。
图 10：低 Re 线性阻力与较高 Re 二次阻力是不同渐近模型；交点位置没有普适数值。
图 11：升力与阻力是同一总流体力在来流坐标中的分量，不是两种独立接触力。
Worked Example 5 — Stokes 区小球终端速度
已知： 半径 \(a=40.0\ \mathrm{\mu m}\)、密度 \(\rho_s=2500\ \mathrm{kg\,m^{-3}}\) 的球在 \(\rho_f=1000\ \mathrm{kg\,m^{-3}}\)、\(\mu=0.100\ \mathrm{Pa\,s}\) 的液体中沉降。
研究对象： 球；向下为正。
模型： 稳态终速、Stokes 阻力、无壁面效应。
终速时重量减浮力等于阻力：
\[\frac43\pi a^3(\rho_s-\rho_f)g=6\pi\mu aV_t,\]
\[V_t=\frac{2a^2(\rho_s-\rho_f)g}{9\mu} =5.23\times10^{-5}\ \mathrm{m\,s^{-1}}.\]
检查
\[\mathrm{Re}_d=\frac{\rho_fV_t(2a)}{\mu}=4.19\times10^{-5}\ll1,\]
所以 Stokes 假设自洽。若算出 \(\mathrm{Re}\) 不小，必须迭代使用适当 \(C_D(\mathrm{Re})\)，不能保留该终速式。
课程综合 / Lesson synthesis
实际流动的建模顺序是：
- 先守质量并声明控制体；
- 指定 Newtonian/非 Newtonian 本构与壁面条件；
- 用 \(\mathrm{Re}\) 判断惯性和黏性相对尺度；
- 选择充分发展管流、边界层、Stokes 流或经验系数模型；
- 用压降、功率与极限数检查结果；
- 若存在转捩、分离或多尺度湍流，承认模型需要实验、关联式或数值求解闭合。
公式—决策表 / Formula-and-decision table
| 任务 | 公式 | 使用条件 |
|---|---|---|
| 简单 Newtonian 剪切 | \(\tau_{xy}=\mu\,du/dy\) | 一维剪切、常 \(\mu\) |
| 运动黏度 | \(\nu=\mu/\rho\) | 定义；不要遗漏密度 |
| 圆管 Poiseuille | \(Q=\pi R^4\Delta p/(8\mu L)\) | 稳恒、充分发展、层流、Newtonian、无滑移 |
| 液压阻力 | \(R_h=8\mu L/(\pi R^4)\) | 线性 Poiseuille 范围 |
| Reynolds 数 | \(\mathrm{Re}=\rho VL/\mu\) | 必须声明特征 \(V,L\) |
| Stokes 球阻力 | \(D=6\pi\mu aV\) | 球、无滑移、\(\mathrm{Re}_d\ll1\)、远离壁面 |
| 系数形式 | \(D=\rho V^2C_DA_{\rm ref}/2\) | \(C_D\) 来自适当 Re/Mach/几何数据 |
模型边界 / Model boundaries
- Poiseuille 定律不适用于入口发展段、明显弯管、脉动强流、非圆管而未改几何因子、湍流或非 Newtonian 流体。
- Reynolds 数相同是动力相似的必要条件之一，不总是充分；自由表面、压缩性和毛细作用会引入其他参数。
- 高 Re 不意味着“黏性为零”；壁面剪切、分离和耗散往往仍由黏性触发。
- \(C_D\) 与 \(C_L\) 必须注明参考面积和实验条件；将其当物体固有常数会产生跨尺度错误。
- 湍流平均方程需要额外闭合；本课不把经验涡黏性写成分子黏度。
误区诊所 / Misconception clinic
- “黏度大就是密度大。” 无必然关系；\(\mu\) 描述剪切动量传递，\(\rho\) 描述单位体积质量。
- “管中心剪切最大，因为速度最大。” 错；剪切取决于梯度。Poiseuille 中中心 \(du/dr=0\)，壁面梯度最大。
- “流量与半径平方成正比。” 截面面积如此，但在固定压差 Poiseuille 中 \(Q\propto R^4\)。
- “Re 超过 2300 必然湍流。” 数值是特定圆管条件下经验范围，转捩受扰动与几何影响。
- “高 Re 可把黏性处处删掉。” 边界层与尾迹中黏性仍不可少。
- “阻力总与速度平方成正比。” Stokes 区为线性；过渡区 \(C_D\) 还随 Re 变。
练习与提取 / Practice and retrieval
- K01： 写出动力黏度与运动黏度的关系和各自 SI 单位。
- K01： 两板间距 \(1.00\ \mathrm{mm}\)，上板速度 \(0.200\ \mathrm{m\,s^{-1}}\)，\(\mu=0.0800\ \mathrm{Pa\,s}\)。求剪应力。
- K02： Poiseuille 条件下半径增加 \(10\%\)，固定 \(\Delta p,\mu,L\)，流量增加多少倍？
- K02： 为什么充分发展圆管层流中 \(u_{\max}=2\bar u\)，而不是等于平均值？
- K02： 两相同流阻串联，每个为 \(R_h\)。求等效阻力与每段压降（总压差 \(\Delta p\)）。
- K03： 水流特征速度 \(0.50\ \mathrm{m\,s^{-1}}\)、长度 \(0.020\ \mathrm m\)、\(\nu=1.0\times10^{-6}\ \mathrm{m^2\,s^{-1}}\)，求 Re。
- K03： 同流体中做 \(1{:}10\) 几何缩尺，仅匹配 Re 时模型速度应为原型的多少倍？
- K04： 解释逆压梯为什么容易使边界层分离。
- K04： 高 Re 外流为何仍会有阻力？至少指出两种机制。
- K05： 半径相同的微球在 Stokes 区速度加倍，阻力如何变？若采用固定 \(C_D\) 二次模型又如何变？
- 综合边界： 你用 Poiseuille 公式计算后得到 \(\mathrm{Re}_D=9000\)。应如何处理，而不是直接报告流量？
参考答案 / Answers
1. \(\nu=\mu/\rho\)；\(\mu\) 的单位是 \(\mathrm{Pa\,s}\)，\(\nu\) 的单位是 \(\mathrm{m^2\,s^{-1}}\)。
2. \(\tau=\mu U/h=(0.0800)(0.200)/(1.00\times10^{-3})=16.0\ \mathrm{Pa}\)。
3. \(Q'/Q=(1.10)^4=1.464\)，约增加 \(46.4\%\)。
4. 抛物线剖面 \(u=u_{\max}[1-(r/R)^2]\) 的面积积分给 \(\bar u=Q/(\pi R^2)=u_{\max}/2\)；局部速度不均匀。
5. \(R_{\rm eq}=2R_h\)；流量相同且阻力相同，每段压降为 \(\Delta p/2\)。
6. \(\mathrm{Re}=VL/\nu=(0.50)(0.020)/(1.0\times10^{-6})=1.0\times10^4\)。
7. \(V_m/V_p=L_p/L_m=10\)，前提是 \(\rho,\mu\) 相同且只要求 Re 相似。
8. 下游压强升高对流动施加减速；近壁流体动量最小，可能先降至零并反向，使流线脱离壁面。
9. 黏性壁面剪切产生 skin-friction drag；边界层分离形成压强不对称和 wake，产生 form drag。
10. Stokes 模型中 \(D\propto V\)，变为两倍；固定 \(C_D\) 二次模型中 \(D\propto V^2\)，变为四倍。二者不能在同一 Re 条件下任意切换。
11. 结果与层流前提不自洽。应改用适合湍流/过渡流的压降关联（含粗糙度），迭代求解并报告模型与不确定性。
Exit ticket
- 从圆管流体圆柱受力平衡推导 \(\tau_{rz}\propto-r\)，再说明这为何导致中心速度梯度为零。
- 给一个“高 Re 但黏性仍决定结果”的实例，并指出黏性集中在哪个区域。
参考资料与来源说明 / References and source notes
- OpenStax, University Physics Volume 1, §14.7 Viscosity and Turbulence — Newtonian 黏性、Poiseuille、Reynolds 数与阻力入门。
- MIT OpenCourseWare 8.09, Chapter 6: Fluid Mechanics — 黏性应力、Navier–Stokes 尺度和管流框架。
- NASA Glenn, Boundary Layer — 无滑移、边界层、层流/湍流与分离。
- NASA Glenn, Aerodynamic Forces — 压力与剪应力积分、升阻力分解。
- NASA Glenn, Factors That Affect Aerodynamics — Reynolds/Mach、几何与流体性质对系数的影响。
- NIST Chemistry WebBook, Thermophysical Properties of Fluid Systems — 流体密度、黏度和其他热物性数据入口。
所有图由本文方程原创绘制，未复制来源插图。来源于 2026-08-20 实际核验；NIST 数据库用于说明物性应在指定 \(T,p\) 下取值，本文数值题所给物性视为题设值。
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| OpenStax §14.7 | 简单剪切、Poiseuille、Re、Stokes 阻力 | K01–K05 | 主要为 Newtonian、层流与经验阻力范围 |
| MIT 8.09 Chapter 6 | 张量应力与尺度分析 | K01, K03 | 连续介质、常黏度近似 |
| NASA Boundary Layer | no-slip、边界层、转捩、分离 | K04 | 临界 Re 随几何与扰动变化 |
| NASA Aerodynamic Forces | 压力/剪切总力与 \(L,D\) | K05 | 参考面积、来流坐标需声明 |
| NASA Factors | 动力相似与系数依赖 | K03, K05 | Re 之外可能需 Mach 等参数 |
| NIST Fluid System | \(\rho,\mu\) 的状态依赖 | K01, model boundaries | 物性需与温压条件一致；例题数值均为题设值 |
复习与拓展
- Lesson 1.8.4 将把 \(-\nabla p+\mu\nabla^2\mathbf v\) 放回局部动量守恒，解释这里各特例的共同来源。
- 边界层、分离和湍流闭合将在更高阶流体力学中展开；不要把本课经验范围误作严格相变点。
- 无量纲化思路会在振动、传热与电磁相似性中重复出现。