[論文レビュー] High-order low-storage explicit Runge-Kutta schemes for equations with quadratic nonlinearities

Pablo D. Mininni, Mininni, Pablo D.+2 more|ArXiv.org|Aug 13, 2008

Computational Fluid Dynamics and Aerodynamics参考文献 21被引用数 3

ひとこと要約

本稿では、補正項を用いてJST法を拡張することで、二次非線形性を有する方程式に対して高次・低記憶型の陽的ルンゲ・クッタスキームを提案する。これにより、第三・第四次精度を維持しつつ、最小限の記憶容量を確保できる。本手法はバーガース方程式および非圧縮性オイラー方程式において検証され、計算コストのわずかな増加で時間精度と安定性が向上することが示された。

ABSTRACT

We show in this paper that third- and fourth-order low storage Runge-Kutta algorithms can be built specifically for quadratic nonlinear operators, at the expense of roughly doubling the time needed for evaluating the temporal derivatives. The resulting algorithms are especially well suited for computational fluid dynamics. Examples are given for the Hénon-Heiles Hamiltonian system and, in one and two space dimensions, for the Burgers equation using both a pseudo-spectral code and a spectral element code, respectively. The scheme is also shown to be practical in three space solving the incompressible Euler equation using a fully parallelized pseudo-spectral code.

研究の動機と目的

• 線形問題に対して高次であるが、二次非線形系に対しては第二次の精度しか達成しないJST法の限界を克服すること。
• 右辺に二次非線形性を含む場合に、JST法に第三・第四次精度を回復させるための補正項を構築すること。
• 元のJSTスキームが持つ低記憶特性を維持し、記憶容量が限られた大規模なシミュレーションへの適用可能性を確保すること。
• 代表的な問題としてバーガース方程式および多次元非圧縮性オイラー方程式に対して、新規スキームの妥当性を検証すること。
• 高解像度シミュレーションにおいて時間誤差が空間誤差を支配する場合、高次時間積分が不可欠であることを示すこと。

提案手法

• 本手法は、低記憶型ルンゲ・クッタスキームであるJST法を出発点とし、正確解のテイラー展開から導出された補正項を追加することで変更する。
• 補正項は、F(u) = L(u) + N(u,u) の形をとる系、すなわちLが線形で、Nが引数に関して対称な二次非線形性である場合に特化して導出される。
• 補正項は、元のJST法が有する二段階記憶構造を維持する形で適用され、低記憶量を確保する。
• 本スキームは、既存の時間積分フレームワークを最小限に変更して、擬スペクトル法およびスペクトル有限要素法のコードに実装された。
• ヘノン＝ハイレス系、1次元および2次元のバーガース方程式、および完全並列化された擬スペクトルコードを用いた3次元非圧縮性オイラー方程式に本手法を適用した。
• 3+スキーム（第三次補正項を含む）は、4回のJST反復で実装可能であり、標準の第三次スキームと比較して著しく低いグローバル誤差を示したため、最も費用効果に優れた実装とされた。

実験結果

リサーチクエスチョン

• RQ1 JST法に補正項を追加することで、記憶要件を増加させることなく、二次非線形系に対して第三・第四次精度を達成できるか？
• RQ2 非線形二次作用素に対してJST法を第二次以上に拡張するために必要な具体的な補正項は何か？
• RQ3 高解像度シミュレーションにおいて、精度と安定性の観点から、補正済みスキームは標準的な低記憶型ルンゲ・クッタ法と比較してどのように性能を発揮するか？
• RQ4 本スキームの向上した時間精度は、たとえば乱流の直接数値シミュレーションのような大規模シミュレーションにおいて実用的に有益であるか？
• RQ5 既存の高次スペクトル法およびスペクトル有限要素コードに、最小限のコード変更で効率的に本補正スキームを統合できるか？

主な発見

• JST法は線形問題に対して正確であるが、二次非線形系に適用すると第三次誤差を生じるため、高次精度を達成するには補正項の導入が不可欠である。
• 導出された補正項は、二次非線形性に対して第三・第四次精度を回復させるとともに、元のJSTスキームの二段階記憶構造をそのままで保持した。
• 3+スキーム（第三次補正項を含む）は、4回のJST反復で実装可能であり、標準の第三次スキームと比較して著しく低いグローバル誤差を示したため、最も費用効果に優れた実装とされた。
• バーガース方程式の数値結果から、スキームの次数が上昇するにつれて切断誤差が減少することが示され、3+スキームは標準の第三次手法を上回ることを確認した。
• 3次元非圧縮性オイラー方程式のシミュレーションでは、新規スキームにより高解像度で長時間にわたる積分が可能となり、低次時間積分で見られる不自然な保存性の誤りを回避できた。
• 本手法は大規模なシミュレーション、特に3次元擬スペクトルコードにおいて実用的かつ効率的であり、ペタスケール規模の乱流直接数値シミュレーションに適している。