差分方程_PDE有限差分方法(11)——常系数对流方程中的常见格式

从今天开始进入双曲型方程的差分方法构造. 主要介绍几个常见的常系数对流方程的构造方法以及CFL方法。下一次我们来复习常系数对流方程差分格式的其他研究手段。

与抛物型方程相比, 双曲线方程缺乏耗散机制, 相应数值困难更为突出.

参考书：
(1) J.W. Thomas - Numerical Partial Differential Equations_ Finite Difference Methods (1995, Springer)
(2) 张强《偏微分方程的有限差分方法》科学出版社，2019年1月版。
(3) K. W. Morton, D. F. Mayers - Numerical solution of partial differential equations (2005, Cambridge University Press)

今天内容：

1 对流方程的特征线

考虑最简单的对流方程

的纯初值问题或周期边值问题, 初值为

这个方程的真解可以用特征线理论求出, 在每条特征线上u都是常值.

特征线就是下面ODE的所有解:

在每条特征线上, 解

满足

利用初值, 可以选择一系列合适的点

并找到通过

的特征线上的数值解, 这个方法就是特征线法.

如果

关于

是Lipshitz函数且关于

连续, 则特征线不会相交.

当

是常数的时候, 由特征线理论, 真解为

特征线为一束平行直线

所以原问题是个行波解, 双曲线方程的解可以由古典解拓展到间断函数(弱解). 所以我们的差分方法应该尽量做到更好描述间断界面.

对于非线性问题, 如果

特征线依然都是直线, 因为

在每条直线上都是常数, 但是这些特征线不平行. 此时, 解可以写成

但是此时特征线可能会相交.

2 基本格式、CFL条件

设

是等距时空网格, 其中

分别是时间步长与空间步长, 双曲线方程的网比是