数值天气预报(4-4-4)--第4节正压原始方程模式的空间差分格式.pptVIP

如果应用恒等式： (4.48) 可以把 (4.47) 式写成： (4.49) 这是平流方程的二次守恒格式，称为半动量格式或松野格式。 N X Y O y x 局地直角坐标系 网格直角坐标系 （这是一种地图投影坐系） 北极点 o a a 二、正压原始方程模式守恒差分方程组的构造 1、正压原始方程模式计算网格区域的选取 第4章 正压原始方程模式 * * * * * * * * * * * * * * * * * * * §4.4 正压原始方程模式的空间差分格式 原始方程模式中含有快波解，其积分的时间步长需取得较小，使计算的非线性不稳定问题较为突出。 通常采用两种方法来解决： 一、在数值积分过程中采用空间平滑； 二、设计守恒的空间差分格式。 本节主要讨论守恒的空间差分格式构造问题。 一个合理的有限差分方程应该保持对应微分方程的物理特性，这样才能保证差分方程的数值解能沿着正确的途径逼近微分方程的解。 通常把按照一定的守恒性质构造的空间差分格式，称为这种性质守恒的空间差分格式。 例如：正压原始方程组，具有全球总能量守恒的积分性质。所以，应当构造一个差分方程组，使其在全球范围内的总能量也保持守恒。这种差分格式称为总能量守恒格式。 一、守恒的空间差分格式的构造 守恒空间差分格式的概念： 来具体讨论守恒的空间差分格式的构造。 现以简单的通量方程和平流方程为例： （4.34） （4.35） 式中Ai，j 为离散网格点上的任一标量函数，上标表示求平均，下标表示进行差分运算,而d 为网格距。同样,在y 坐标方向上也有类似的差分算符。? 为了方便,引进一些常用的有限差分算符(X坐标方向): （4.36） a、一次守恒的通量格式 1、一次守恒格式 对某一封闭区域或一有周期边界条件的区域 S 积分 上式表明：对微分方程来说，物理量F 在S 区域内的平均值是守恒的。 如果将通量的微分方程改写成差分方程，在S 区域内对所有格点求和，也能得到与(4.37)式相类似的结果，我们称这种差分格式为一次守恒格式。这是因为它保持了相应的微分方程对F 本身的积分性质。 (4.37) 如采用中央差格式将(4.34)式写成： (4.38) 此时，对图2所示的网格区域S 内的各网格点求和：考虑用边界条件： 图4.2 网格区域示意图 a、 b、 如果网格距d 取得很小,边界 条件b 则相当于区域S边界上的法向风速为零，即区域S是封闭的。 1 2 3 4 式中，?S=d 2 为面积元。 可以得到： (4.39) 这表明(4.39)式与(4.37)式具有相同的性质，证明了采用简单的中央差格式，能使区域S 内物理量F 的总值（平均值）守恒。所以(4.38)式是通量方程(4.34)的一次守恒格式。? 1 2 3 4 (4. 39) 应用一次守恒的通量格式，可以证明以下差分格式 是平流方程的一次守恒格式。 （4.42） b、一次守恒的平流格式 对某一封闭区域或一有周期边界条件的区域S 积分 新的问题： 实际的数值计算表明，如果不采用附加的平滑运算，一次守恒差分格式的计算稳定性能较差，原因在于物理量F 的平均值守恒并不能保证其绝对值不无限增长。 为此，需要进一步来设计能使物理量 F2 守恒的所谓的二次守恒差分格式。? 如果将通量方程（4.34）的有限差分方程乘以F并对区域S 内的所有网格点求和也能得到与上式相类似的结果，则因它保持了相应微分方程对 F2 的积分性质，称这种差分格式为二次守恒格式。 2、二次守恒格式 a、二次守恒的通量格式 乘以F （4.43） 对区域S 积分： 是通量方程（4.34）的二次守恒格式。 下面来证明形式如下的差分格式 （4.44） 将上式乘以F，并对区域S 内的各格点求和，应用前面给出的边界条件b，整理后可得： （4.45） 显然,(4.45)式与(4.43)式相似，它保持了相应的微分方程对物理量F2 的积分性质。因此，差分格式(4.44)式是通量方程(4.34)的二次守恒格式。 写成： 同样，可以证明差分格式 （4.47） 是平流方程