跳转至

第 8 章 PDE 问题的长波渐近（Long-Wave Asymptotics for PDE Problems）

作者

同 Ch 1-7，Thomas Witelski（Duke University，数学系）和 Mark Bowen（Waseda University，国际理工学中心）合著。本章是 Part II（Solution Techniques）的第 4 章——Part II 的最后一章。

本章在全书中的定位：Ch 7 的边界层理论从 ODE 扩展到 PDE——在"细长矩形"（slender rectangle）上用长波近似（\(\varepsilon = H/L \ll 1\)）把 2D Laplace 方程 \(\varepsilon^2 u_{xx} + u_{yy} = 0\) 约化为 1D ODE \(u_{yy} = 0\)，再用边界层方法处理两个端点。这是润滑理论（lubrication theory）的数学基础，也是 Ch 12 末尾"流体动力学应用"的核心。

内容概述

本章回答一个核心问题：在细长（长宽比 \(\varepsilon \ll 1\)）的二维区域上，Laplace 方程如何用长波近似 + 边界层系统化求解？

章节结构： 1. §8.1 经典分离变量解：复习 Laplace 方程的 Fourier 级数解（特征函数 \(\sin(n\pi X/L)\), \(\sinh(n\pi Y/L)\)） 2. §8.2 细长矩形上的 Dirichlet 问题：\(\varepsilon^2 u_{xx} + u_{yy} = 0\) 的"长波近似"（outer solution）+ 端点边界层（inner solution） 3. §8.3 绝缘导线问题：Dirichlet 边界条件换为 Neumann（无通量），引入可解性条件（solvability condition）——边界层反过来决定外解 4. §8.4 非均匀绝缘导线：边界 \(Y = H f(X/L)\) 是非常数（domain perturbation），导出更复杂的可解性方程 (8.37) 5. §8.5 进一步方向：润滑理论、slender body asymptotics、约化维数模型 6. §8.6 习题：包括KdV 方程的推导（8.5）和多孔介质方程的推导（8.6）——Ch 5 主题的 PDE 来源

核心思想：在细长区域上，长波方向（\(X\) 沿细长方向）的二阶导数被 \(\varepsilon^2\) 抑制，使 2D Laplace 方程约化为 1D ODE——这是降维（dimensional reduction）的典型例子。端点边界层修正"假定的边界条件不能两端都满足"。

前置知识：Ch 6（扰动方法）、Ch 7（边界层 + 匹配）、Fourier 级数（§8.1 复习）、Laplace 方程基础。

核心方程与概念

1. 经典分离变量解（§8.1, eq. 8.1-8.9）

Laplace 方程： $\(\nabla^2 U \equiv \frac{\partial^2 U}{\partial X^2} + \frac{\partial^2 U}{\partial Y^2} = 0 \tag{8.1}\)$

矩形域上的 Dirichlet 问题（8.3）：\(0 \le X \le L, 0 \le Y \le H\)，三边齐次 \(U = 0\)，顶边 \(U(X, H) = F(X)\)。

分离变量解（8.8-8.9）： $\(U(X, Y) = \sum_{n=1}^\infty c_n \sinh\left(\frac{n\pi Y}{L}\right) \sin\left(\frac{n\pi X}{L}\right)\)$ $\(c_n = \frac{2}{L \sinh(n\pi H/L)} \int_0^L F(\tilde X) \sin\left(\frac{n\pi \tilde X}{L}\right) d\tilde X\)$

注意（8.9 末尾）：对任何 \(H/L\) 都成立，但无穷级数形式不简洁。长波极限 \(H/L \to 0\) 给出更简洁的 ODE 解。

2. 细长矩形的 Dirichlet 问题（§8.2, eq. 8.10-8.20）

无量纲化（8.11-8.12）：\(X = Lx, Y = Hy, U = \bar U u\)，得 $\(\varepsilon^2 u_{xx} + u_{yy} = 0, \quad \varepsilon = H/L \tag{8.12a}\)$ - 顶/底：\(u(x, 0) = 0, u(x, 1) = f(x)\) - 左/右：\(u(0, y) = g_0(y), u(1, y) = g_1(y)\)

