跳转至

第 12 章 应用流体力学建模（Modelling in Applied Fluid Dynamics）

作者

同 Ch 1-11，Thomas Witelski（Duke University，数学系）和 Mark Bowen（Waseda University，国际理工学中心）合著。本章是 Part III（Case Studies）的第 1 章——也是全书的最后一章。

本章在全书中的定位：综合应用 Ch 8（长波近似）、Ch 9（MMTS）、Ch 10（快/慢系统）、Ch 11（矩方法）等技术解决两个工程流体力学问题——(1) 空气轴承滑块（硬盘读头）和 (2) 楔形几何中的小溪流（泡沫排水）。本章不引入新方法，而是展示"建模者如何把前 11 章的工具组合使用"。

内容概述

章节结构： 1. §12.1 润滑理论：从 Navier-Stokes 方程（12.2-12.3）出发，用 Ch 8 的长波近似 + Ch 4 的无量纲化，导出Reynolds 方程（12.18） $\(\partial_t (\rho h) = \frac{1}{3} \partial_x \left( \rho h^3 \frac{\partial p}{\partial x} \right)\)$ 这描述薄液膜在平板上的扩散。 2. §12.2 空气轴承滑块（12.19-12.29）：应用润滑理论到硬盘读头——读头与旋转磁盘之间有 ~\(1 \mu m\) 的空气间隙。引入轴承数 \(\Lambda = 12 \mu \bar L \bar U / (P_{atm} \bar H^2)\)（12.23），\(\Lambda\) 大意味着对流主导，\(\delta = 1/\Lambda\) 是小参数。压力方程 (12.24) 是奇异扰动（边界层在 \(x_* = 1\) 滑块后缘）。读头的稳定性分析导出延迟微分方程（12.28），由 \(a(t - 2x)\) 的"延迟"项产生——这是稳定性分析的极限。 3. §12.3 楔形几何中的小溪流（12.30-12.45）：在垂直楔形中的黏性液体，分析两种情形： - §12.3.1 毛细吸液（imbibition）：毛细作用 vs 重力 → 相似解 (12.36) \(A = t^{-2/5} f(\eta)\)，\(\eta = x/t^{2/5}\) - §12.3.2 排液（draining）：重力 + 毛细作用 → 出现激波（12.41），由 Rankine-Hugoniot 条件（Ch 2）确定位置

核心思想：本章是Ch 8-11 的工程应用展示——所有方法在润滑理论、空气轴承、楔形流中协同工作。前 11 章是"工具箱"，本章是"用工具箱解决工程问题"。

前置知识：Ch 1-11 全部内容、Ch 8 长波近似、Ch 2 激波与特征线、Ch 5 自相似解、Ch 10 慢流形。

核心方程与概念

1. 润滑理论（§12.1, eq. 12.1-12.18）

Navier-Stokes 方程（12.3）：质量守恒 (12.3a) + 动量守恒 (12.3b, c) 在 2D \((X, Y)\)。

无量纲化（12.4-12.5）：用长宽比 \(\varepsilon = H/L\) 做小参数。关键选择： - 时间尺度 \(T = L/U\)（对流时间） - 垂直速度 \(V = \varepsilon U\)（来自质量守恒 \(\partial_x U + \partial_y V = 0\)）

Stokes 极限（12.6-12.9）：减少 Reynolds 数 \(\text{Re} = \varepsilon^2 \bar \rho \bar U \bar L / \mu \ll 1\) → 惯性项 \(\ll\) 黏性项 → Stokes 方程 \(0 = -\nabla P + \mu \nabla^2 V\)（12.9）。

润滑压力尺度（12.8, 12.10）：\(\bar P = \mu \bar U / (\varepsilon^2 \bar L)\)（黏性 vs 压力梯度平衡）。

Stokes 方程在长波极限（12.10）： $\(\partial_x p_0 = \partial_y^2 u_0, \quad \partial_y p_0 = 0\)$ → \(p_0 = p_0(x, t)\)（与 \(y\) 无关），\(u_0 = \frac{1}{2} \partial_x p_0 \cdot y^2 + C_1 y + C_2\)（抛物线）

