紧致差分格式报告.ppt

五点紧致差分格式的原理与应用,汇报学生：吴广智指导老师：孙建安教授,一，基本原理介绍二，实际应用：数值求解RLW方程1.方程介绍2.数值格式的建立3.算例与数值结果的对比三，未来方向1.紧致差分格式的隐式应用12.紧致差分格式的隐式应用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方程孤波解),注:,守恒量精确值为3.97995,表2误差范数比较(RLW方程孤波解),守恒量精确值为3.97995,注:,表3误差范数比较(RLW方程孤波解),注：,守恒量精确值为2.10941,三，未来方向,1，紧致差分格式的隐式应用1,以Kdv方程为例简要说明紧致差分格式的隐式应用，方程形式如下：,（14）,方程（14）的隐式格式可写为：,上式可改写为：,（15）,对（15）式要做类似（13）式的处理，但是对于n+1时间层的一阶导数项与三阶导数项应配备不同的系数，因此引入新变量（15）式可改写为：,（16）,使用（1,5）式对（16）式做类似（13）式的处理：,（17）,由于非线性项的限制（17）式无法直接求解，对（1）式做如下变换。,（18）,（17）式与（18）式联合可以确定一组线性方程组，求解可得,2，紧致差分格式的隐式应用2紧致差分格式的公式（1,2,3,4,5,6,7,8）可以用矩阵乘法来表示：,所以可做如下变换：,这样也可实现对隐式格式的求解。根据现阶段计算实验显示，紧致差分格式