[論文レビュー] Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation

Hendrik Ranocha, Jan Glaubitz+2 more|arXiv (Cornell University)|Jun 3, 2016

Medical Imaging Techniques and Applications参考文献 13被引用数 6

ひとこと要約

この論文は、双曲型保存則に対する補正手順による再構成（CPR）スキームの安定性を、和分ごとに部分（SBP）および同時近似項（SAT）フレームワーク内での人工拡산を導入することで向上させる。明示的オイラー時間離散化を用いて完全離散的安定性を証明し、十分に小さな時間ステップに対して安定性を保証するための適応的拡散戦略を提案する。

ABSTRACT

The correction procedure via reconstruction (CPR, also known as flux reconstruction) is a framework of high order semidiscretisations used for the numerical solution of hyperbolic conservation laws. Using a reformulation of these schemes relying on summation-by-parts (SBP) operators and simultaneous approximation terms (SATs), artificial dissipation / spectral viscosity operators are investigated in this first part of a series. Semidiscrete stability results for linear advection and Burgers' equation as model problems are extended to fully discrete stability by an explicit Euler method. As second part of this series, Glaubitz, Ranocha, Öffner, and Sonar (Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering, 2016) investigate connections to modal filters and their application instead of artificial dissipation.

研究の動機と目的

• 双曲型保存則に対する高次CPRスキームにおける非線形安定性の欠如を解消すること。
• 半離散的SBP-SAT安定性結果を明示的時間積分を用いた完全離散スキームに拡張すること。
• 明示的時間積分において、十分に小さな時間ステップに対して安定性を維持することができる適応的人工拡散戦略を開発すること。
• SBP作用素を用いたCPRスキームにおける人工拡散の理論的基盤を確立すること。
• 第2部のシリーズの基盤を築くこと。第2部では、モーダルフィルタリングとの関係を検討する。

提案手法

• 証明可能な安定性を可能にするために、CPRスキームを和分ごとに部分（SBP）作用素および同時近似項（SAT）を用いて再定式化する。
• 連続的設定で人工拡散作用素を導入し、その後SBP作用素を用いて離散化することで、保存則および安定性を維持する。
• 半離散的設定において、安定性を維持し、部分積分を模倣するための拡散作用素の特別な表現を用いる。
• 完全離散スキームを達成するために、明示的オイラー時間離散化を適用する。
• 時間ステップサイズに応じて人工拡散の強度を調整するための適応的アルゴリズムを提案する。
• 離散的固有値を連続的ルジャンドル多項式の対応物と比較することで、安定性特性を検証する。

実験結果

リサーチクエスチョン

• RQ1 SBP-SATフレームワークを用いて、CPRスキームに人工拡散を効果的に組み込むことができるか。これにより完全離散的安定性が保証されるか。
• RQ2 人工拡散の強度の選択が、明示的時間積分を用いた完全離散CPRスキームの安定性にどのように影響するか。
• RQ3 拡散作用素の離散的固有値と連続的ルジャンドル多項式スペクトルとの関係は何か。
• RQ4 明示的CPRスキームにおいて、十分に小さな時間ステップに対して安定性を維持できる適応的拡散戦略を設計できるか。
• RQ5 SBPに基づく人工拡散作用素は、バーガース方程式のような非線形問題において、保存則および安定性をどのように維持するか。

主な発見

• 提案された人工拡散作用素は、CPRスキームにおけるSBP-SATフレームワークで保存則および安定性を維持する。
• 線形移流およびバーガース方程式に対する半離散的安定性が、明示的オイラー法を用いた完全離散スキームに拡張された。
• 拡散作用素の離散的固有値は、連続的ルジャンドル多項式スペクトルとよく一致しており、作用素の設計が妥当であることが検証された。
• 適応的拡散戦略により、時間ステップが十分に小さい場合に完全離散的安定性が保証された。
• 陰的時間積分を必要とせず、効率的な高次シミュレーションが可能になった。
• 理論的枠組みは、第2部のシリーズにおけるモーダルフィルタリングとの関係を検討する基盤を提供した。