跳转至

第 2 章 输运方程（Transport Equations）

作者

同第 1 章，由 Thomas Witelski（Duke University，数学系）和 Mark Bowen（Waseda University，国际理工学中心）合著。本章是 Part I 的第 2 章，与第 1 章形成"时间维度上的 ODE"与"时空维度上的 PDE"的衔接：从 Ch 1 描述单点状态的时间演化，扩展到描述密度场随时间和空间的演化。

内容概述

本章回答两个核心问题：(i) 如何从守恒原理（质量、动量、粒子数）严格推导出连续介质的演化方程（输运 PDE）？ (ii) 给定一个一阶 PDE（线性/半线性/拟线性），如何用特征线法（method of characteristics）把 PDE 转化为可解的常微分方程组？

章节结构： 1. §2.1 Reynolds 输运定理：从"被动粒子被流场携带"出发，区分 Euler（场）描述和 Lagrangian（粒子）描述，给出"随体导数"（material derivative）公式。这是连续介质力学的核心数学工具。 2. §2.2 守恒律推导：把"每个物质 blob 上某物理量的总和不随时间改变"（守恒原理）转化为"局部 PDE 形式"。关键技术是 du Bois-Reymond 引理（积分恒等式 → 微分恒等式）。应用：质量守恒 → 连续性方程；动量守恒 → Cauchy 动量方程。 3. §2.3 线性对流方程 \(\rho_t + c \rho_x = 0\)：最简单的输运模型，引入行波解 \(\rho(x, t) = P(x - ct)\)、色散关系 \(\omega = \omega(k)\)、相速度 \(c_p(k) = \omega/k\) 等基本概念。 4. §2.4 线性对流方程组：通过特征值问题 \(M^T \mathbf{u} = \lambda \mathbf{u}\) 把耦合 PDE 解耦为独立行波。如果 \(M^T\) 的所有特征值为实数且有完整特征向量组，系统称为双曲型（hyperbolic）。 5. §2.5 特征线法：半线性 PDE \(p_t + c(x,t) p_x = r(x,t,p)\) 通过 Lagrangian 形式转化为 ODE 组（2.34a），给出一族参数化曲线（特征线）。 6. §2.6 拟线性方程的激波：当特征速度依赖于解本身（\(c = c(p)\)），特征线可能相交，导致多值解（物理上不可接受）。通过插入激波（shock, 移动间断点）并由 Rankine-Hugoniot 激波速度关系（eq. 2.48）确定其位置，恢复单值性。例子：无黏 Burgers 方程 \(p_t + p p_x = 0\)（2.39）的等面积律（equal-area rule）。

前置知识：Ch 1 全部内容（ODE、相线分析），单变量 PDE 基础，Fourier 级数（用于理解色散关系），Leibniz 积分法则。

核心方程与概念

1. Euler / Lagrangian 描述的转换（§2.1, eq. 2.1–2.5）

被动粒子的运动方程： $\(\frac{d\mathbf{X}}{dt} = \mathbf{v}(\mathbf{X}, t), \quad \mathbf{X}(0) = \mathbf{A} \tag{2.1b}\)$

随体导数（material derivative, eq. 2.5）：对跟随粒子的标量场 \(f_{part}(t) \equiv F(t; A) = f(\mathbf{X}(t; A), t)\)： $\(\frac{df_{part}}{dt} = \frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla f \equiv \frac{Df}{Dt}\)$

这是 Euler 描述下的"流体质点所感受到的 \(f\) 的总变化率"，包含局地变化 \(\partial f/\partial t\) 和对流变化 \(\mathbf{v} \cdot \nabla f\)。

Reynolds 输运定理（eq. 2.9）：对运动物质 blob 的积分： $\(\frac{d}{dt} \int_{D(t)} f \, dV = \int_{D(t)} \left[ \frac{\partial f}{\partial t} + \nabla \cdot (f \mathbf{v}) \right] dV\)$

这是把"对运动区域的积分"转化为"对固定区域的积分"的标准工具。Leibniz 法则在一维情形给出 (2.8)： $\(\frac{d}{dt} \int_{a(t)}^{b(t)} f(x, t) dx = \int_a^b \frac{\partial f}{\partial t} dx + f(b, t)\frac{db}{dt} - f(a, t)\frac{da}{dt}\)$

