计算地球流体力学：第4讲 有限差分法稳定性分析.pptx

4.1 Fourier 方法 ( 或Von Neumann方法 )4.2 直接方法(或矩阵方法)4.3 Hirt 启示法4.4 能量分析法(或范数法)4.5 数值算例;4.1 Fourier 方法 ( 或Von Neumann方法 )4.2 直接方法(或矩阵方法)4.3 Hirt 启示法4.4 能量分析法(或范数法)4.5 数值算例;Fourier 方法 判别准则 算例;Fourier 方法 判别准则 算例;以一维线性平流方程初值问题;Fourier 方法;扩充离散函数定义域;Fourier 积分变换;由此得出：;由此得出：;上式中因子 称为增长因子(或增幅因子，放大因子等)。;上式中因子 称为增长因子(或增幅因子，放大因子等)。;上式中因子 称为增长因子(或增幅因子，放大因子等)。;Fourier 方法 判别准则 算例;上式中因子 称为增长因子(或增幅因子，放大因子等)。;上式中因子 称为增长因子(或增幅因子，放大因子等)。;Von Neumann 条件;正规矩阵;判别准则;判别准则;判别准则;Fourier 方法 判别准则 算例;例 4.1 考虑逼近一维线性平流方程;（2）时间向前差，空间向后差格式：;（2）时间向前差，空间向后差格式：;（2）时间向前差，空间向后差格式：;（2）时间向前差，空间向后差格式：;所以，各式的增长因子分别为;所以，当 时，有 ，即 Von Neumann 条件满足，格式是稳定的。;所以，当 时，有 ，即 Von Neumann 条件满足，格式是稳定的。;（3）时间向前差，空间中央差格式：;例 4.2 考虑逼近一维线性平流方程;例 4.2 考虑逼近一维线性平流方程;例 4.2 考虑逼近一维线性平流方程; 这是一个三层格式，讨论这种类型的差分格式的稳定性，一般先化成与其等价的二层差分方程组：;如果令 ，则上述方程组可以写成：;如果令 ，则上述方程组可以写成：;如果令 ，则上述方程组可以写成：;所以，增长矩阵为;下面求增长矩阵的特征值;下面求增长矩阵的特征值;当 时，为方便起见，取 ，则增长矩阵;当 时，为方便起见，取 ，则增长矩阵;4.1 Fourier 方法 ( 或Von Neumann方法 )4.2 直接方法(或矩阵方法)4.3 Hirt 启示法4.4 能量分析法(或范数法)4.5 数值算例;直接方法(或矩阵方法);直接方法举例;采用??列显式差分格式来逼近;先把差分方程(4.2.3a)写成;先把差分方程(4.2.3a)写成;其中 ，可以把(4.2.4)写成向量形式，即;其中;关键是求出S 的特征值。设 和 分别为S 的特征值和特征向量。;关键是求出S 的特征值。设 和 分别为S 的特征值和特征向量。;写成分量的形式：;其中 为矩阵S 的元素，所以 。设(4.2.6)的解具有形式：;解为：;显然 ，;如果采用的是隐式差分格式;如果采用的是隐式差分格式;如果采用的是隐式差分格式;可以把(4.2.7)写成向量形式：;4.1 Fourier 方法 ( 或Von Neumann方法 )4.2 直接方法(或矩阵方法)4.3 Hirt 启示法4.4 能量分析法(或范数法)4.5 数值算例;Hirt 启示法;Hirt 启示法举例;在格点 处进行 Taylor 级数展开，有;在格点 处进行 Taylor 级数展开，有;即;即;在格点 处，由差分方程得到;由平流方程可以得到;略去高阶误差项，得到第一微分近似;当 时，方程化为原来的平流方程。否则为抛物型方程，要有意义，必须;当 时，方程化为原来的平流方程。否则为抛物型方程，要有意义，必须;当 时，方程化为原来的平流方程。否则为抛物型方程，要有意义，必须;Hirt 启示法举例;4.1 Fourier 方法 ( 或Von Neumann方法 )4.2 直接方法(或矩阵方法)4.3 Hirt 启示法4.4 能量分析法(或范数法)4.5 数值算例;能量分析法(或范数法);离散空间及其范数;非负算子;主要稳定性定理; 定理 4.6 若 为非负算子，则当 时，差分格式(4.4.2)是无条件稳定的