FV混合格式在二维非结构网格中的推广.pdfVIP

第28卷 第2期 计 算 物 理 Vo1．28，No．2 2011年 3月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Mar．，2011 文章编号：1001-246X(2011J02-01~8．1I 基于静动态重构的高阶DG／FV混合格式 在二维非结构 网格中的推广 张来平 ， 刘 伟 ， 贺立新 ， 赫 新 ， 邓小刚 ’ (1．空气动力学国家重点实验室，四川 绵阳 621000；2．中国空气动力研究与发展中心，四川 绵阳 621000) 摘 要 ：通过比较 间断 Galerkin有限元方法 (DGM)和有限体积方法 (FVM)，提出 “静态重构”和 “动态重 构”的概念，进一步建立基于静动态 “混合重构”算法的三阶DG／FV混合格式．在DG／FV混合格式中，单元 平均值和一阶导数由DGM方法 “动态重构”，二阶导数利用FVM方法 “静态重构”；在此基础上，构造高阶多 项式插值函数 ，得到三阶精度 的DG／．FV混合格式．将 DG／FV混合格式推广应用于二维非结构网格 ，求解二维 标量方程和 Euler方程．典型算例数值实验表明，DG／FV混合格式达到了理想 中的三阶精度． 关键词 ：DG有限元法 ；有限体积法 ；Taylor基 ；混合格式；非结构网格 中图分类号 ：035 文献标识码 ：A 0 引言 近年来，高阶精度计算方法越来越受到CFD工作者的高度重视 ¨ ．高阶精度格式具有较小的数值色 散和耗散，适合求解多尺度流动 问题，如湍流、气动声学 问题等．在众多的高阶精度计算方法 中，间断 Galerkin有限元方法(DGM)备受关注．间断有限元方法最早 由Reed和Hill为解决中子运输方程提出 ．在 此基础上 ，特别是近年来 ，出现了种类丰富多样 的DG方法 ，DGM保持了传统有 限元方法 (FEM)和有限体 积法 (FVM)的优点，融入了高分辨率有限差分方法 (FDM)和有 限体积方法 中如数值通量、Riemann间断分 解、TVD和限制器等思想 ；同时它又是一种积分形式的计算方法 ，避免了有限差分方法 (FDM)对计算网格的 限制 ，可以直接应用于非结构网格或混合 网格．与其它数值方法相 比，DGM具有突出的优势．然而 ，DGM 的 计算量和存储量 巨大，其在复杂外形的大型数值模拟方面仍有许多不足，大大限制了其工程应用． 相比DGM，传统的二阶FVM 的计算量和存储量均要小很多，但是 FVM方法提高精度需要扩展模板，如 对于三维问题四阶精度的有限体积格式 ，至少需要 20个单元 的信息．由于非结构 网格的数据存储方式是随 机的，使得搜索临近单元费时费力 ，同时扩展模板也导致边界条件的处理不便．因此，综合 FVM和DGM 的优 势 ，建立二者混合的算法是一种必然也是 自然的选择．目前 ，混合算法主要有 以下三种形式 ：1)控制方程的 无粘和粘性通量分别采用 FVM和DGM计算的混合方法 ；2)基于计算 区域分解 的混合方法 ；3)基于多 项式混合重构的混合方法，如Dumbser等人提出的PP 格式 ，HermitWENO格式 等． 在以往的工作中 ，我们沿着第三种混合算法的思路 ，通过 比较 FVM和 DGM的差异 ，提出了 “静态重 构”和 “动态重构”的概念 ，在统一的框架下对 FVM和 DGM进行表述 ，进而 自然地构造 出一类高阶的DG有 限元／有限体积混合格式(简称 DG／FV混合格式)，并在一维线性和非线性标量方程的计算 中得到 了较好的 应用．数值实验表明DG／FV混合格式能够达到理想的精度 ，有些格式甚至出现超收敛的现象． 本文旨在将 DG／FV混合格式推广应用于二维非结构网格的计算 ，求解二维标量方程和 Euler方程．其 基本思想是利用基于Taylor基的DGM计算单元上的低阶导数(“动态重构”)，利用FVM重构单元的高阶导 数 (“静态重构”)，通过 “混合重构”实现在二维非结构 网格上的高阶DG／FV混合格式．通过二维标量方程 收稿 日期 ：2010—03一Ol；修回 日期 ：2010—07—16 基金项 目：国家重点基础研究发展计划 (973)(项 目编号 ：2009CB723802)及 国家 自然科学