[論文レビュー] A fully adaptive explicit stabilized integrator for advection-diffusion-reaction problems

Ibrahim Almuslimani|arXiv (Cornell University)|Jan 25, 2022

Numerical methods for differential equations参考文献 26被引用数 5

ひとこと要約

本稿では、対流拡散反応方程式に適した完全に適応可能な2次精度の明示的安定化ルンゲ・クッタ・チェビシェフ法、ARKCを提案する。この手法は2種類のチェビシェフ多項式を用いて、時間刻み幅、段数、減衰パラメータを動的に調整する安定性領域を構築し、特に高Peclet数において、高価な対流項の評価回数を最小限に抑えつつ、高い精度と安定性を維持する。既存手法に比べ顕著に優れた性能を示す。

ABSTRACT

A novel second order family of explicit stabilized Runge-Kutta-Chebyshev methods for advection-diffusion-reaction equations is introduced. The new methods outperform existing schemes for relatively high Peclet number due to their favorable stability properties and explicitly available coefficients. The construction of the new schemes is based on stabilization using second kind Chebyshev polynomials first used in the construction of the stochastic integrator SK-ROCK. An adaptive algorithm to implement the new scheme is proposed. This algorithm is able to automatically select the suitable step size, number of stages, and damping parameter at each integration step. Numerical experiments that illustrate the efficiency of the new algorithm are presented.

研究の動機と目的

• 本稿の目的は、標準的な明示的スキームが不安定であり、陰的スキームが計算コストが高いため、高Peclet数を有する対流拡散反応方程式を効率的に統合することの課題に対処することである。
• PIROCK や IMPRKC といった既存の安定化手法の限界、すなわち減衰パラメータが固定されていること、または対流／反応項の評価回数が多すぎるという点を克服することを目的としている。
• 局所誤差とPeclet数に基づき、各ステップで時間刻み幅、段数、減衰パラメータを自動的に選択できる完全適応的積分法の開発を目的としている。
• 関数評価回数、特に高価な非剛性反応項の評価を最小限に抑えながら、高い精度と安定性を維持することを目的としている。
• 科学計算およびPDEシミュレーションにおける剛性のある対流拡散反応問題のための、強固で効率的かつ柔軟な時間積分スキームの提供を目的としている。

提案手法

• 本手法は、虚軸に沿って広い安定性領域を有する2種類のチェビシェフ多項式を用いて導出された2次精度の明示的安定化ルンゲ・クッタ・チェビシェフスキームに基づいている。
• SDEのSK-ROCKにインspiredされた、1種類目と2種類目のチェビシェフ多項式を組み合わせた新規な安定化技術を採用し、負の実軸に沿った長さを損なわずに安定性領域の幅を拡大している。
• 局所誤差推定器に基づき、各積分ステップで段数、時間刻み幅、減衰パラメータを動的に選択する適応アルゴリズム（アルゴリズム4.4）を採用している。
• 安定性領域は、ODE系の固有値分布に合わせて動的に形状を変化させ、Peclet数に依存するため、最適な時間刻み幅選択が可能である。
• 任意の段数に対して係数を明示的に計算可能であり、効率的な実装と適応的制御を可能としている。
• 半離散化された対流拡散反応PDEに適用され、拡散項（ROCK2に類似した安定化を用いて明示的に処理）と対流／反応項（各段で評価）が分離されている。

実験結果

リサーチクエスチョン

• RQ1 大きなPeclet数を有する対流拡散反応問題に対して、高い安定性と精度を維持しつつ、2次精度の明示的安定化ルンゲ・クッタ・チェビシェフ法をどのように構築できるか？
• RQ2 安定性と精度を維持しつつ、効率性を最大化するための、減衰パラメータ、段数、時間刻み幅の最適な適応的制御法は何か？
• RQ3 半離散化から生じるODE系の固有値分布の変化に合わせて、安定性領域を動的に調整できる完全適応的積分法を設計できるか？
• RQ4 PIROCK、IMPRKC、PRKCといった既存の安定化手法と比較して、新しい適応的スキームの関数評価回数と誤差制御における性能はどの程度か？
• RQ5 本手法では係数が明示的に利用可能であるが、閉形式の係数式を持たないPIROCKなどの手法と比較して、その明示的可視性が適応性と計算コストの低減にどの程度寄与するか？

主な発見

• Peclet数10の場合、ARKCはわずか18ステップ、拡散項の評価が3593回で済んだが、PIROCKは28ステップ、4306回の評価を要した。これは、極めて優れた効率性を示している。
• 許容誤差10−5の条件下で、ARKCは最大誤差3.3×10−7を達成し、237段と2132回の対流項評価で、PIROCKおよびIMPRKCを両方の観点で上回った。
• 可変Peclet数を有するバーガース方程式に対して、ARKCはIMPRKCに比べて対流項評価回数を最大70％、PIROCKに比べて50％削減したが、精度は同等であった。
• ARKCの適応的減衰により、変化するPeclet数に対しても広い安定性領域を維持でき、固定パrameterスキームに比べてより大きな時間刻み幅と少ない段数を実現した。
• テストされたすべての状態において、ARKCはIMPRKCおよびPIROCKに比べて、対流項の評価回数が著しく少なかった。これは、時間ステップ数が少ないためであり、拡散項の評価回数はほぼ同等であった。
• 本手法の適応性と係数の明示的計算により、特に高Pe領域において、動的安定性領域制御のおかげでPIROCKをすべてのPeclet数領域で上回った。