网格与高精度差分计算问题.pdfVIP

网格与高精度差分计算问题 张涵信 (中国空气动力研究与发展中心 国家计算流体力学实验室) 摘耍：本文研究TNS方程差分求解时来流雷诺数、计算格式精度和计算网格之间的关系，给出 了判定空同三个方向上的粘性贡献在给定雷诺数、格式精度和网格下是否eB够正确计入 的定量准则，指出在NS方程的二阶差分方法的数值模拟中，由于物面法向采用了压缩冈 格技术，物面附近的同格间距很小，该方向上的粘性贡献可被计人，但是如果流向和周围 的网格较粗，相应的差分方程中的粘性贡献仅以截断误差相同的量级被计人，因此在精 度上等于仍是求解略去流向和周围粘性项的薄层近似方程．文中指出，高阶精度的差分 计算格式。可以避免对网格要求苛刻的困难。本文进一步讨论了建立高阶精度格式的问 题，提出了建立高阶精度格式应该满足的三个原则：稳定性或耗数控制原则、熵增或色散 控制原则，以及激渡关系控制原则．根据这些原则，论文倡导发展混合高阶精度格式，即 在连续区，采用高阶精度格式，在激渡区，采用能近似代表激渡关系的NND格式． 关键词：NS方程，网格判则，高精度格式，格式建立原则 一、引言 随着计算机的快速发展和数值方法能力的迅速提高，复杂外形绕流的NS方程 数值模拟的研究工作越来越多．在这些计算模拟工作中，有的计算采用数十万个网格 点，有的采用数百万个甚至上千万个网格点。分析这些研究给出的结果大致可得到以 下结论：第一，当用高阶精度差分格式求解时，所需网格点数比低阶格式要少；第二， 当用二阶精度格式求解时，对于高Reynolds流动，在物体的迎风区，粗网格的结果和 实验值尚能接近，但对于背风、拐角等具有分离、旋涡和强烈压缩或膨胀的流动区域， 其粗网格的结果和实测值相差甚远[1]；第三，细网格可给出与粗网格完全不同的旋涡 和分离结构乜]。这些情况给我们提出如下问题： 1．当用二阶格式求解复杂外形绕流的NS方程时，应该如何选择网格点数才能 获得满意的结果?有无准则可以采用? z．粗网格的数值模拟，在多大程度和多大范围内能反应粘性的效应?其机理和规 律是什么7 3．对于高阶精度的差分格式，为什么网格点的数目可以放松? 4．对于有激渡的流场．如何建立高阶精度的差分格式7 ． ’ 本文目的就是探讨以上问题。包括引言本文共分四节。第二节讨论网格数与差 分格式精度的关系，这里我们提出了确定网格数的规则，并阐明了网格点数控制粘性 影响的机制。第三节讨论了建立高阶格式的原则，这些原则可帮助我们选择和分析各 种高阶格式．第四节是简单的结论。 二，网格点数和格式精度的关系 大家知道，三维流动无量纲化的Ns方程可写成： ) (1)… 警+篆+豸+警—一I：．(a。F,+等+百81-1,Re盘’缸’却’玉 L口’西‘玉。 这里为了简单，我们略去了无量纲化的方法和方程(1)中各项和各个符号意义的说明 差分格式求解方程(1)，与m阶差分方程等价的修正方程是： 甜．甜．靠，aH 1，积。，箔。，aH。、 i十磊十万十i2瓦‘面十万十虿’ (2) +O(nz“，△，，△矿，…) 以上的小量。(2)式可进一步被写成： 型+鬈+雾+i81-18t ‘知‘西‘＆ mE血4缸’Re，△， ＆7 (3) +O(zXz“，△，，＆”…) ● 选择(3)中的a，p和丫，使其满足 a一馐