2. 守恒律与 du Bois-Reymond 引理（§2.2, eq. 2.10–2.18）

质量守恒 → 连续性方程（1D, 3D）： $\(\frac{\partial \rho}{\partial t} + \frac{\partial (\rho v)}{\partial x} = 0, \quad \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \tag{2.12, 2.13}\)$

带源/汇项的守恒律（eq. 2.16）： $\(\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = R\)$

Cauchy 动量方程（eq. 2.18）： $\(\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v}) = \nabla \cdot \boldsymbol{\sigma}\)$

其中 \(\boldsymbol{\sigma}\) 是应力张量，叠加本构关系（constitutive relation）就得到 Navier-Stokes 方程（流体）或弹性方程（固体）。这是连续介质力学的核心方程组。

du Bois-Reymond 引理（2.11）：若对所有子域 \([a,b]\) 都有 \(\int_a^b g(x) dx = 0\)，则 \(g(x) \equiv 0\)（假设 \(g\) 光滑）。这是把积分形式（弱形式，blob-wise）转化为微分形式（强形式，pointwise）的关键引理。

3. 线性对流方程（§2.3, eq. 2.19–2.22）

\[\frac{\partial \rho}{\partial t} + c \frac{\partial \rho}{\partial x} = 0, \quad \rho(x, t) = P(x - ct)\]

所有解都是恒定形状的行波。色散关系（dispersion relation）描述波数 \(k\) 和频率 \(\omega\) 的关系。

对五阶色散方程（2.22）\(\partial_t + c_1 \partial_x + c_3 \partial_x^3 + c_5 \partial_x^5 = 0\)： $\(\omega(k) = c_1 k - c_3 k^3 + c_5 k^5\)$

  • 相速度 \(c_p(k) = \omega(k)/k\)：单色波传播速度
  • 群速度 \(c_g(k) = d\omega/dk\)：波包能量传播速度（习题 2.3 推导）
  • 若 $c_p(k) = $ const，方程非色散（无色散 dispersionless）；否则色散，不同波数分量以不同速度传播，初始脉冲随时间分散
  • 若 \(\omega\) 含虚部，波随时间耗散（dissipation），如扩散方程 \(\rho_t = \rho_{xx}\) 的 \(\omega = -ik^2\)（虚部为 \(-k^2\)）

4. 线性对流方程组的对角化（§2.4, eq. 2.23–2.30）

两方程系统 \(\mathbf{p}_t + M \mathbf{p}_x = 0\)（eq. 2.29）通过特征值问题 \(M^T \mathbf{u} = \lambda \mathbf{u}\)（2.30）解耦。如果 \(M^T\) 有实特征值和完整特征向量组，系统是双曲型（hyperbolic）： $\(\mathbf{p}(x, t) = \sum_k c_k \mathbf{v}_k P_k(x - \lambda_k t)\)$

5. 特征线法（§2.5, eq. 2.31–2.36）

半线性 PDE \(p_t + c(x,t) p_x = r(x,t,p)\) 转化为常微分方程组： $\(\frac{dX}{dt} = c(X, t), \quad \frac{dP}{dt} = r(X, t, P) \tag{2.34a}\)$ $\(X(0) = A, \quad P(0) = f(A) \tag{2.34b}\)$

\((X(t), P(t))\) 在 \((x, t)\) 平面上的曲线是特征线。先解 \(X(t; A)\)，再代回 \(P\) 方程。这是 PDE → ODE 的标准降维方法。

例（eq. 2.35）：\(p_t + 2x p_x = x p^2\)，\(p(x,0) = 2 + \sin x\)。特征线 \(X = A e^{2t}\)，\(P = \frac{1}{C - \frac{1}{2} A e^{2t}}\)，最终解 \(p(x, t) = \left[ \frac{x}{2}(e^{-2t} - 1) + \frac{1}{2 + \sin(x e^{-2t})} \right]^{-1}\)。

6. 拟线性方程与激波（§2.6, eq. 2.37–2.50）

拟线性 PDE \(p_t + c(x, t, p) p_x = r(x, t, p)\) 的特征 ODE 耦合（2.38），特征线可相交（特征相交 = 多值解 = 物理不合理）。

