求解双曲守恒律方程的WENO型熵相容格式.pdfVIP

第 34卷 第 4期 爆 炸 与 冲 击 Vo1．34。No．4 2014年 7月 EXPLOSION AND SHOCK WAVES Ju1．，2014 文章编号 ：i001—1455(2014)04—0501—07 求解双 曲守恒律方程的WENO型熵相容格式 程晓晗 ，封建湖 ，聂玉峰 (1．西北工业大学应用数学系，陕西 西安 710072； 2．长安大学理学院，陕西 西安 710064) 摘要 ：通过在单元交界面处进行高阶WENO重构，得到了一种求解双曲型守恒律方程的WENO型熵相 容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟 ，结果表明，该格式具有高精度、基本无振 荡性等特点。 关键词 ：流体力学；wEN0型熵相容格式；WENO重构 ；双 曲守恒律 中图分类号 ：0354；0242 国标学科代码 ：1302564 文献标志码 ：A 许多物理模型都采用双曲型守恒律方程描述 ，如气体动力学 中的Euler方程、浅水动力学 中的浅水 波方程和等离子物理学中的磁流体方程等 。在一维情况下，有如下形式 ： 8 + 一0 (1) t 。 az 、一 式中：．zC-R，t0，H一 ( “， ) ，是守恒型 向量 ，，(z1)是通量。求解式 (1)的困难是 ，即使初始值充分 光滑，解也会在有 限时刻 出现间断。当解 出现间断时，古典解 的理论不再适用 ，为此 Lax提出了弱解的 概念 。然而弱解并不唯一 ，这就需要外加的标准确定满足相关物理特性的唯一解。定义熵对 (己，，F)，其 中U(H)是ll的一个凸函数 ，即Hessian矩阵L厂(H)对称正定 ，而 F(“)满足 (F ())：==(u (11))f (H)。 如果弱解 U满足熵稳定条件 ：8U(“)／ot+3F(H)／az≤0，则 “是唯一的且具有物理意义。一些数值格式 不满足熵稳定条件 ，如 Roe格式，往往造成数值解产生膨胀激波等非物理现象口]。为满足熵稳定 ，长期 以来的做法是添加数值黏性 ，但应该添加多少数值黏性则依赖所求解 的具体问题 ，且包含经验性常数 。 为研究熵稳定和数值黏性之间的关系 ，E．Tadmor_2定义了一类熵守恒格式 ，并证明了一个三点格式只 需含有 比熵守恒格式更多的黏性则是熵稳定 的。熵稳定格式 的提出，为获得具有物理意义且数值上稳 定的解提供了一种有效的方法 。P．L．Roe[3提出在熵守恒格式基础上加入 Roe格式的数值黏性项 ，使 新格式在保证熵稳定的同时保持 了Roe格式对激波的 良好捕捉效果 ，但该格式是一阶精度的。F．Is— mail等_4发展 了熵相容格式。对 比熵稳定格式，熵相容格式进一步控制了激波处 的熵增量 ，可有效避 免膨胀激波以及间断处的伪振荡等现象。 本文 中，进一步研究熵相容格式，通过在单元交界面处进行高阶重构，得到一类高精度 的加权基本 无振荡(weightedessentiallynon—oscillatory，WENO)型熵相容格式 。最后 ，通过一些数值算例验证所 构造的新格式的性能。 1 WENO型熵相容格式 1．1 一维标量守恒律方程 1．1．1 熵守恒格式 设给定熵对 (U，F)，定义熵变量 73—3U()／0u、熵势 (())一 (“)f()一F(“)，满足 f()一 * 收稿 日期 ：2012一¨一22；修 回 日期 ：2013-06—03 基金项 目：国家 自然科学基金项 目；中央高校基本科研业务费专项项 目(CHD2010JC060)； 西北工业大学博士论文创新基金项 目(CX201426) 作者简介 ：程晓晗(1987一 )，男，博士研究生 。 5O2 爆 炸 与 冲 击