Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

技術要約：流体・構造連成問題における仮想領域フレームワークを用いた非整合格子上の求積誤差評価

問題設定
本論文は、分布ラグランジュ乗数法に基づく仮想領域定式化を用いた、流体・構造連成（FSI）問題の数値離散化を取り扱う。このフレームワークでは、流体と固体は、完全に独立した非整合な2つのメッシュを用いて離散化される。固体は、基準領域 $B$ 上のラグランジュ的な枠組みで記述され、固定されたオイラー的な流体領域 $\Omega$ 内の現在の位置へと写像される。流体の速度と固体の変形との間の運動学的結合は、固体基準領域上で定義されたラグランジュ乗数 $\lambda$ を通じて強制される。

計算上の重要な課題は、結合（界面）行列のアセンブリにおいて発生する。この行列は、流体の基底関数（固体写像によって合成されたもの）と固体の基底関数の積を、固体要素上で積分する必要がある。これらの関数は非整合な格子上に定義されているため、厳密な積分を行うには、合成された複合求積規則を構築するために、流体メッシュと写像された固体メッシュの幾何学的な交差を計算する必要がある。このプロセスは計算コストが高く、特に結合項は各タイムステップでアセンブリされなければならない。あるいは、代替的な手法として、メッシュの交差を考慮せずに標準的な求積規則を用いて固体要素上で直接積分する方法があるが、これには求積誤差が生じる。本論文は、この近似的な積分戦略に関連する理論的妥当性と誤差評価について調査している。

手法
著者らは、定常状態のストークスに基づくFSIモデルを分析しているが、これらの結果はより一般的な非線形ケースに対しても代表的なものであると述べている。変分形式は、流体の速度 $u$、圧力 $p$、固体の変形写像 $X$、およびラグランジュ乗数 $\lambda$ を求め、サドルポイント系を満たすものである。

本研究は、結合 $c(\cdot, \cdot)$ を定義するために使用される2つの具体的な選択肢に焦点を当てている：

1. $c_0$ (双対ペアリング): $H^1(B)$ の双対としての乗数空間に基づく、$L^2(B)$ 内部積に基づくもの。
2. $c_1$ (スカラー積): $H^1(B)$ 自体を乗数空間とする、$H^1(B)$ 内部積に基づくもの。

手法は以下のステップで進行する：

• 離散化: すべての変数に対して有限要素空間を導入する。流体にはBercovier–Pironneau要素（$P1$-iso-$P2/P1$）を用い、固体の写像と乗数には連続な1次要素を用いる。
• 近似積分: 著者らは、合成された流体基底関数が固体メッシュに対して不連続な区分的多項式である事実を無視して、固体メッシュ上の標準的な求積規則を用いて結合項を計算する、離散的な双線形形式を定義する。
• 抽象的誤差解析: inf-sup条件の枠組みに従い、著者らは、厳密な離散解と、不正確な積分を用いて得られた解との差を抑える一般的な誤差評価（定理2）を導出する。この境界は、求積誤差汎関数に依存する。
• 求積誤差の推定: 核心となる技術的貢献は、結合項の具体的な誤差評価を証明することである。著者らは、合成された関数 $v_h(X)$ が連続かつ固体要素上で区分的に線形であることを利用し、特定の分数ソボレフ空間（$H^{1+s}$）における正則性特性を利用する。彼らは、誤差を抑えるために補間理論と逆不等式を用いている。
• 安定性の証明: 結合項が近似的に構成される場合でも、特定のメッシュサイズに関する条件が満たされる限り、離散問題が適切に設定されている（離散inf-sup条件を満たす）ことを証明する。

主な貢献と結果

1. 近似結合による適切設定: 本論文は、結合行列が近似的な求積を用いてアセンブリされる場合でも、メッシュの交差を必要とせず、離散問題が安定かつ適切に設定されていることを証明している。
2. 求積誤差評価:
   • $L^2$ ベースの結合 ($c_0$) について、求積誤差は $O(h_B^{3/2} |\log h_B|)$ で抑えられる。著者らは、固体メッシュサイズ $h_B \to 0$ となる場合、近似的な手法が最適な収束率を与えることを示している。
   • $H^1$ ベースの結合 ($c_1$) について、勾配項があるため誤差解析はより複雑である。誤差境界には、メッシュサイズの比に依存する項が含まれる。具体的には、最適な収束には $h_B \to 0$ だけでなく、$h_B/h_\Omega \to 0$（ここで $h_\Omega$ は流体メッシュサイズ）という条件も必要となる。もし $h_B/h_\Omega$ が一定に保たれる場合、手法は最適性を失う可能性がある。
3. 逆不等式: 誤差評価、特に双対ノルムのケースにおけるラグランジュ乗数のノルムを扱うために、特定の逆不等式を導出し、利用している。
4. 数値検証: 理論的評価を検証するために、4つの数値テストが提示されている：
   • テスト 1 & 2: 異なるメッシュ細分化戦略（一定比率 vs 減少する比率）の下での $c_0$ と $c_1$ の収束率を比較する。結果は、$c_0$ が近似積分に対して頑健である一方で、$c_1$ は精度を維持するために $h_B/h_\Omega \to 0$ を必要とすることを裏付けている。
   • テスト 3: 界面を横切る速度勾配のジャンプ（FSIに典型的）を持つ解を使用し、解の特異性に起因する期待される劣最適な収束率 $O(h^{1/2})$ を確認する。
   • テスト 4: 高次要素（$P2/P1$）を調査し、$c_0$ は依然として頑健であるが、$c_1$ は最適な収束率を達成するために厳密な積分または特定のメッシュ比を依然として必要とすることを示す。

意義と主張
本論文は、計算効率のためにしばしば採用される、近似積分を仮想領域FSI法で使用することに対する厳密な理論的正当化を提供すると主張している。著者らは以下を実証している：

• 「不正確な」アプローチは、数値スキームの安定性を損なわない。
• $L^2$ ベースの結合 ($c_0$) は、精度を維持するために流体メッシュと固体メッシュを特定の比率で細分化する必要がないため、近似積分に特に適している。
• $H^1$ ベースの結合 ($c_1$) は、異なる安定特性を提供する可能性があるものの、求積誤差が解の誤差を支配することを避けるために、より厳しいメッシュ要件（$h_B/h_\Omega \to 0$）を課す。

著者らは、これらの理論的結果が、厳密な（メッシュ交差を行う）アセンブリ技術と近似的な（非交差の）アセンブリ技術のどちらを選択すべきかについて、特定の結合形式とメッシュ構成に基づいた判断を行うための強固な数学的基礎を提供していることを確認した。