无黏 Burgers 方程（eq. 2.39）：\(p_t + p p_x = 0\)，\(p(x,0) = \max(0, 1 - |x|)\)。特征线 \(X = Pt + A\)（P 沿特征线守恒）： - 区间 \(A \in [-1, 0]\)：膨胀扇（expansion fan / 稀疏波），特征线发散，斜率 \(\partial p / \partial x\) 随时间减小 - 区间 \(A \in [0, 1]\)：压缩波（compressive wave），特征线汇聚，\(t = 1\) 时形成折断（overturning），\(t > 1\) 时多值

激波修复：插入移动间断 \(x = x_s(t)\)，分离 \(p^-(x, t)\)（左侧）和 \(p^+(x, t)\)（右侧）。Rankine-Hugoniot 激波速度关系（eq. 2.48）： $\(\frac{dx_s}{dt} = \frac{q(p^+(x_s, t)) - q(p^-(x_s, t))}{p^+(x_s, t) - p^-(x_s, t)}\)$

其中 \(q(p)\) 是通量函数（flux）。对 Burgers 方程 \(q(p) = p^2/2\)，且 \(p^+ = 0\)、\(p^-\) 已知时给出 \(\frac{dx_s}{dt} = \frac{1 + x_s}{2(1+t)}\)，\(x_s(1) = 1\)，解 \(x_s(t) = \sqrt{2(1+t)} - 1\)。

等面积律（equal-area rule, Fig. 2.4）：激波位置不仅使解单值，还保证两侧被切除的多值区域面积相等——这实际上是质量守恒的几何表达。

关键结论

  1. Reynolds 输运定理 + du Bois-Reymond 引理 是从"blob-wise 守恒原理"推导"pointwise PDE"的两步法：先把 blob 积分用 Reynolds 转化为 Euler 积分（产生散度项），再用 du Bois-Reymond 把"对所有 blob 成立"转化为"被积函数恒为零"（产生局部 PDE）。这套两步法是连续介质力学、流体动力学、量子力学（Born 几率守恒）共同的语言。
  2. 行波 + 色散关系是分析线性 PDE 的核心工具：\(\omega = \omega(k)\) 同时给出传播（实部）和耗散（虚部）信息。无色散（\(\omega = ck\)）时所有波同步传播，色散时不同波数分离，扩散时纯虚部（驻波）。这种分类对 KdV 方程（Ch 8, 9）、水波（习题 2.5）、光纤孤子（习题 2.4）都成立。
  3. 特征线法 = PDE → ODE：对半线性方程，特征 ODE 天然解耦（\(X\) 不依赖于 \(P\)），因此特征线不相交，经典解唯一存在较长时间。但对拟线性方程，特征 ODE 耦合，特征线可在有限时间内相交（blow-up），此时需要引入激波。激波的形成时间 \(t^*\) 可以通过"两相邻特征线相交的最小时间"严格计算（习题 2.12 给出了公式）。
  4. Rankine-Hugoniot 关系是激波动力学的"运动方程"：\(dx_s/dt = \Delta q / \Delta p\)。对 Burgers 方程，激波速度 \(dx_s/dt = (p^+ + p^-)/2\)（激波速度 = 两侧值的算术平均），这是 Riemann 问题解的核心。
  5. Burgers 方程 \(p_t + p p_x = 0\) 是"流速场 \(p\) 在自身流场中平流"的数学模型：如果把 \(p\) 解释为流体速度，这就是无黏、无压、纯对流的 1D 流体方程。它的"折断"行为就是激波形成（shock formation）的最小例子。Ch 4（无量纲化）和 Ch 5（自相似解）会再次出现。
  6. Boltzmann 输运方程的微观-宏观桥梁：§2.1 的 Euler-Lagrangian 转换本质上是把"个体追踪"（individual-based models）和"场论描述"（continuum）联系起来。Ch 11（Reduced Models）会通过矩方法（method of moments）把 6D 相空间（位置 + 速度）的 Boltzmann 方程约化为 1D 扩散方程——这是更深层的"由 Lagrangian 到 Eulerian"约化。

