基于非均匀有理B样条的间断Galerkin算法通用参数化边界修正.pdfVIP

中国工程热物理学会 热机气动热力学 学术会议论文 编号：092094 基于非均匀有理B 样条的 间断Galerkin 算法通用参数化边界修正 朱晓峰 卢义 袁新 清华大学热科学与动力工程教育部重点实验室，北京市 100084 Tel: 010 Email: sxf03@mails.tsinghua.edu.cn 摘 要：本文采用隐式间断Galerkin 方法对C3X 叶型Euler 流场进行了求解。为改善低阶边界对高阶 间断Galerkin 精度的影响，采用4 阶3 次非均匀有理B 样条（NURBS ）对几何边界进行统一描述，避 免了高阶网格边界之间的不连续性，给出了处理高阶精度算法所要求的高阶边界的通用的处理方法， 其计算量和存储量的要求相比高阶网格大大减少。计算结果表明，本文所作的参数化边界修正方法有 助于间断Galerkin 方法使用通用非结构网格进行高精度的计算。 关键词：间断有限元，高阶曲线边界，参数化修正，NURBS ，非结构 0 前言 复杂几何形状的高阶精度计算一直是CFD 计算中的难点，间断有限元法(DGM)综合 了有限差分法(FDM) 和有限体积法(FVM) 的优点，使用局部任意阶重构，突破了 ENO(Essential Non-Oscillatory)和WENO(Weighted Essential Non-Oscillatory)方法模板选 择的限制，能够使用非结构网格、非匹配网格以及重叠网格进行任意阶精度的计算[1]。 边界条件的准确处理对于高精度CFD 程序至关重要，当数值程序的边界与几何边界 不相吻合时会引起误差，从而有可能影响主区域内的求解。Bassi 和Rebay 指出，即便使 用非常细密的一阶网格，二阶或大于二阶精度的DGM 程序计算都会产生非物理熵增， 必须至少使用二阶精度的网格[2]。 1 间断有限元方法 守恒型的二维Euler 方程表述为如下形式： ?u +??F(u) 0 (1) ?t 守恒变量u 和组成通量F(u) 的笛卡儿分量f (u) 、g(u) 分别为： 基金项目：国家自然科学基金资助项目）；国家重点基础研究发展计划（973 计划）资助项 目（2007CB210108 ） ρ ?ρu ? ?ρv ? ? ? ? ? ? 2 ? ? ? ρu ?ρu +p ? ?ρuv ? u ? ?； ( ) ； ( ) (2) ? ? f u ? ? g u ? 2 ? ρv ?ρuv ? ?ρv +p ? ? ? ? ? ? ? ρE ρuH ρvH ? ? ? ? ? ? 上式中ρ 为流体密度，u 和v 为速度分量，p 为压强，E 和H 分别为总能和总焓，总 焓