求解双曲型守恒律的半离散中心差分格式.pdfVIP

第22卷 第4期 应 用 力 学 学 报 Vol_22 No．4 2005年9月 CHINESEJOURNAL OFAPPLIED MECHANICS Dee．2005 文章编号：1000-4939(2005)04-0536-05 求解双曲型守恒律的半离散中心差分格式 陈建忠 封建湖 史忠科 胡彦梅 (西北工业大学 西安 710072)‘(长安大学 西安 710064) 摘要：给出了一种求解双曲型守恒律的三阶半离散中心差分格式。该格式以一种推广的三阶重构 为基础，同时考虑了波传播的局部速度。格式的构造方法是利用重构，先计算非一致交错网格上的 均值，再将该网格均值投影回原来的非交错网格，得到新的全离散 中心差分格式，该格式有半离散 形式。本文半离散格式保持了中心差分格式简单的优点，即不需用Riemann解算器，避免了进行特 征解耦。它具有守恒形式，数值通量满足相容性条件。数值试验结果表明该格式是高精度、高分辨 率的。 关键词：双曲型守恒律；中心差分格式；半离散；重构 中图分类号：035 文献标识码： A 器，从而避免了进行复杂且耗时的特征分解过程。中 1 引 言 心格式又具有迎风格式的高分辨率。本文基于文献 [4～6]中的基本方法，利用文献[7]的三阶重构， 本文讨论如下双曲型守恒律的求解问题 给出了一种新的三阶半离散中心格式。数值试验结 + ) =0 (1) 果表明该格式是高精度、高分辨率的。 式(1)在流体力学和空气动力学等领域有着重要应 用。解决好式 (1)的求解问题，具有重要意义。 2 半离散的中心格式 近年来，用中心差分格式(简称中心格式)求解 式(1)成为研究的热点问题之一，国外许多学者作 根据图1所示，下面介绍文献 [4～6]中构造半 了大量的研究工作。Nessyahu和Tadmor…提出二阶 离散中心格式的基本方法 ，在此基础上再具体给出 中心格式，Liu和Tadmor2J，Levy等人 将该格式 本文格式的推导过程。假定 t=t时刻的网格均值 进行了更高阶推广。由于采用了交错网格，这些全离 {-n}已知，重构分片多项式”(，￡)=∑ ()xj()， 散格式以计算交错网格均值为基础。Kurganov等 是第 个网格 = [一下1，Xj+牛]上的特征函数， 人 提出了半离散的中心格式，其主要思想是利 用波传播的局部速度，更精确地估计黎曼扇的宽度。 g(x)是定义在 上的多项式。然后估计网格边界 + 处波传播的局部速度，它的上界定义为 与全离散格式相比，半离散中心格式的数值耗散更 3， 小，无需计算交错网格均值，方法也更易于推广到对 0士= ma．x P( ()) (2) 流扩散方程。与传统的迎风格式如Godunov格式和 ‘ (}·寺) MUSCL格式相比较 ，中心格式的主要优点是简单： 这里P(A)表示矩阵A的谱半径，工牛=p。(f+}，t) 不需用精确 Riemann解算器或近似 Riemann解算 和 }= tl,(+下1，t)分别表示V(x，t)在 +÷右边 ·基金项 目：国家 自然科学基金重点项 目 来稿 日期 ：2004-02-22 修回日期 ：2004