N-S方程求解技术报告.ppt

* * 流过平板的超音速流动 物理模型 超音速平板绕流的类型 10.3.3 有限差分方程 向量形式的控制方程可以表示成： 利用泰勒级数展开，每个网点（i，j）处的流动变量都可以 推进到下一个时间步 其中，U还是表示流场变量，而且在时刻t是已知的，我们定义： 为了求 ，使求解能够推进，采用下列步骤： 由t时刻的已知的流动，用向前空间差分计算控制方程的右 端项，从而得到 的值 用步骤1）的结果，可以得到（t+△t）时刻流场变量的预 估值 合并步骤1），2），预估值的表达式为： 3）将上预估值代入控制方程右端并利用向后空间差分计算， 可以得到预估值的时间导数 4）将上得到的时间导数代入方程（10-6），得到（t+△t）时 刻U的校正量（具有二阶精度），合并步骤3），4），有： 重复1）至4 ）步骤，直到流场变量达到一个定常值，就是所求 的定常解。 为了保持二阶精度，E中关于x的导数项用于 方向相 反的差分离散，而关于y的导数项用中心差分近似。同样的， F中关于y的导数项用于 方向相反的差分离散，而F 中关于x的导数项用中心差分近似。 每个预估步和校正步结束后，由向量U可以得到原始变量 当