双曲守恒律方程ENO优化.pdfVIP

第34 卷第1期 中 　国 　科 　学 　技 　术 　大 　学 　学 　报 Vol. 34 ,No. 1 2 0 0 4 年 2 月 JOURNAL OF UNIVERSITY OF SCIENCE AND TECHNOLOGY OF CHINA Feb. 2 0 0 4 文章编号 (2004) 双曲守恒律方程 WENO 格式的优化方法 徐振礼 ,邱建贤 ,刘儒勋 ( 中国科学技术大学数学系,安徽合肥 230026) 摘要 :Weighted Essentially NonOscillatroy (WENO) 是求解双曲守恒律方程的高精度 高分辨率数值格式. 论文讨论了双曲守恒律方程 WENO 格式的一些优化策略 ,减少 了非线性权的计算次数和特征分解的次数 ,通过数值算例证明了这些策略的可行性 , 并比较了优缺点. 关键词 :WENO 格式 ;双曲守恒律方程 ;RungeKutta 方法 中图分类号:V211. 1 　　　文献标识码 :A AMS Subject classif ication( 2000) : 65M06 ,65M99 ,35L65 1 0 　引言 在流体力学和工程计算中 ,双曲守恒律方程的数值格式的设计具有很重要的地位 ,尤其 是常常在求解区域中会遇到一些物理量有很大的梯度变化或者间断的问题 ,如水跃问题 ,气 动力学上的激波问题等 ,模拟这一类问题是非常困难的. 自上世纪90 年代以来 ,基于 ENO ( Essentially NonOscillatory) [4 ,5 ] 发展起来的 WENO (Weighted ENO) 格式一直是人们研究的热点. WENO 格式具有高精度 ,高分辨率的良好性 质 ,因此特别适合间断问题的处理 , 比如激波湍流问题, 航空声学问题等. Liu , Osher 和 Chan[12 ]首先构造了一维空间上的三阶有限体积 WENO 格式 ,Jiang ,Shu , Hu 和 Qiu 等许多 人对发展 WENO 做了许多工作[1 ,3 ,6 ,7 ,14 ,16 ] . WENO 是通过低阶数值流通量的加权平均得 到高阶数值流通量 ,其中的非线性权值通过数值解的局部光滑性来确定 ,这部分计算占用了 大量的 CPU 工作时间 ,使得计算效率很低. Levy 、Puppo 和 Russo [10 ] 首先基于 TVDRunge Kutta 时间离散方法 ,将在每一时间步 ,只作一次非线性权计算应用到中心 WENO 中 ,并证 明是可行的[8 ,9 ,13 ] ,本文在迎风 WENO 中也证明了它的可行性 ,并讨论了在求解双曲守恒 律方程组时 ,在一个时间步中只进行一次特征分解的可行性. 1 　WENO 格式优化方法讨论 首先以一维标量方程 ( ) ( ) ut + f u x = 0 1 收稿日期 ( ) 　基金项目 : 国家自然科学基金 和中国科学技术大学火灾科学国家重点实验室创新基金 　作者简介 :徐振礼 ,男 ,1978 年生 ,硕士生. Email :xuzl @mail. ustc. edu. cn © 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 30 中国科学技术大学学报 　　　　　　　　　　　　　　　　　第34 卷 为例介绍有限差分 WENO 格式的重构思想[7 ] .