挑战和开放性问题

  1. 从"弱形式"（积分）到"强形式"（微分）的合法性条件：du Bois-Reymond 引理需要 \(g\) 光滑，但激波解本身就是间断的——这构成一个根本性矛盾。本章的处理是"先做 PDE 推导（合法），再处理激波作为修正（启发式）"。严格处理需要在 BV 函数空间或分布意义下重新定义"守恒律解"，这一理论由 Lions、Lax 等人在 1970-80 年代建立（参考 Lax [X6]）。本章没有触及。
  2. 激波厚度问题：激波在物理上是有限厚度的薄层（如空气中 ~10 nm，金属中 ~1 μm），但无黏 Burgers 方程把它当作几何间断。这一矛盾由黏性 Burgers 方程（习题 2.14，\(p_t + p p_x = \epsilon^2 p_{xx}\)，\(\epsilon \to 0\)）的极限解弥合——该极限正好给出 Rankine-Hugoniot 速度。这暗示激波是一种"奇异极限"（singular limit），物理黏性 \(\epsilon\) 在激波内部不能忽略。这与 Ch 7（边界层理论）密切相关。
  3. 特征线相交后：本章只处理单个激波的形成。多激波（multi-shock）相互作用（\(p^+, p^-\) 都是时间函数）、激波-激波碰撞、激波-稀疏波相互作用是 Riemann 问题的完整图像，本章未涉及。这一理论在 1950-60 年代由 Lax、Oleinik 等人发展。
  4. 守恒律的离散化：特征线法是连续的精确解法，但数值求解（有限差分、有限体积）需要单调格式（monotone schemes）才能保持"熵条件"（entropy condition）——即物理正确的激波。Lax-Wendroff 定理保证单调收敛格式收敛到弱解，但可能收敛到非物理弱解（如膨胀激波）。熵条件（Lax 1967）从一族弱解中筛出唯一的物理解。
  5. §2.4 末尾的"双曲型"定义不完整：作者说"如果 \(M^T\) 的所有特征值为实数且有完整特征向量组"是双曲型，但没给出 "强双曲" vs. "弱双曲" 的区分（取决于特征向量是否构成完备基）。这一区分对线性稳定性和数值稳定性都至关重要。
  6. 特征线法对 2D 问题的扩展：习题 2.8（2D 放射性物质扩散）和习题 2.15（一般一阶 PDE 的 Lagrange-Charpit 方法）展示了特征线法在 2D 情形下的工作方式。但对多维守恒律（如 2D 浅水方程、3D 欧拉方程），特征线变成特征锥（characteristic cone）或双特征曲面（bi-characteristic surface），几何分析更复杂。本章未触及。
  7. "被动" vs. "主动"输运：作者在 §2.1 区分了被动输运（粒子的存在不影响流场）和主动输运（如 Ch 1 末尾提到的"structured population models"——粒子的运动依赖于密度）。后者需要耦合 ODE/PDE 系统，本章只处理被动情形。Ch 11（化学反应-扩散-对流耦合）会再次出现这一区分。

个人反思与批判性分析

