时空切比雪夫伪谱方法求解Burgers方程.docVIP

时空切比雪夫伪谱方法求解Burgers方程 刘家春 哈尔滨华德学院，通识教育学院，哈尔滨，150025 关键词: 切比雪夫伪谱方法; Burgers方程; 高精度. 摘要：本文提出一个新的方法求解一维Burgers方程，使用时空切比雪夫伪谱方法求解Burgers方程. 首先使用切比雪夫伪谱方法进行一维空间离散，得到一个常微分方程组，然后再使用切比雪夫伪谱方法求解常微分方程组，最后比较数值解和精确解，数值结果证实该方法有很高的精度和稳定性. 0 引言 本文考虑下面Burgers方程[1]的初边值问题 (1) 其中为常数，,和为已知函数. 本文首先使用切比雪夫伪谱方法[2,3,4,6]离散空间导数，然后再使用切比雪夫伪谱方法求解离散后的常微分方程组. 1 切比雪夫伪谱方法 由区间的Chebyshev-Gauss-Lobatto节点, 做下面映射 则为区间的Chebyshev-Gauss-Lobatto节点. 由Lagrange插值多项式，有 , 其中为二阶微分矩阵， 为Lagrange插值基函数. 对于一阶非线性常方程组，其一般形式为： 初始条件为. 令 , , 有 . 取的区间的Chebyshev-Gauss-Lobatto节点 of , 则化为矩阵形式 . 因此 , 其中 ,, ，, ，. 转化为向量的等价形式为 , (2) 其中为 Kronecker积, 且 , , 利用Newton-Raphson法[5] 求解(2), 便得原非线性常微分方程组解在离散点的值. 取Chebyshev-Gauss-Lobatto节点 of . 在,上有 , 其中 . 则有 (3) 初值 . (4) 我们使用切比雪夫伪谱方法求解一阶非线性常微分方程组(3)-(4)，进而得出原问题的近似解. 2 数值结果 实例 [1] 考虑如下初边值问题 精确解为 . 如下数值结果给出了本文方法的优越性. Fig. 1. 时的网格图 Fig. 2. 时的数值解 Fig. 3. 时的绝对误差曲面图 Fig. 4. 时的误差性能分析图 3 结论 本文使用时空切比雪夫伪谱方法，通过数值实例，和性能分析图，可以说明该方法具有高阶精度和稳定性且在时间和空间上都具有指数阶收敛性（谱精度）. 本文的方法也可以求解其他时空发展方程和刚性常微分方程. 参考文献 Peng Y J, Yan T, Zhang S M, Wang B E. MOL Numerical Solution Method for Burgers Equation[J]. Journal of Xi’an University of Technology, 2004,20(3), 276-279. Shen J, Tang T. Spectral and high-order methods with applications[M]. Beijing: Science Press, 2006. Guo B Y. Spectral methods and their applications[M]. Hong Kong: World Scientific, 1998. Canuto C, Hussaini M Y, Quarteroni A, et al. Spectral methods in fluid dynamics[M]. Berlin: Springer-Verlag, 1987. Deu?hard P. Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms[M], Springer-Verlag, NewYork, 2011. Wang Z Q, Guo B Y. Legendre-Gauss-Radau Collocation Method for Solving Initial Value Problems of first Order Ordinary D