解Euler方程组的一类二阶精度无振荡隐式有限差分格式＋一t－.PDF

维普资讯 J990年第4期 西 北 大 学 学 报 NO．4 1990 第20鞋 JOURNAL OFNORTHWESTUNIVERS兀．Y Vo1．20 解 Euler方程组的一类二阶精度 无振荡隐式有限差分格式 一， ’ 董海涛 (数学研究所) t 摘 要 奉丈构造了一类二阶精度无振荡隐式有限差分格式采计算双曲型守恒律的弱 解。这些格式由两步构成。第一步用线化一阶守恒格式，第二步用钱化二阶无振 荡格式。由这两步构成的格式具有绝对稳定性。无振荡性．二阶精度及较好 的守 恒性．数值实验也说明了格式的这些性质 ． 关键词：敏波；奎变差不增：无振荡；臆式有限差分格式． 0 引 言 1983年，Hartcn发表了关于双曲型守恒律的显式TVD格式。此格式的解对激波具’ 有商分辨率。本文用隐式TVD格式求解 Eulcr方程组。采取了特殊的线化方法 。既保持 TVD性，又有守恒性。还不失绝对稳定性。 双曲型守恒律方程组的初值问题 r ．十，( =0。 ‘ (1) ， 0)= )。 一∞ 。。． ．1)是m维向量I酒量 为m维向量函数．如果矩阵 ( 一 的所有特征值口1。… 。是实数且右特征向盘 。…，五 的集合是完备的．则称方程组0)是双曲型． 1 格 式 逼近(1)的守恒臆式为 + (：：，一t-”1，：)一 本文1990年5月31日收到． 维普资讯 一 14一 西北大学学报(自然科学翮 1990年 一~1” f… k+1， … ， l+i／2 - 先从标量方程开始，构造一阶格式 ，再修正成二阶格式 ．取方程 (1)中 ／n=I．( = At／Ax)． 格式(2)可写成形式 ’ 一 c 。，：A，+．，：v + c[一 A，一。， v = - (3) -- 1／2 越 A1+1／2= ，+1一V，，格式(3)为TVD的充分条件为 0≤c 1／2≤c∞ ·