[论文解读] Spectral algorithms for reaction-diffusion equations

Richard V. Craster, Roberto Sassi|arXiv (Cornell University)|Oct 17, 2018

Advanced Mathematical Modeling in Engineering参考文献 24被引用 13

一句话总结

本论文提出了一套用于在MATLAB和Fortran 77中求解一维和二维反应-扩散方程的谱算法，采用谱方法结合指数时间差分和傅里叶谱空间离散化。该方法有效消除了扩散项带来的刚性，使得能够采用较大的显式时间步长并实现高阶精度，显著优于标准方法在速度和效率方面的表现，同时在基准问题上保持了高精度。

ABSTRACT

A collection of codes (in MATLAB & Fortran 77), and examples, for solving reaction-diffusion equations in one and two space dimensions is presented. In areas of the mathematical community spectral methods are used to remove the stiffness associated with the diffusive terms in a reaction-diffusion model allowing explicit high order timestepping to be used. This is particularly valuable for two (and higher) space dimension problems. Our aim here is to provide codes, together with examples, to allow practioners to easily utilize, understand and implement these ideas; we incorporate recent theoretical advances such as exponential time differencing methods and provide timings and error comparisons with other more standard approaches. The examples are chosen from the literature to illustrate points and queries that naturally arise.

研究动机与目标

• 为求解一维和二维反应-扩散方程，提供实用且可重用的MATLAB和Fortran 77谱代码。
• 展示通过积分因子和指数时间差分法的谱方法在克服扩散项带来的刚性方面的有效性。
• 在标准方法计算成本过高的二维问题中，实现高阶显式时间推进。
• 从精度、速度和实现复杂度方面，将谱方法与传统格式进行比较。
• 通过数学生物学、化学和物理学中的实例，展示该方法的通用性。

提出的方法

• 通过快速傅里叶变换（FFTs）实现傅里叶谱配点法进行空间离散化，以达到谱精度。
• 应用积分因子以消除扩散项带来的刚性，将系统转化为非刚性常微分方程系统。
• 采用指数时间差分（ETD）结合显式龙格-库塔时间积分方法，实现高阶时间精度。
• 结合近期理论进展，如ETD格式和傅里叶插值函数，以提高分辨率。
• 在FFTs中引入零填充，以对解进行重采样，提升图形可视化效果，同时不损失谱精度。
• 与交替方向隐式（ADI）方法进行比较，以评估精度与计算成本之间的权衡。

实验结果

研究问题

• RQ1 谱方法结合积分因子和指数时间差分法是否能够实现二维反应-扩散方程的稳定、高阶、显式时间推进？
• RQ2 对于二维反应-扩散问题，谱算法在精度和速度上与标准有限差分或有限元方法相比如何？
• RQ3 与传统方法相比，谱方法在二维模拟中能在多大程度上降低计算成本？
• RQ4 该方法在Fisher方程和Gray-Scott模型等基准问题上的表现如何？
• RQ5 通过FFT零填充实现的解重采样对视觉质量与精度有何影响？

主要发现

• 结合积分因子和指数时间差分法的谱方法可实现稳定、高阶的显式时间推进，支持大时间步长，显著快于标准方法。
• 在50×50区域的256×256网格上，当Δt = 10^-2时，谱方法的相对误差达到10^-6量级，其速度和精度均优于ADI方法。
• 该方法将原本需在超级计算机上耗时数小时的二维问题，缩短至在普通PC上仅需数分钟，展现出极高的计算效率。
• 通过零填充实现的傅里叶插值函数提高了视觉分辨率，同时不损失谱精度，解决了粗网格绘图中的问题。
• 该谱方法在一维问题中与自适应有限差分格式相比具有竞争力，且通常更快、更易于实现。
• ADI方法虽无条件稳定，但需采用更小的时间步长，且因具有O(Δt²)的精度和矩阵求逆开销，对反应-扩散问题效率较低。