Nusselt 速度剖面（12.11-12.13）：无滑移 \(u_0(0) = 0\) → \(C_2 = 0\)；自由面 \(\partial_y u_0 |_{y = h} = 0\) → \(C_1 = -\partial_x p_0 \cdot h\) → $\(u_0(x, y, t) = \frac{1}{2} \partial_x p_0 (y^2 - 2 h y)\)$

Reynolds 方程（12.18）：通过 (12.16) 的质量守恒 + 速度剖面 (12.13) 积分 → $\(\partial_t (\rho h) = \frac{1}{3} \partial_x \left( \rho h^3 \frac{\partial p}{\partial x} \right)\)$

这是润滑理论的核心方程——描述薄液膜在压力梯度作用下的扩散。

2. 空气轴承滑块（§12.2, eq. 12.19-12.29）

几何（12.19, Fig. 12.2）：\(H(X, T) = A(T) + S(X)\)，读头有特定形状 \(S(X)\)，垂直位置 \(A(T)\)（前缘高度）。磁盘在 \(Y = 0\) 以速度 \(\bar U\) 移动。

润滑方程（12.22）： $\(\partial_T (\tilde \rho H) + \frac{\bar U}{2} \partial_X (\tilde \rho H) = \frac{1}{12 \mu} \partial_X \left( \tilde \rho H^3 \frac{\partial P}{\partial X} \right)\)$

注意：与 (12.18) 相比，对流项（来自 Couette 流动 \(\bar U(1 - Y/H)\)）出现——磁盘的移动驱动了对流。

无量纲化（12.23）： $\(\partial_t (\rho h) + \frac{1}{2} \partial_x (\rho h) = \frac{1}{\Lambda} \partial_x \left( \rho h^3 \frac{\partial p}{\partial x} \right)\)$

轴承数（bearing number, 12.23）： $\(\Lambda = \frac{12 \mu \bar L \bar U}{P_{atm} \bar H^2}\)$

这是对流（Couette）vs 扩散（Poiseuille）的比值。小间隙 + 大速度 + 长滑块 → \(\Lambda\) 大 → \(\delta = 1/\Lambda\) 小。

理想气体律（12.24）：\(\rho = k p\)（等温），Reynolds 方程变为 $\(\partial_t (ph) + \frac{1}{2} \partial_x (ph) = \delta \partial_x \left( p h^3 \frac{\partial p}{\partial x} \right)\)$

奇异扰动（12.24-12.27）：\(\delta \to 0\) 极限下，外解 \(p_0 h\) 满足线性平流方程 \(\partial_t (p_0 h) + \frac{1}{2} \partial_x (p_0 h) = 0\)，解为行波 \(p_0(x, t) h(x, t) = f(t - 2x)\)。应用 \(p(0, t) = 1, h(0, t) = a(t)\) → \(f = a\) → $\(p_0(x, t) = \frac{a(t - 2x)}{a(t) + s(x)}\)$

边界层（12.24 末段）：在 \(x_* = 1\)（滑块后缘）有边界层——这是奇异扰动边界值问题（Ch 7 末段预告的"triple deck"在润滑问题中常见）。

稳定性分析（12.28-12.29）：用力平衡——升力 \(F_L = \int_0^1 (p - 1) dx\) = 重力 \(F_S\) + 惯性 \(m \ddot a\)。关键：\(a(t - 2x)\) 是延迟项（\(x\) 从 0 到 1 → \(a\) 的"过去"值）——这是延迟微分方程（DDE, delay differential equation）。

线性化（12.29）：\(a = \bar a + \sigma e^{\lambda t}\) → 特征方程是 \(\lambda\) 的超越方程（含 \(e^{-2\lambda x}\) 的积分）——这比 ODE 的代数特征方程复杂得多。稳定性由这个超越方程的根决定。

工程意义：空气轴承的稳定性分析需要解 DDE 特征问题——这是设计硬盘驱动器时必须验证的核心条件。作者在 [109] 给出了详细分析。

3. 楔形小溪流（§12.3, eq. 12.30-12.45）

几何（12.30, Fig. 12.3）：在 \(X\) 方向上、\(Y, Z\) 横截面的楔形（三角形横截面）中的液体。Plateau 边界（泡沫排水中的"三相接触"）的物理模型。

无量纲化（12.30b）：\(\bar V = \varepsilon^2 \gamma / \mu\) 是由表面张力 \(\gamma\) 主导的速度尺度。