外解（8.13-8.14）：\(u \sim u_0 + \varepsilon^2 u_1 + \varepsilon^4 u_2 + \ldots\) - \(O(\varepsilon^0)\)：\(u_{0yy} = 0, u_0(x, 0) = 0, u_0(x, 1) = f(x)\) → \(u_0 = f(x) y\) - \(O(\varepsilon^2)\)：\(u_{1yy} = -u_{0xx} = -f''(x), u_1(x, 0) = 0, u_1(x, 1) = 0\) → \(u_1 = \frac{1}{6} f''(x) (y - y^3)\) - ...

问题：外解不满足 \(u(0, y) = g_0(y), u(1, y) = g_1(y)\)（因为 \(u_0 = f(x) y\) 一般不匹配 \(g_0, g_1\)）。需要边界层。

内解（8.15-8.18）：设 \(X = x/\varepsilon\)（在 \(x = 0\) 附近的"内尺度"）： - \(\alpha = 1\)（主导平衡）→ 内 ODE 是 \(U_{XX} + U_{yy} = 0\)（回到 Laplace 方程！） - 边界条件：\(U_0(0, y) = g_0(y), U_0(X \to \infty, y) = f(0) y\) - 提取"边界层修正" \(V(X, y) = U_0(X, y) - f(0) y\)，得 \(V\) 满足"建筑块"问题 (8.18)： \(V_{XX} + V_{yy} = 0, V(X, 0) = V(X, 1) = 0, V(0, y) = g_0(y) - f(0) y, V(\infty) = 0\)

分离变量解（8.19）： $\(V(X, y) = \sum_{n=1}^\infty c_n e^{-n\pi X} \sin(n\pi y), \quad c_n = 2 \int_0^1 [g_0(y) - f(0) y] \sin(n\pi y) dy\)$

复合解（8.20）： $\(u_0^{unif}(x, y) = u_0(x, y) + V^L(X^L, y) + V^R(X^R, y), \quad X^L = x/\varepsilon, X^R = (x-1)/\varepsilon\)$

注意：内解中的 \(e^{-n\pi X}\) 因子在 \(X \to \infty\) 时衰减——匹配自然成立（没有"不可匹配"的指数发散）。

3. 绝缘导线问题（§8.3, eq. 8.21-8.31）

Neumann 边界条件（顶/底无通量）： $\(u_y(x, 0) = u_y(x, 1) = 0 \tag{8.21b}\)$ $\(u(0, y) = g_0(y), u(1, y) = g_1(y) \tag{8.21c}\)$

外解（8.22-8.23）：\(u_0(x, y) = C_2(x)\)（与 \(y\) 无关！）

可解性条件（8.24-8.25）：从 \(O(\varepsilon^2)\) 得 \(u_1 = -\frac{1}{2} C_2''(x) y^2 + C_3(x) y + C_4(x)\)，\(u_1|_y = 0\) → \(C_3 = 0\)，\(u_1|_y = 1\) → \(C_2''(x) = 0\)。

这是"边界层反过来决定外解"的范例： - \(C_2(x) = D_1 x + D_2\)，两个常数未定 - 必须从边界层（8.26-8.30）确定 \(D_1, D_2\)

边界层（8.28）：\(V_{XX} + V_{yy} = 0, V_y(X, 0) = V_y(X, 1) = 0, V(0, y) = g_0(y) - D_2, V(\infty) = 0\)。

分离变量解（8.29）： $\(V(X, y) = \sum_{n=0}^\infty c_n e^{-n\pi X} \cos(n\pi y)\)$

关键约束（8.30）：\(c_0 = 0\) 来自 \(V(\infty) = 0\) 约束（\(n = 0\) 项不退化），即 $\(D_2 = C_2(0) = \int_0^1 g_0(y) dy\)$

类似地，\(C_2(1) = \int_0^1 g_1(y) dy\)。最终外解（8.31）： $\(u_0(x, y) = x \left(\int_0^1 [g_1(y) - g_0(y)] dy\right) + \int_0^1 g_0(y) dy\)$

信息流方向（8.3 末段）： - Dirichlet（8.12）：外解独立确定 → 边界层匹配外解 - Neumann（8.21）：边界层决定外解的"等效边界条件"（常数 = \(g\) 的平均值）

这是边界层理论的"两种流向"，对实际工程应用非常重要。

4. 非均匀绝缘导线（§8.4, eq. 8.32-8.37）

