(完整版)4-迎风格式.docVIP

PAGE1

PAGE14

§4.迎风格式

4.1迎风原则

先考虑线化的模型方程-对流方程

其中为常数。

空间导数用向后差分近似

空间导数用向前差分近似

向后差分格式

向前差分格式

迎风原则的示意图

4.2迎风格式

向后差分格式和向前差分格式都可以看成是对中心差分格式的修改

（用于）

（用于）

不仅如此，经过这样的改写，我们发现它们还可以统一起来，写成

如果再经过进一步的改写

并重新组合格式中的各项，就得到

记，，则上述格式可写成

注意到，无论还是，都有

，

而且

于是我们可以对上面得到的格式做出新的解释：

将分解成了“正”和“负”两个部分

，，

对流项也相应地分裂成了两项。然后根据迎风原则，对这两项分别用向后差分和向前差分来近似。这样得到的差分格式

称为迎风格式。

4.3泰勒展开法

显式格式要计算的是网格点处的精确解，将它按自变量做泰勒展开，可将它与时刻的解联系起来，

但是由原方程，有

代入上面的泰勒展开式，得

如果用中心差分作为一阶和二阶空间导数的差分近似，即

就得到

它也可以写成

这就是Lax-Wendroff格式，具有二阶精度。

上面介绍的推导方法，称为泰勒展开法。这种方法利用了精确解的泰勒展开，还利用了原方程本身，但是并没有用到方程的双曲型特性，所以适用于所有非定常问题（含一阶时间导数的方程），当然也可以用来推导迎风格式。

4.4二阶迎风格式

下面用泰勒展开法推导迎风格式。先假设对流方程中的。与前面完全一样，将精确解按自变量做泰勒展开

由原方程，

代入上面的泰勒展开式，得

将上式截断到一阶导数项，就是

如果用向后差分作为一阶空间导数的近似，就得到

即

这就是后差格式，即时的一阶迎风格式。

如果截断到二阶导数项，就是

此时，一阶和二阶空间导数都用后差近似，即

就得到差分格式

即

这就是二阶迎风格式，文献中也称为Warming-Beam格式。它还可以写成

为了和一阶迎风格式进行对比，将上式改写成

引入“开关”变量，就可以将一阶和二阶的迎风格式统一写成

显然，就是一阶迎风格式，就是二阶迎风格式。

当时，我们可以给出与上面平行的一系列结果。此时只需利用

就得到差分格式

再利用开关函数将它和一阶格式写在一起，就是

为了将两种情形进一步统一，还是需要把分解成

（，）

然后根据迎风原则，就可以写出

4.5守恒型迎风格式

对于守恒律方程组

迎风格式也可以推广过来，方法仍然是泰勒展开法。

将按自变量做泰勒展开

对于守恒律方程，记，则由原方程，有

代入上面的泰勒展开式，得

此时，按照守恒格式的要求，需直接对通量的空间导数使用差分近似。这里暂且不论是否符合迎风原则，就用向后差分，可得

式中的，。

注意到，所以有

于是得到具有二阶精度的迎风格式

当然也可以将它与一阶迎风格式统一起来，就是

为了运用迎风原则，必须设法将通量分解成，这种处理手段就是通量分裂。分裂之后，对使用向后差分，对使用向前差分，并记，，就有

上式中，对以及，省略了代表时间步的上标。同时，为了更大的灵活性，我们允许开关函数随网格变化，不再统一取值。

定义数值通量

可将迎风格式写成

在上述格式中，若将所有开关函数关闭，就成为

这就相当于通量分裂之后，将守恒律方程组写成

然后根据守恒原则，直接对通量进行差分近似。并根据迎风原则，对使用向后差分，对使用向前差分。

进一步定义总的数值通量

，

就能够把迎风格式写成守恒格式的一般形式

1

一月

二月

三月

产品名称

数量

金额

利润

产品名称

数量

金额

利润

产品名称

数量

金额

利润

合计

合计

合计

四月

五月

六月

产品名称

数量

金额

利润

产品名称

数量

金额

利润

产品名称

数量

金额

利润

合计

合计

合计

下午13：00—17：00

B．实行不定时工作制的员工，在保证完成甲方工作任务情况下，经公司同意,可自行安排工作和休息时间。

3．1．2打卡制度

3.1.2.1公司实行上、下班指纹录入打卡制度。全体员工都必须自觉遵守工作时间，实行不定时工作制的员工不必打卡。

3.1.2.2打卡次数：一日两次，即早上上班打卡一次，下午下班打卡一次。

3.1.2.3打卡时间：打卡时间为上班到岗时间和下班离岗时间；

3.1.2.4因公外出不能打卡：因公外出不能打卡应填写《外勤登记表》,注明外出日期、事由、外勤起止时间。因公外出需事先申请，如因特殊情况不能事先申请，应在事毕到岗当日完成申请、审批手续，否则按旷工处理。因停电、卡钟（工卡）故障未打卡的员工，上班前、下班后