[論文レビュー] Velocity‑vorticity‑pressure formulation for the Oseen problem with variable viscosity

Verónica Anaya, Rubén Caraballo+3 more|arXiv (Cornell University)|Apr 20, 2022

Advanced Numerical Methods in Computational Mathematics参考文献 57被引用数 14

ひとこと要約

本稿では、可変粘性を伴うOseen問題に対して、速度・渇度・圧力の混合有限要素法を提案し、最小二乗安定化を用いて非定数粘性および対流項を扱う。この手法は、任意のストークス適合速度-圧力ペアおよび任意の次数の渇度空間に対してinf-sup安定性を保証し、適応メッシュ・リファインメント下でも最適収束率と頑健な性能を達成する。数値実験により、複雑な幾何形状および可変粘性場を含む状況でも検証された。

ABSTRACT

We propose and analyse an augmented mixed finite element method for the Oseen equations written in terms of velocity, vorticity, and pressure with non-constant viscosity and homogeneous Dirichlet boundary condition for the velocity. The weak formulation includes least-squares terms arising from the constitutive equation and from the incompressibility condition, and we show that it satisfies the hypotheses of the Babuška-Brezzi theory. Repeating the arguments of the continuous analysis, the stability and solvability of the discrete problem are established. The method is suited for any Stokes inf-sup stable finite element pair for velocity and pressure, while for vorticity any generic discrete space (of arbitrary order) can be used. A priori and a posteriori error estimates are derived using two specific families of discrete subspaces. Finally, we provide a set of numerical tests illustrating the behaviour of the scheme, verifying the theoretical convergence rates, and showing the performance of the adaptive algorithm guided by residual a posteriori error estimation

研究の動機と目的

• 非定数粘性のため、標準的な渇度定式化が破綻する状況下でも、安定かつ収束する混合有限要素法をOseen問題に構築すること。
• 従来、Brinkman方程式およびストークス方程式に用いられてきた拡張混合定式化を、対流項および可変粘性を有するOseen問題へ拡張すること。
• 粘性が空間的に変化する場合でも、Babuška-Brezzi理論に基づき、離散問題の安定性と可解性を保証すること。
• 2種類の有限要素部分空間族について、事前および事後誤差見積もりを導出し、適応的リファインメント戦略を可能にすること。

提案手法

• 粘性が非定数である場合のOseen問題を、速度、渇度、圧力の変数として定式化し、構成方程式および非圧縮性条件からの最小二乗項を導入して変分定式化を安定化する。
• 強い形から導かれる安定化項を有する拡張混合有限要素法を採用し、混合法におけるBabuška-Brezzi理論に適合する。
• 渇度方程式を運動量保存の回転から導出し、ひずみ率テンソル ε(u) = 1/2(∇u + ∇u^T) を用いて可変粘性に適応した弱形式を採用する。
• ガラーキン最小二乗法を用いてシステムを安定化し、特に可変粘性および対流項 (β·∇)u に起因する課題に対処する。
• 残差に基づく事後誤差推定器を用いて適応的メッシュリファインメントを駆動し、再入角や境界層など誤差が大きい領域に注目する。
• 誤差指標に従い、ブルー・グリーンリファインメント戦略を採用し、特異性や急勾配が存在する場合でも最適収束を達成する。

実験結果

リサーチクエスチョン

• RQ1 可変粘性を伴うOseen問題に対して、速度、渇度、圧力を主変数として用いる安定かつ収束する混合有限要素法を構築可能か？
• RQ2 非定数粘性下では恒等式 curl(curl v) = −Δv + ∇(div v) が直接適用できない場合に、標準的な渇度に基づく定式化をどのように拡張できるか？
• RQ3 粘性が空間的に変化し、かつ対流項が存在する状況下でも、inf-sup安定性と最適収束を保証する安定化技術は何か？
• RQ4 事後誤差見積もりが、再入角を持つL字型領域のような特異領域において、最適収束率を回復するための効果的な適応的メッシュリファインメントをどのように駆動できるか？
• RQ5 臨床的応用に近い、可変粘性を有する患者特異的降主动脈アーチモデルのような現実的で複雑な幾何形状において、この手法はどのように性能を発揮するか？

主な発見

• 提案された拡張混合有限要素法は、任意のストークス安定な速度-圧力ペアおよび任意の次数の渇度近似に対して、inf-sup条件を満たし、離散問題の一意解を保証する。
• 事前誤差見積もりにより、すべての変数で最適収束率が確認された：速度はH1ノルム、渇度はL2ノルム、圧力はL2ノルムにおいて、構造的および非構造的メッシュ上でも実証済み。
• 事後誤差推定器の効果的インデックスは1に近く（約1.12–1.16）、信頼性の高い誤差推定と効果的な適応的リファインメントを示している。
• 再入角を持つL字型領域では、均一リファインメントが特異性のため最適収束率を達成できないのに対し、適応的リファインメントにより最適収束率が回復された。
• 3次元降主動脈アーチシミュレーションでは、境界層や渇度構造といった複雑な流れの特徴を的確に捉え、適応的リファインメントにより下降大動脈および壁付近に要素が集中し、臨床的血行動態応用における頑健性を確認した。