半拉格朗日平方守恒计算格式的研究-大气科学和地球流体力学数值.PDF

第25卷 第6期 大 气 科 学 25,No.6 2001年 11月 ChineseJournalofAtmosphericSciences NVoo1r,2001 半拉格朗日平方守恒计算格式的研究* 陈嘉滨 季仲贞 (中国科学院大气物理研究所 北京 10002引 (中国科学院大气物理研究所大气科学和地球流体 力学数值模拟实验室.北京 100029) 摘 要 作者研究了 种新的显式平方守恒的半拉格朗日计算格式 该格式是欧拉空间平 方守恒格式在半拉格朗日空间的推广和发展，并采用了守恒插值方法 此外，还给出这种格 式在一维原始方程上的应用口 关健词二半拉格朗日计算格式;完全能量守恒二守恒插值 前言 目前半拉格朗日计算格式广泛应用于数值天气预报业务，不仅提高了计算时效而且 改进了计算精度。但存在的问题是原来偏微分方程的整体守恒性得不到保证，长时间数 值积分会带来计算误差 世界上使用的守恒格式，是指在时间保持连续情况下的守恒， 典型格式是Arakawa的格式[I1在时间离散情况下 这样的守恒格式却不能保证物理 量的守恒 并且其守恒性问题仅局限在欧拉空间讨论。目前，国外研究半拉格朗日空间 的守恒格式尚不多见 本研究是应用欧拉空间显式完全平方守恒格式2[[.31和守恒插值方 法41.构造一种新的平方守恒半拉格朗日计算格式。本文将这种新格式应用到一维原始 方程上.并得到初步的计算结果 2 一维浅水波方程的分析 欧拉空间完全平方守恒格式的关键，是偏微分方程右端项的微分算子的离散空间算 子具有反对称性质，它们的平方变量 (如能量)可写为辐散形式 在拉格朗 日空间，偏 微分方程也应具有此性质 该格式才能推广 对一维重力波方程 ) l ( ) 2 ( 其能量方程为 旦厂{{一,+‘p(2、一 乒(urp2、-P,u2+,e2互 (3) dr\ 2 / Cx \ 2 J 2 CX 2000-0--2收到.2000-09-28收到修改稿 国家重4基础研究发展规划项 目G1998040900和国家 自然科学基金资助项 同资助 大 气 科 学 25卷 其中:、是速度，W是位势。显然，(3)式右端项不具有辐散形式。因此，方程 (1)和 (2)的右端项很难构造出具有反对称性质的离散空间算子 引人辐散的拉格朗日表示a[l =止兰(6x) (4} ecuxJx dr “ 其中6二表示随流体运动的无限小的长度间距。于是，在拉格朗日空间一维重力波方程 (1)、(2)可改写为 dt}(_JSx)_一嵘T}71 ) 5 ( dr}(}dx)_一疡2cx ) 6 ( 其 中 pl二40. L= lix 方程 (5)和 (6)的总能量方程为 d从2+衬、_dfu`}p十Wz、_、_ 、__户2_WZ、_ ;八42、 dr}一一万一