Stokes 极限（12.31）： $\(\partial_x p = \partial_y^2 u + \partial_z^2 u, \quad \partial_y p = 0, \quad \partial_z p = 0\)$ → \(p\) 仅依赖于 \(x\)。

几何关系（12.31 末段 + 12.32 末段）：\(P \propto A^{-1/2}\)（\(A\) 是横截面积）——这是Young-Laplace 方程 \(P = \gamma (\kappa_1 + \kappa_2)\) 的结果，对楔形横截面，\(A\) 越小 → 曲率越大 → 毛细压差越大。

质量守恒（12.32-12.33）：通过横截面积分 + 速度积分 → 演化方程 $\(\partial_t A = \partial_x \left( A^{1/2} \partial_x A \right) \tag{12.33}\)$

这是多孔介质方程（Ch 5.21）的变种！用 \(A^{1/2}\) 代替 \(A^3\)——物理上对应不同的扩散机制**（毛细驱动 vs 重力驱动）。

自相似解（12.36-12.37）：\(A = t^{-2/5} f(\eta), \eta = x/t^{2/5}\)，\(f\) 满足 $\(-\frac{2}{5}[f + \eta f'] = [f^{1/2} f']'\)$ 解为 \(f(\eta) = \frac{1}{100} (\eta_0^2 - \eta^2)^2_+\)（有限支撑）。

这是 Ch 5 自相似解的工程实现——同样的标度分析，应用到不同 PDE。

§12.3.1 毛细吸液（imbibition）（12.38-12.40）：垂直楔形，重力 vs 毛细 → 演化方程 (12.38) 是拟线性对流： $\(\partial_t A - 2A \partial_x A = \partial_x (A^{1/2} \partial_x A)\)$

大时间结构（Fig. 12.4）： - 外解（重力主导，稳态，12.39）：\(A = 4/(x + x_0)^2, x_0 = 4/M\)（\(M\) 是总质量） - 内解（毛细主导，自相似，12.40）：\(A = t^{-2/3} g(\mu), \mu = x/t^{1/3}\)，接触线 \(g \sim (\mu_0 - \mu)^2\) 在 \(\mu \to \mu_0\) - 匹配：\(\mu \to 0+\) → \(g \sim 4/\mu^2\)（与外解 (12.39) 匹配）

§12.3.2 排液（draining）（12.41-12.45）：重力 + 毛细同向 → 演化方程 (12.41) 是 $\(\partial_t A + 2A \partial_x A = \partial_x (A^{1/2} \partial_x A)\)$

激波出现（12.43-12.45）：外解（重力主导）有自相似行波 \(A = t^{-1/2} h(v), v = x/t^{1/2}\)，\(h(v) = \frac{1}{2} v\)（12.44），但有限质量条件要求 \(A\) 在 \(v = v_0 = 2 \sqrt M\) 处突然归零（Rankine-Hugoniot 条件，Ch 2 末段）。

内解（激波附近，12.45）：\(\Phi(\xi) = \frac{v_0}{2} \tanh^2(\frac{1}{2} \sqrt{v_0/2} \xi)\) 是激波层的精确解——与 Burgers 方程的激波结构（Ch 2.6 末段"黏性 Burgers"）形式相同。

这是 Ch 2 激波 + Ch 5 自相似 + Ch 11 矩方法的完美综合。

关键结论

  1. 润滑理论 = "长波 + 低 Reynolds 数"双重简化（12.6-12.18）：从完整的 Navier-Stokes (12.3) 出发，用 \(\varepsilon = H/L \ll 1\)（长波）和 \(\text{Re} = \varepsilon^2 \bar \rho \bar U \bar L / \mu \ll 1\)（低 Reynolds 数）同时得到 Stokes 方程（12.9） + Reynolds 方程（12.18）。双重简化让薄液膜问题从 4D 不可压缩 Navier-Stokes 降为 1D 非线性扩散。
  2. Nusselt 速度剖面（12.13）：\(u_0 = \frac{1}{2} \partial_x p_0 (y^2 - 2 h y)\) 是抛物线（Poiseuille 流）——但与无滑移边界 \(u_0(0) = 0\) 和自由面 \(u_0'(h) = 0\) 同时兼容。这是 Ch 8 长波近似的"线性 + 二次"组合。
  3. 空气轴承的"轴承数"\(\Lambda\) 是 Péclet 数（12.23）：\(\Lambda = 12 \mu \bar L \bar U / (P_{atm} \bar H^2)\) 量纲分析上 = $U L / (\mu/P) = $ 对流/扩散。小间隙 + 大速度 + 长滑块 → \(\Lambda\) 大 → 对流主导 → 边界层在滑块后缘。这种边界层位置与 Ch 7 ODE 例子 (7.23) 末段"激波在右"是同一现象。
  4. 延迟微分方程（DDE）的出现（12.28-12.29）：空气轴承的力平衡是 DDE——\(a(t - 2x)\) 的"延迟"项来自压力信号从滑块前缘传到后缘的时间（\(T_{delay} = 2L/U\)）。这种延迟让稳定性分析需要解超越特征方程（含 \(e^{-2\lambda x}\) 的积分）——比 ODE 复杂得多。这暗示了"代数特征值 → 超越特征值"的推广。
  5. §12.3 小溪流的"非线性扩散方程"是 Ch 5 多孔介质方程的变种（12.33）：\(\partial_t A = \partial_x (A^{1/2} \partial_x A)\) vs Ch 5.21 的 \(\partial_t u = \partial_x (u^3 \partial_x u)\)。物理上：\(A^{1/2}\) 对应毛细驱动（Plateau 边界曲率 ~ \(1/A^{1/2}\)），\(u^3\) 对应Darcy 流动（多孔介质渗透率 ~ \(1/u^3\)）。数学结构相同：都是有限支撑自相似解。
  6. §12.3.2 激波 = "重力 + 毛细"协同效应（12.41-12.45）：当重力和毛细同向驱动液体时，对流项 \(\partial_t A + 2A \partial_x A\) 与 Ch 2 Burgers 方程的拟线性对流同形。特征线相交（Ch 2.6 末段）→ 激波形成 → 激波速度由 Rankine-Hugoniot 条件确定 → 激波层结构由黏性扩散确定（与 Burgers 激波同结构）。
  7. §12.3.1 vs. §12.3.2 的对比（12.38 vs 12.41）：对流项的符号不同 → \(\partial_t A \mp 2A \partial_x A = \partial_x (A^{1/2} \partial_x A)\) → 两种动力学：
  8. 毛细吸液（\(-\) 符号）：特征线发散（膨胀扇），无激波
  9. 排液（\(+\) 符号）：特征线收敛（压缩），有激波

这是 Ch 2 末段 Burgers 方程的"对流方向决定激波"在毛细作用 + 重力系统中的实现。 8. 本章是"全书工具箱"的工程实现：(a) Ch 8 长波 → 润滑方程 (12.18)；(b) Ch 4 无量纲化 → 轴承数 (12.23)；(c) Ch 5 自相似解 → 楔形流 (12.36)；(d) Ch 2 特征线 + 激波 → 排液激波 (12.41-12.45)；(e) Ch 7 边界层 → 滑块后缘 (12.24 末段)；(f) Ch 1 相线分析 → DDE 稳定性 (12.28-12.29)。所有 Ch 1-11 的方法在 Ch 12 中协同工作**。

挑战和开放性问题

  1. 润滑理论的"高阶修正"（§12.1）：作者只给出首阶 Stokes 极限——\(O(\varepsilon^2)\) 修正让 \(u_0\) 包含水平扩散 \(\varepsilon^2 u_{0xx}\)（12.8a 的右边第二项）。这导致 Reynolds 方程 (12.18) 的 \(O(\varepsilon^2)\) 修正——长波展开实际应该是 \(O(\varepsilon^2)\) 修正的偏微分方程而非首阶抛物线。严格处理需要 Ch 8 的"长波展开"系统化推导。
  2. 空气轴承的"惯性效应"（12.6-12.8）：作者在 (12.6) 引入减少 Reynolds 数 \(\text{Re} = \varepsilon^2 \bar \rho \bar U \bar L / \mu\)，严格说这个 \(\text{Re} \ll 1\) 假设是过分严格。实际工程中空气轴承的 \(\text{Re}\) 在 \(0.1\) 到 \(10\) 之间——惯性项不能忽略。更精确的模型需要包含惯性项（\(\text{Re} \cdot \partial_t u\)），这让方程变为非线性。作者没展开——但这是工程上的核心问题。
  3. DDE 稳定性分析的"超越特征方程"（12.28-12.29）：作者只说"the governing equation for \(a(t)\)... can be related to a delay-differential equation and yields a more challenging transcendental equation for the linear stability analysis"。严格解需要DDE 数值方法（dde23 in MATLAB, BIFTOOL）——这是工程分析的实际工具。作者引用 [109] ——详细解在该文献。
  4. §12.3 的"Plateau 边界近似"（§12.3 注释 4, 12.31 末段）：作者明确说"the sides of a Plateau border are actually circular arcs and there is no free-surface"。真实的 Plateau 边界有曲率，不是简单的楔形。严格的 Plateau 边界分析需要更复杂的几何 + 表面张力平衡（Young-Laplace 方程 \(P_{in} - P_{out} = \gamma (\kappa_1 + \kappa_2)\)）。作者用了"楔形"近似——这是 Plateau 边界的线性化几何。
  5. §12.3.1 vs 12.3.2 的"接触线奇异性"（12.40 末段）：作者给出内解 \(g \sim (\mu_0 - \mu)^2\) 在 \(\mu \to \mu_0\)（接触线附近）——这是零接触角的"自然"边界条件。实际流体有有限接触角（与固体表面的浸润性）——接触线是流动问题的经典困难（"moving contact line problem"），需要滑移模型或扩散界面。作者没展开。
  6. §12.3.2 排液激波的"内解精确解"（12.45）：\(\Phi(\xi) = \frac{v_0}{2} \tanh^2(\frac{1}{2} \sqrt{v_0/2} \xi)\) 是精确的——因为 Burgers 方程的激波结构是 \(\tanh\)。对其他非线性扩散（如 \(\partial_x(A^n \partial_x A)\)）激波结构不是 \(\tanh\)——需要数值积分。这是 Burgers 方程的特殊性。
  7. 缺失主题：弹性润滑（EHL）。真实的轴承（如汽车发动机、关节假体）涉及弹性基底（金属在压力下变形）。这导致 Reynolds 方程需要耦合弹性方程——弹性流体动力学润滑（EHL）。这超出 SUMS 范围——是工程流体力学的高级主题。
  8. 润滑理论的"现代应用"（§12.1 导言）：作者列举"air bearing slider, oil lubricating parts of engine"——但现代润滑理论应用于软物质（泪液膜、粘液运输、油墨打印）、生物流变学（关节液、血液流变）、微流体（MEMS 器件）。这些应用的物理与经典轴承不同（非牛顿流体、表面张力主导、可压缩气体）。作者的"基础润滑"是入门——实际工程需要扩展。

个人反思与批判性分析

本章是全书的"应用展示"——综合 Ch 1-11 工具解决两个工程问题。这种"工具箱 + 应用"的结构是科学方法论的标准模式：先学方法，再用方法解决具体问题。几个值得反思的点：

  1. 润滑理论 = "Ch 8 + Ch 4 + Ch 3" 的工程化。§12.1 几乎只用 Ch 8 长波近似和Ch 4 无量纲化就完成了——这暗示Ch 8 的"长波近似"是所有薄液膜问题的核心技术。学生学完 Ch 8 应该能独立推导润滑理论——这是一个能力测试。
  2. 空气轴承 = "Ch 2 + Ch 7 + Ch 8 + Ch 9 + Ch 11" 的综合。§12.2 涉及：
  3. Ch 8 长波 → 润滑方程 (12.22)
  4. Ch 7 边界层 → 滑块后缘 (12.24 末段)
  5. Ch 2 行波解 → \(p_0(x, t) = f(t - 2x)\)
  6. Ch 11 矩方法 → 升力 \(F_L = \int_0^1 (p - 1) dx\)
  7. DDE 稳定性 → 超越特征方程 (12.29)

这是 Ch 1-11 的"工程集成"——一个工程问题往往需要多个数学方法协同。

  1. §12.3 的"激波 + 自相似 + 有限支撑"三重奏。§12.3.2 排液有：
  2. Ch 2 激波（12.41, 12.45）
  3. Ch 5 自相似（12.43）
  4. Ch 5 有限支撑（12.37）
  5. Ch 7 激波层匹配（12.45 是内解）

这种"所有方法在同一个问题中"是数学建模的极致——把一个看似简单的"小溪流"问题展现为多尺度、多方法的交响乐。

  1. 延迟微分方程（DDE）的出现（12.28-12.29）是"章末彩蛋"——前 11 章的所有方法都是 ODE/PDE 框架，DDE 出现在最后一章。这是对读者的提示："前 11 章的方法不能解决所有问题——有些问题需要超越 ODE/PDE 的工具"。DDE 在工程中无处不在（控制系统、通信网络、生物振荡）——作者把这一深刻洞察放在最后一章，是对全书结构的隐喻。

  2. 缺失主题：3D 润滑。真实润滑问题通常是 3D——球轴承、关节、齿轮。2D 长波近似（12.1）是入门。3D 完整问题需要椭圆型 Reynolds 方程——这是计算流体动力学（CFD）研究的核心。作者没展开——但 12.1 的 2D 推导是基础。

  3. "压力场 + 几何形状"的双向耦合。§12.2 空气轴承：压力场 \(p(x, t)\) 决定升力 → 决定读头位置 \(a(t)\) → 决定几何 \(H = A + S\) → 又决定压力场。这是一个"双向耦合"问题。§12.3 楔形流：压力场 \(P \propto A^{-1/2}\)（Young-Laplace）决定速度场 → 决定 \(A(x, t)\) 演化 → 又决定 \(P\) 场。这是"自由边界问题"——边界本身是变量。这种双向耦合是流体力学的核心困难——也是现代数值方法（如 ALE, Arbitrary Lagrangian-Eulerian）需要解决的核心问题。

  4. "工程常数"vs"普适结构"。8/945（Taylor 弥散）、\(\frac{1}{2}\)（压缩波速度，12.44）、\(\frac{4}{M}\)（吸液常数，12.39）——这些数字看起来"工程化"但它们有"普适结构"。例如：\(8/945\) 是 Poiseuille 流的特征（不同流动几何有不同的"magic number"）；\(\frac{1}{2}\) 是 Burgers 方程的特征（特征线速度）；\(\frac{4}{M}\) 是质量守恒 + 几何的结果。这些常数告诉读者"每个具体问题有自己的"特征数""，但结构是普适的。

  5. 从润滑理论到"软物质"。§12.1 的 Reynolds 方程适用于牛顿流体。现代润滑（关节液、泪液膜、粘液）涉及非牛顿流体（黏弹性、剪切稀化、屈服应力）。非牛顿修正让 Reynolds 方程变为 $\(\partial_t h = \frac{1}{\mu(h, \dot\gamma)} \partial_x \left( h^3 \partial_x p \right)\)$ 其中 \(\mu\) 依赖于剪切率 \(\dot \gamma\)。这是"流变学"的核心——但SUMS 不涉及。

重要参考文献

  • [X1] G.K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, 1967. (流体力学经典)
  • [X2] L.G. Leal, Advanced Transport Phenomena, Cambridge University Press, 2007. (高等输运现象)
  • [X3] O. Reynolds, "On the Theory of Lubrication and Its Application to Mr. Beauchamp Tower's Experiments, Including an Experimental Determination of the Viscosity of Olive Oil," Philosophical Transactions of the Royal Society of London 177 (1886) 157–234. (Reynolds 润滑理论原始论文)
  • [X4] G.I. Taylor, "Dispersion of Soluble Matter in Solvent Flowing Slowly through a Tube," Proceedings of the Royal Society A 219 (1953) 186–203. (Taylor 弥散原始论文)
  • [X5] A.E. Hughes, Lubrication, Cambridge University Press, 1990. (润滑工程)
  • [X6] Bernard J. Hamrock, Steven R. Schmid, and Bo O. Jacobson, Fundamentals of Fluid Film Lubrication, Marcel Dekker, 2004. (流体膜润滑基础)
  • [X7] Daniel T. Valentine, Fundamentals of Incompressible Flow, World Scientific, 2017. (不可压缩流体力学)
  • [X8] Stephen Childress, An Introduction to Theoretical Fluid Mechanics, Courant Lecture Notes, 2009. (理论流体力学)
  • [X9] Avner Friedman, Free Boundary Problems in Fluid Dynamics, Springer, 2018. (流体自由边界问题)
  • [X10] L.M. Milne-Thomson, Theoretical Hydrodynamics, Dover, 1996. (理论流体力学经典)