特征线法
一阶偏微分方程看起来吓人,因为未知量同时依赖好几个变量,但有一个优美的招法能驯服它:找出穿过区域的那些特殊曲线,偏微分方程沿之坍缩成一个常微分方程。沿其中一条曲线前行,方程就只是告诉你解随之如何变化。这些特权路径就是特征曲线,沿它们行进就是特征线法。
取一个一阶方程,如 a(x, y) u_x + b(x, y) u_y = c(x, y, u)。左边是 u 沿向量 (a, b) 的方向导数,所以沿任一切向为 (a, b) 的曲线——即满足 dx/dt = a、dy/dt = b 的曲线——u 的变化由 du/dt = c 支配。单个偏微分方程已变成一组常微分方程:在 (x, y) 平面里解出特征曲线,再沿每一条积分 du/dt = c,把解从数据携带向前。从几何上看,解曲面由这些特征曲线编织而成,每一条都从初始曲线的一点出发、携带它的值。对最简单的情形——输运方程 u_t + c u_x = 0——特征是直线 x - c t = 常数,u 沿每条简单地保持常数——一个以速度 c 向右刚性平移的轮廓。
特征线法是一阶偏微分方程以及一般双曲问题的天然工具:交通流、气体动力学、几何光学,以及任何守恒律。它最诚实的一课是非线性方程会发生什么:当特征速度依赖于 u 本身时,特征会相交,解曲面折叠,光滑解不复存在——激波形成。这个方法不仅解简单情形,它还精确预言经典解必将在何处、何时崩溃,这是任何其他初等方法都给不了你的信息。
对输运方程 u_t + c u_x = 0、初值 u(x, 0) = f(x),特征是 x - c t = 常数,u 沿每条保持常数,于是 u(x, t) = f(x - c t):初始轮廓只是以速度 c 向右滑动,形状从不改变。
沿每条特征,偏微分方程都是平凡的;把这些直线缝在一起就给出行波解。
对非线性一阶偏微分方程,特征可能相交;相交之处不存在单值的光滑解,激波形成——这个方法诚实地标示出崩溃,而非隐藏它。