[論文レビュー] Stability and error analysis of a splitting method using Robin-Robin coupling applied to a fluid-structure interaction problem

Erik Burman, R. Durst+1 more|arXiv (Cornell University)|Nov 15, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 19被引用数 5

ひとこと要約

本稿では、流体構造連成（FSI）問題に対するロビン-ロビン分割法の安定性および誤差解析を提示している。明示的で緩い結合スキームを用い、流体-構造界面でロビン伝達条件を適用する。主な貢献は、ロビンパラメータ $\lambda$ が $O(1)$ のとき、空間離散化に依存しない時間 $T$ におけるグローバル誤差が $O(\sqrt{T\Delta t})$ であることを示す厳密な誤差推定であり、$h^{-1}$ スケーリングの安定化項を必要とせず収束を証明する。

ABSTRACT

We analyze a splitting method for a canonical fluid structure interaction problem. The splittling method uses a Robin-Robin boundary condition, explicit strategy. We prove the method is stable and, furthermore, we provide an error estimate that shows the error at the final time $T$ is $O(\sqrt{TΔt})$ where $Δt$ is the time step.

研究の動機と目的

• 追加質量効果に起因する古典的な緩い結合スキームの不安定性に対処すること。
• 数値的に安定であることが観察されたが理論的根拠の欠如していた、追加の圧力安定化を要しないロビン-ロビンカップリング法の分析。
• 時間離散化された分割スキームに対する厳密な誤差推定を確立し、空間離散化に依存しない収束を示すこと。
• $h^{-1}$ スケーリングの安定化項を回避することで、時間ステップサイズの制限や反復補正の必要性を排除すること。

提案手法

• 本手法は、各時間ステップで流体および固体のPDEを順次解く時間分割戦略を採用し、界面でロビン-ロビン伝達条件を適用する。
• ロビン-ロビンカップリングは、法線方向速度と応力の界面フラックスを用い、パラメータ $\lambda$ がカップリング強度を制御する。
• 解析では、局所PDE解の十分な正則性を仮定し、空間離散化を避けて時間離散化誤差に焦点を当てる。
• 重要な技術として、修正されたエネルギーノルムを用いたエネルギー推定と、G. ベイカーの手法にインspiredされた方法により、誤差境界の指数的増大を回避する。
• 打切り誤差は、正確な解と分割スキームを比較することで導出され、安定性と組み合わせてグローバル誤差を評価する。
• 誤差蓄積を制御するために、パラメータ $\delta = \Delta t / (2T)$ をヤングの不等式に用いる。

実験結果

リサーチクエスチョン

• RQ1 流体構造連成問題において、追加の圧力安定化を要しないロビン-ロビン緩い結合スキームは安定か？
• RQ2 空間離散化に依存しない時間離散化されたロビン-ロビン分割法に対して、厳密な誤差推定を導出できるか？
• RQ3 誤差は時間とともに指数的に増大するのか、それとも $\sqrt{T\Delta t}$ のように増大するのか？ そしてエネルギー法を用いてこれを証明できるか？
• RQ4 ロビンパラメータ $\lambda = O(1)$ のとき、$O(\sqrt{T\Delta t})$ の誤差境界を達成できるか？ これにより $h^{-1}$ スケーリングの安定化パラメータの必要性が回避されるか？
• RQ5 反復補正が不要な状況下でも、時間離散化スキームの収束が保証されるか？

主な発見

• ロビン-ロビン分割スキームのグローバル誤差は、適切なエネルギーノルムにおいて $O(\sqrt{T\Delta t})$ で有界であり、定数は空間離散化に依存しない。
• 誤差推定は $\sqrt{T}$ の増大を示し、指数的増大よりも顕著な改善である。
• 解析により、圧力安定化や $h^{-1}$ スケーリングのロビンパラメータを要せず、時間離散化スキームの安定性と収束性が証明された。
• パラメータ $\lambda = O(1)$ で収束が達成され、非常に小さな時間ステップや反復補正手順の必要性が排除された。
• 誤差境界は、G. ベイカーのアプローチにインspiredされた新規エネルギー技術を用いて導出され、指数的増大を回避した。
• この結果は、長時間シミュレーションに対して本手法が頑健であり、効率的な空間離散化選択を可能にする。