数值天气预报(3-5-5)--第5节非线性计算稳定性.pptVIP

NANJING UNIVERSITY OF INFORMATION SCIENCE TECHNOLOGY 第3章 数值计算方案 大气科学学院 陈海山 * NANJING UNIVERSITY OF INFORMATION SCIENCE TECHNOLOGY 第3章 数值计算方案 大气科学学院 陈海山 * CH3 数值计算方案 3.1 差分方法和差分格式 3.2 差分格式的基本性质 3.3 时间积分方法和积分格式 3.4 有限差分格式的误差分析 3.5 非线性计算稳定性 主要内容 3. 5 非线性计算稳定性 前面讨论了线性偏微分方程差分格式的计算稳定性问题，而数值天气预报用的往往是一组非线性偏微分方程。 对于非线性偏微分方程，线性偏微分方程的稳定性条件，只能给出其计算稳定性的必要条件，即使满足这一条件，也可能会因为差分方程的边界条件和非线性项的不正确表示而产生计算的不稳定现象，我们把这种由于非线性作用而产生的不稳定，称为非线性不稳定。 本节将通过一个计算实例的分析，来说明非线性计算不稳定的产生，并进一步分析其产生的原因。 一、计算实例： 对一维平流方程： （3.131） （3.132） 第一种形式 第二种形式 上述方程可构造不同的差分方程，下面给出两种方案： 其中： 〖1〗 〖2〗 （3.133） （3.134） 给定两种不同的初值： 〖1〗 （3.135） 〖2〗 （3.136） 此时差分格式是满足线性计算稳定性条件的。 假设 根据前面有关线性稳定性的讨论，其线性稳定性判据为： 对两种差分格式，分别用以上给定的初值进行计算。为了考察二者的稳定性，对每一步的格点总动能： 进行了计算，考察其总动能随时间的变化。 ?图3.4 差分格式总动能随时间的变化曲线 （a对应初值1，b对应初值2；实线-方案1，虚线-方案2） 以上分析表明： ??? 在满足线性稳定性的条件下，非线性差分方程仍然会产生不稳定现象，而且不稳定具有突变的特点。 ????非线性计算不稳定的产生和差分格式有关，初值对计算不稳定具有显著的影响。 以上简单说明了非线性计算不稳定的产生，下面进一步分析其产生的原因。? 有关非线性计算不稳定产生的可能原因，Philips（菲利普斯）曾指出，由于有限网格不能正确地分解短波，从而造成混淆现象，而混淆误差可能会引起计算的不稳定。 差分方法是用有限的自由度系统来代替原来的连续介质系统的。而有限格点上的函数值只能分解有限的波数，其最短波长为 ，对于非线性作用产生的波长小于 的波动，网格系统不能正确地将它表示出来，而把它错误地表示成为某种波长大于 的波，从而产生了误差，我们把这种波的误差称之为混淆误差。 二、混淆误差-非线性计算不稳定产生的可能原因 我们把长度为 的空间域分为I等分，即 ， 为空间步长，可得I+1个空间格点。函数u在格点上的值可表示为有限级数的形式： （3.137） 由于 ，上式也可写成： 其中为k波数，而 为网格系统可以分辨的最大波数。 混淆误差是产生： 假设 中包含两个波：波数分别为 和 。 由于非线性作用 上式表明：由于波数为 和 的两个谐波的非线性相互作用，有限网格不能正确分解波数为 的波，而是把它错误地表示为 S 的分量。 假设 如果引进S 此时 可得到 （3.138） 使原有的波数为 的波的能量发生了 变化，这是由于采用了差分近似而产生的现象，并不反 映真实的物理过程。 取 I=12 则有限网格可以分解的最大波数为I/2=6 设 波数 （实线） 只能表示为波数为 图3.5 混淆误差示意图（其中实线表示波数为 的波分量，而虚线表示波数为 S 的波分量） 假设， ，则 ，此时二倍格距波长的波与四倍格距波长的波发生相互作用，有能量反馈到波长为四倍格距的波中去。 如果这种过程不断重复，则四倍格距的波能量就会迅速增长，这种虚假的计算现象可能会使得计算结果迅速增长而产生计算的稳定，称之为非线性计算不稳定。 可以证明：只有对波长很短的波才发生能量反馈；非线性不稳定的产生主要是由于短波能量的虚假增长所造成的。 NANJI