8TVD格式.doc

§8. TVD格式 8.1 守恒格式 将上一节的统一格式改写成 注意到 Lax格式 迎风格式 Lax-Wendroff格式 于是可以定义函数 Lax格式 迎风格式 Lax-Wendroff格式 （注）通量差分格式里介绍过的熵修正函数，可以看成是函数 的第四种形式，所以当时也用了相同的记号。 利用函数 ，统一格式可写成 如果再记 ，则 将格式写成这种形式，是为了便于推广到非线性的情形。 对于非线性标量方程， 是变量，定义数值通量 式中 就有守恒格式 8.2 显式TVD格式 充分条件一：显式格式 具有TVD性质的充分条件是 ， ， 【证明】根据全变差的定义， 如果判据成立，则上式中各项系数都是非负的，可以从绝对值里面提出来，得到， 注意到，对于无限求和，求和指标可以调整，就有 所以显式格式是TVD的。 将上面给出的判据应用于守恒格式 利用 的表达式，就有 由此可以看出（注意 的下标要从 调整为 ）， ， 所以，函数 ， 于是刚才给出的判据成为 ， 即 这就是显式TVD格式的判据。这一结果表明： 显式TVD格式的下限就迎风格式，具有最少的格式粘性。 Lax-wendroff 格式没有落在上述范围之内，所以它不是TVD格式。 熵修正函数满足显式TVD格式的判据。 因此，使用熵修正函数不会破坏格式的TVD性质。 在显式TVD格式的判据中， 、（也就是 ）是实质性的，而 （即 ）只不过是保证显式格式稳定性的CFL条件。 8.3 隐式TVD格式 充分条件二：对于 ，隐式格式 具有TVD性质的充分条件是 ， ， 【证明】将隐式格式改写成 所以 即 将上式取绝对值并无限求和，就是 将上式的左右两端分别记作 和 。先来看左边。利用不等式 有 如果 ， ，则因为 ，上式中各项系数都是非负的。类似于显式格式的证明过程，可以推出， 另一方面，在判据成立的条件下，可以推出 ， ， 于是，右端 中各项系数也都是非负的。还是类似于显式格式，有 （注）当 时，直接有 。 综合起来，就得到 所以隐式格式是TVD的。 将守恒格式推广成隐式格式 其中数值通量 的表达式与 类似，只是将上标换成了。 这一隐式格式具有TVD性质的充分条件就是 （ ） 对于显式格式（），上式就是前面给出的条件 对于完全隐式格式（），函数 的值没有上限，只需 对于平均隐式格式（），就是 8.5 TVD格式的设计 上一节关于格式粘性的讨论，是针对含激波的流场，而且仅限于激波附近。对于光滑的流场，以及含激波流场中远离激波的区域，我们还是希望数值格式能够有较高的精度。 为此，重点考察具有二阶精度、不带格式粘性的Lax-Wendroff格式 和格式粘性最少的迎风格式 它们的右端项之差为 所以，在Lax-Wendroff格式右端加上这样一项，就能得到迎风格式。由于这一项的作用类似于格式粘性，但不是格式本身就有的，而是人为加上的，所以称为人工粘性（人为粘性）。 反过来，如果从迎风格式右端减去这样一项，就能得到Lax-Wendroff格式。这种做法相当于给格式加上了一个负的人工粘性，故称为反扩散方法。 首先能想到的，就是引入一个开关函数，在这两个格式之间来回切换。在激波附近打开开关，或将开关调到一个合适的值，使格式具有一定的格式粘性，对抑制振荡来说刚好够用。而在远离激波的区域，关掉开关，或将开关调成一个小量，又可以使格式粘性（几乎）消失，把格式的精度恢复到二阶或接近二阶。这就是构造TVD格式的基本思路。为了实现这一想法，通常的做法是把限制器函数作为系数与人工粘性项相乘，或者说是用限制器来代替开关函数。 下面就以标量方程为例，介绍几种显式TVD格式，然后再推广到守恒律方程组。这些格式都是守恒形式的， 因此只需给出格式的数值通量 （将下标减一，就可得到 ）。 为简洁起见，下面省略代表时间的上标 。 8.6 通量修正 — Harten的TVD格式 Harten是TVD格式的开创者，他给出了第一个具有TVD性质的二阶格式，当然是非线性格式。尽管后来的研究表明，TVD格式只有二阶精度，而且并不是处处都能有二阶精度，在极值点附近，格式的精度会降为一阶的。 Harten从一阶精度的迎风格式入手，其第一修正方程是 记 ，则 ，而上式成为 Harten的想法是：既然已经知道会产生格式粘性 ，就事先打好埋伏，将迎风格式应用于方程 此时的第一修正方程就应该是 这里， 对应着 的格式粘性，而 则是与 对应的格式粘性，并且 。于是就有 这表明，对于原方程来说，这一做法将得到一个二阶格式。而且，迎风格式是TVD的，因此这个格式也应该是TVD的。 Harten将通量 修改成 ，因此称为通量修正法。由