学习目标 / Learning objectives
- 区分流线、迹线、流管与稳恒流，并从速度场求流线或粒子轨迹。
- 从控制体质量账推导积分与微分连续性方程，说明 \(A_1v_1=A_2v_2\) 的额外假设。
- 从无黏 Euler 方程沿流线推导 Bernoulli 常数，解释各项的单位和能量意义。
- 建模 Torricelli 出流、Venturi/Pitot 测速与压强积分产生的升力。
- 用黏性、非稳恒、可压缩、旋转性、泵功与空化等条件审核理想流模型。
先修知识 / Prerequisites
范围、约定与建模契约 / Scope, conventions, and modeling contract
- 采用固定地面惯性系，位置 \(\mathbf x\)，速度场 \(\mathbf v(\mathbf x,t)\)；粒子标签为 \(\mathbf a\)，轨迹为 \(\mathbf X(\mathbf a,t)\)。
- 控制面外法向为 \(\mathbf n\)，\(\mathbf v\cdot\mathbf n>0\) 表示流出。
- \(Q=\int_A\mathbf v\cdot\mathbf n\,dA\) 是体积流率；\(\dot m=\int_A\rho\mathbf v\cdot\mathbf n\,dA\) 是有符号质量流率。只有截面速度近似均匀时才写 \(Q\approx Av\)。
- “理想流体”在本课指无黏；许多例题还另加不可压缩、稳恒、一维截面平均等假设。它们不是同义词。
- 统一取 \(z\) 向上，\(\mathbf g=-g\mathbf e_z\)；SI 单位制。
符号表 / Notation
| 符号 | 定义 | SI 单位 |
|---|---|---|
| \(\mathbf v(\mathbf x,t)\) | Eulerian 速度场 | \(\mathrm{m\,s^{-1}}\) |
| \(\mathbf X(\mathbf a,t)\) | 带标签流体粒子的轨迹 | m |
| \(Q\) | 体积流率 | \(\mathrm{m^3\,s^{-1}}\) |
| \(\dot m\) | 质量流率 | \(\mathrm{kg\,s^{-1}}\) |
| \(H=p/(\rho g)+v^2/(2g)+z\) | Bernoulli 水头 | m |
| \(p_0\) | 停滞压强 | Pa |
| \(M=v/c\) | Mach 数 | 1 |
Knowledge 1 — Streamlines, pathlines, and steady flow / 流线、迹线与稳恒流
Knowledge ID: P01-U08-L02-K01
1.1 两类曲线回答不同问题
在固定时刻 \(t_0\)，流线处处与瞬时速度相切。参数 \(s\) 下
\[\frac{d\mathbf r}{ds}\parallel\mathbf v(\mathbf r,t_0),\]
二维时可写
\[\frac{dy}{dx}=\frac{v_y(x,y,t_0)}{v_x(x,y,t_0)}.\]
迹线是同一物质粒子的时序位置：
\[\frac{d\mathbf X}{dt}=\mathbf v(\mathbf X,t), \qquad \mathbf X(\mathbf a,t_0)=\mathbf a.\]
若 \(\partial\mathbf v/\partial t=\mathbf0\)，流动称为稳恒（steady），此时连续释放的染料线、流线和迹线重合；非稳恒流中它们通常不同。
图 1：左图冻结整个速度场，右图追踪一个粒子；非稳恒时不能把瞬时流线当作粒子未来轨迹。
1.2 流管
由一束流线围成的侧壁称为流管。因为速度切于侧壁，理想数学流管的侧壁没有质量穿越；质量只经端面进出。
图 2：流管变窄并不自动推出速度增大；还需质量守恒、稳恒与密度关系。
Knowledge 2 — Mass conservation and continuity equation / 质量守恒与连续性方程
Knowledge ID: P01-U08-L02-K02
2.1 控制体质量账
对固定控制体 \(V\)，内部质量变化率加净流出率为零：
\[\frac{d}{dt}\int_V\rho\,dV +\oint_{\partial V}\rho\mathbf v\cdot\mathbf n\,dA=0.\]
对任意固定小体积用散度定理，可得局部连续性方程
\[\boxed{\frac{\partial\rho}{\partial t}+\nabla\cdot(\rho\mathbf v)=0}.\]
若流体质点密度保持不变，即 \(D\rho/Dt=0\)，则 \(\nabla\cdot\mathbf v=0\)。这称为不可压缩流动；它比“某时刻密度场均匀”更严格。
图 3：外法向约定让流入通量为负、流出为正；稳恒只令储存项为零，并不令每个端面的通量为零。
2.2 一维截面平均形式
稳恒、单入口单出口时
\[\dot m=\rho_1A_1\bar v_1=\rho_2A_2\bar v_2.\]
再加不可压缩条件 \(\rho_1=\rho_2\) 才有
\[Q=A_1\bar v_1=A_2\bar v_2.\]
实际管流的速度剖面并不均匀，\(\bar v=Q/A\) 是面积平均速度；不要把中心线最大速度直接代入 \(Av\)。
图 4：在不可压缩稳恒一维模型中，较小截面对应较大平均速度；箭长按定性比例绘制。
Worked Example 1 — 分支管的质量守恒
已知： 水以 \(Q_1=6.00\times10^{-3}\ \mathrm{m^3\,s^{-1}}\) 进入三通，一支流出管直径 \(d_2=30.0\ \mathrm{mm}\)、平均速度 \(v_2=4.00\ \mathrm{m\,s^{-1}}\)；另一支直径 \(d_3=40.0\ \mathrm{mm}\)。
研究对象： 包住三通的固定控制体；流入取负、流出取正。
模型： 水不可压缩、稳恒、截面平均一维。
\[Q_2=A_2v_2=\frac{\pi(0.0300)^2}{4}(4.00) =2.83\times10^{-3}\ \mathrm{m^3\,s^{-1}},\]
\[Q_3=Q_1-Q_2=3.17\times10^{-3}\ \mathrm{m^3\,s^{-1}},\]
\[v_3=\frac{Q_3}{\pi d_3^2/4}=2.52\ \mathrm{m\,s^{-1}}.\]
三位有效数字。\(Q_2+Q_3=Q_1\)，质量账闭合；若算得 \(Q_3<0\)，就表示实际流向与假定相反。
Knowledge 3 — Bernoulli equation and energy interpretation / 伯努利方程及能量解释
Knowledge ID: P01-U08-L02-K03
3.1 从 Euler 方程沿流线积分
稳恒、无黏、密度恒定且只有重力体力时，Euler 方程为
\[\rho(\mathbf v\cdot\nabla)\mathbf v=-\nabla p+\rho\mathbf g.\]
沿流线取弧长 \(s\)，利用 \(\mathbf v\cdot d\mathbf r=v\,ds\) 与切向加速度 \(v\,dv/ds\)：
\[v\,dv+\frac{dp}{\rho}+g\,dz=0.\]
沿同一流线从 1 到 2 积分：
\[\boxed{p+\frac12\rho v^2+\rho gz=C_{\rm streamline}},\]
或以单位重量能量写成水头
\[\frac{p}{\rho g}+\frac{v^2}{2g}+z=H.\]
三项分别是压强能、动能与重力势能的单位体积形式，单位都是 Pa；不是三种“压强的定义”。若流动无旋，Bernoulli 常数可在整个连通流场相同；一般旋转流中只保证沿同一流线相同。
图 5：Bernoulli 方程比较同一流线上的两个状态点；曲线本身不是管壁。
图 6：同高、同一 Bernoulli 流线中速度增大伴随静压下降；这不是任意流动都满足的因果口号。
Worked Example 2 — 向上喷嘴中的压强变化
已知： 水的密度取 \(\rho=1000\ \mathrm{kg\,m^{-3}}\)；下截面 1 的速度 \(v_1=2.00\ \mathrm{m\,s^{-1}}\)、表压 \(p_{1g}=180\ \mathrm{kPa}\)；上截面 2 高 \(5.00\ \mathrm m\)，管径为下截面的 \(0.600\)，忽略损失。
研究对象： 连接两截面的稳恒流管；\(z\) 向上。
模型： 不可压缩、无黏、一维截面平均、同一流线族。
连续性给出
\[v_2=v_1\frac{A_1}{A_2} =v_1\frac{d_1^2}{d_2^2} =2.00\left(\frac{1}{0.600}\right)^2 =5.56\ \mathrm{m\,s^{-1}}.\]
以表压使用 Bernoulli（两端大气基准相同）：
\[p_{2g}=p_{1g}+\frac12\rho(v_1^2-v_2^2)-\rho g(z_2-z_1) =118\ \mathrm{kPa}.\]
压强因加速和升高都下降。若计算得到低于液体蒸气压，理想连续液柱模型必须因空化而停止使用。
Knowledge 4 — Torricelli, Venturi, and lift applications / 托里拆利、文丘里与升力应用
Knowledge ID: P01-U08-L02-K04
4.1 Torricelli 出流
大水箱自由面与小孔均接同一大气，自由面速度远小于孔口速度，二者高差为 \(h\)。Bernoulli 给出
\[v_{\rm ideal}=\sqrt{2gh}.\]
真实孔口因收缩与黏性损失常写 \(Q=C_dA\sqrt{2gh}\)，\(C_d<1\) 由实验给出；不能把 \(C_d\) 假装由理想模型推导。
图 7：理想 Torricelli 速度来自自由面到孔口的压强、动能和高度交换；孔口射流随后是抛体运动。
4.2 Venturi 与 Pitot
水平 Venturi 管中，连续性与 Bernoulli 联立：
\[Q=A_1A_2\sqrt{\frac{2(p_1-p_2)}{\rho(A_1^2-A_2^2)}}.\]
Pitot 探头的迎流停滞点速度降为零，同高时
\[p_0-p=\frac12\rho v^2, \qquad v=\sqrt{\frac{2(p_0-p)}{\rho}}.\]
这些公式均需要相应的理想条件；工程仪表以校准系数修正探头、边界层和可压缩效应。
图 8：水平收缩段中平均速度升高、静压降低；测压柱高度体现静压而不是“总能量损失”。
图 9：简化 Pitot–static 管以迎风总压孔连到中央停滞通道，并以与外壁齐平的侧孔连到彼此隔离的静压通道；停滞腔速度近零，但并非整个流场停住。
4.3 升力不是“等时到达”
对物体表面，理想流动的压力贡献为
\[\mathbf F_p=-\oint_S p\mathbf n\,dA.\]
升力和阻力只是该合力相对于来流方向的分量。翼型周围的速度与压强场由边界条件、环流和起动过程共同决定；上下分开的两个流体粒子没有必须同时到达尾缘的定律。实际升力还含黏性剪力，并受分离和失速控制。
图 10：升力来自完整表面应力积分；红箭定性表示压力，不宣称按数值比例。
Worked Example 3 — 大水箱孔口射程
已知： 小孔位于自由面下 \(h=1.80\ \mathrm m\)，孔口离地 \(y=1.25\ \mathrm m\)；自由面和射流均接大气。
参考系： 地面系，出流初速度水平。
模型： 大水箱、理想无黏出流，空气阻力忽略。
\[v=\sqrt{2gh}=5.94\ \mathrm{m\,s^{-1}}, \qquad t=\sqrt{\frac{2y}{g}}=0.505\ \mathrm s,\]
\[x=vt=3.00\ \mathrm m.\]
有趣的是 \(x=2\sqrt{hy}\)，显式的 \(g\) 消去；但真实 \(C_d<1\) 会缩短射程。
Worked Example 4 — Venturi 流量
已知： 水平管 \(d_1=60.0\ \mathrm{mm}\)、喉部 \(d_2=30.0\ \mathrm{mm}\)，测得 \(p_1-p_2=18.0\ \mathrm{kPa}\)；水密度 \(998\ \mathrm{kg\,m^{-3}}\)。
研究对象： 两测压截面间流管。
模型： 稳恒、不可压缩、无黏，截面平均速度代表一维流动。
面积比 \(A_1/A_2=4\)，故 \(v_2=4v_1\)。Bernoulli 给出
\[p_1-p_2=\frac12\rho(v_2^2-v_1^2) =\frac{15}{2}\rho v_1^2,\]
\[v_1=1.55\ \mathrm{m\,s^{-1}}, \qquad Q=A_1v_1=4.38\times10^{-3}\ \mathrm{m^3\,s^{-1}}.\]
流量为 \(4.38\ \mathrm{L\,s^{-1}}\)。实际装置需乘流量系数，理想结果通常是相同压差下的偏高估计。
Worked Example 5 — 低速 Pitot 测速
已知： 空气密度 \(1.18\ \mathrm{kg\,m^{-3}}\)，停滞压强比静压高 \(420\ \mathrm{Pa}\)。
模型： 低 Mach、同高、无黏近似，测点不在强边界层中。
\[v=\sqrt{\frac{2(420)}{1.18}}=26.7\ \mathrm{m\,s^{-1}}.\]
若声速约 \(343\ \mathrm{m\,s^{-1}}\)，\(M\approx0.078\)，不可压缩近似合理。若 \(M\) 接近 \(0.3\) 或更高，应使用可压缩停滞关系而非本式。
Knowledge 5 — Limits of ideal-flow models / 理想流动模型的边界
Knowledge ID: P01-U08-L02-K05
5.1 使用 Bernoulli 前的审核
经典形式 \(p+\rho v^2/2+\rho gz=C\) 需要逐项检查：
- 稳恒吗？ 非稳恒势流会多出 \(\partial\phi/\partial t\) 项。
- 无黏近似可靠吗？ 边界层、长管损失、尾迹和分离区通常不可靠。
- 密度可视为常量吗？ 气体在较大压差或 \(M\gtrsim0.3\) 时需检查可压缩性。
- 比较的是同一流线吗？ 有旋流中不同流线可有不同常数。
- 有泵、涡轮或热交换吗？ 需加入轴功、损失或热力学焓项。
- 压强是否低于蒸气压？ 若是，可能空化，单相模型失效。
连续性方程本身是质量守恒，不因黏性存在而失效；失效的是把它简化成 \(Av=\mathrm{constant}\) 的那组假设。
图 11：公式选择流程。先守质量，再判断能量模型；不能因为题目出现“流体”就自动套 Bernoulli。
课程综合 / Lesson synthesis
理想流动的核心不是三条孤立公式，而是三个层次：
- 几何层： \(\mathbf v(\mathbf x,t)\) 定义流线，粒子积分定义迹线。
- 守恒层： 连续性核算质量；Bernoulli 在额外条件下核算沿流线的机械能。
- 诊断层： 压差、流量和速度可相互推断，但必须用边界条件和模型边界闭合问题。
公式—决策表 / Formula-and-decision table
| 目的 | 关系 | 成立条件 |
|---|---|---|
| 求流线 | \(dy/dx=v_y/v_x\) | 固定同一时刻的二维速度场 |
| 固定控制体质量账 | \(d\int_V\rho dV/dt+\oint\rho\mathbf v\cdot\mathbf n\,dA=0\) | 普适连续介质质量守恒 |
| 局部质量守恒 | \(\partial_t\rho+\nabla\cdot(\rho\mathbf v)=0\) | 场足够光滑；激波处用弱形式 |
| 一维不可压缩流量 | \(A_1\bar v_1=A_2\bar v_2\) | 稳恒、单流管、密度相同 |
| 经典 Bernoulli | \(p+\rho v^2/2+\rho gz=C\) | 稳恒、无黏、不可压缩、沿流线、无轴功 |
| Torricelli | \(v=\sqrt{2gh}\) | 大水箱、小孔、同外压、无损失 |
| Pitot | \(p_0-p=\rho v^2/2\) | 低速不可压缩、同高、同一流线 |
模型边界 / Model boundaries
- 无黏外流可在边界层外很好地近似压强场，却不能预测壁面摩擦或失速。
- 管长很长时，即使速度不大，累积黏性损失也可主导压降。
- 流线弯曲时存在法向压强梯度；只用“快则低压”可能漏掉曲率动力学。
- 对气体，\(M<0.3\) 只是常用经验筛查，不是严格定理；还应检查密度变化。
- 流经激波时总机械能形式需用可压缩热力学守恒，经典不可压缩 Bernoulli 不跨激波成立。
误区诊所 / Misconception clinic
- “窄处速度一定大。” 只有在指定质量流率、稳恒且密度关系确定时成立。
- “Bernoulli 说速度大导致压强小。” 它是沿特定条件下流线的能量关系；共同原因是整体动力学与边界条件。
- “流线就是粒子轨迹。” 只在稳恒流中重合。
- “升力来自上下空气必须同时到尾缘。” 不存在等时到达定律；应从速度场、压力/剪应力积分与环流边界条件分析。
- “连续性在黏性流中不成立。” 错；质量守恒仍成立，变的是速度剖面和能量损失。
- “停滞压强就是容器内静水压强。” 停滞压强是流体沿适当减速过程达到零速度时的压强。
练习与提取 / Practice and retrieval
- K01： 给出稳恒流的数学定义，并说明它是否要求每个流体粒子静止。
- K01： 速度场 \(\mathbf v=(ax)\mathbf e_x-(ay)\mathbf e_y\)。求通过点 \((x_0,y_0)\) 的流线方程。
- K02： 直径从 \(8.00\ \mathrm{cm}\) 缩至 \(4.00\ \mathrm{cm}\) 的水管中，上游平均速度 \(1.20\ \mathrm{m\,s^{-1}}\)。求下游平均速度。
- K02： 判断场 \(\mathbf v=(2x)\mathbf e_x-(2y)\mathbf e_y\) 是否满足二维不可压缩连续性。
- K03： 水平理想流中，速度从 \(2.0\) 增至 \(6.0\ \mathrm{m\,s^{-1}}\)。求静压变化 \(p_2-p_1\)，取 \(\rho=1000\ \mathrm{kg\,m^{-3}}\)。
- K03： 写出 Bernoulli 常数能从“一条流线一个值”推广到“全流场同一值”的附加条件。
- K04： 水箱孔口位于水面下 \(0.800\ \mathrm m\)，求理想出流速度。
- K04： 空气密度 \(1.20\ \mathrm{kg\,m^{-3}}\)，Pitot 压差 \(150\ \mathrm{Pa}\)，求低速近似速度。
- K04/K05： 为什么不能仅凭翼面上方“路程更长”就断言速度更大？
- K05： 对长粗糙管、正在启动的泵送水流、低速开阔水槽、跨声速喷管四种情景，分别指出经典 Bernoulli 的主要风险。
- 综合： 某水平收缩管测得喉部压强接近水的蒸气压。继续减小喉径时，为什么理想模型给出的流量趋势可能失效？
参考答案 / Answers
1. \(\partial\mathbf v/\partial t=\mathbf0\)（在所选参考系内）；粒子仍可运动，只是每个固定位置的场量不随时间变。
2. \(dy/dx=(-ay)/(ax)=-y/x\)，积分得 \(xy=x_0y_0\)。
3. 面积缩为 \(1/4\)，所以 \(v_2=4v_1=4.80\ \mathrm{m\,s^{-1}}\)。
4. \(\nabla\cdot\mathbf v=\partial_x(2x)+\partial_y(-2y)=0\)，满足不可压缩局部条件。
5. \[p_2-p_1=\frac12\rho(v_1^2-v_2^2)=-1.60\times10^4\ \mathrm{Pa}=-16.0\ \mathrm{kPa}.\]
6. 流动还需无旋，即 \(\nabla\times\mathbf v=\mathbf0\)，并位于同一连通区域。
7. \(v=\sqrt{2gh}=3.96\ \mathrm{m\,s^{-1}}\)。
8. \(v=\sqrt{2\Delta p/\rho}=15.8\ \mathrm{m\,s^{-1}}\)。
9. 路程差本身不给时间或速度；上下粒子没有同时到达尾缘的约束。速度场由控制方程、环流与边界条件决定。
10. 长粗糙管：黏性水头损失；启动流：非稳恒项；低速开阔槽：若缓变且损失小可局部使用，但自由面/非均匀速度需谨慎；跨声速喷管：强可压缩性与可能的激波。
11. 低压可触发空化，出现两相、蒸气泡与额外损失；固定密度单相 Bernoulli 和几何流量关系不再闭合。
Exit ticket
- 写出固定控制体质量守恒，并说明从它得到 \(A_1v_1=A_2v_2\) 必须增加哪三个假设。
- 给定两个流动截面，先不用公式作答：列出你在决定能否使用经典 Bernoulli 前会检查的五项条件。
参考资料与来源说明 / References and source notes
- OpenStax, University Physics Volume 1, §14.5 Fluid Dynamics — 稳恒流、流线、体积流率与连续性。
- OpenStax, University Physics Volume 1, §14.6 Bernoulli's Equation — Bernoulli 能量解释、Torricelli 与应用。
- MIT OpenCourseWare 8.09, Chapter 6: Fluid Mechanics — 控制体连续性、Euler 描述与流体动力学方程。
- NASA Glenn, Bernoulli's Equation — Bernoulli 推导假设与航空应用边界。
- NASA Glenn, Pitot-Static Tube Speedometer — 静压、总压与低速测速。
- NASA Glenn, Bernoulli and Newton — 升力解释及“等时到达”错误。
- NASA Glenn, Aerodynamic Forces — 表面压力与剪应力积分、升阻力分解。
所有图依据方程原创，未复制来源插图。资料于 2026-08-20 实际核验；NASA 页面使用英美工程单位的示例时，本文全部转换为 SI，并采用 \(z\) 向上的统一符号。
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| OpenStax §14.5 | 流线、稳恒流、\(Q=Av\) | K01, K02 | 本文明确截面平均与不可压缩条件 |
| OpenStax §14.6 | Bernoulli、Torricelli、Venturi | K03, K04 | 稳恒无黏不可压缩 |
| MIT 8.09 Chapter 6 | 积分/局部连续性与 Euler 框架 | K02, K03 | 连续介质与光滑场；激波处需弱形式 |
| NASA Bernoulli | 方程假设检查 | K03, K05 | 沿流线、无轴功与低损失 |
| NASA Pitot-static | 停滞压强测速 | K04 | 低速不可压缩近似 |
| NASA Bernoulli and Newton; Aerodynamic Forces | 升力与表面应力积分 | K04, K05 | 拒绝等时到达解释；实际流含黏性 |
复习与拓展
- Lesson 1.8.3 将把壁面无滑移、剪切应力和耗散加入本课的质量账。
- Lesson 1.8.4 将用物质导数与 Navier–Stokes 把 Bernoulli 作为受限积分结果重新定位。
- 分析力学中“沿路径积分得到常数”的思想，会与保守力和势函数再次联系。