数值天气预报(3-4-6)--差分格式的误差分析.pptVIP

NANJING UNIVERSITY OF INFORMATION SCIENCE TECHNOLOGY 第3章 数值计算方案 * CH3 数值计算方案 3.1 差分方法和差分格式 3.2 差分格式的基本性质 3.3 时间积分方法和积分格式 3.4 有限差分格式的误差分析 3.5 非线性计算稳定性 主要内容 3.4 有限差分格式的误差分析 差分可以精确代替微分； 差分方程可以精确代替微分方程； 从而保证数值方法可以得到足够准确的解； 0时 实际工作中，不可能任意缩小步长（由于计算量过大等原因），实际计算是在有限的网格下进行的，一定程度的误差是不可避免的。 本节将分析在有限网格区域下差分格式引起的常见误差。 一、计算解问题 二、频率误差 三、波数误差 四、相速度、群速度误差 计算解是使用三个时间层积分格式所引起的一个普遍问题。以平流方程为例，分析采用三层差分格式时，计算解的产生及其特点。 一、计算解问题 对上式的空间导数取中央差格式，则有： 对于一维平流方程： （3.98） （3.99） 假定u具有 波动解形式， 相当于振幅： 令 时间取中央差： （3.80） 可以得到关于增幅因子的二次方程： 其对应的两个根为： （3.81） 对于中央差时间积分格式，其稳定性判据为： 由于 故当 时，格式是稳定的。 尽管格式满足稳定性判据，其增幅因子对应两个根， 可见，与 相联系的波，具有实际的物理意义，称之为物理波型，而与 相联系的波则不具有任何物理意义，是在计算过程中产生的虚假波型，称之为计算解。 物理解： 计算解： 其中a，b为常数。上式说明中央差分格式的数值解为两个波动的叠加，其第一个波与真解性质相近，称之为物理解，而第二个波是由于采用了中央差格式而产生的虚假波型，称其为计算解。 而计算所得的数值解为： （3.82） 因此，在设计差分格式时，总是使得该差分格式下计算解与物理解振幅的比尽量小，实际工作中通常采用提高网格分辩能力来抑制和减小计算解带来的影响。 以惯性振荡为例: 其解为： ，其中 为惯性振荡频率。 二、时间的截断误差（频率误差） 以下通过数值方法求上述问题的振荡频率的差分解（数值解），通过与 比较，说明差分格式所引起的频率误差。 差分方程形为： 设其数值解为： （ 为数值解频率或差分解频率）。 代入上式，可得： （a）中央差格式（显式格式） （b）采用梯形格式（隐形格式） 差分方程： 设其数值解为： （ 为数值解频率/差分解频率，下标表示隐式） 代入上式有： 展开之后有： 令以上等式两边虚部、实部系数相等，则有： 利用半角公式： 故有： 即为差分解频率。 可见用时间差商来代替时间微商时会引起频率误差。 图3.2时间差分格式的频率比较 对于某要素场，假设其可以表示为形如 的波动则有： 其中k为波数，或称之为微分解的波数。 三、空间的截断误差（波数误差） 令 ，对上式用中央差格式展开为： 为差分解波数或计算波数 为了反映空间差商代替空间微商时引起的波数误差，定义 来表示空间差分格式的精确度。 R与 k 和 有关； 可见，差分近似式的精度随k或 的减小而增大，这也就是说对于波长较短的波，其产生的波数误差较大；而波长较长的波，则差分方程可以比较精确地表示其空间微商，波数误差很小，精度较高。 另外：为了提高差分近似的精确度，减小波数误差，可以构造高阶差分格式。 例如对于一阶偏微商的四阶精度的差分格式： 其计算解波数为： 对于任何波长，四阶差分格式均比二阶差分格式精度高。因此，在实际工作中，缩小格距或采用较高阶的差分格式均能提高计算的精确度。 NANJING UNIVERSITY OF INFORMATION SCIENCE TECHNOLOGY 第3章 数值计算方案 *