低耗散非结构网格多维限制器艾邦成.pdfVIP

网络出版时间：2017-09-29 15:39:18 网络出版地址：/kcms/detail/11.2011.O41539.013.html 低耗散非结构网格多维限制器 1 2 艾邦成, 张亮, 陈智 （中国航天空气动力技术研究院，北京，100074） 摘 要:基于极值原理分析了非结构网格多维限制器的构造准则，通过引入约束位置和约束模板的定义，将非结构网 格多维限制器构造方法统一化，并以此为基础提出了一种新型非结构网格多维限制器。通过典型算例计算与已有限 制器进行了比较。结果表明:对于包含激波、接触间断等复杂波系干扰流动，新型多维限制器在保证数值稳定的条 件下具有更低的数值耗散，对于复杂波系结构具有更高的流动分辨率。 关键词: 非结构网格; 多维; 限制器; 极值原理; 低耗散 中图分类号: O362 文献标识码: A 0 引言 随着计算流体力学和高性能计算的不断发展，数值模拟技术在各类工程问题中得到了广泛的应用[1] [2] 。由于实际外形往往较为复杂，高质量结构网格生成困难，所以非结构网格逐渐成为数值模拟技术发 展的重要趋势。在非结构网格技术中，通量格式、时间格式以及加速收敛等相关技术一般均可直接利用 已有的结构网格研究成果。然而由于非结构网格具有特殊的几何特征，在结构网格领域发展得相对成熟 的限制器技术难以直接拓展到非结构网格框架。近年来，发展与非结构网格几何特性相匹配的多维限制 器技术一直是非结构网格技术发展的重要领域。 [3-6] 限制器概念的引入源于上世纪七、八十年代高分辨数值格式的发展 ，其基本思想是通过对耗散 较高的低阶格式添加“限制的”负耗散以实现数值格式在保证稳定的前提下最大程度降低耗散水平，提高 格式精度。在Harten[7]提出TVD概念以后，Sweby[8]通过理论分析将限制器构造技术统一化，建立了经典 的TVD 限制器构造准则，成为了目前不同限制器研究的重要理论基础。van Leer 、minmod 、superbee及 van Albada等在内的多种经典限制器均可纳入这一理论框架之下，其之间的差异性仅体现在不同的耗散 水平以及限制函数的光滑性。虽然这一理论框架仅适用于一维条件，但结构网格具有天然的维度分解特 性，在每一维度上分别应用该限制器技术可有效实现多维流动的稳定准确模拟。因此，TVD限制器在结 构网格数值模拟技术中取得了巨大的成功。 而对于非结构网格，由于网格的无序性导致并不存在显著的维度分解方向，一维TVD限制器技术无 法得以应用，所以需要建立严格的多维限制器技术以保障数值模拟的精度和稳定性。Goodman和LeVeque 收稿日期: ; 修回日期: 1作者简介: 艾邦成（1972-），博士，主要从事高超声速气动热和热防护研究，联系电话：010，E-mail： aimen011@163.com 2通讯作者：张亮（1982-），博士，主要从事高超声速数值模拟技术研究，联系电话：010，