本章是全书从 ODE（PDE-free）到 PDE 的桥梁，处理得很扎实。Reynolds 输运定理的引入非常自然——先讲 Euler/Langrangian 的等价（§2.1），再讲守恒律的"blob-wise 表述 → pointwise PDE"（§2.2），最后讲如何求解（§2.5–2.6）。这种"建模 → 求解"的递进是本书 Part I-IV 反复出现的节奏。但从数学严谨性角度看，本章有几个可议之处：

  1. du Bois-Reymond 引理的证明完全省略。作者在 §2.2 直接陈述"若积分对所有子域为零，则被积函数为零"，没有给出任何证明或动机。这对读者的"知识信任"要求较高——读者必须相信这是一个"显然的"分析引理（实际上其证明需要变分法基本引理和一些紧性论证）。一个简单的反证（如果 \(g\) 在某点 \(x_0\) 不为零，则在 \(x_0\) 附近取 \(g\) 主导的小区间，积分不为零）会让这引理更"可触"。建议在脚注中加一行证明梗概。

  2. Rankine-Hugoniot 关系的推导"工程化"。作者在 §2.6 通过 Leibniz 法则 + 重组 + 减项得到 (2.48)，这是一种"凑"的过程——读者很难第一遍就看懂。更"几何化"的推导方式：(a) 把激波视作移动边界，(b) 对激波两侧的控制体分别应用守恒律，(c) 减去得到"质量守恒"对激波速度的约束。这种推导更符合直觉，也更容易记忆。作者也可以在脚注给出等价推导。

  3. 激波物理与数值方法的"分割"。本章完全没有提到数值激波捕捉（shock capturing）——这是 CFD 教科书（如 LeVeque [X4]）的核心主题。读者如果只读本书，会形成"激波 = 解析构造"的印象，但实际工程计算中激波 = 数值伪振荡的斗争是 CFD 60 年的主线。建议在 §2.6 末尾或 §2.7 提一句"激波的数值处理是另一主题（参考 LeVeque）"，作为衔接。

  4. 色散关系只在 §2.3 简短出现，但重要性贯穿全书。色散关系 \(\omega = \omega(k)\) 在 KdV 方程（Ch 8, 9）、水波（习题 2.5）、孤子理论（习题 2.4, BBM 方程）、量子力学（de Broglie 关系）、光纤通信（非线性 Schrödinger）都出现。本章用 1 页（§2.3 + §2.4 末尾）介绍这一概念，密度不够。建议给出色散关系的几何解释（\(\omega\) vs. \(k\) 图）、群速度的几何构造（\(\omega(k)\) 在 \((k, \omega)\) 平面上的切线斜率），这样读者能在 Ch 8 看到 KdV 时不至于完全陌生。

  5. Euler 描述 vs. Lagrangian 描述的"等价性"是有条件的。作者在 §2.1 强调"两种描述等价"，但严格来说只有当流场是光滑且不发生特征相交时才等价。一旦发生激波或湍流涡量场中的"粒子轨迹纠缠"（stochastic particle paths），Lagrangian 描述就失去严格意义——这是湍流理论的起点（Onsager 猜想，2010 年代由 Isett 等人证明）。本章的措辞在物理上无害（只处理被动粒子和平滑流），但对湍流学生而言是误导。

  6. §2.6 末的"等面积律"与 Ch 5 自相似解的衔接。作者在 §2.6 末尾说"Ch 5 会用相似解重新得到 (2.50)"——这是一个重要的内部预告。Ch 5 引入自相似变量 \(\eta = x/t\) 后，Burgers 方程的激波解（等面积律）会自然出现，但那时是从相似变换的角度重新推演，视角不同。读者如果跳过 Ch 5 直接做习题 2.13（不同初值下的 Burgers 激波），会错失这种"两种方法殊途同归"的洞察。

  7. 缺失主题：交通流模型（Lighthill-Whitham-Richards 模型）。§2.6 末尾的"激波 = 移动交通堵塞" 是个完美的伏笔，但作者没有展开交通流模型 \(h_t + (h v(h))_x = 0\)，其中 \(h\) 是车流密度，\(v(h)\) 是密度依赖的速度。这是宏观交通流理论的核心方程，也是输运方程在社会科学中的成功应用。如果作者在 §2.7 "Further Directions"中给一段引文（如 Whitham [X8] 的 Chapter 11），会很有用。

重要参考文献

  • [X1] Rutherford Aris, Vectors, Tensors, and the Basic Equations of Fluid Mechanics, Dover, 1989. (Reynolds 输运定理的原始教材)
  • [X2] George B. Whitham, Linear and Nonlinear Waves, Wiley, 1974. (色散关系、行波、激波理论的经典教材)
  • [X3] Joel Smoller, Shock Waves and Reaction-Diffusion Equations, Springer, 1994. (激波理论严格数学处理，包括熵条件)
  • [X4] Randall J. LeVeque, Numerical Methods for Conservation Laws, Birkhäuser, 1992. (激波数值捕捉、有限体积方法)
  • [X5] Peter D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM, 1973. (双曲守恒律严格理论)
  • [X6] Peter D. Lax, "Hyperbolic Conservation Laws in Continuum Physics," Bulletin of the AMS 46 (2009) 1–10. (Lax 对自己工作的综述)
  • [X7] M.J. Lighthill and G.B. Whitham, "On Kinematic Waves. II. A Theory of Traffic Flow on Long Crowded Roads," Proceedings of the Royal Society A 229 (1955) 317–345. (LWR 交通流模型原始论文)
  • [X8] Constantine M. Dafermos, Hyperbolic Conservation Laws in Continuum Physics, Springer, 2016. (现代严格处理，包括熵-弱解唯一性)
  • [X9] Lawrence C. Evans, Partial Differential Equations, AMS, 2010. (du Bois-Reymond 引理和 PDE 弱解的严格数学背景)
  • [X10] Michael E. Gurtin, An Introduction to Continuum Mechanics, Academic Press, 1981. (连续介质力学的标准教材，覆盖 Cauchy 动量方程)