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

第4章 正压原始方程模式 * * * * * * * * * * * §4.5 正压原始方程模式的时间积分格式 在数值预报模式的空间差分格式确定之后，通常根据数值天气预报模式种类、预报或模拟的天气系统的特点来选择较为适当的时间积分格式。 前面曾对普遍的时间积分格式进行过较为详细的讨论。本节将给出针对正压原始方程模式常用的时间积分格式，进行进一步讨论。 由已知时刻 t 或 t-?t 的变量值，求 t+?t 时刻的变量值。主要有中央差格式和欧拉--后差格式。 一、显式格式 式中c为一常数，时间中央差格式为： 1、中央差格式 仍以一维平流方程为例： （4.84） （4.85） 或 前面曾指出，采用中央差时间积分方案，存在计算解的问题。为了抑制计算解对结果的影响，通常采用 “三步法”的时间积分方案，来减小初始时刻计算解的振幅，并达到抑制计算解对最终结果的影响。 （4.86） “三步法”的具体做法是：时间积分起步时，先用 ?t/2向前差；再用?t /2中央差；以后一直用?t 做中央差 对于非线性方程，由于不同解之间的非线性相互作用，当物理解的振幅发生改变时，通常计算解的振幅也会发生变化。随着时间积分的延长，计算解的振幅会变得越来越大，在实际计算中通常采用时间平滑来抑制计算解的增长。其平滑公式为： (4.87) 式中S为时间平滑系数。如果取S=1/2，则可以阻尼高频振荡，可以完全滤去周期为2?t的计算解。 欧拉--后差格式由欧拉格式和后差格式结合而成，以平流方程为例，先用欧拉格式： 2、欧拉--后差格式 （4.88） 计算出 作为第一近似值,然后用时间后差格式: （4.89） 可见，这种格式每积分一步相当于进行一次迭代，所以又称之为迭代格式。 在以上两式中消去 ,可以把它写成如下形式： (4.90) 上式表明，欧拉--后差格式相当于在时间前差 格式上加一个平滑系数为 的空间平滑运算 ，对于快速移动的高频波，系数很大；而对于移动 较慢的低频波而言，该系数则很小。 欧拉--后差格式具有如下特点：1、计算很稳定；2、计算简单；3、不存在计算解的问题；4、具有选择性衰减；但是，5、计算精度较低；6、有线性计算稳定性条件的限制，积分的时间步长不能取得太长，因而计算耗时多。 由于以上特点，通常很少在数值积分过程中，单独长时间使用这种格式，而经常是与其他格式结合使用。 对于正压原始方程模式的时间积分方案，可在时间积分的最初阶段，采用欧拉--后差格式，来阻尼由于初始资料不平衡而引起的高频振荡；积分一段时间后，再改用“三步法”起步的中央差格式进行积分。 上面介绍的两种格式，预报方程的右端量均是已知量，可直接计算出预报要素，这两种格式都是显式格式。 尽管这两种格式计算简单，积分一步所耗的时间较少，但由于这两种格式受线性计算稳定性条件的限制，时间步长不能取得太长，因而整个时间积分过程计算步数仍然很多。 第4章 正压原始方程模式 * * * * * * * * * * *