新型间断Galerkin激波捕捉方法的研究.pdfVIP

第十五届全国计算流体力学会议论文集 新型间断Galerkin激波捕捉方法研究 于剑，姜振华，阎超 (北京航空航天大学航空科学与工程学院，北京100191) 摘要：开展了适用于间断Galerkin有限元方法的典型激波捕捉格式研究，包括HWENOfHermite weighted essentially 重构方法，并将其用于高阶DG格式限制过程中。发展的方法充分利用了vailNeumann neighborhood可获 得的信息以保证DG模板的紧致性。同时，提出了一种简单而有效的方法增加单侧模板的光滑性；谱人工粘 性方法方面，发展了一种基于高阶导数的人工粘性方法，用于抑制间断附近所产生的非物理振荡。为了保证 人工粘性模型中高阶导数的计算精度，提出了一种高阶导数的计算方法。数值实验表明，本文发展的方法能 够高分辨率地捕捉间断，同时对光滑区域的污染较小，有较高的计算精度。 关键词：计算流体力学；间断Galerkin有限元方法；激波捕捉；Hermite weighted non．oscillatory： essentially 人工粘性 0引言 【4l首先将满足TⅦ性质的minmod限制器引入 RKDG方法，但是该方法存在一个常数需要调节， 并且限制器作用于一阶分量，在抑制虚假振荡的同 可压缩流动的高精度数值模拟仍然是一项颇 时也具有TVD格式的缺点，即可能将格式的精度 具挑战性的工作。典型的TVD格式能够达到较高 限制在一阶精度，尤其对于大于二阶精度的情况。 的分辨率，也是目前最为常用的数值方法。但是理 论分析和计算实践都表明，在流场的极值点附近， 限制器【5】，即按照从高阶分量到低阶分量的顺序进 这类格式的精度往往会降到一阶，从而带来较大的 行选择性限制。该方法放弃了TVD性质，能较好 数值误差。况且TVD格式的理论是针对标量守恒 地保持格式的高阶精度，但是向高阶精度和非结构 律发展的，推广到非线性系统后，这些格式的数学 网格的扩展较为困难。Qiu等人提出在包含间断等 特性便不能保证了【11。WENO格式克服了TVD格 大梯度区域采用WENO格式来代替TvB限制器进 式的缺点，在有效抑制非物理振荡的同时，保证了 行重构计算【6】。为了减小WENO格式的模板，他们 格式的高阶精度。然而，WENO格式在非结构网格 又引入HWENO格式【7】来代替上述的WENO格式。 上的设计复杂，为了增强稳定性，所需网格模板的 抑制虚假振荡的另外一种方法就是引入人工 规模通常远大于精度对模板单元数的要求，这给边 粘性项。结合DG的人工粘性方法目前也发展了很 界的处理和编程都带来了较大的困难[21。间断 Galerkin有限元方法(discontinuou