非结构六面体网格下间断Galerkin有限元方法的限制器及间断探测器研究.pdfVIP

第十六届全国计算流体力学会议 非结构六面体网格下间断Galerkin有限元方法的限制器 及间断探测器研究 1 2 2 龚小权 ,马明生 ,邓有奇 （1. 西北工业大学，陕西西安，710072 ） （2. 中国空气动力研究与发展中心,四川绵阳，621000 ） 摘要：随着对精细的流场结构、准确的气动数据以及网格质量、生成速度要求越来越高，基于非结构网格的高精度 计算方法得到快速发展。本文针对建立在大拉伸比的非结构六面体网格的间断Galerkin有限元方法(DGM)，对方法 中限制器及间断探测器的使用进行了初步研究。首先采用等熵涡问题验证了发展的方法的精度；然后在NACA0012 翼型上分析比较了Barth–Jespersen 限制器、Venkatakrishnan 限制器以及Hermite WENO 限制器三种限制器在DGM 中 的表现；针对方法对限制器敏感的特点，配合基于密度的间断探测器来侦测“ 问题单元” ，对“ 问题单元”采用上面的 三种限制器，同时给出“ 问题单元”位置。NACA0012翼型的跨音速无粘流动数值模拟试验表明这三种限制器能抑制 激波附近的震荡；限制器配合使用间断探测器具有更高的计算精度和更好的激波捕捉能力。最后数值模拟三维 ONERA M6 机翼来考察了Venkatakrishnan 限制器在DGM 中的表现。结果表明该类斜率限制器能够捕捉复杂流场中 的激波，但是全流场采用该限制器可能会对整个流场的计算精度有较大影响。 关键词：间断Galerkin有限元方法，限制器，间断探测器，六面体网格 此后该方法迅速发展，被广泛用于计算空气 0 引言 动力学、计算声学、计算电磁学。与其它高 [1][2] 精度格式一样，DGM 在间断附近也会产生数 间断 Galerkin 有限元方法 (Disconti 值振荡，导致非物理解例如出现负压和负密 nuous Galerkin finit element Method,DGM)在求 度从而导致计算中断。目前通用的做法是通 解守恒方程组问题中得到广泛应用，近年来 过激波捕捉方法或者构建限制器来抑制激波 许多学者对 DGM 开展了大量的研究工作。 附近的震荡。前者需要加入一些与网格和计 间断 Galerkin 有限元方法结合了有限体积法 算问题都相关的系数，因此非常不方便。故 和有限元法的优点。作为一种有限元方法， 通常采用限制器来抑制振荡。大多数有限体 它采用有限元法中基函数思想来构造单元内 积法中的斜率限制器都能用于 DGM 中，例 的变量分布，容易实现高精度。同时它又具 [4] 如 Barth –Jespersen 限 制 器 ， 有有限体积法积分形式的计算形式且融入了 Venkatakrishnan[5] 限制器。Chi-Wang Shu 和 J. 有限体