差分精度对水平变量配置结果的影响研究.pdf

差分精度对水平变量配置结果的影响* 田福厚 刘宇迪 总参大气环境研究所 北京 100029 摘 要 本文采用一种推导频散关系的通用方法 研究在Arakawa AE网格上采用二阶中央差或四阶紧致差分 格式描述惯性重力波产生误差的情况 分可分辨和不可分辨两种情况从频率和群速两个方面进行了对比分 析 结果表明不论在可分辨或不可分辨 只有Arakawa C网格当采用差分精度较高的紧致差分格式时 模拟 惯性重力波的精度才会提高 而其它的四种网格 Arakawa AB D E 当差分精度提高以后 模拟惯性 重力波的精度反而下降 关键词 中央差 惯性重力波 紧致格式 差分精度 水平网格 1 引言 1 2 关于变量配置的水平研究早在上个世纪60年代Winninghoff就有过论述 后来刘宇迪 9按照Winninghoff的思路对到目前为止出现的所有水平网格从模拟惯性重力波和Rossby波 的频率和群速等角度进行了广泛而深入地研究 所得结论较前人有所发展 并且还设计了 两种新的垂直跳点网格和几种三维网格 同时对垂直网格和三维网格的计算频散性也进行 了深入地研究 但他们前面的研究是基于二阶中央差格式的 如果差分精度再提高 即采 用4阶的格式如何呢 是否会改变二阶精度时所得的结论 为此 本文采用目前较为流行的 1016 紧致差分格式 进行研究 不过本文的研究思路较前人有所突破 不再象以前那样每种网 格写出一个差分方程进而导出其频散方程和群速方程 而是抽象出一个通用的差分格式 这个格式适用于任何差分精度和任何网格类型 也适用于微分情况 然后分可分辨和不可 分辨两种情况 以及分别采用二阶中央差和四阶紧致差分格式的差分精度 用图示的方法 比较了Arakawa AE网格描述惯性重力波的频率和群速方面产生的误差情况 下面首先看一下推导惯性重力波频散关系的通用方法时采用的一些算子 2 一维空间中的算子 如果网格点x i i N (N为网格点总数)上的变量f 用f i 表示 则 算子S 表示 网格的中间点x =(x +x )/2处 变量f 的近似值f 0 i 1 2 i i 1 i 1 2 算子S1 2 表示 网格点xi 1 2 处变量f 一阶导数f 的近似值f i 1 2 / 算子S 表示 网格点x 处变量f 一阶导数f的近似值f 1 i i 为比较不同的差分精度下变量配置的情况 本文选择了两种精度的差分格式 一种是 一般的二阶中央差格式 用C2表示 它是数值模式中最常用的格式 第二种是四阶紧致差 分格式 用C4表示 这两种格式的S S 和S 见表1 0 1 2 1 表1两种格式C2和C4的算子S0 S1 2 和S1 算子 具体形式 C2 S0 f i 1 2 (f i f i 1 ) / 2