一维ENO格式以及一维WENO格式在双曲守恒方程数值计算中的应用.pdfVIP

20O9年第 1期 桂林航天工业高等专科学校学报 (总第 53期) JOURNALOFGUILINCOLLEGEOFAEROSPACETECHNOLOGY 数学研究与应用 一 维ENO格式以及一维WENO格式在 双曲守恒方程数值计算 中的应用 周立新 陈凯 (桂林航天工业高等专科学校 计算机系，广西 桂林 541004) 摘 要 论文研究了高阶精度本质无振荡(ENO)及其高阶精度加权本质无振荡(wEN0)格式在双曲型守恒方程中的 应用 。并应用五阶WENO空间离散格式和三阶 TVDRunge—kutta时间离散格式对一维 SodRiemann问题 进行 了数值模拟。数值结果表 明该格式具有高精度性和本质无振荡性 。 关键词 高精度；双曲型守恒方程；ENO；WENO；TVDRunge--kutta方法 中图分类号：O241．3 文献标志码 ：A 文章编号：l0。9—1o33(20O9)o1一o1o8—04 目前，计算流体力学难点问题的研究中，高精度 、高分 ()，它在单元 内部是 (z)的k阶精度逼近 ： 辨率格式是一个热点研究领域 ，引起 了广大计算数学工作 A (z)一 (z)+O(z~x )，z∈Ii，一1，2，… ，N (5) 者极大的关注和研究热情 。 特别地 ，函数 (z)在单元边界上有 ： 双 曲守恒律方程一维 问题 的一阶单调格式具有较强的 j专一 (z+专)，二{一 t(zt一号)， 1，2，…，N 稳定性并且可以正确地收敛到熵解 ，但是在许多实际应用 (6) 中表现出太强的耗散性从而得不到高分辨率的效果。为了 它是 阶精度 的： 克服 TVD格式 的弱 点，出现 了ENO格式 。ENO格式放 ÷一 (z ÷)+0(△ )，三{一 (z一一i1)+ 宽了对数值解总变差不能增长 的限制 ，允许总变差微小增 O(△ )，i一 1，2，… ，N (7) 加 ，并使格式达到一致高阶精度且基本无振荡 。WEN0格 给定单元 f。和精度 后 ，我们首先选择一个模板 ： 式 由Liu、Osher等人在 ENO格式构造思想的基础上提 出 S()一 {I一 ，… ，+s) (8) 的，其核心思想是在 ENO重构过程中将所有的基架都引 它包括单元 j，它左边 的r个单元和右边 的s个单元 ， 入插值过程 ，对于不同的基架根据各 自的光滑程度分配不 并且 r+s+1一k，P()是唯一 的至 多 k—1次的多项式 。 同的权系数来优化组合 ，跨越 间断的基架 的权系数近似为 它在 S()上 的单元均值与 (z)一致 。 0，从而提高离散格式在光滑 区域 的计算精度 ，达到高分辨 率 的 目的。 一 — h 川 ，… s (9) 1 ENO格式