一种求解euler方程的间断galerkin方法的改进研究.pdfVIP

126 LHD 2011 年度夏季学术研讨会 一种求解Euler 方程的间断Galerkin 方法的改 进研究 1 1 2 刘洪伟 ，樊菁 ，徐昆 （1 中国科学院力学研究所高温气体动力学国家重点实验室 （筹），北京海淀区 100190） （2 香港科技大学数学系，香港九龙清水湾） 摘要 本文将对一种数值求解Euler 方程的间断Galerkin 方法提出改进的策略。 该间断Galerkin 方法 由Tang 和Warnecke 提出（Computers Fluids 34(2005), 375-398 ），数值试验已表明，此方法的3 阶及以上格式会有相当程度的精度损失。针对此问题，本文提出了一种改进的办法，可以使原方法 的精度不会随着格式阶数的增加而有所损失，从而确保了该间断 Galerkin 方法的高精度特性。对于 一维的情况，本文也研究了不同限制器对于此方法求解含间断问题的影响，给出了相关的建议。 关键词 RKDG 方法，气体动理学格式，数值通量，收敛精度 通量耗散性太大所导致，据此我们构造了一种 引 言 新的耗散性比较小的通量用于此方法，数值结 间断 Galerkin 方法是近年来发展比较迅速 果表明，格式的精度不会再随着阶数的增加而 的一种数值计算方法。这种方法有一些比较明 下降。 显的特点，例如，单元的计算模板只包括与其 1 基于气体动理学的间断Galerkin 方法 相邻的计算单元，这将有利于边界的处理、非 结构网格的应用和程序的并行等；可以达到高 这里我们首先来简要回顾一下由 Tang 和 精度；具有比较好的激波捕捉的特性等。间断 Warnecke 提出的基于气体动理学的求解 Euler Galerkin 方法的这些特点使其越来越受到相关 方程的间断 Galerkin 方法。简单起见，我们以 学者的关注和认可，目前它不仅被应用于计算 一维的Euler 方程为例，它可写成如下形式   流体力学领域，而且已被拓展到了很多其他领 Ut  F(U) x 0 。 （1） [2] 对于给定的网格I [x ,x ] ，我们可将 域的数值计算和模拟 。 i i1/ 2 i1/ 2 [3] 数值解近似为 气体动理学方法 是基于气体动理学理论 来构造求解气体运动方程数值格式的方法，目 k (l ) (l ) U (x ,t) U (t) (x), x I 。（2 ） h i i i 前已被应用于求解 Euler 方程、Navier-Stokes l 0 方程、稀薄气体流动控制方程和磁流体动力学 (l) I