紧致差分格式说课.ppt

五点紧致差分格式的原理与应用 汇报学生：吴广智 指导老师：孙建安教授 一，基本原理介绍 二，实际应用：数值求解RLW方程 1.方程介绍 2.数值格式的建立 3.算例与数值结果的对比 三，未来方向 1.紧致差分格式的隐式应用1 2.紧致差分格式的隐式应用2 四，致谢 主要内容 一，基本原理介绍 　紧致有限差分方法是使用函数值的某种线性组合来表示该函数导数值的线性组合的一类差分方法, 该方法能有效增加差分格式的精度与稳定性。其基本原理如下： 对于函数 ，将自变量区间n等分, 插入n+1个节点 , 相邻节点间距为h, 则函数一阶导数的对称紧致差分格式为： 内点六阶精度： （1） 近边界点四阶精度： （2） （3） （4） 同样的，函数二阶导数的紧致差分格式可表示为： 内点六阶精度： （5） （6） 虽然一阶导数和二阶导数的差分格式精度的理论值一样，但是实际一阶导数的数值解更接近精确值。 其次，由于紧致差分格式是使用函数在离散点上的函数值计算该函数在单个离散点处的一阶或高阶导数值的线性组合，所以导数值的求出需要求解线性方程组，所以该方法很难用于求解时间空间混合导数项，或用于建立求解某一方程的隐式格式。 （7） （8） 二，实际应用：数值求解RLW方程 1，方程介绍 RLW方程是由Peregrine提出的一类非线性演化方程, 是描述许多物理现象(如浅水波、等离子体声波等)的一种非常好的模型, 尤其是在研究非线性色散波方面起了非常重要的作用, 因此对其数值解法的研究十分重要. 众多数值方法都曾经用于求解RLW方程, 例如五次和二次B样条Petrov-Galerkin有限元法、伽辽金线性有限元法、二次B样条集中伽辽金有限元法等。 其具体形式如下： （9） 2.数值格式的建立 对时间差分采用四阶龙格库塔方法，对空间差分采用 五点紧致差分格式，由于混合导数项的存在龙格库塔方法 与紧致差分方法都无法直接作用与方程，采用如下方法处 理。 引入变量 将方程（9）改写为 （10） （11） 对方程（11）使用四阶龙格库塔方法: （12） 其中 表示 在第n时间层的值， 由于 已知, 由式（1）和式（5）可解得 进而可得到 与 的值。下面以 为例说明 的求解方法。 记 由（10）式可得 利用（5） 式可得： （13） 由于 已知求解（13）式确定的线性方程组即可求得 ，利用（1）式可得到 ，继而得到： 同样的可以求得 ，代入（12）式可求出 再将（13）式中的 替换为 重新求解（13）式 即可得到 ，到此完成求解。 3.算例数值结果的对比 考虑如下初始条件的RLW方程 其对应的方程的精确解为 其中 ，实际计算时取 由于算例为孤波解，所以边界处近似满足一阶导数值为0.因此, 为了简化边界点与近边界点的处理方式, 实际计算时在求解区间左右端点的外侧分别外插了四个节点，且该四个节点处函数值为0，这样只使用公式（1）（5）即可。 为了方便比较引入误差范数与守恒量，定义为： 表1给出了本文算法与几种其他算法求解RLW方程孤波解所得到的数值结果在时间t=20时的误差范数与守恒量对比. 表2给出本文算法所得的数值结果与伽辽金有限元方法与余弦微分求积所得的数值结果在不同时刻的误差范数与守恒量对比. 表3给出在振幅不同的情况下本文算法所得的数值结果与余弦微分求积法得到的的数值结果在不同时刻的误差范数与守恒量对比. 表1 误差范数与守恒量比较(RLW方程孤波解) 数值方法 L2 L∞ I 紧致差分 0.599×10-5 0.214×10-5 3.97995 有限差分 0.796×10-3 0.281×10-3 3.97996 CDQM 0.209×10-3 0.078×10-3 3.98002 线性伽辽金 0.511×10-3 0.198×10-3 3.98206 伽辽金有限元 0.219×10-3 0.086×10-3 3.97988 分裂法 196.1×10-3 67.35×10-3 4.41219 二次B样条有限元法 0.227×10-3 0.081×10-3 3.97986 B样条配置法 0.53