[論文レビュー] High-order generalized-$\alpha$ methods

Quanling Deng, Pouria Behnoudfar+1 more|arXiv (Cornell University)|Feb 14, 2019

Numerical methods for differential equations参考文献 8被引用数 4

ひとこと要約

本稿では、第2次精度の一般化-$\alpha$法を高次離散化展開を用いて拡張することで、無条件安定性を維持したままユーザーが制御可能な高周波数減衰を有する第3次精度の一般化-$\alpha$時間積分法を提案する。この手法は、元のスキームの単一ステップ構造を保ちつつ第3次精度を達成する。数値的減衰はパラメータ $\alpha_m$ と $\alpha_f$ を用いて制御される。

ABSTRACT

The generalized-$\\alpha$ method encompasses a wide range of time integrators. The method possesses high-frequency dissipation while minimizing unwanted low-frequency dissipation and the numerical dissipation can be controlled by the user. The method is unconditionally stable and is of second-order accuracy in time. We extend the second-order generalized-$\\alpha$ method to third-order in time while the numerical dissipation can be controlled in a similar fashion. We establish that the third-order method is unconditionally stable. We discuss a possible path to the generalization to higher order schemes. All these high-order schemes can be easily implemented into programs that already contain the second-order generalized-$\\alpha$ method.

研究の動機と目的

• 第2次精度の一般化-$\alpha$法を時間的に第3次精度に拡張するが、その主要特徴を維持する。
• 新たな高次スクリプトにおいても無条件安定性を維持する。
• 一般化-$\alpha$法の特徴である高周波数数値的減衰のユーザー制御を保持する。
• 既存の第2次一般化-$\alpha$コードベースへのスムーズな統合を可能にする。
• より高次スクリプトへの一般化のためのフレームワークを提供する。

提案手法

• 一般化-$\alpha$法を再定式化し、$U_{n}^{\alpha_f}$ と $V_{n}^{\alpha_m}$ を線形サブステップではなく、高次テイラー展開として扱う。
• $U_{n+\alpha_f}$ および $V_{n+\alpha_m}$ のテイラー展開における高次項を表すために、追加のパラメータ $\alpha_f$ と $\alpha_m$ を導入し、第3次精度を達成する。
• テイラー級数解析を用いて、最適な精度を得るための自由パラメータ $\alpha_m$ と $\alpha_f$ を導出する。
• 増幅行列を構築し、その固有値を分析することで、無条件安定性領域を確立する。
• 高周波数固有値 $\rho_\infty$ を $[0,1]$ のユーザー定義値に設定することで、数値的減衰を制御する。
• 記号的および数値的解析を用いて、無条件安定性を検証し、安定性および所望の減衰を実現するパラメータ領域を特定する。

実験結果

リサーチクエスチョン

• RQ1 一般化-$\alpha$法を、無条件安定性を損なわず、時間的に第3次精度に拡張することは可能か？
• RQ2 高次時間積分法において、ユーザー制御可能な高周波数減衰をどのように維持できるか？
• RQ3 どのパラメータ選択 $\alpha_m$ と $\alpha_f$ が第3次精度と無条件安定性を達成するか？
• RQ4 得られたスキームは、既存の第2次一般化-$\alpha$コードに対して最小限の変更で実装可能か？
• RQ5 第3次一般化-$\alpha$法の安定性領域は何か？第2次バージョンと比較するとどうなるか？

主な発見

• 提案された第3次一般化-$\alpha$法は、解および速度項の高次テイラー展開を用いることで、時間的に形式的に第3次精度を達成する。
• すべての $\alpha_m \geq \frac{7}{12}$ および $\frac{1}{2} \leq \alpha_f \leq \alpha_m - \frac{1}{12}$ に対して、この手法は無条件安定であることが示され、記号的および数値的に確認された。
• 高周波数における数値的減衰はパラメータ $\rho_\infty$ によって制御され、$0 \leq \rho_\infty \leq 1$ に対して $\alpha_m = \frac{13 + 20\rho_\infty - 5\rho_\infty^2}{12(\rho_\infty + 1)^2}$ および $\alpha_f = \frac{1 + 3\rho_\infty}{2(\rho_\infty + 1)^2}$ が成り立つ。
• 増幅行列の固有値は、無条件安定性領域全域で絶対値が1以下に制限されており、安定性が確認された。
• 数値的プロットにより、第3次スキームは第2次バージョンよりも大きな無条件安定性領域を有することが示された。
• この手法は単一ステップ・陰的スキームのままであり、既存の第2次一般化-$\alpha$ソルバーへの統合が、最小限のコード変更で容易に可能である。