[論文レビュー] Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering

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

Medical Imaging Techniques and Applications参考文献 12被引用数 5

ひとこと要約

本稿では、双曲型保存則に対する補正手順による再構成（CPR）スキームの安定性を、和分ごとに部分和（SBP）フレームワーク内でのモーダルフィルタリングを導入することで向上させる。1次精度までにおいてモーダルフィルタリングとスペクトル粘性が等価であることを確立し、適応的フィルタ強度戦略を提案するとともに、明示的エーラー時間積分を用いた完全離散スキームの証明可能な安定性を示し、従来の人工拡散法で見られる時間刻み制限を克服する。

ABSTRACT

A recently introduced framework of semidiscretisations for hyperbolic conservation laws known as correction procedure via reconstruction (CPR, also known as flux reconstruction) is considered in the extended setting of summation-by-parts (SBP) operators using simultaneous approximation terms (SATs). This reformulation can yield stable semidiscretisations for linear advection and Burgers' equation as model problems. In order to enhance these properties, modal filters are introduced to this framework. As a second part of a series, the results of Ranocha, Glaubitz, Öffner, and Sonar ("Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation", 2016) concerning artificial dissipation / spectral viscosity are extended, yielding fully discrete stable schemes. Additionally, a new adaptive strategy to compute the filter strength is introduced and different possible applications of modal filters are compared both theoretically and numerically.

研究の動機と目的

• 双曲型保存則に対する補正手順による再構成（CPR）スキームの安定性を、SBP-SATフレームワーク内でのモーダルフィルタリングを用いて向上させること。
• 第I部で得られた人工拡散の安定性結果をモーダルフィルタリングに拡張し、完全離散スキームの安定性を保証すること。
• ユーザー入力の削減と耐障害性の向上を図るため、フィルタ強度を自動的に選択する新しい適応戦略を開発すること。
• 異なるモーダルフィルタリング適用戦略（解のフィルタリング、時間微分のフィルタリング、および作用素分割法）を比較・分析すること。
• SBP-CPRスキームの文脈において、モーダルフィルタリングとスペクトル粘性との理論的および数値的等価性を確立すること。

提案手法

• 線形移流方程式およびバーガース方程式に対して、同時近似項（SATs）を用いた和分ごとに部分和（SBP）作用素を用いてCPR法を再定式化し、保存的かつ安定な半離散化を保証する。
• 解または時間微分に対して指数関数的フィルタを用いてモーダルフィルタリングを適用し、フィルタ行列を $\underline{\underline{F}}$ で表す。
• 作用素分割法を用いて、モーダルフィルタリングとスペクトル粘性の等価性を分析し、時間微分のフィルタリングが $\underline{\underline{M}}\underline{\underline{F}}^{-1}$ が誘導するノルムの下で安定な半離散化に対応することを示す。
• $(\Delta t)^2$ の誤差項を補償するための新しい適応的フィルタ強度戦略を導入し、明示的エーラー時間積分を用いた証明可能な安定性を可能にする。
• 時間離散化には明示的エーラースキームを用い、フィルタド時間微分を $\partial_t \underline{u} \mapsto \partial_t \underline{\underline{F}} \underline{u}$ とし、$\underline{u}_{+} = \underline{u} + \Delta t \, \partial_t \underline{\underline{F}} \underline{u}$ を得る。
• 数値実験では、解のフィルタリング、時間微分のフィルタリング、および作用素分割法のアプローチを比較し、エネルギーノルムを用いて安定性を評価する。

実験結果

リサーチクエスチョン

• RQ1 モーダルフィルタリングは、SBP-CPRスキームに効果的に適用可能であり、証明可能な安定性を持つ完全離散スキームを達成できるか？
• RQ2 安定性および精度の観点から、モーダルフィルタリングは人工拡散およびスペクトル粘性と比べてどのように異なるか？
• RQ3 解のフィルタリングと時間微分のフィルタリングの違いが、数値的安定性および振動制御に与える影響は何か？
• RQ4 時間刻みと解の挙動に応じて自動的にフィルタ強度を調整できる適応的戦略を設計できるか？
• RQ5 提案されたフィルタリング戦略は、人工拡散法で観察された時間刻み制限を解消するか？

主な発見

• モーダルフィルタリングは、時間に関して1次精度までにおいてスペクトル粘性と理論的に等価であり、人工拡散の代替として安定な手法を提供する。
• 適応的フィルタ強度戦略は、$(\Delta t)^2$ の誤差項を効果的に補償し、明示的エーラー時間積分を用いた証明可能な安定な完全離散スキームの実現を可能にする。
• 時間微分のフィルタリングは、$\underline{\underline{M}}\underline{\underline{F}}^{-1}$ が誘導するノルムの下で安定な半離散化をもたらすが、全時間ステップにおける逆フィルタの適用がこの安定性を破壊するため、解のフィルタリングに比べて劣る。
• 数値実験では、解のフィルタリングが振動を低減させ、エネルギーの振るまいを改善するのに対し、時間微分のフィルタリングはエネルギーノルムの振動を引き起こし、理論的予想を確認する。
• 不連続な初期データを有する線形移流問題では、スキームは安定かつ正確に保たれるが、衝撃形成を伴う非線形問題では振動が持続する。
• 提案されたフレームワークにより、SBP-CPRスキームの非線形問題への適用範囲が拡張され、特に適応的フィルタリング機構のおかげで耐障害性が向上する。