双曲型方程的差分法.doc

双曲型方程的有限差分法 §0 预备知识 0.1双曲型方程的常见类型： （1）、一阶线性双曲型方程 （2）、一阶常系数线性双曲型方程组 其中u为未知函数向量，A为p阶常数方阵。 （3）、二阶线性双曲型方程（波动方程） 一维 a(x)为正值函数 二维 三维 （4）、对流扩散方程 等等。 这些方程的定解条件，可以是仅有初始条件，也可以是初始条件与边界条件的混合。如对波动方程（一维），有 （1）、初值问题 （2）、混合问题 第一类： 第二类：边界条件改为： 第三类：边界条件改为： 波动方程及其特征 线性双曲型方程的最简模型：波动方程初值问题 （1） 下面讨论它的特征和解析解。 由二阶偏微分方程理论，上述方程的特征方程为 或 . 分解为： ， 由此可定出两个方向，称为特征方向。 解得： 是XT平面上的两族直线，称为特征。用u沿特征的偏导数来表示u沿x，t 同理，有 由方程（1），得 由初值条件， 得 解此常微分方程，可得解：（称为Dalembert公式） 由此公式可看出：u在点的值仅依赖于初始函数在区间上的值。（此区间就称为点的依存域） 过作第一特征（斜率为正），过作第二特征（斜率为负），相交于，所得三角形域为此依存域的决定域，依存域上的初值决定了三角形域内u的值。反之，由作两条特征，与x轴（）相交截出的闭区间就是依存域。如下图。 依存域 ================================================================== 从解的表达式还可看出，对x轴上任一点，它所能影响到的所有的集合是以为顶点，两条特征为边的角形域，称为的影响域。如上图。 双曲型方程这种对初值的局部依赖关系和特征关系是其他方程所没有的。因此初值函数的一些性质（如间断等）也会沿特征线传播，从而使解不具有光滑性。在构造双曲型方程的差分格式时，应考虑这些特性。 §1 一阶线性双曲型方程组 一般形式：n个未知函数 n个方程： 其中 ，，都是已知光滑函数。 称为一阶线性偏微分方程组。 记，， 则有 记B可逆。不是一般性，上面方程又可化为 （2） 若A有n个实的互异特征值，对某个，称（2）是在点的（狭义）双曲型方程组。若对G中每点都成立，称（2）是在G的双曲型方程组。 注：更一般地，若A的n个实特征值固有n个线性无关的特征向量时（或A可对角化时），称（2）是双曲型方程组。 若A和C与u有关，称（2）为拟线性偏微分方程（双曲型）。 若A与u无关，C与u有关，称为半线性偏微分方程（双曲型）。若A与u无关，C与u无关或线性地依赖于u，称（2）为线性偏微分方程组。 1.1一阶线性双曲型方程初值问题 简单介绍。现对一个方程，但a可正可负。 （3） 由于，所以沿特征线的变化率为0，即沿特征线的值为常数。过任意一点的特征线 与的交点为，由于沿此特征线的值为常数，故有 由点的任意性，可知初值问题（3）的解为 迎风格式（Upwind Scheme） 由气体力学的含义，a(x)表示气流速度，故得此名。 现考虑一个方程， 设a可正可负。 设分别为时间和空间步长。对空间偏导数用不同的离散化方法（向后、向前、中心差分），可得： 1） 称左偏心格式 2） 称右偏心格式 3） 称中心差分格式 1），2）的误差阶为，3）的误差阶为。 收敛性：设。如（P.127）时，格式1）和3）： 时，满足Courant条件； 格式2）不能满足Courant条件； 当时，格式2）和3）：时，满足Courant条件； 格式1）不能满足Courant条件。 稳定性：由分离变量法，可得结论如下： 格式1）稳定的充要条件：当时，； 格式1）只能用于的情况。 格式2）稳定的充要条件：当时，； 格式2）只能用于的情况 格式3）绝对不稳定，不能用。 2. Lax-Friedrichs格式 在上述格式3）中如用代替，即得Lax-Friedrichs格式（也简称为Lax格式）： 误差阶为。收敛的Courant条件和稳定的充要条件均为。 利用特征线构造差分格式 a) Lax-Wendroff格式 利用特征线来构造（构造过程略）。 . 截断误差阶为。稳定的充要条件为。 利用特征线还可以构造一些其他的格式，如 b) Beam-Warming格式 它的稳定的充要条件为，从而可以加大时间步长,减少计算工作量。 跳蛙格式(Leap-frog Scheme