第六讲TVD格式简介.docVIP

第六讲 TVD格式简介 （一）保单调差分格式与单调差分格式 设迁移方程 （11-30a） 现考虑纯初值问题，设有初值 （11-30b） 则其解析解为 （11-31） 从方程（11-30a）的精确解（11-31）中可以看出，如果初始函数是的单调函数，则时的解也是的单调函数。于是自然想到，对于差分格式如果初值是单调的非增函数（或非减），希望差分方程的解也能保持非增（或非减）。 对于非线性迁移方程 （11-32a） （11-32b） 定义：求解（11-32）的点守恒型格式的一般形式可写为 式中，。 如果是单调函数而也是单调函数，则称该格式为保单调的。 如果满足关系式 ， （11-33） 则称该格式为单调差分格式。 注意单调差分格式与保单调差分格式不同： 在定义保单调差分格式时只谈到了在时间层上的量与时间层上量之间的关系，而未直接涉及到格式本身结构，在定义单调差分格式时涉及到了格式的结构，要求式（11-33）得到满足。 Godunov证明了如下定理。 定理：差分格式 （11-34） 是单调差分格式的充要条件是所有。Godunov在这里讨论的是线性方程的差分格式。对线性情况单调差分格式和保单调差分格式是一致的，即保单调格式也是单调格式。Godunov还证明了，常系数单调差分格式的截断误差是一阶的。在一般情况下（非线性方程）单调差分格式一定是保单调的，但保单调差分格式不一定是单调格式。 考察迁移方程 （） （1）Euler后差格式 （） 可写成 可见，Euler后差格式是单调格式（满足稳定性条件）； （2）Lax-Wendroff格式 可写成： 可见Lax-Wendroff格式是非单调格式； 格式的单调性质对于避免在激波附近产生非物理振荡具有决定作用。 需要指出的是真正严格意义上的单调格式只能是一阶精度格式。 （二）TVD格式与单调格式的关系 具有下列特性的格式被称为TVD（Total Variation Diminishing，总变差递减）格式或TVNI（Total Variation Nonincreasing，总变差非增）格式。 1．对标量非线性方程及常系数双曲型方程组，格式的解是总变差递减的； 2．与守恒律和熵不等式是相容的。 对于偏微分方程数值解，总变差定义为 （11-35） 下面来讨论非线性迁移方程 （11-32a） （11-32b） 解的总变差变化情况。该方程的解沿特征线（其斜率为） 是常数，因为此时，但却是变量。 设从时刻到时刻特征线没有相交（如图所示）。在时刻有任意划分，可在时刻找到相应的划分，使得，。由此可知，在特征线相交之前方程（11-32）的解对于的总变差为常数。 现假设特征线和于时刻相交（如图所示），此时时刻点之左方的总变差与时刻点之左方的总变差相等，和点之右方的总变差相等，而 因而当有激波时总变差是非增的。 由上面简单讨论中可知，方程（11-32）的解随时间的推移总变差是非增的。它一定满足以下两点： （1）方程解在随时间推移的增长中不可能产生新的局部极值点； （2）解的局部最小值在时间推移过程中不再减小，而局部最大值不再增加。 因此，如果在初始时刻解是单调的，在以后的时刻将保持这种单调性。可见，总变差递减与格式的单调性相关。 Harten证明了下述结论： 1．任何单调格式都是TVD格式； 2．任何TVD格式都是保单调格式。 由于任何线性的保单调格式必为单调格式，而单调格式只有一阶精度，因此高精度的TVD格式就只能是非线性的保单调格式。 （三）TVD格式的充分条件 以常用的三点（）显式格式为例。对于非线性迁移方程（11-32a） 显式守恒型差分格式可写为 （11-37） 式中，，为数值通量，它们的不同取值代表了不同的显式差分格式。例如 取=，式（11-37）就是Euler后差格式； 取 式（11-37）就是Lax格式。 为了判断格式是否TVD（总变差递减），一般要将守恒型差分格式写成增量形式，即将式（11-37）写成： （11-38） 式中： ； ； Harden证明，只要满足下列关系（设） ， ， （11-39） 三点显式差分格式（11-38）一定是TVD的。 证明如下： 在点上列出格式（11-38） （11-40） 式中 令（11-40）-（11-38），得到 整理得： 两边取绝对值，并考虑关系式（11-39），得到 （11-41） 在式（1