谱方法和高阶时间离散方法及其应用.pdfVIP

摘 要 在偏微分方程的数学理论里，我们很少能够给出偏微分方程的解析解；在应 用科学里，我们需要数值求解偏微分方程．许多包含低阶的非线性项与高阶的 线性项，与时间相关的半线性抛物型方程(抛物型是指线性部分光滑，半线性意 程，Hodgkin—Huxley方程，Korteweg—de 精确的数值解，我们需要在空问，时间方向用高阶的数值逼近方法．然而许多数 值计算在空间上采用有限差分，有限元等低阶方法离散，由于非线性与刚性带来 的困难，在时间上只有二阶精度． 本文采用谱方法与高阶时间离散方法，研究实际问题的偏微分方程模 型，包括Burgers方程，Korteweg—de 型偏微分方程，如Helmholtz方程，双调和方程等，当方程在几何规则区域上有 光滑解时，相比有限差分方法，有限元方法的低阶精度，采用谱方法离散，仅用 作基函数展开，得到带状代数方程组，求解后谱方法能达到指数阶收敛的谱精度 『83，84，661．对于无界区域上的问题，有限差分，有限元常截断为有限区域，采 用人工边界条件，然而谱方法不用将无界区间截断为有限区问，可以直接采用 解无穷区间上的问题『94]．对于周期问题j采用F0urier谱方法，无论一维，高维， 都容易求解．谱方法的不足在于对方程解的光滑性要求高，大多数应用于理想 的数值例子中，在具体实际问题里，不能普遍应用． 方法并没有像有限差分普及，国内没得到广泛的关注．经典谱方法分成三类： 谱Galerkin方法，谱配点法和谱tau方法．因为编程简单，容易实现，能解变系 数，非线性的偏微分方程，谱配点法在国内得到最多的应用，而谱配点法所用 的谱微分矩阵是满的，不能利用代数线性方程组的快速解法．谱Galerkin方法 谱方法与高阶时间离散方法及应用 如有限元方法一样采用变分形式，国内研究主要是集中于理论误差分析，很少 用于具体的应用计算．其原因在于谱Galerkin方法要将物理空问的值变换到谱 空间，在谱空间计算后，再变换回到物理空间，在编程实现上不如谱配点法容 法，受到普遍的欢迎．本文在第二章采用Trefethen著作的格式，对周期问题的 程的Legendre—Galerkin谱方法编程实现．我们采用理论一算法一代码一算例的形 式，编程模块化，看过后就能理解算法，读懂代码．特别是Legendre—Galerkin谱 方法的代码．可以作为学习谱Galerkin方法编程的模扳．更多的基本MATLAB 代码，可以参考沈捷等人的专著f911． ofLines 对于半线性抛物型偏微分方程的数值求解，我们常采用Method (MoL)方法『821，首先离散空问变量，得到在时间上耦合的常微分方程组 (ODEs)，然后采用合适的时间离散格式求解．常用的ODEs高阶时间离散方法 大致分成两类：第一类由本性高阶离散格式组成，如Runge—Kutta(RK)方法，线 性多步法等；第二类是通过Richardson外插或者延迟校正(DeferredCorrection) 来加速低阶格式的收敛．对于非刚性问题．存在极有效的方法．对于刚性问题， 隐式Runge—Kutta方法有良好的稳定性，但高阶精度需要很大的代价；隐式多步 法有高阶的收敛，但稳定性相对弱．因此，有人采用Runge—Kutta方法达到五阶， 如要求更高阶精度，再采用外插法加速．然而，外插需要计算在更细网格上的一 系列数值解，尽管有效但运算量很大．目前，对于含有高阶空间变量导数的偏微 分方程，特别是这些高阶空间变量导数前面有不小的系数时，经过谱方法离散 空间变量后，生成在时间上耦合的常微分方程组(0DEs)，采用显式格式，局部 时间离散，为保持数值格式稳定，对时间步长有严格的约束．因此，有必要采用 Deferred 隐式稳定的数值格式．当前研究的高阶时间格式有半隐式的Spectral Correction(SDC)方法，AdditiveRunge—Kutta(ARK)方法，以及指数时间差分 Exponential 间离散格式，如广泛应用于Hamiltonian动力学模型的分数步方法，来自于常微 分方程理论的积分因子(IF)方法，指数时间差分(ETD)方法，重点介绍了半隐