双曲守恒律的几种新数值方法的比较研究.pdfVIP

维普资讯 200]年 6月 应用数学与计算数学学报 第 15卷 第 1期 June．2001 CoMM oN APPL M 盯 H．AND COMPUT Vol15 N01 双 曲守恒 律 的几 种 新 数 值 方 法 的 比较 研 究 窦 红 汪 继 文 ，0 (1．安徽 大学 教学 系，合 肥 230039) (2 中国科学 技术 太 学数 学 系 ，告 吧 230026，3 安暂 大 学计 算机 科 学 与工 程系】 摘 要 ：本 文就 一维 线 性 双 曲方 程 的光 滑和 间断 两 种 初 值 问题 的求 解 ，对 双 曲 守 恒律 的三种 新 数 值 方 法 ，即 ，W ENO 方挂 ，间 断 Galetkin方 法 和 全局 复 台方 法 ， 进 行 了数值 比较 实验 ，在 精度 、计 算速度 等方面 的比较 上 ，对这三个 方法有 了 个 较详 细 的 了解 ，得 到 了 些 有 用 的结论 ． 关 键 词 ：职 曲守 恒 律 ， wENo 方挂 ．间 断 Ga~rkin方 法 ，全局 复 舍 方 法 ． 1．引 言 双 曲守恒 律 问题 的数值 求 解 方 法 的研 究 是计 算 数 学、计 算 流体 力 学 的重 要 研 究 方 向之一 ．由于 非线性 双 曲守 恒方 程 (组)存 在 间 断解 ，这给 方 程 的求 解 带 来非常 大 的匿难 ．人 们遇 到 的 匿难 是 ：如 果采 用低 阶 方法 ，计算稳 定 ，不 会产 生 振 荡 ，但 数值 解 在应 该 出现 间 断 的地 方 会有 较 大 的耗 散 ，会 抹 平 激 波 ．如果 采 用高 阶 方 法 ，在解 的光 滑部 分 有很 好 的逼 近 ，但数 值 解 在激 波 附近 会 出现 虚 假 的非物 理振 荡 ，甚 至振荡 使得计 算无 法进 行．因此 ，寻求好 的数值 方法 去求解 双 曲守恒律 问题 ，使得 既能高精 度 地 逼 近真解 ，求 出锐 利 的间断，叉不 产 生非 物 理振 荡，一直 是人们 努 力 的方 向．几十 年 来 ，人们 为此 发展 了各种各样 的理 论和 方法 ，如人 工粘 性方法 ，自调节混 台方法 ，反扩 散方 法 ，变 差减 步 滤波 器， 以及 著 名 的TVD方 法 和恨 制 器方 法 ．但 这 些方 法 都有 些不 尽 如人 意 的方 面 ． 最 近十 余年 ，特 别是近 五 、六年 ，出现 并 发 展 了几个 值得 注 意 的新 方 法 ．一 个是本质无振荡 (ENO)和加权 ENO(WENO)方法 l【，2，3】，一个是 间断 Galerkin方法 [4，5】，这 两个 方 法 结合在 时间离散上 使 用 新发展 的高 阶 TVDRunge-Kutta方 法 ， 取 得 了较 已有 方 法有 较 大改进 的数 值结 果 ．另 外 ，98年 底 ，出现 了一种 新 的全 局复合方法 6【，7】，在 Euler方程 、浅水波方程 的求解方面取 得 了较好 的数值实验 结果 这几 种新 方 法各 自展 现 了 良好 的结果 ，但 是 ，它们 比较 起 来 有 些什 么 特 点和 优 缺 点，人们 却不 十分 清 楚 ， 本 文 就 线 性 双 曲方程 的光 滑 和 间 断 两 种 初 值 问题 ，对 WENO 方 法 、间 断 Galerkin方法 和全 局 复合 方 法进 行 了数值 比较实验 ，在 精度 、计 算速度 等 方面 的 比较 上 ，对 这 三个 方法 有 了一 个较 详 细 的 了解 ，得 到 了一些 有用 的结论 ． 本 文 2000年 t0月 17 日收 到 本文得 到 国家 自然科 学基金 的支持 维普资讯 应 用 数学 与 计 算 数 学 学 报 15卷