双曲守恒律方程WENO格式的优化方法-core.pdf

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