辛算法在双曲型偏微分方程中的应用.pdfVIP

摘要 摘 要 研究不可积系统主要依赖于数值计算，传统数值算法不可避免地带有人为耗 散性的缺陷，故在稳定性和保结构方面均具有良好优势的辛算法越来越被人们关 注。 二阶辛算法虽然形式简单，但却非常重要，常被用来构造高阶辛算法。 Yoshida 提出通过辛映射将低阶辛算法通过特定的系数组合成高阶辛算法。Li 和 Wu 通过将二阶辛算法的三阶误差项消去构造了适用势函数为二次型函数的新型 四阶力梯度辛算法。这种辛算法形式简单且在数值精度上也显示良好优越性。 本文主要是以波动方程和非线性Sine-Gordon 方程为例研究辛算法在双曲型 偏微分方程的应用。一方面，通过对波动方程的空间方向进行数值离散得到有限 维的 Hamilton 系统。然后利用辛算法对该有限维系统进行数值求解，从而得到 对原波动方程而言的半离散化辛算法，并对由不同空间差分格式和不同辛算法构 造得到的几类不同精度半离散化辛算法的稳定性加以讨论。此外，注意到波动方 程数值离散后得到的 Hamilton 系统的势函数是二次型的，本文将新型的三步四 阶辛算法推广运用，构造了适用于波动双曲型方程的高阶半离散数值辛算法。另 一方面，文章探讨了非线性Sine-Gordon 方程的复合构造多辛算法，并与同阶的 半离散化辛算法进行了数值比较，得到半离散辛算法和多辛算法在数值精度及计 算效率上各自的优越性。 关键词：Hamilton 系统；辛几何算法；多辛算法 I 万方数据 ABSTRACT ABSTRACT Numerical calculations have been mainly used on the research of non-integrable systems ．However ，the traditional numerical algorithms appear the deviation of their trajectry which are caused by their artificial dissipation ．Therefore ，symplectic algorithms take more and more attraction because of the advantage on the stability and symplectic structure ． The second-order symplectic method ，as the simplest of the explicit symplectic integrators，often be used to construct high order symplectic integrators ．Yoshida propose a method to construct the higher order integrators through by the combination of lower order integrators ．Li and Wu construct two new three stage fourth-order symplectic integrators which have the similar form with the second-order integrators and