误差对大气数值模式设计的影响研究.pdfVIP

误差对大气数值模式设计的影响 廖洞贤 (国家气象中ti，北京100081) 一、一般的考虑 把一个偏微分方程化作代数方程或代数方程组时，会引人两类误差:一是由于取近似 引起的，另一类是进行离散引起的。虽然近似和离散的精度已为人所熟知，但误差在方程 中的作用，特别是它们对方程各项的影响却很少有人讨论。为了说明问题，设大气模式方 程组可以写成 daFt,·}i,`Fv二，“二‘’，，，” (1) 其中F,是预报变量，凡是对aF;ia有‘贡献的项。 设(1)中除第一项外对各项取近似，则 aF; +ZFyi+名Eyi二。 (2) at 其中Fy二Fy,+Ey,。对上方程用差商代替徽商，则 餐·岑iFv}z·耳。“·丫EyzEz;二。 (3) 其中SF;/at,Fq,洛是aFiat和凡，的差商;Ei2和Eo,各是8F/和与2的截断误差。令 max(EIyom=E;,,min(介，)，二F,; max(Eq,,Eiyz,E,2)，二;，‘·(FS,;二。，=Ftu 则不等式 瓦，F;, (4) 和 E;凡12 (5) 必须满足，否则在方程(2)中考虑F,1项或在(3)中考虑Fn2项没有意义;而且，在(2)中考 虑所有Fy，中不大于E;，的项，或在(3)中考虑所有Fy,2中不大于的项也没有意义。 二、取近似的情况 1.取互力加速度 ‘为常数 在一般情况下这样处理影响不大，但并不是在任何情况下都不大。 8二So(1+0.005302sinzp一0.000078i内9)(1一a,z) (6) al二3.14x10-m- 利用上式，定义E二S一9.8。一?，可以得到图1(图略)。 ·29 · 从图中可以看出:在10000M高度，E可以达到5x10-z二一，，一般也有10-2ms2，这 和中一Y运动中aw/at等的量级相当，显然，不考虑8随，和:的变化，要算好aw/at等 是困难的。 2.用地球平均半径a代替: 根据计算，在对流层下部，在大尺度运动中。水平运动方程各项误差中最大者为气压 梯度力项的误差，约10-6.-1，比二/a项为大，和2f2cosg。相当。故考虑这两项没有 意义。 三、离散的情况 设气压P可以表示为 P二PH+P, (7) 其中Pu是静力部分，P,是其偏差。再设气a随z的分布是 T 一几， 0鉴 ; ‘rZ 一 一 T T (S) T,一几r, Zr ‘z. 则静力方程成为 了 卫 . . ， 0 ‘:蕊zr 1 益 ‘ ﹂ 一 瞥 L 旗(’一”“ z, ‘z 仍 其中下标: