多守恒差分格式的构造及其数值试验 - 大气科学和地球流体力学数值.PDF

简 报 第48 卷 第7 期 2003 年4 月 多守恒差分格式的构造及其数值试验 王 斌 季仲贞 (中国科学院大气物理研究所, 大气科学和地球流体力学数值模拟国家重点实验室, 北京100029. E-mail: wab@) 摘要 针对全球正压浅水波方程, 通过选取合理的表示形式, 并对半离散方程的空间差分部分加以巧 妙的修正, 使得其守恒性由原来的3 个增加到4 个, 最终构造出具有4 个守恒性的差分格式, 并在数值 检验中取得好的模拟效果. 关键词 多守恒性 差分格式 正压浅水波方程 总能量守恒是全球大气动力体系(在不考虑外源 总绝对涡度守恒的瞬时差分格式; Lilly[5]构造了具有 和摩擦的前提下) 的一个重要特征, 隐式和显式平方 总拟能和总绝对涡度守恒的瞬时差分格式等. 那么, 守恒格式就是基于这一重要的守恒性构造出来的. 能否构造出具有4 个或5 个守恒性的差分格式呢? 这 但除了该守恒以外, 全球大气动力系统还有诸如拟 是很具有挑战性的一个问题, 也是本文试图解决的 能守恒、质量守恒等很多其他的物理守恒性. 关于这 问题. [1] 一点, 曾庆存 已作过详细论证, 并强调在近似求解 1 方程的表示形式及其守恒性 中绝不能轻易将这些重要的基本守恒量加以破坏. 球坐标下的正压浅水波方程可以有多种表示形 为了在数值模拟中能很好地保持大气体系的这些重 式, 如下是其中的一种: 要的基本特性, 就必须构造多守恒差分格式. 本文选 u 1 E hv 0, 取有代表性的正压浅水波方程组作为对象, 研究多 t a cosq l 守恒差分格式的构造方法. v 1 E hu 0, (1) 正压浅水波方程组具有5 个基本的守恒性: 3 个 t a q j 1 一次守恒性和两个二次守恒性. 3 个一次守恒性分 t a cosq l uj q vj cosq 0, 别是总质量守恒、总绝对角动量守恒和总位涡守 其中, q 和l 分别为纬度和经度, a 为地球半径, u, v 和 恒, 它们都有重要的物理意义, 这在文献[1,2] 中分 j 分别为纬向风、经向风和位势高度, E 和h 为系统的 别作了阐述. 两个二次守恒性, 即总能量守恒性和 局地能量和绝对涡度. 总位涡拟能( 位涡平方) 守恒性, 在数值模拟中比起 1 2 2 E (u v ) j, 一次守恒性来说更具有实用意义. 首先, 它们是数 2