1.时间推进法简介.docVIP

1．时间推进法简介 对于定常流动，假定是二维的，那么它的控制方程为： 但是，可以证明，上述方程在超音流和亚音流中呈现不同的特征，难以用统一的算法求解。因此，即使模拟定常流场，我们也是求解非定常Euler方程 （7-1） 在求解过程中，检验，如果满足 （为小量） 就认为，就是原方程的数值解。 这时，方程（7-1）实际上退化为定常方程。 2．关于N-S方程的数值求解 仍以二维流动为例，N-S方程可写成（无量纲形式） （7-2） 由于，我们遇到的大部分流动都是高雷诺数的（），N-S方程中起重要作用的是对流项（输运项）F、G。所以用于求解Euler方程的大多数格式，只需加上对粘性应力和导热项（R、S）的处理就可用于求解N-S方程，并且这种处理是很简单的，因为粘性应力和热传导均为二阶项，具有耗散特性，只需对它们采用中心差分即可。 一． 一维Euler方程 式中 ， 补充关系式 MC格式可写成： （7-3）预测 （7-3） 校正 1．精度分析 将预测步代入校正步，得 整理得： （7-4） 利用Taylor展开式，得 （7-5） （7-6） 将式（7-5）和（7-6）代入式（7-4） 整理得： （7-7） 由Euler方程得：