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

谱方法与高阶时间离散方法及应用 如有限元方法一样采用变分形式，国内研究主要是集中于理论误差分析，很少 用于具体的应用计算．其原因在于谱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)方法，重点介绍了半隐 Deferred 式的Spectral Correction(SISDC)方法，以及采用GMRESKrylov子 Deferred 空间方法加速的Krylov Correction(KDC)方法． 本文主要是应用谱方法及高阶时间离散方法求解实际问题中的偏微分方程 模型．在第四章，我们研究浅水波模型的Korteweg—deVries方程，反应一扩散系 摘 要 程，以及数学物理的非线性Schr6dinger方程．这里我们重点研究描述结晶固体 反相边界运动的Allen—Cahn方程[5】 魏“：^，[△u一去，(札)】 与描述固体复杂相分离与晶体粗化现象的Cahn—Hilliard方程『201 1 魏¨=吖△f一△1f+去，(“)]． ‘ 。‘ S‘ 这两个方程在材料科学里作为相跃迁，模式形成模型，有广泛的应用与数值研 究。『30，29，3l，113]．当￡很小时，构造高阶，有效，无条件能量稳定的数值格 Deferred 们构造稳定的半隐式Spectral Correction(SSISDC)方法[65]．SSISDC 方法是无条件稳定，在每一次迭代过程中只要解简单的线性方程组，在每一个 时问步仅要几次迭代就能达到高阶时间精度．大量的数值结果说明SSISDC方 f102，106，