群速度直接控制四阶迎风紧致格式Core.PDFVIP

( ) A 　辑 　　　　　　 　 中　国　科 　学 　　　 第 31 卷 　第 6 期 SCIENCE IN CHINA ( Series A) 　 2001 年 6 月 群速度直接控制四阶迎风紧致格式 马延文① 　傅德薰② ( 中国科学院力学研究所 ①高温气体动力学开放实验室 ; ②非线性力学国家重点实验室 , 北京 100080) 摘要 　　为求解多尺度复杂流动在有限网格点的情况下高精度方法是可供选择的方 法之一 ,其中紧致型格式数值解的分辨率更好. 用迎风型紧致格式计算激波时 ,数值 解中仍有数值振荡产生 ,这是由于对应于不同波数之群速度不均一所致. 采用群速 度直接控制方法重构紧致型格式以提高捕捉激波的能力. 该方法简单 ,精度高 ,网格 基架点少. 算例表明 ,这一新的方法是可行的. 关键词 　　紧致型格式 　高分辨率激波解 　群速度控制 紧致型高精度格式受到人们的关注. 在同样逼近精度条件下 ,紧致型格式更高的分辨率 , 且所需网格点较少. 对于求解多尺度复杂流动 , 要求数值方法能很好地捕捉到高波分量 ,要 求高波量的耗散误差及色散误差都比较小. 在计算资源有限的情况下 ,采用高精度格式是有 益的[ 1～5 ] . 采用现有的一些格式 ,在计算激波时数值解中将有数值振荡产生. 为克服数值振 荡已发展了很多有效方法 ,其中包括 TVD ,NND 及 ENO 等类型的格式[6～12 ] . 这些格式已成功 地用于求解实际问题. 这些方法的发展都是从简单的一维问题出发进行研究 ,并在实践中得 到了验证 ,但紧致型对激波有高分辨率的格式并不多见. 现有的一些高分辨率格式多是从数学观点出发而构造的. 在这些方法中没有很好地研究 数值解中产生振荡的原因 ,也不是针对数值振荡产生的直接原因而对格式进行改造的. 文献 [ 11] 中对个别具体格式通过对应数值解的群速度 ,分析了数值解的行为. 文献[ 5 ] 中参照数值 ( ) ( ) ( ) 解的群速度把格式分为 3 类 :快型 FST , 慢型 SLW 和混合型 MXD . 对于慢型格式 ,数值振 荡出现于激波之后 ,对于快型格式数值振荡出现于激波之前 ,对于混合型格式数值振荡出现于 激波之两侧 ,但主要出现于激波之前. 文 中将 引入一个函数用 以控制数值解 中的群速度 (group velocity control ,简称 GVC) ,使得在激波之前格式为 SLW 型的 ,在激波之后主要表现为 FST 型的 , 以此来达到控制数值解中非物理振荡的目的. 1 　依赖于自由参数的紧致差分 为简便起见 ,考虑如下模型方程及对应的半离散化方程 : 　　收稿 ( ) ( ) 　　 国家自然科学基金 批准号 : 、科学技术部“九五”预选及“九七三”基金 批准号 : 1999032805 资助项 目 　 Email : fud @cc5 . imech . ac . cn 第 6 期 马延文等 : 群速度直接控制四阶迎风紧致格式