数值天气预报(4-2-5)--第2节正压原始方程组的线性计算稳定性.pptVIP

第4章 正压原始方程模式 * * * * * * * * * §4.2 正压原始方程组的线性计算稳定性 在前一章中，已经讨论过线性平流方程的计算稳定性问题。它是确定线性方程时间显式格式积分步长的重要依据；对于非线性方程，一般也是根据其线性计算稳定性的要求，来确定时间显式格式积分步长的。 本节将进一步讨论正压原始方程组的线性计算稳定性问题并得出其判据，以便为数值求解正压原始方程模式时，提供一个选取时间积分步长的依据。 式中 分别为纬向平均的风速分量和平均位势高度， 为扰动量。上式略去了科氏力项。 如果令 ， 为重力外波波速，上述第一式除 ；第二式除 则可得： 正压原始方程组的线性计算稳定性的分析 为了讨论的方便，仅考虑一维运动并进行线性化: 线性化后 运动方程 连续方程 第一式除 ;第二式除 ,则上两式可写成为： （4.24） （4.25） 将上述两式改写为差分方程，其时间、空间均采用中央差分格式，可得如下差分方程： 其中d 和?t分别为空间步长和时间步长。 式中 将上两式相加，并令： （4.26） 可得： （4.26） 设差分方程的波动解为： （4.27） 式中 为F 的常值波振幅，G 为振幅放大因子，而 k 为纬向波数。 请同学们自己思考 与动力气象学中波动的形式解有什么类似的关系？ 下面用冯-纽曼方法来讨论并得出其线性稳定性判据。 将波动解代入差分方程： （4.28） （4.26） 将上式除以 （4.28） 可求解出： （4.29） 此时，差分方程的解和相应的差分格式均是计算稳定的。 在满足 情况下： （4.31） 请同学们自己思考并证明当 解和相应的差分格式均是计算不稳定的. 第4章 正压原始方程模式 * * * * * * * * *