半隐式半拉格朗日平方守恒计算格式的构造.pdfVIP

第28 卷 第4 期 大 气 科 学 Vol28 　No 4 2004 年 7 月 Chinese Journal of Atmospheric Sciences July 2004 半隐式半拉格朗日平方守恒计算 格式的构造 1) 2) 陈嘉滨 　季仲贞 1) ( 中国科学院大气物理研究所, 北京 100029) 2) ( 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室, 北京 100029) 摘 　要 　　在显式半拉格朗日完全平方守恒格式基础上 , 构造出半隐式半拉格朗日完全平 方守恒计算格式 , 它继承了显式半拉格朗日完全平方守恒格式的优点 , 并突破计算不稳定 柯朗条件对时间步长的约束 , 使时间步长大为增长。此外 , 还给出这种新的计算格式在一 维原始方程上的应用。 关键词 : 半隐式半拉格朗日格式 ; 平方守恒 文章编号 (2004)　　中图分类号　P435 　　文献标识码 　A 1 　引言 半拉格朗日平流格式 , 已被广泛用于天气和气候预报和摸拟, 这是由于它与欧拉格 式相比在节省计算工作量和计算精度方面有明显优越性。关于半拉格朗日平流格式和 有关问题 , 早在 1991 年 Staniforth 和 Cote 就有一个很好的评论[1 ] 。然而 , 与欧拉格式不 一样 , 半拉格朗日平流格式普遍存在的问题是 , 由于内插造成物理量, 如总质量、总能 量等缺乏守恒性 , 这给长时间数值积分带来误差。针对这一守恒性问题 , 10 多年来 , 有 不少人致力于这方面的研究[2～6 ] , 他们研究半拉格朗日质量守恒格式的构造 , 并应用到 二维和三维问题。主要想法是 : 对于一个流体微团 (用坐标格点表示为格网元) 从出发 ( ) 点移动到标准网格点上 , 在质量守恒条件下计算这个微团 格网元 在出发点和到达点 的质量。不同作者对标准网格内质量分布有不同假设 , 例如逐段常定、逐段线性、逐段 ( ) ( ) 抛物线 PPM 、逐段立方曲线 PCM 等。 本文构造的守恒格式的思路和前面叙述的不同之处在于 , 本文中的控制方程引入 耗散算子 , 在每步计算中将多余的物理量 , 如质量、能量等耗散掉 , 使之总物理量保持 守恒。 在文献 [7 ] 中 , 欧拉空间显式完全平方守恒格式[8 ] 已被发展和推广到半拉格朗日 空间 , 并放松对 B 算子为正定假定的要求 , 使格式更为灵活。 但需指出的是 , 这种半拉格朗日平方守恒计算格式 , 仅具有理论意义。因为 , 在此 格式中 , 对产生重力波的项 , 例如气压梯度项和辐散 (合) 项 , 仍然取显式格式 , 这就 收到 ,收到修改稿 中国科学院知识创新工程项目 ZKCX2SW210 、国家重点基础研究发展规划项目 G1999032801 和 2001BA603B4 共同资助 528 　　 大 　　气 　　科 　　学 28 卷 　 限制了时间步长的拉长。 在本文中 , 我们将研究半拉格朗日平方守恒的半隐式格式。这就可以突破重力波对 时间步长的约束 , 增大时间步长 , 节省计算机时间 , 使格