求解双曲型守恒律的五阶松弛格式.pdfVIP

第 26卷 第 4期 空 气 动 力 学 学 报 VoL26。No．4 20o8年 12月 ACTA AERODYNAM ICA SINICA Dec．．2OO8 文章编号：O258—1825(2008)O4一O5O8一O5 求解双 曲型守恒律的五阶松弛格式 陈建忠，史忠科 ，封建湖。，胡彦梅 (1_西北工业大学，陕西 西安 710o72；2．长安大学理学院，陕西 西安 710064) 摘 要 ：给 出了一种求解一维双曲型守恒律 的五阶松弛格式。该格式 以五阶wEN0重构和显隐式 Runge-Kutta 方法为基础。本文格式保持了松弛格式简单的优点，即不用 Riemann解算器和计算非线性通量函数的雅可比矩 阵。用该格式对一维 Eu1er方程进行了数值试验 ，并与三阶和 四阶松弛格式的计算结果进行 了比较 ，结果表明本文 的格式具有更低 的数值耗散和更高的分辨率 。 关键词 ：双 曲型守恒律 ；松弛格式；wEN0重构 中图分类号：V2l1．3 文献标识码 ：A 数值试验 ，并与三阶和四阶松弛格式的计算结果进行 0 引 言 了比较，结果表明本文格式具有更低的数值耗散和更 高的分辨率。 松弛格式是 由Jin和 Xin[]l995年提 出的一种 高精度、高分辨率有限差分方法，其主要思想是将对 1 五阶松弛格式的构造 守恒律的求解转化为求解松弛方程组 ，当松弛率趋于 考虑一维双曲型守恒律 零并满足一定的限制条件时，松弛方程组的解趋于守 ll+，(n) 一 O， H∈R (1) 恒律的解 。松弛格式的优点是不用求解 Riemann问 利用松弛方法[1]，对守恒律式(1)的求解可转化为求 题和计算非线性通量函数的雅可 比矩阵，因此特别适 解松弛方程组 合复杂模型 问题 的数值模拟。一阶和二阶松弛格 式L1在 Hamilton_Jacobi方程 、动力学方程、浅水方程 a + 一 o， ￡ 。az 一 的计算和反应流及两相流 的模拟等方面得到应用。 2O04年，SeaIdL2基于三阶紧凑中心 wENO重构[3]， +A 一÷(— )) (2) 构造了一种具有三阶精度的松弛格式，但应用该格式 其中 ’∈R”，A=diag{n，n，…，口}是正的对角矩阵， 计算浅水方程和气动力学Euler方程组时，数值粘性 rO是松弛率。这种转化的主要优点在于利用松弛 还是偏大。为降低数值耗散 ，提高分辨率，文献[4]提 方程组(2)的线性特征场和刚性低阶源项的特殊结构 出了一种四阶松弛格式，对一维问题 ，该格式以四阶 构造简单 、有效的数值方法。 中心WENO重构[5]为基础 ，对二维问题，用逐维计 松弛方程组 (2)有两个特征变量 算的方法将 四阶中心 WEN0重构进行了推广。作 士A ／H (3) 为对文献[1，2，4]中格式的高阶推广，本文给出了一 当r一0时，若下列的子特征条件(subcharacteristic 种具有五阶精度的松弛格式，将守恒律转化为松弛方 condition)成立 。 程组，并用