非线性发展方程的非守恒格式的计算稳定性问题.pdfVIP

第28 卷 第4 期 大 气 科 学 Vol28 　No 4 2004 年 7 月 Chinese Journal of Atmospheric Sciences July 2004 非线性发展方程的非守恒格式的 计算稳定性问题 林万涛　谢正辉 ( 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室, 北京　100029) 摘 　要 　　针对非线性发展方程的非守恒格式 , 以一维浅水波方程为例, 对非守恒格式的 计算稳定性进行了研究分析 , 探讨了非线性发展方程的非守恒格式与初值的关系。理论分 析和数值试验表明 , 在格式结构已经确定的情况下 , 非守恒格式的计算稳定性主要由初值 的形式所决定。 关键词 : 非线性发展方程 ; 非守恒格式 ; 计算稳定性 ; 初值 文章编号 (2004)　　中图分类号　P413 　文献标识码 　A 1 　引言 近代研究表明 , 大气和海洋的运动、发展与变化遵从非线性规律。描述大气和海 洋运动的方程是一类比较复杂的非线性偏微分方程 , 通常称之为非线性发展方程。基 于描述大气和海洋运动的非线性发展方程本身的复杂性 , 很难求得解析解 , 由此就必 须将其进行离散化后数值求解。 大气和海洋系统短期运动的主要特点是能量守恒。对这类问题进行数值求解时使 用最多的是有限差分格式 , 因此 , 能否设计出长时间计算稳定的差分格式是一个关键 性问题。正是基于这种情况 , 曾庆存等[1 ,2 ] 和季仲贞[3 ,4 ] 系统地研究了绝热或无耗散情 况下非线性发展方程的计算稳定性问题 , 探讨了产生非线性计算不稳定的原因 , 首先 构造了计算稳定的隐式完全平方守恒差分格式。之后 , 季仲贞等[5 ,6 ] 和王斌等[7～10 ] 又 设计发展了计算稳定的显式完全平方守恒差分格式。最近 , 陈嘉滨和季仲贞[11 ] 研究了 保持总能量守恒的半拉格朗日平方守恒格式。季仲贞等[12 ] 对总能量守恒和辛几何算法 之间的关系也进行了讨论。林万涛等[13 ] 分析了大气和海洋短期运动与守恒及非守恒格 式的关系 , 进一步证实了对大气和海洋系统的短期运动问题用平方守恒格式进行数值 求解是稳定的 , 而用非守恒格式则是不稳定的。所以用平方守恒格式解决这类问题有 更多的优势。 大气和海洋的中、长期运动一般情况下能量是不守恒的。由于大气和海洋运动本 ( ) 质上的非线性 , 其中、长期运动状态依赖于基态 初始场 。基态的改变会导致大气和 海洋运动的改变[14 ] 。在中、长期数值天气预报和海流数值模拟中, 对非线性大气海洋 收到 ,收到修改稿 中国科学院重大创新项目 KZCX2208 、中国科学院百人计划项目“气候与植被相互作用”和国家自然科 学基金资助项共同资助 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　4 期 林万涛等 : 非线性发展方程的非守恒格式的计算稳定性问题 511 动力学方程组同样需要使用有限差分格式进行数值求解。因此 ,