《一种基于间断 Galerkin 格式的新型限制器设计.》.pdfVIP

第 32卷 第 4期 空 气 动 力 学 学 报 Vol_32。No．4 2014年 O8月 ACTA AERoDYNAM ICA SINICA Aug．，2O14 文章编号 ：0258-1825(2014)04—0462—06 一 种基于间断 Galerkin格式的新型限制器设计 郭永恒，杨 永，刘逸鸥 (西北工业大学 翼型 叶栅空气动力学 国防科技重点实验室 ，陕西 西安 710072) 摘 要：使用高阶问断Galerkin格式求解守恒律方程组时，激波附近的Gibbs效应容易导致非物理解的产生。为 抑制这一现象，必须构造合理的限制器对数值解进行处理。 目前间断 Galerkin格式 中的限制器多源于有限体积 法 ，在非结构网格上只对低阶导数项进行限制 ，对高阶导数项则很难给 出普适判据 。文章对 间断 Galerkin解进行 广义 Fourier展开 ，实现不 同频域范围内的谱分解 ；在新 的模板坐标系下描述各阶方 向导数 的变化规律 ；结合当前 单元和相邻单元的信息，分层限制各阶方向导数 ，实现对非物理解的抑制。通过求解 Euler方程，对二维Riemann 问题 、翼型周围的亚、跨声速流动 问题、前台阶问题 以及超燃冲压发动机 内流场激波反射 问题进行数值模拟，检验 了新型限制器的可靠性 以及 向高阶格式推广 的可行性 。 关键词 ：间断 Galerkin格式；限制器 ；广义 Fourier展开 ；模板坐标系；方 向导数 中图分类号 ：V211．3 文献标识码 ：A doi：10．7638／kqdlxxb-2012．0157 0 引 言 Galerkin格式的具体形式 ，把握 “局部解的方向导数” 航空航天事业的迅速发展 ，对计算流体力学提出 这一信息，并在新 的模板坐标系下对其进行分层 限 了越来越高 的要求 ，特别是在研究复杂流动问题方 制 ，有效避免了非物理解 的产生，在保持高阶精度不 面，传统的有限体积法 已经很难适应新形势的需要 。 变的同时，使计算过程得以Jl~,N进行。 于是，许多新型的高精度算法应运而生 ，其 中，间断 1 间断 Galerkin格式的建立 Galerkin格式口的研究与应用在最近几年逐渐成为 一 个热点问题 ，并在计算流体力学领域取得了一系列 直角坐标系下二维 Euler方程的守恒形式为 令人瞩 目的结果。然而 ，当处理间断解问题时，Gibbs + 0下F(Q)+ OG(Q)一o (1) 一 3t 。 a-z ’ av 、一 效应很容易导致非物理解 的产生，使 间断 Galerkin 格式的计算过程 出现中断。为了避免这种现象的发 首先 ，把整个求解 区域 n划分成三角形非结构 生 ，人们在 间断 Galerkin格式中应用大量的限制器 ， 单元的集合 ，即 ===Un ，其 中NE为网格单元的个 如 Shu的斜率限制器l2]、Hoteit改进的VanLeer限 数 ；同时，在第 k号单元上 ，设局部有限元解为 制器 ¨3]、Tu和Aliabadi改进 的VanA