[論文レビュー] High Order Implicit-Explicit General Linear Methods with Optimized Stability Regions

Hong Zhang, Adrian Sandu+1 more|arXiv (Cornell University)|Jul 9, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 20被引用数 3

ひとこと要約

本稿では、剛性部分微分方程式を高い精度と効率で解くために、高次の陰・陽的一般線形法（IMEX-GLM）を提案する。具体的には、第四および第五次のDIMSIMスキームを用いる。安定性領域の最適化と高い段階順序を活用することで、従来のIMEXルンゲ＝クッタ法で一般的に生じる次数低下を回避し、より大きな時間刻みを許容し、移流拡散や圧縮性流れなどのマルチスケール問題において優れた性能を発揮する。

ABSTRACT

In the numerical solution of partial differential equations using a method-of-lines approach, the availability of high order spatial discretization schemes motivates the development of sophisticated high order time integration methods. For multiphysics problems with both stiff and non-stiff terms implicit-explicit (IMEX) time stepping methods attempt to combine the lower cost advantage of explicit schemes with the favorable stability properties of implicit schemes. Existing high order IMEX Runge Kutta or linear multistep methods, however, suffer from accuracy or stability reduction. This work shows that IMEX general linear methods (GLMs) are competitive alternatives to classic IMEX schemes for large problems arising in practice. High order IMEX-GLMs are constructed in the framework developed by the authors [34]. The stability regions of the new schemes are optimized numerically. The resulting IMEX-GLMs have similar stability properties as IMEX Runge-Kutta methods, but they do not suffer from order reduction, and are superior in terms of accuracy and efficiency. Numerical experiments with two and three dimensional test problems illustrate the potential of the new schemes to speed up complex applications.

研究の動機と目的

• 現代のPDEソルバーにおける高次空間離散化と整合する高次時間積分法の必要性に対応する。
• 既存のIMEXルンゲ＝クッタ法および線形多段法には、剛性問題や複雑な境界条件において次数低下を生じるという限界があるため、それを克服する。
• 完全な次数精度を維持し、安定性と効率性を向上させる高次IMEX一般線形法（IMEX-GLM）を開発する。
• 数値的最適化により、新しいスキームの安定性領域を最適化し、剛性系に対して頑健であることを保証する。
• 大規模かつ現実世界の応用において、標準的なIMEX-RKスキームと比較して優れた性能を示すことを実証する。

提案手法

• ルンゲ＝クッタ法と線形多段法を一般化するDIMSIM（対角的陰的多段統合法）フレームワークに基づき、第四および第五次のIMEX-GLMを構築する。
• L安定性を満たす陰的成分と、制約付き安定性領域を最大化するように最適化された陽的成分を用いる。特に、剛性固有値が集中する虚軸をカバーするよう配慮する。
• 数値的最適化手法を用いて安定性領域を形状化し、特に虚軸付近に固有値を有する剛性問題に対して頑健性を向上させる。
• 次数低下を防ぐために高い段階順序を確保する。これは、IMEX-RK法が剛性または複雑な境界値問題に適用された際に顕著に現れる主要な問題である。
• 有限差分法および非連続ガラーキン法から得られる空間離散化された常微分方程式系に、時間積分スキームを適用するための方法論的ラインズフレームワークに統合する。
• 2次元および3次元問題における数値実験を通じて、アレン＝チャン方程式、バーガーズ方程式、圧縮性Euler方程式を対象に、精度、安定性、計算コストを比較して検証する。

実験結果

リサーチクエスチョン

• RQ1 従来のIMEX-RK法が失敗するような剛性系においても、次数低下を回避できる高次IMEX-GLMを構築できるか？
• RQ2 剛性成分と複雑な固有値分布を有する問題において、最適化されたIMEX-GLMの安定性特性はIMEX-RK法と比べてどのように異なるか？
• RQ3 高次空間離散化を伴う大規模シミュレーションにおいて、新しいIMEX-GLMはどれほど精度と効率を向上させるか？
• RQ4 IMEX-GLMの安定性領域を、虚軸など複素平面の重要な領域を効果的にカバーするように最適化できるか？特に、固有値が虚軸に集中する問題に対して性能向上が見られるか？
• RQ5 3次元圧縮性流れや2次元移流拡散系といった現実世界の問題において、新しいスキームは既存のIMEX-RK法と比較して実際の性能をどのように発揮するか？

主な発見

• 提案された第四および第五次のIMEX-DIMSIMスキームは、従来のIMEX-RK法とは異なり、剛性問題や複雑な境界条件下でも次数低下を示さない。
• IMEX-DIMSIMはIMEX-RK法と同等の最大安定時間刻みを達成しているが、広い時間刻み範囲で完全な次数精度を維持している。
• 2次元および3次元のテスト問題（上昇気泡、アレン＝チャン方程式など）において、ワーク・プレシジョン図においてIMEX-RK法を上回る効率性を示した。
• 3次元圧縮性Euler問題では、IMEX-DIMSIM4は約10⁻⁷の誤差でワーク・プレシジョンのプラトーに達するが、IMEX-RK4は次数低下を示し、大きな時間刻みで誤差が高くなる。
• 3次元上昇気泡テストにおいて、IMEX-DIMSIM4とIMEX-RK4の両者で最大安定時間刻みは約1.0秒であり、安定性特性が同等であることが確認された。
• IMEX-DIMSIMの安定性領域は、固有値が虚軸に集中する問題に対して特に重要であるため、虚軸の大部分をカバーするように最適化されており、2次元上昇気泡ケースの固有値プロットによって裏付けられた。