总能量守恒格式的构造及检验-空气动力学学报.pdfVIP

第 卷 第 期 空 气 动 力 学 学 报 ， ! # $% ’! (% ’# 年 月 !#! !$%’()!*+! ,+)+! ， !! )! *+, ’ !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） !-./).!- !! #/#)0/0 总能量守恒格式的构造及检验 季仲贞，赵 颖，王 斌，林万涛 （中国科学院大气物理研究所大气科学和流体力学数值模拟国家重点实验室，北京 )!1 ） 摘 要：由于大气海洋方程大都带有拟线性形式的对（平）流项，本文依据一个有关格式构造 的定理，提出一个构造总能量守恒格式的新方法，不仅使以往构造的许多平方守恒格式成其 为特例，而且为大气海洋数值模拟提供了灵活地选择计算格式的机会。文中给出的算例表明 所构造格式的守恒性很好。 关键词：大气海洋方程；总能量守恒；数值算法；计算稳定性 中图分类号： 文献标识码： 234)’! 5 - 引 言 在忽略摩擦和外源强迫的条件下，大气可视为一个保守系统。把地球大气作为一个 整体，研究其整体性质是十分重要的。通过这些整体性质，既可以简单而清晰地揭示出外 界因子对整个大气运动的基本影响，又可以推出大气运动内部各物理量之间的最一般的 关系，从而揭示出有关大气运动的最基本的性质。在文献［］中，曾庆存详细地论证了全 ) 球大气运动在不考虑强迫和摩擦的情况下，具有若干重要的整体性质，即其总能量、总质 量、总涡度、总拟能及总角动量是守恒的。然而，当人们作具体的气候数值模拟或数值天 气预报时，不可避免地要采用近似的计算方法。根据当代计算方法研究的一条不成文的 ［］ ! 基本法则 问题原型的基本特征在离散化后，应尽可能地得到保持 ，否则，就会给我们的 / 物理系统引入虚假的外来影响，或者使最基本的物理定律遭到破坏，那时，就不可能期望 得到较准确的气候模拟和天气预报结果了。 正是基于以上的考虑，尤其是针对大气海洋方程非线性和复杂性特征，为了保证其数 值计算能长时间进行下去，计算稳定性、计算准确性和计算省时性成了三个需同时关注的 重要问题。大量的理论分析和数值实验结果表明，保持总能量守恒是整个问题的一个重 要关键。过去的实践，已使我们这里形成了一套总能量守恒的算法，针对不同类型的问题 相应地构造了各自保持总能量守恒的差分格式（在线性情形和若干简单情况就是所谓的 ［ ］ 4/. 平方守恒格式） 。 下面，我们从另一角度出发，给出构造总能量守恒格式的一个新方法。 收稿日期： ； 修订日期： !!/#/! !!/0/4 ’ 基金项目：国家重点基础研究发展规划项目（编号6)1114!.)），国家杰出青年科学基金项目（编号#1.!-)1 ） 和国家自然科学基金（编号#113-! ）共同资助。 万方数据 作者简介：季仲贞（ ），研究员，方向：计算地球流体力学 )1#/ 7 空 气 动 力 学 学 报 第 卷