微可压缩模型(SCM)预处理技术深入的研究.pdfVIP

2005’全国流体力学青年研讨会论文集 微可压缩模型(SCM)预处理技术深入研究奉 邓小刚丛成华毛枚良向大平 中国空气动力研究与发展中心，四川绵阳621000 摘要：文献[1]提出的微可压缩模型(Sa曲是求解低马赫数流动的～种有效馍型，S蕊在求解过程中 不必满足速度散度为零的条件，可直接应用时间推进方法求解得到不可压缩N心方程的解。为了加速 收敛，文献[5]将预处理技术引入到SOI中，并在二维计算中得到了应用。_本文继续深入研究了预处理 技术，将其应用推广至B复杂二维獭和三维情形中，i疆过对NACA4412嚣型、椭球以及硒机嚣的数值 棱}jj；l，—方面进—步展示了S例的健壮性，表明S叫适1分于低马赫数流动的数值模拟：另—旒越■ 步验证了预处理技术在微可压缩模型中的作用，解决了简单三维低马赫数流动的收敛问题，并提高了 求解的准确性和精度。这为SOl应用于工程实际创造了—定的条件。 0引言 在文献[5]中提到，低马赫数流动求解的三种方式：基于密度的方法，基于压力的方法，基于渐近 思想的扰动方法，都面临着求解精度的问题。文献[1]对低马赫数流动进行了深入探索，在分析了低马 赫数理想气体的流动特性后，提出了微可压缩模型(SCM)，这种方法的最大贡献是将控制方程转换成了 双曲型方程，使整个求解过程保持在严}，框架内，不必满足速度散度为零的条件，可直接胡时间推 进方法求解得到不可压缩N-S方程的解。同时求解可压缩方程的几乎—切成熟方法都可以应用在S铡 的数值求解中。 S(]Vl提出后，以邓小刚为首的工作组对S蕊的数值模拟方法、特点以及应用做了很多探索“剐，结果 显示微可压缩模型具有较高的计算精度和效率，更适合刊氐马赫数流动计算，但也要看到，S锄的收敛 速度依然不能让^满意。我们知道，收敛速度是求解低马赫数流动时各种方法都面临的一个问题，主 要原因是对流速度和声速之间存在很大的差异。为了提尚挠龟嫩，文献[5]将预处理技术引入到S例 中，比较好地克服了控制方程的刚性问题，提高了S嘣的收敛速度，同时也提高了求解的准确性和精 度。 本文在文献[5]的基础上，继续深入研究了预处理技术，将其应用推广到复杂二维流动以及三维流 动的计算中，通过对NACA4412翼型、椭球以及M6机翼的数值模拟，—方面来验证SCM的健壮性，另 —方面进—步验证预处理技术在微可压缩模型中的作用。 1控制方程的预处理 文献[1]通过理论分析建立了微可压缩模型张量形式的控制方程，在此基础E推导出了曲线坐标系 下矢量形式的控制方程删，以非守恒变量．Sf=[p，“，1，，川7’为自变量对守恒变量形式的方程进行转换 有： 罢+N-Ⅳ罢+N一1删要+N～CNo。s：N。。YIS (1) af d二 d卵 d 其中： ‘基金项目：国家杰出青年科学基金(N0；国家自然科学基金(No；创新研究群体科学基金 (No 6 邓小刚、丛成华、毛枚良等： 微可压缩模型(SCM)预处理技术深入研究 Ⅳ=警’Ⅳ_1=嘉，Q-[愀pv，耐，彳=丽OF，缀黼慨 设P为预处理矩阵，对(1)式进行预处理得到： 罢+胛1斛罢+科～BNo。s+PⅣ～CNo，s：斛一’VIS (2) Ut U与 uT| U与 实际应用中，已经发展了—些相对成熟的预处理矩阵㈨∞，本文采用Weiss啪的预处理矩阵，它是 非奇异正定的，保证了控制概统的双曲性质，满足对称性、稀疏性等，形式简单，便于实际应用， 其具体形式为： 了1 0 0 0 D2 P=O 1 O 0其中∥=Min[1，Max(Ma，0．5Mr)](3) O 0 1 O 0 O O 1 其中^痧为来流马赫数，Ma为当地马赫数。 这样可以得到预处理后的Jacobi矩阵的特征