变宽导线：\(0 \le Y \le H f(X/L)\)，\(f(x)\) 给定。无通量条件（8.32e）变为 $\(-\frac{H}{L} f'(X/L) \frac{\partial U}{\partial X}\bigg|_{Y = Hf(X/L)} + \frac{\partial U}{\partial Y}\bigg|_{Y = Hf(X/L)} = 0\)$

无量纲化（8.33）： $\(u_y(x, f(x)) - \varepsilon^2 f'(x) u_x(x, f(x)) = 0 \tag{8.33d}\)$

外解（8.34）：同前 \(u_0 = C_2(x)\)。

\(O(\varepsilon^2)\) 兼容性条件（8.35-8.37）： $\(u_1 = -\frac{1}{2} C_2''(x) y^2 + C_3(x) y + C_4(x)\)$ $\(u_1|_y = f(x) = f'(x) C_2'(x)\)$ $\(\Rightarrow \frac{d}{dx} \left[f(x) \frac{dC_2}{dx}\right] = 0 \tag{8.37}\)$

这是一个"加权 ODE"——\(f(x)\) 是变系数。边界条件仍然由边界层决定 \(C_2(0) = \int_0^{f(0)} g_0(y) dy\)（注意积分上限是 \(f(0)\)，因为边界层区域在 \(0 \le y \le f(0)\)）。

这是 domain perturbation（域扰动）问题——边界不是直线。Laplace 方程的分离变量解不再可用，但长波近似仍然有效。

5. 与 Ch 5-7 的对比

章节 主题 关键技术
Ch 5 自相似解 标度对称性 + ODE 约化
Ch 6 扰动方法 主导平衡 + 重新缩放
Ch 7 边界层（ODE） 外 + 内 + 匹配
Ch 8 长波近似（PDE） 细长 + 降维 + 边界层

Ch 8 是 Ch 7 的"高维推广"——把"一维细长"推广到"二维细长"（slender rectangle）。

关键结论

  1. 长波近似 = 二维 Laplace 方程的"降维"：\(\varepsilon^2 u_{xx} + u_{yy} = 0\) 在 \(\varepsilon \to 0\) 极限下退化为 \(u_{yy} = 0\)（1D ODE）。这是物理学中"薄层近似"的数学基础——润滑理论（lubrication theory）、薄膜流动、河口流动、晶圆涂层、泪液膜都基于此。
  2. Dirichlet 边界条件 vs. Neumann 边界条件的信息流方向：
  3. Dirichlet：外解 \(u_0 = f(x) y\) 完全由 \(f(x)\) 决定 → 边界层"修正"端点条件
  4. Neumann：外解 \(u_0 = C_2(x)\)（独立于 \(y\)）→ 边界层决定 \(C_2\) 的端点值（"等效边界条件"）

这是工程上重要的区分——许多物理问题（电流、热流、扩散）使用 Neumann 条件，外解的有效性依赖边界层分析。 3. 可解性条件（solvability condition）= 高阶问题的相容性约束：对 Neumann 问题，\(O(\varepsilon^2)\) 问题只有在 \(C_2''(x) = 0\) 时有解。这就是"Fredholm alternative"的微缩版——线性算子有非平凡解当且仅当"右端项与零空间正交"。作者没有用这个名字，但本质相同。 4. 绝缘导线的"等效电路"理解（8.31 末尾）： - \(u_0(x, y) = x \cdot \bar g + g_0^{avg}\)，其中 \(\bar g = \int_0^1 (g_1 - g_0) dy\) 是平均电压差、\(g_0^{avg} = \int_0^1 g_0 dy\) 是端点 0 的平均电压 - 电流（\(\partial u / \partial x\)）沿 \(x\) 方向是常数——这是欧姆定律的极限形式 - 电阻 \(R = 1/\bar g\)（无量纲形式）——长波近似自动导出电路模型！ 5. Domain perturbation（域扰动, 8.32-8.37）： - 当边界 \(Y = Hf(X/L)\) 是 \(X\) 的函数时，Laplace 方程 + 边界 \(f\) 给出变系数 ODE（8.37） - 变系数反映"导线宽度变化" - 物理上：细窄处电阻大、电流小；细宽处电阻小、电流大——等价于"变截面导线" 6. §8.5 末段的"降维模型"哲学：把 2D PDE 约化为 1D ODE，降维的代价是丢失"横向模式"（\(y\) 方向的高阶 Fourier 模式被边界层吸收）。但端点信息通过边界层传递到外解——这是"降维 + 边界层"的核心价值。 7. Ch 5 多孔介质方程的 PDE 来源（习题 8.6）：从 Navier-Stokes 出发做长波近似（\(h_t + \partial_x \int_0^h u dy = 0\)）给出多孔介质方程 \(h_t = \frac{1}{3} \partial_x (h^3 \partial_x h)\)。Ch 5 的自相似解有严格的 Navier-Stokes 长波推导。 8. KdV 方程的 PDE 来源（习题 8.5）：从"自由表面 + 不可压缩流"做长波近似 → KdV 方程 \(F_t + 6 F F_x + F_{xxx} = 0\)（Ch 8.5 末段预告，Ch 9 完成）。这是孤子物理的核心。

挑战和开放性问题

  1. 多边界层 / 多内区（§8.2-8.3）：作者处理了两个端点边界层（\(x = 0, 1\)），但每个边界层的"内部" 也有可能"再嵌套"——例如"内-内层"（inner-inner layer）。这种"层中层"在 §7.4 末段出现过（triple deck）。在 PDE 情形下，两维空间让"层中层"更复杂。
  2. 非均匀截面的可解性条件（8.37）：\(\frac{d}{dx}[f(x) dC_2/dx] = 0\) 看起来简单，但当 \(f(x)\) 变化剧烈（如 \(f \to 0\) 在某点）时，\(dC_2/dx\) 出现奇点（电流密度 → ∞）。这是工程上的"瓶颈效应"——导线最细处最易烧毁。严格处理需要"边界层 + 边界层"的多尺度分析。
  3. 大长宽比极限的有效性（8.12a）：当 \(\varepsilon = H/L \to 0\) 时，细长近似渐近精确。但对 \(\varepsilon = O(1)\)（如 \(H \sim L\)），外解 \(u_0 = f(x) y\) 是错误的——应该回到完整 Laplace 方程。作者没有给出"何时近似失效"的判据。一个常用判据是比较 \(\varepsilon\) 和 \(f''(x) H^2\)——如果 \(f''(x) H^2 > \varepsilon\)，\(O(\varepsilon^2)\) 项不能忽略。
  4. §8.3 的"信息流方向"是物理洞察。作者在 §8.3 末段明确指出 Dirichlet 和 Neumann 问题的"信息流"是相反的。这在工程上有重要含义：当你设计一个 PDE 模型时，先识别"哪个方向主导"（哪条边界是"主"），再选择主边界条件的种类（Dirichlet / Neumann / Robin）。作者没有进一步展开，但这是一个值得思考的建模哲学。
  5. 缺失主题：圆柱绕流的细长近似（高雷诺数流体力学）：Ch 8 的 Laplace 方程是"势流"极限。真实流体有黏性（边界层在表面，Oseen 方程）和惯性（对流项）。Oseen 方程的细长近似是 Prandtl 的开创性工作（1904），引出了边界层理论（Ch 7 末段引用的 [76]）。Ch 12 末尾"流体动力学应用"会涉及。
  6. 第 8.5 末段的"长波极限"失效情形：当边界 \(f(x)\) 出现"尖角"（如 \(f(0) = 0\)），\(x = 0\) 附近的边界层可能"塌缩"为奇点（corner singularity）。这在 PDE 中产生"\(r^{2/3}\) 奇异性"——经典结果是满足 Dirichlet 边界条件的 Laplace 方程在 \(90°\) 角有 \(r^{2/3}\) 梯度奇异性。作者没提。
  7. 第 8.3 末段的"等效边界条件"在多孔介质中很重要：多孔介质中的 Darcy 流动（渗流）用 Laplace 方程 + 边界条件描述。当边界是 Neumann 条件（通量给定），"等效边界条件"对应注入/抽取井的流量。这是地下水文学、油藏工程的标准工具。Ch 11 矩方法可以推广到多井问题。
  8. 润滑理论的应用范围（§8.5 末段）：作者提到"machine parts"（机械零件之间的油膜）、"tear films on the eye"（泪液膜）。这是软物质物理和生物流体力学的核心。现代润滑理论包括弹性润滑（elastohydrodynamic lubrication, EHL）—— 当润滑层很薄时，基底弹性形变会影响压力分布。Ch 12 末尾"rivulet in a wedge"是另一种 EHL 简化。

个人反思与批判性分析

本章是 Ch 7 边界层理论的"高维推广"——但结构上反而比 Ch 7 简洁：因为 Laplace 方程 + 长波近似的组合在结构上自然降维，匹配过程也对称。几个值得反思的点：

  1. "长波近似"是物理学中的"宏观极限"。作者在 §8.5 末段列举"oil between machine parts, tear films on the eye"——这些是"细长"几何的真实物理例子。长波近似的核心思想：当某个方向远小于其他方向时，垂直方向上的"细节"（如 Laplace 方程的 Fourier 模式）被边界层吸收，只留下"沿细长方向"的平均行为。这是降维物理的数学基础。

  2. "信息流方向"是工程上的"决策"工具。§8.3 末段的"信息流方向对比"是Ch 8 最深刻的洞察之一：Dirichlet 问题的边界条件直接决定外解；Neumann 问题的边界条件通过边界层决定外解的有效边界条件。这意味着：当一个工程问题同时含 Dirichlet 和 Neumann 条件时，"主"信息是哪条？——通常是 Dirichlet（直接传递到外解）。

  3. §8.3 末段的"等效电路"理解是把数学问题转化为工程直觉的好例子。\(u_0(x, y) = x \bar g + g_0^{avg}\) 直接就是"线性电压分布"——这是均匀导线的欧姆定律。长波近似自动给出"降维的等效电路"——这是非常漂亮的物理洞察。

  4. 作者没有展开"匹配的高阶修正"。§8.2 末段给出了首阶匹配，高阶（\(O(\varepsilon^2)\)）需要重新做一次匹配——确定内解中 \(c_n\) 的高阶修正。这与 Ch 7 末段相同——本书的"高阶匹配"留作习题，没有完全展开。

  5. "边界层 + 外解"是 Ch 7-8 的核心模式：

  6. 外解：在"大区域"上构造
  7. 内解：在"小区域"（端点）上构造
  8. 匹配：两者在重叠区域一致
  9. 复合：两者相加减去重叠

Ch 8 把这一模式在 PDE 上完整实施——展示了多尺度方法的 PDE 版本。Ch 9-10 的弱非线性、快/慢系统沿用同一思想（不过变量是时间，不是空间）。

  1. "Domain perturbation"（§8.4）是建模中的常见技术：当边界不是直线时，直接求解往往困难。作者的方法是"在原问题基础上做域扰动"——把"变边界"看作"直边界 + 小扰动 \(\varepsilon^2 f'(x) u_x\)"。这是工程上常用的简化。

  2. Ch 8 与 Ch 5-6 的衔接：Ch 5 给出多孔介质方程的精确解（自相似），Ch 8 习题 8.6 给出它的 PDE 来源（Navier-Stokes + 长波近似）。这是"逆向追溯"——精确解 + 它的 PDE 来源 + 长波近似三者形成一个完整的知识闭环。这是本书的精妙结构——习题不仅是练习，也是全书知识的"再连接"。

  3. Ch 8 与 Ch 9-10 的衔接：Ch 8 处理长波 + 边界层（空间多尺度），Ch 9 处理多时间尺度（时间多尺度）。Ch 10 处理快/慢系统（时间多尺度的极端情形）。三章构成"多尺度方法"的完整体系：空间多尺度（Ch 7-8）+ 时间多尺度（Ch 9-10）。

重要参考文献

  • [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] G. Bayada and M. Chambat, "The Transition between the Stokes Equations and the Reynolds Equation: A Mathematical Proof," Applied Mathematics and Optimization 14 (1986) 73–93. (Stokes-Reynolds 转换的严格证明)
  • [X4] A. Oron, S.H. Davis, and S.G. Bankoff, "Long-Scale Evolution of Thin Liquid Films," Reviews of Modern Physics 69 (1997) 931–980. (薄膜长波演化)
  • [X5] J. Eggers and M.A. Fontelos, Singularities: Formation, Structure, and Propagation, Cambridge University Press, 2015. (奇点的形成、结构和传播)
  • [X6] Stephen H. Davis, Theory of Solidification, Cambridge University Press, 2001. (凝固理论中的长波近似)
  • [X7] D.J. Acheson, Elementary Fluid Dynamics, Oxford University Press, 1990. (流体力学入门)
  • [X8] L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Pergamon, 1987. (流体力学的金标准)
  • [X9] A.J. McCarty and H. Soo, Thin Film Coating Processes, World Scientific, 2009. (薄膜涂层过程)
  • [X10] J.H. Snoeijer, "Free-Edge Thin Film Flows," Annual Review of Fluid Mechanics 53 (2021) 369–393. (自由边缘薄膜流的现代综述)