跳转至

第 5 章 微分方程的自相似尺度解（Self-Similar Scaling Solutions of Differential Equations）

作者

同 Ch 1-4，Thomas Witelski（Duke University，数学系）和 Mark Bowen（Waseda University，国际理工学中心）合著。本章是 Part II（Solution Techniques）的开篇——四章"求解方法"的第一章。

本章在全书中的定位：Ch 4 无量纲化的"高级形式"——求保持方程形式不变的标度变换，并把 PDE 降维为 ODE。Ch 4 是"用尺度减少参数"，Ch 5 是"用尺度构造解"。这是从"分析"到"求解"的过渡。

内容概述

本章回答一个核心问题：什么样的 PDE 具有"标度不变性"（scale invariance），如何利用这种不变性把 PDE 约化为 ODE，并给出精确解（exact solution）？

章节结构： 1. §5.1 寻找标度不变对称性：通过引入 \(u = U \tilde u, x = L \tilde x, t = T \tilde t\) 寻找使方程形式保持不变的 \((U, L, T)\) 关系。例：inviscid Burgers 方程 \(u_t + u u_x = 0\) 要求 \(UT/L = 1, UL = 1\) → \(L = T^{1/2}, U = T^{-1/2}\)。 2. §5.2 确定相似解的形式：通过要求 \(u^a x^b t^c\) 是标度不变的（Π 定理的"代数版本"），构造相似变量 \(\eta = x t^{-1/2}\) 和相似函数 \(f = t^{1/2} u\)。这给出 \(u(x, t) = t^{-1/2} f(x t^{-1/2})\) 的形式。 3. §5.3 求解相似函数：代入 PDE 给出 ODE（5.12），求解得到 \(f = \eta\) 即 \(u = x/t\)（膨胀扇 / rarefaction fan），与 Ch 2 的特征线解一致。 4. §5.4 进一步评注：与 Ch 4 的"固定尺度"对比——Ch 4 给出确定的尺度，Ch 5 给出一族解。 5. §5.5 热方程的相似解： - §5.5.1 Cauchy / 源型解（5.17）：\(u = (M/\sqrt{4\pi t}) e^{-x^2/(4t)}\) - §5.5.2 Boltzmann 解（5.20）：\(u = 1 - \text{erf}(x/\sqrt{4t})\)（半无界区间上的扩散边界层） 6. §5.6 非线性扩散方程（多孔介质方程）：\(\partial_t u = \partial_x (u^3 \partial_x u)\)，有限质量自相似解给出 \(u \sim (C - x^2/10 t^{2/5})^{1/3} / t^{1/5}\)，有限支撑（compact support）的"展平水滴"。

核心结论：如果 PDE 在标度变换下形式不变（标度对称性），则存在一族标度不变解——参数由 \(T\) 自由选择。把 \(u = t^{-\alpha} f(\eta = x t^{-\beta})\) 代入 PDE，PDE → ODE，\(t\) 和 \(\eta\) 分离。这是降低维数的标准技巧。

前置知识：Ch 4（无量纲化和 Buckingham Π 定理）、Ch 2（特征线法、激波）、单变量 ODE 求解（分离变量、积分）。

核心方程与概念

1. 标度对称性（§5.1, eq. 5.1-5.6）

标度变换：\(u(x, t) = U \tilde u(\tilde x, \tilde t), x = L \tilde x, t = T \tilde t\)，\((U, L, T)\) 是待定正常数。

标度不变性：若变换后的方程与原方程形式相同（且至少一个 \(U, L, T\) 自由），则原方程具有标度对称性（scaling symmetry）。

Burgers 例子（5.3-5.5）：\(u_t + u u_x = 0, \int_0^\infty u dx = 4\)。代入标度变换得 \(U T / L = 1\)（来自对流项）和 \(UL = 1\)（来自积分），解出 \(L = T^{1/2}, U = T^{-1/2}\)。一族解： $\(u(x, t) = T^{-1/2} \tilde u(T^{-1/2} x, T^{-1} t) \tag{5.6}\)$

2. 相似变量与相似函数（§5.2, eq. 5.7-5.11）

标度不变量：\(\Pi = u^a x^b t^c\)。在标度变换下独立于 \(T\) 的条件是 \([U^a L^b T^c] = T^0\)。

