[論文レビュー] Spectral algorithms for reaction-diffusion equations

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

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

ひとこと要約

本稿では、1次元および2次元の空間において、指数時間差分法とフーリエスペクトル空間離散化を用いたスペクトル法により反応拡散方程式を解くための、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.

研究の動機と目的

• 1次元および2次元の反応拡散方程式を解くための実用的で再利用可能なスペクトルコード（MATLABおよびFortran 77）を提供すること。
• 積分因子と指数時間差分法を用いたスペクトル法が、拡散項に起因する剛性を克服する有効性を示すこと。
• 標準的手法が計算的に非現実的となる2次元問題において、高次の明示的時間積分を可能にすること。
• 精度、速度、実装の複雑さの観点から、スペクトル法と従来のスキームを比較すること。
• 数学的生物学、化学、物理学の分野における例を通じて、本手法の汎用性を示すこと。

提案手法

• 高速フーリエ変換（FFTs）を用いたフーリエスペクトルコロケーションによる空間離散化により、スペクトル精度を達成する。
• 拡散項に起因する剛性を除去するために積分因子を適用し、系を剛性のない常微分方程式系に変換する。
• 高次の時間精度を達成するために、指数時間差分法（ETD）と明示的ルンゲ＝クッタ時間積分を組み合わせる。
• 最近の理論的進展（ETDスキームやフーリエ補間関数による補間）を組み込み、解の分解能を向上させる。
• スペクトル精度を損なわずにグラフィカルな可視化を向上させるために、FFTにゼロパディングを適用して解をアップサンプリングする。
• 精度と計算コストのトレードオフを評価するために、交替方向陰解法（ADI）手法と比較する。

実験結果

リサーチクエスチョン

• RQ1 積分因子と指数時間差分法を用いたスペクトル法は、2次元の反応拡散方程式に対して安定的で高次の明示的時間積分を可能にするか？
• RQ2 2次元の反応拡散問題において、スペクトルアルゴリズムは、標準的な有限差分法や有限要素法に比べて精度と速度で優れているか？
• RQ3 従来のアプローチに比べて、スペクトル法は2次元シミュレーションにおける計算コストをどの程度削減できるか？
• RQ4 フィッシャー方程式やグレイ＝スコットモデルといったベンチマーク問題において、本手法の性能はいかがであるか？
• RQ5 FFTのゼロパディングによる解のアップサンプリングは、視覚的品質と精度にどのような影響を与えるか？

主な発見

• 積分因子と指数時間差分法を用いたスペクトル法により、安定的で高次の明示的時間積分が可能となり、大規模な時間ステップが可能となり、標準的手法よりも著しく高速なシミュレーションが実現された。
• 50×50領域における256×256グリッドでは、Δt = 10^-2 の条件下で相対誤差が約10^-6のオーダーに抑えられ、ADI法に比して速度と精度の両面で優れた性能を示した。
• 2次元問題において、従来のスーパーコンピュータを用いた数時間のシミュレーションを、標準PCで数分に短縮するという高い計算効率を示した。
• ゼロパディングを用いたフーリエ補間関数による補間により、スペクトル精度を損なわずに視覚的分解能が向上し、粗いグリッドでのプロットに起因する問題が解消された。
• 1次元問題では、適応的有限差分スキームと同等の性能を示し、しばしばより高速かつ実装が容易であった。
• ADI法は無条件安定であるが、時間ステップが小さくならざるを得ず、O(Δt²)の精度と行列逆行列計算のオーバーヘッドにより、反応拡散問題に対しては非効率であった。