高精度迎风偏斜格式的比较与分析.pdfVIP

第 3 1 卷 第 2 期 大 　气 　科 　学 Vol3 1 　No2 2007 年 3 月 Chinese J ournal of At mo sp heric Sciences Mar . 2007 高精度迎风偏斜格式的比较与分析 冯涛1 , 2 　李建平1 1 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室 ,北京 　100029 2 中国科学院研究生院 ,北京 　100039 摘 　要 　利用一种具有任意阶精度的一般显式有限差分公式构造出高精度迎风偏斜格式 , 并利用 Fourier 分析法 评估了这些迎风偏斜格式的耗散误差与频散误差 。结果表明 , 偶阶精度格式的数值相速度快于实际相速度 , 而奇 阶精度格式的数值相速度慢于实际相速度 。并且 , 偶阶精度格式的耗散误差与频散误差低于相邻的奇阶精度格 式 。为了检验这些格式的计算性能 , 在一维问题上进行了应用 。首先 , 考虑恒定风场条件下的一维平流试验 。主 要选择两种不同的初始条件来评价数值格式的精度 , 这两种试验问题分别是高斯函数 、方波函数 。试验结果表 明 , 随着数值格式精度的提高 , 数值格式的误差逐渐减小 。而对于高于六阶精度的格式来说 , 改进的程度并不是 很大 。其次 , 应用各阶格式到具有两种不同初始条件的无粘 Bur ger s 方程 。数值结果表明 , 随着数值格式阶数的 增加 , 数值结果也得到了明显改进 。而对于高于六阶精度的格式来说 , 进一步的变化并不明显 。总之 , 在兼顾效 率与精度条件下六阶迎风偏斜格式是最好的。 关键词 　显式有限差分公式 　迎风偏斜格式 　平流方程 文章编号 　1006 9895 (2007) 02 0245 09 　　　中图分类号 　P435 　　　文献标识码 　A A Comparison and Analysis of High Order UpwindBiased Schemes F EN G Tao 1 ,2 and L I J ianPing1 1 S t ate Key L aborat ory of N umerical M odel i ng f or A t mosp heric S ciences an d Geop hy s ical Fl ui d Dy nam ics , I ns t it ute of A t mos p heric P hy s ics , Chi nese A ca demy of S ciences , B eij i ng 　100029 2 Gra d uate Uni vers ity of t he Chi nese A ca demy of S ciences , B eij i ng 　100049 Abstract 　High or der up windbia sed schemes ba sed on t he general exp licit difference formulas wit h arbit rary or der accuracy f rom Li (2005) are develop ed The di ssip ation and di sp er sion error s of t hese up windbiased schemes are as sessed by mean s of t he Fourier analysi s The result s show t hat t he numerical p ha se sp