相似变量（\(a = 0\)）：\(\eta = x t^c\)，对 Burgers 例子 \(c = -1/2\) → \(\eta = x t^{-1/2}\)。

相似函数（\(a = 1\)）：\(f = t^c u\)，对 Burgers 例子 \(c = 1/2\) → \(f = t^{1/2} u\)。

Π 定理类比（§5.2 末段）：\(F(\Pi_1, \Pi_2) = 0 \Rightarrow F(\eta, f) = 0 \Rightarrow f = f(\eta)\)。

最终形式（5.11）： $\(u(x, t) = t^{-1/2} f(\eta = x t^{-1/2})\)$

3. Burgers 方程的相似解（§5.3, eq. 5.12-5.14）

代入 (5.11) 到 (5.3) 得 ODE（5.12）： $\(f + \eta f' - 2 f f' = 0, \quad f(0) = 0, \quad \int_0^\infty f d\eta = 4\)$

等价于 \((\eta f - f^2)' = 0\)，积分给出 \(\eta f - f^2 = C\)。\(f(0) = 0\) 要求 \(C = 0\)，故 \(f(\eta) = \eta\) 或 \(f(\eta) = 0\)。非平凡解： $\(u(x, t) = \frac{x}{t} \tag{5.14}\)$

这是膨胀扇（rarefaction fan），与 Ch 2 特征线法在初始条件 \(p(x, 0) = \max(0, 1 - |x|)\) 下的解的右半部分一致。

4. 热方程的相似解（§5.5, eq. 5.15-5.20）

热方程：\(\partial_t u = \partial_x^2 u\)。标度要求 \(L = T^{1/2}\)（扩散时间尺度，\(T \sim L^2/D\)），\(U\) 由边界条件决定。

Cauchy / 源型解（5.17，全空间 + 有限质量 \(M\)）： $\(u(x, t) = \frac{M}{\sqrt{4\pi t}} e^{-x^2/(4t)}\)$

这是高斯函数，宽度 \(\sim \sqrt{t}\)，高度 \(\sim 1/\sqrt{t}\)，保持总质量 \(M\)。

Boltzmann 解（5.20，半无界区间 + Dirichlet \(u(0,t) = 1, u(\infty) = 0\)）： $\(u(x, t) = 1 - \text{erf}\left(\frac{x}{\sqrt{4t}}\right) = 1 - \frac{2}{\sqrt{\pi}} \int_0^{x/\sqrt{4t}} e^{-s^2} ds\)$

这是误差函数型扩散边界层，半宽度 \(\sim \sqrt{t}\)（扩散前沿），与 \(u = 1\) 渐近的常数边界条件匹配。

关键差异： - Cauchy 解质量守恒（高度衰减、宽度增长） - Boltzmann 解在 \(x = 0\) 处固定为 1，扩散前沿是误差函数

长时间行为（5.5.1 末段）：Cauchy 解在 \(t \to \infty\) 时是任何合理初始条件的吸引子——只要给定质量、中心、方差，长时间后解都趋近于 Cauchy 解。这是自相似吸引子（self-similar attractor）的概念。

5. 多孔介质方程（§5.6, eq. 5.21-5.23）

多孔介质方程（非线性扩散）： $\(\partial_t u = \partial_x \left( u^3 \partial_x u \right), \quad u_x(0, t) = 0, \quad \int_{-\infty}^{\infty} u dx = M \tag{5.21}\)$

物理意义：极黏性流体（蜂蜜）在多孔介质/平板上的扩散。\(u\) 是流体层厚度。

标度要求：\(U^3 = L^2/T\)（来自非线性项）+ \(UL = 1\)（质量守恒） → \(L = T^{1/5}, U = T^{-1/5}\)。

相似解（5.22-5.23）： $\(u(x, t) = \frac{1}{t^{1/5}} \left[ C_2 - \frac{x^2}{10 t^{2/5}} \right]^{1/3} \quad \text{for } |x| < x_*(t) = \eta_* t^{1/5}\)$ $\(u(x, t) = 0 \quad \text{for } |x| \ge x_*(t)\)$

有限支撑（compact support）特性：解在 \(x = \pm x_*(t)\) 处突然变为零——这是激波的扩散类比（接触线 / contact line）。这种"有限支撑解"是非线性守恒律的典型行为。

与 Ch 8 长波、KdV 方程的连接：多孔介质方程的 \(u^3\) 系数是非线性的（非线性扩散），与 KdV 方程的非线性对流项互补。多孔介质方程 ≈ 非线性扩散守恒律，Burgers 方程 ≈ 非线性对流守恒律。两者都是非线性 PDE 求解的"原型"。

关键结论

  1. 标度不变性 = 解的可参数化一族：如果 PDE 在 \((U, L, T)\) 变换下形式不变，则任何解 \(u(x, t)\) 都对应一族解 \(u_T(x, t) = T^{-\alpha} f(x T^{-\beta} / t^\gamma)\)。这一族是 PDE 的精确解子集——不需要积分常数。Ch 5 的核心方法就是找到这种"自相似族"。
  2. 相似解 = PDE 降维为 ODE：把 \(u = t^{-\alpha} f(x t^{-\beta})\) 代入 PDE，所有 \(t\)-依赖分离，给出关于 \(f(\eta)\) 的 ODE。这一步把2 个独立变量 \((x, t)\) 减少到1 个 (\(\eta\))——这是 PDE 求解的"超能力"。
  3. 自相似解作为长时间吸引子：热方程的 Cauchy 解 (5.17) 在 \(t \to \infty\) 时是任何具有相同质量、中心、方差的初始条件的吸引子（§5.5.1 末段）。这是自相似解的"普遍性"（universality）——长时间行为只依赖于宏观守恒量（mass、center、variance），不依赖于初始形状。这是 Ch 11（Reduced Models for PDE）矩方法的伏笔。
  4. 有限支撑解：多孔介质方程的相似解 (5.23) 在 \(x = \pm x_*(t)\) 处突然变为零——这是非线性的结果（线性热方程永远不会有有限支撑）。这种"扩散激波"在物理上是接触线（contact line）——水滴的边缘。Ch 12 末尾的"Plateau 边界"问题也涉及这种有限支撑。
  5. 相似变量 \(\eta = x / \sqrt{t}\) 的几何意义：它把所有时间 \(t > 0\) 的解映射到单一的 \(\eta\) 轴。这意味着解在 \((x, t)\) 平面上的行为由 \(f(\eta)\) 完全决定——不同的 \(t\) 只是 \(f\) 的重新缩放。这是"自相似"的几何含义。
  6. Boltzmann 解 vs. Cauchy 解：同样对热方程，不同的边界条件给出不同的相似解——Boltzmann 解（半无界 + Dirichlet）和 Cauchy 解（全空间 + 有限质量）。这说明相似解的"分类"依赖于边界条件的选择。
  7. 自相似解 = 渐近分析的工具：作者在 §5.4 末尾说"解可能在某些极限下趋近于相似解"——称为渐近自相似解。这是 Ch 6-7 的核心：当某个无量纲参数 \(\Pi\) 是小量时，精确解趋近于某个相似解。Ch 5 的相似解给出了目标函数 \(f\)，Ch 6-7 的扰动方法给出如何"接近"这个目标。

挑战和开放性问题

  1. 相似解的"存在性"非平凡：作者默认相似解存在，然后推导它满足的 ODE。严格的存在性需要：(a) 标度对称性确实存在（代数条件已检查），(b) ODE 边界条件匹配（\(t \to 0\) 和 \(t \to \infty\) 行为），(c) ODE 解的唯一性。对多孔介质方程，ODE (5.22) 是二阶非线性 ODE + 一个移动边界 \(x_*(t)\)，需要巧妙处理（Stellar 1980s）。作者在 §5.6 没展开这个分析。
  2. 多参数相似解（多标度族）：本章只处理单参数相似解（\(T\) 自由）。多参数族（如 \(T_1, T_2\) 自由）需要更复杂的标度分析。例如 Navier-Stokes 方程的 Kolmogorov 湍流解含外尺度和内尺度（Kolmogorov 微尺度）——两个独立参数。本章没有处理。
  3. Burgers 相似解 \(u = x/t\) 在 \(t \to 0\) 时的行为：当 \(t \to 0\) 时 \(u = x/t\) 在固定 \(x \neq 0\) 处发散。这意味着这种相似解只在 \(t \to \infty\) 有效（远离奇点）。这是"渐近自相似"的典型——\(t \to 0\) 的行为由初始条件决定，\(t \to \infty\) 的行为由相似解决定。
  4. §5.4 末尾的"非比例标度"：作者提到"模化或去掉某些项"可以创造标度不变性。这是奇异扰动方法（singular perturbation）的伏笔——例如去掉 Burgers 方程的 \(\kappa u_{xx}\) 项（无黏极限）时，方程变为非扩散的，反而具有更简单的标度结构。这种"减少物理内容获得数学结构"的方法是建模的常见技巧。
  5. 多孔介质方程的 \(\eta_*\) 决定（5.22-5.23）：作者把 \(\eta_*\) 当作常数，但实际需要总质量守恒 \(\int_{-\eta_*}^{\eta_*} f d\eta = M\) 联立求解。这是移动边界问题（Stefan 问题）的一种简化，严格处理需要 Stefan 条件（界面能量平衡）——多孔介质方程省略了这一条件。Ch 12 末尾的"Plateau 边界"是类似问题。
  6. §5.5.1 末段的"中心"和"方差"作为吸引子参数：作者提到 Cauchy 解在 \(t \to \infty\) 时的吸引子由（质量、中心、方差）决定。这隐含热方程的自相似解空间是三维的（三个守恒量），而不是一维的（一个尺度参数 \(T\)）。对非线性扩散，守恒量数目可能不同——例如多孔介质方程的动量（一阶矩）守恒但方差（二阶矩）不守恒（因为解有有限支撑）。这种"非线性动理学"理论由 De Masi, Presutti 等人在 1990s 系统化。
  7. 缺失主题：群论观点。Ch 5 的"标度对称性"实质上是一维李群（\(T \in \mathbb{R}^+\)）作用在解空间上。现代 PDE 理论（Ovsiannikov, Olver, Bluman）用完整的李群（包括平移、旋转、Galilean 变换、标度变换）生成无穷多守恒律。本章没用李群语言，但 §5.5 末尾"heat equation has several classes of symmetries"是这种观点的入口。

个人反思与批判性分析

本章是全书数学结构最简洁的一章——核心思想只有一句话："如果 PDE 在标度变换下不变，则存在一族自相似解"。但应用范围却非常广：热方程、Burgers 方程、多孔介质方程都适用。这种"思想简单、应用广泛"是优秀数学方法的标志。

从建模哲学角度看，本章有几个值得反思之处：

  1. 自相似解的"特殊性" vs. "普遍性"。作者在 §5.4 强调"即使问题不严格标度不变，解也可能在某些极限下趋近于自相似解"。这意味着自相似解是"渐近行为"而不是"严格行为"——一个真实物理系统可能永远不会精确满足相似性，但在 \(t \to \infty\) 时趋近于自相似解。这是相似解的"普遍性"——长时间行为由宏观守恒量决定，与初始细节无关。这种思想在重整化群（Kadanoff, Wilson）和临界现象（统计物理）中是核心。

  2. 相似解 = "降维的代价"是"解的家族"。把 PDE 降为 ODE 看似"得到解"，但实际上失去了某些自由度——解不再依赖于所有初始条件，而只依赖于标度不变的组合。例：Burgers 相似解 \(u = x/t\) 没有"参数"——它是一个特定解，而不是一族解的全体。这意味着相似解丢失了"通解"中的一部分信息。在物理上，这意味着真实系统的细节（如初始条件）会在长时间演化中被"磨平"，只剩下宏观守恒量决定的行为。这是"熵增"的数学版本。

  3. "自相似" vs. "行波" 的区别。Ch 2 的行波（travelling wave）\(\rho(x, t) = P(x - ct)\) 和 Ch 5 的自相似 \(u = t^{-1/2} f(x t^{-1/2})\) 都是 PDE 的特殊解。区别：

  4. 行波：形状恒定，平移速度恒定（无标度变换）
  5. 自相似：形状变化（标度因子 \(t^{-\alpha}\)），扩展速度 \(\sim t^{-\beta}\)

在 \((x, t)\) 平面上，行波是直线（斜率 \(c\)），自相似是曲线（\(x \sim t^\beta\)）。多孔介质方程的接触线 \(x_*(t) = \eta_* t^{1/5}\) 就是这种"自相似推进"——它不保持恒定速度，而是 \(dx_*/dt \sim t^{-4/5}\)（减速扩散）。

  1. §5.5.1 末段的"中心"和"方差"是 Ch 11 矩方法的伏笔。作者提到 Cauchy 解作为 \(t \to \infty\) 吸引子由（质量、中心、方差）三个守恒量决定。这正是 Ch 11（Reduced Models for PDE）矩方法的核心——用有限个矩（mass、center、variance）代替整个解，得到常微分方程组（ODE）作为 PDE 的"约化模型"。这是 Ch 5 思想的延伸：自相似解是一种"约化"（用 ODE 代替 PDE），矩方法是另一种"约化"（用有限维 ODE 代替 PDE）。

  2. 多孔介质方程的"水滴模型"。\(u(x, t) \sim (C - x^2/10 t^{2/5})^{1/3}\) 描述了一个有限大小的液体团在多孔介质上的扩散。这与 Ch 12 末尾的"Plateau 边界"（液膜边缘）问题在物理上同源——液体边缘的扩散。Ch 12 给出的是液膜在大气中的流动，Ch 5 给出的是液膜在多孔介质中的扩散。两个问题的解有相似的数学结构（\(u \sim (C - x^2/A t^B)^{1/n}\)）。

  3. 缺少对"自相似性破缺"的讨论。本章全部例子都成功构造了自相似解。但真实物理系统中自相似性常被打破——如粘性指进（viscous fingering）现象，流体的有效扩散系数依赖于局部梯度，导致非线性反馈破坏标度不变性。作者没提这些"反例"，这会让读者误以为自相似解是普适的。这是教学上的简化，但也使学生失去对"反例"（如湍流的反常标度律）的敏感度。

  4. 与 Ch 8 长波、KdV 方程的衔接。Ch 5 的非线性扩散（多孔介质方程）与 Ch 8 的非线性色散（KdV 方程）形成对照：

  5. 多孔介质：\(u_t = (u^3 u_x)_x\)（非线性扩散） → 扩散前沿
  6. KdV：\(u_t + 6 u u_x + u_{xxx} = 0\)（非线性色散） → 孤子

两者都是非线性 PDE，但耗散 vs. 色散导致完全不同的解结构。Ch 5 给出扩散类，Ch 8 出色散类——这是 Part II 对非线性 PDE 的两种处理范式。

  1. §5.4 末段的"渐近自相似"是 Ch 6-7 的入口。作者说"在某些极限下解趋近于相似解"——这是匹配渐近展开（matched asymptotic expansion, Ch 7）的核心。Ch 5 的相似解提供了"内层"或"外层"的渐近行为，Ch 6-7 的扰动方法给出"中间层"的修正。

重要参考文献

  • [X1] G.I. Barenblatt, Scaling, Self-Similarity, and Intermediate Asymptotics, Cambridge University Press, 1996. (尺度分析与自相似性的标准专著)
  • [X2] G.I. Barenblatt, Similarity, Self-Similarity, and Intermediate Asymptotics, Plenum, 1979. (前书的精简版)
  • [X3] L.V. Ovsyannikov, Group Analysis of Differential Equations, Academic Press, 1982. (微分方程的群分析)
  • [X4] Peter J. Olver, Applications of Lie Groups to Differential Equations, Springer, 1993. (李群方法的标准教材)
  • [X5] George W. Bluman and Stephen Kumei, Symmetries and Differential Equations, Springer, 1989. (对称性方法)
  • [X6] Daniel Zwillinger, Handbook of Differential Equations, Academic Press, 1997. (ODE/PDE 求解手册)
  • [X7] John Ockendon, Sam Howison, Andrew Lacey, and Alexander Movchan, Applied Partial Differential Equations, Oxford University Press, 2003. (应用 PDE 教材，覆盖自相似解)
  • [X8] Lawrence C. Evans, Partial Differential Equations, AMS, 2010. (PDE 严格理论，包括 Stefan 问题)
  • [X9] David G. Duffy, Transform Methods for Solving Partial Differential Equations, Chapman & Hall, 2004. (变换方法，包括 Boltzmann 变换)
  • [X10] Lawrence S. Schulman, Techniques and Applications of Path Integration, Wiley, 1981. (路径积分与扩散方程的 Feynman-Kac 公式)