3交错网格紧致差分格式-于欣博士.docVIP

不可压缩流体力学问题关于压力计算的 一个改进算法 于欣 （ 高温气体动力学国家重点实验室,中国科学院力学研究所，北京海淀区 100190） 其中 （1.5） 其中V=(u,v)，设 则(2.1)可以写为 以下近似是四阶的 其中 对离散格式用同样的符号则(2.3)(2.4)为连续方程(1.2)的四阶离散格式。 边界上的三阶离散格式为（2.3）with （j=1，2，…,M） 类似于(2.6)可以给出 3 交错网格紧致差分格式 考虑二维非定常粘性不可压缩Navier-Stokes方程(1.1)(1.2)中的 我们考虑(1.1)(1.2)(3.1)(3.2)的差分格式 其中 ， （1）中的一阶导数的差分离散为 （2）中的二阶导数的差分离散为 （3）压力的一阶导数的差分离散 （4）连续方程中速度一阶导数的差分离散 （本文主要采用另外的离散格式，见第二部分） （5）非线性项中的里的用以下插值公式 其中二阶导数的近似可用(3.6)(3.12)(3.13)的结果。 下面给出边界上的差分逼近。我们主要给出x=0上的差分方程： 其中: 其中 其中 其中可以由(3.5)1,(3.10)和 得到。 以上差分格式中的u,v,p的右上角的撇表示导数对应的差分方程中的未知量。 4 非定常问题时间方向离散的Runger-Kutta方法 这样f 就是V的函数。 求解的四阶Ruger-Kutta方法为： 其中: 边界上常用 （4.2）在中间层（n+0.5时间层）精度只有一阶。我们给出如下改进，精度比(4.2)提高一阶 5 方腔问题数值计算 考虑在顶部受剪切力驱动的二维方腔流动。控制方程为(1.1)(1.2)(3.1)(3.2),计算区域为0x1, 0y1,边界上的速度： 上边界(y=1)上速度u=1,v=0,其他边界上速度均为零。 FVC格式连续方程用第二部分的离散公式。 [19]中的数值结果证实了FVC格式的高精度。 取Re=100, N=M=16. 图1是用N-S方程计算到t=T（T=0.25, 0.5625, 1.1875, 2, 3, 5, 8）时压力置零(初始化)且对应的非线性部分不随t改变（用已经算出的速度计算非线性项）继续计算的压力的误差曲线，Re=100, N=M=16。压力的误差ErrP是压力与网格增加（x，y方向均加倍）后计算出的压力在时间较大时的稳定结果的差（所有点差的平方和 除以压力的网格点数 再开方。网格增加后（x，y方向均加倍）到原来网格上压力插值计算公式精度：内部四阶，边界上3阶）。压力计算结果均在最后加一个常数使得平均值为零。 图1中t=T（1.5）时误差的一个‘尖’是由于压力置零导致的。压力置零是为了模拟在 t = T 时刻开始计算压力时不知道压力值（只知道速度值）的情况。 我们也试算了非线性项一部分不变的情况，即中只是前面的对应的项不变，见图2。从图中可以看出，计算结果比非线性项全部不变的情况有所改善。 从图1和图2可以看出，从t=5开始用线性方程计算（图中T=5的曲线）与非线性方程计算结果误差相当接近。对定常问题，在一定程度上说明这样的改进算法是成功的，既节省了计算量（将非线性问题转化为线性问题）又保证了计算精度。  图1.压力的误差的曲线 图2. 非线性项部分不变时压力的误差的曲线 本文内容更新： /fvcp 作者：于欣，yu@, , 010 参考文献: 18．于欣,解不可压缩流体力学问题的降阶法 II. 空间Vh的基函数和误差估计计算物理3(2), 217-226(1986)，/yuxin/paper/div2.htm 19. 于欣，不可压缩/FVC 高温气体动力学国家重点实验室（筹）2012年夏季研讨会 LHDsummer-2012 2012年8月 长春 - 1 -