第十四讲 平方守恒差分格式-协调耗散算子课件.pptx

第十四讲 平方守恒差分格式-协调耗散算子第十三讲回顾平方守恒格式 变步长格式定步长格式完全平方守恒隐式瞬时平方守恒？显式Y第十三讲回顾-显式平方守恒格式定步长显式平方守恒格式：变步长显式平方守恒格式：第十三讲回顾-显式平方守恒格式（定步长）（变步长）其中14 协调耗散算子协调耗散算子关于平方守恒格式的基本定理正压浅水波方程的显式平方守恒格式14 协调耗散算子协调耗散算子关于平方守恒格式的基本定理正压浅水波方程的显式平方守恒格式耗散算子是两个关键环节之一显式平方守恒格式的构造有两个关键的环节：在空间离散方面，反对称离散算子的构造十分关键；在时间离散方面，从前面的讨论，大家看到了，耗散算子的构造很关键。协调耗散算子所谓协调耗散算子，是这样一种耗散算子，它既起着耗散作用以保障格式的平方守恒性，同时又能减少因时间离散而造成的截断误差以保持格式的高精度，比人工耗散算子协调、自然。协调耗散算子协调耗散算子的数学定义：如果耗散算子B满足：或则称B为m阶协调（时间）耗散算子。耗散性条件对于半离散方程： ，若L反对称，则当 时， 恒成立；而当 时， 成立的条件是 。其中B的定义如下：耗散性证明注：注：利用L的反对称性以及下列关系式：可以验证：耗散性证明忽略二阶及以上小量，即：则利用上式便可直证明B的耗散性（负定性）。协调耗散算子的迭加原理设 均为m阶协调耗散算子，则它们的算术平均：亦为m阶协调耗散算子，其中 为整数。协调耗散算子的最低阶原理设 分别为m1,m2,?,mp阶协调耗散算子，则它们的算术平均：为q阶协调耗散算子，其中 。协调耗散算子的加速原理设 均为m阶协调耗散算子，且具有不同的余项，则它们的某种组合：可成为m+p-1阶协调耗散算子。协调耗散算子的加速原理m阶协调耗散算子 具有不同的余项的含意为：且协调耗散算子的加速原理 的组合为m+p-1阶协调耗散算子的充要条件为：一阶协调耗散算子的构造利用二阶Runge-Kutta方法构造一阶协调耗散算子：其中一阶协调耗散算子的证明先证明其耗散性算子。利用算子L的反对称性可得：一阶协调耗散算子的证明这样再由 得：一阶协调耗散算子的证明由于 ，则 。用数学上的??-??语言描述，即为：???0, ???0, 使得当0≤????时，有 ，即 。取 , 则???00，使得当0≤????0时，有 , 因此 ，即B具有耗散性。一阶协调耗散算子的证明再证明其精度。若LF是关于F的光滑函数，古典Runge-Kutta法已有证明。对于地球流体力学方程，空间微商算子￡ 一般是关于F、 和 的充分光滑函数，即 ,若用中央差对其离散, 则离散算子可表示成一阶协调耗散算子的证明为了证明简练，这里只考虑一维情形，即只考虑x方向： ，对于二维情形的证明，完全类似，其中g(p, q)为充分光滑函数。一阶协调耗散算子的证明利用半离散方程可得：一阶协调耗散算子的证明另一方面，一阶协调耗散算子的证明因此，便有：当 时，得到：一阶协调耗散算子的证明最后有：根据协调耗散算子的定义，上述B算子为一阶协调耗散算子。同样可证后面的二阶和三阶协调耗散算子，不再推导。二阶协调耗散算子的构造利用三阶Runge-Kutta方法构造二阶协调耗散算子：其中三阶协调耗散算子的构造利用四阶Runge-Kutta方法构造三阶协调耗散算子：其中m阶线性协调耗散算子的构造假设L为线性算子，则对任意正整数m，若取： ，其中则由 定义的B算子为m阶协调耗散算子。14 协调耗散算子协调耗散算子关于平方守恒格式的基本定理正压浅水波方程的显式平方守恒格式非负算子定义定义14.1：若算子A满足：则称A为非负算子（当等号成立时，即为广义反对称算子；当大于号成立时，即为正定算子）。非负算子举例定步长显式平方守恒格式的空间离散算子A为一非负算子，其中L为反对称算子，B为耗散算子；变步长显式平方守恒格式的空间离散算子A为一非负算子，其中算子L和B性质同上。稳定性定义定义14.2：对于求解发展方程 的差分格式： ,若当 ??足够小时，上述格式的解满足：则称该格式是计算稳定的，其中 。稳定性定理定理14.1：若A 为非负算子，则当 时差分格式：是无条件计算稳定的。稳定性定理证明事实上，用 对差分格式两端同时进行内积，得到（利用A的非负性）：稳定性定理证明而稳定性定理证明代入不等式稳定性定理证明得到：稳定性定理证明由于 ，故有：即