迎风型weno格式 - read.docVIP

迎风型WENO格式 考虑一维标量双曲守恒型方程s。它的半离散守恒型格式如下 (1) 其中L是空间离散算子。 下面我们将通量进行矢通量分裂为： (2) 在此的分裂规则是：，可以采用多种分裂形式，例如采用简单的Lax-Friedrichs分裂格式： 则有 (3) 对于时间项，我们可以使用Runge-Kutta法来提高精度。三阶TVD型Runge-Kutta法（简写为RK-3）为： 四阶非TVD型Runge-Kutta法（简写为RK-4）为： 下面我们假定，讨论的计算，并且略去上标“+”；当时，与关于对称。 简要介绍一下迎风型WENO格式的重构，详细过程请参考Jiang和Shu的论文。其基本思想通过线性组合低阶通量得到高阶通量，这些系数叫做线性权重。WENO格式是基于ENO格式构造出来的，有关ENO格式的构造这里就不再赘述，可以参考Shu和Osher的文章。根据r阶精度的ENO格式的选取模板（即插值区域）思想，考虑迎风效应，选取初始模板为，经过r-1次扩充，可能的模板有r个： r阶ENO格式是从上述模板中选取一个最光滑的模板，在其上构造插值多项式来逼近f在附近的数值通量，其数学表达式为： (4) 其中是一个插值多项式，它满足下面的关系式： (5) a2l为已知常数，可以从Shu和Osher的文章查到，例如：。 现以三阶的ENO格式为例说明之。此时被选模板有三个： 它们对应的数值通量分别为：。首先我们来确定，根据式(4)可以得到 其中是在上构造的二次插值多项式，用Newton插值方式表示如下： 于是我们可以得到 所以 (6) 同理可得和 所以 (7) 所以 (8) ENO的思想就是在若干个模板里选出一个最光滑的，并由它求出数值通量。在包含间断的模板里，我们非常希望采用这种方法，因为它能排除包含间断模板中的不是很精确甚至完全不精确的信息。然而在光滑的模板里面，这种做法却不是我们所希望的，因为此时所有可能的插值模板都提供着同样精确的信息，从而浪费了其他r－1个模板。事实上对于r阶ENO格式而言，有r个被选模板，如果将它们合并在一起将形成一个包含2r－1个格点的插值区域，按照式子(4)可以给出网格界面处2r－1阶的近似： (9) 下面，我们依旧以三阶ENO为例，它的三个模板一共包含五个点，所以可以构造一个四次插值多项式： 则数值通量函数可以构造为： 求得界面处的数值通量： (10) 它正好是五阶上游中心格式的数值通量，但是它不能排除包含间断的模板所带来的不精确信息，用它所计算的通量在间断附近会产生非物理的振荡。 为此，构造WENO格式以解决此问题。WENO格式的数值通量表达式为： (11) 其中由(4)确定，是赋予相应模板的权重，问题的关键是如何选择。为了获得基本无振荡的性质，应该使与相应模板的光滑性联系起来，使得包含间断的模板被赋予几乎为零的权重。同时要求在光滑的模板中权重的分布将使得式(11)与式(9)的通量相接近。其实可以组合r－1阶数值通量来求得2r－1阶数值通量： (12) 通过简单的代数运算可以确定这些系数，并且它们满足。 对于三阶ENO格式而言，把式(6)(7)(8)和(10)带入(12)就可以确定这些系数，它们的值为： (13) 由式(11)和(12)等号两边相加，可得 因为，如果要求，那么上式就可以改写为： (14) 上式右端第一项是的2r-1阶近似，其精度为。因为式(14)第二项的求和符号内的两项的乘积的第二项的精度是，所以如果再要求权重系数满足，那么式(14)右端第二项的精度在光滑区域将达到。也就是说如果的选择符合下述三个条件，那么式(11)的精度在光滑区域将有；同时此数值通量将有ENO基本无振荡的性质，此时称为最佳权重。 ① 组合系数符合； ② 在光滑区域内组合系数满足； ③ 任何包含间断的模板都被赋予几乎为零的权重，即。 Jiang和Shu选取如下符合上面三个条件的系数： (15) 其中，是一个小量，一般取10-6，用来避免分母为0；p是大于等于2的正整数，它的作用是放大模板的不光滑性，使得包含间断的模板的权重更小；ISk是第k个模板的光滑性度量。Jiang和Shu采用如下方法度量模板的光滑性： (16) 对于三阶ENO而言，将前面的带入上式，有 (17) 当r=3时，我们验证一下是否符合最佳权重条件。首先，由上述各式构造的符合。其次，分析的精度，将式(17)进行Taylor展开，我们可以得到： 如果，有；将进行Taylor展开，有，把它带入的表达式，并利用Tayl