全矩阵无积分间断有限元方法实现.pdfVIP

全矩阵无积分间断有限元方法实现# 谭勤学，任静，蒋洪德** （清华大学热能工程系） 5 摘要：本文主要针对间断有限元方法计算量及存储量开销巨大的问题，提出一种节约计算量 及存储量的实现方法。使用该方法，在相同自由度条件下，间断有限元方法计算速度及存储 量与普通有限体积方法相当。另外结合有限体积隐式求解方法，实现间断有限元方法精度多 重网格，进一步有效加快间断有限元方法收敛速度，改善收敛特性。 关键词：间断有限元，精度多重网格，无积分, 矩阵运算 10 中图分类号：V231.1 Quadrature-free implementation of Discontinuous Galerkin Method Using Matrix Product Formula Tan Qinxue, Ren Jing, Jiang Hongde 15 (Tsinghua University, Beijing 100084) Abstract: In this paper, a discontinuous galerkin only include matrix operation formulation that avoids the use of discrete quadrature formulas is described. Using this method, the memory requirements and computing speed are almost equal to the finite volume method When degrees of freedom of the same. A p-multigrid discontinuous Galerkin method that combine explicit DGM and implicit FVM is 20 presented for the solution of thecompressible Euler equations on unstructured grids. And the computation can converge much better and faster with this p-multigrid method. Key words: Discontinuous Galerkin Method, p-multigrid, Quadrature Free, Matrix Operation 0 引言 25 间断有限元方法(DGM, Discontinuous Galerkin Method)最早是由 Reed 和 Hill 在 1973 年 [1]提出的，用于中子输运方程的求解。该方法结合了有限体积方法(FVM)和有限元方法(FEM) 的基本思想，特别易于处理复杂边界及边值问题，同时具有灵活处理间断的能力，且可以通 过适当选取基函数，提高单元插值多项式的次数来实现任意高阶精度，另外具有容易实现自 适应、并行计算等优势[2] 。由于这些优异性质，在欧洲探讨用于未来工业设计高阶CFD 程 [3] 30 序的 ADIGMA 计划中 ，间断有限元方法被认为是最有前途的算法之一。 近年来 DGM 得到了非常大的发展，其中以 CockburnShu 提出的龙格库塔间断有限元 方法尤为突出，得到了广泛的应用。Cockburn 和 Shu[3-7]引入龙格库塔DGM 求解 Euler 方程。 Bass