[論文レビュー] On a finite element approximation of the Stokes problem under leak or slip boundary conditions of friction type

Takahito Kashiwabara|arXiv (Cornell University)|Dec 22, 2010

Contact Mechanics and Variational Inequalities参考文献 18被引用数 5

ひとこと要約

本稿では、摩擦型のスリップまたはリーク境界条件を伴うストークス方程式に対する有限要素法を提示する。変分不等式の定式化により、解の存在と一意性を保証する。正則性条件下で収束速度が $O(h^2)$ まで達成され、数値解法としてUzawa型反復アルゴリズムを提案し、数値例により最適収束が確認された。

ABSTRACT

A finite element approximation of the Stokes equations under a certain nonlinear boundary condition, namely, the slip or leak boundary condition of friction type, is considered. We propose an approximate problem formulated by a variational inequality, prove an existence and uniqueness result, present an error estimate, and discuss a numerical realization using an iterative Uzawa-type method. Several numerical examples are provided to support our theoretical results.

研究の動機と目的

• 非線形摩擦型境界条件の下でのストークス問題に対する包括的な有限要素理論の構築。
• スリップまたはリーク境界条件の摩擦型を伴うストークス方程式の有限要素近似における存在性、一意性、誤差推定の確立。
• 離散変分不等式問題の数値解法を目的としたUzawa型反復アルゴリズムの設計および解析。
• さまざまなパrameter設定下での数値実験を通じて理論的収束速度の妥当性の検証。
• リーク境界条件における圧力の非一意性とその誤差推定への影響の取り扱い。

提案手法

• スリップおよびリーク境界条件の両方を、摩擦型の変分不等式によって連続問題として定式化する。
• 2次元多角形領域において、P2/P1有限要素ペア（テイラー＝フッド要素）を用いて問題を離散化する。
• 変分不等式理論を用いて、SBCFおよびLBCFの両ケースにおいて離散解の存在と一意性を証明する。
• 誤差推定を導出：$H^{1+\epsilon}$ 正則性下で $O(h^{\text{min}{\epsilon,1/4\u007d})$、追加の接線速度の符号条件があると $O(h^{\text{min}\u007b\epsilon,1\u007d})$ に改善される。
• 離散変分不等式を解くためのUzawa型反復アルゴリズムを提案し、反復列が解に収束することを証明する。
• LBCFにおける圧力の非一意性に対処するため、基準点を用いて圧力定数を固定し、誤差比較を可能にする。

実験結果

リサーチクエスチョン

• RQ1 摩擦型スリップまたはリーク境界条件を伴うストークス方程式に対して、有限要素法を厳密に定式化・解析することは可能か？
• RQ2 このような非線形境界条件下での有限要素近似の収束速度はいかほどか？
• RQ3 離散問題を効率的かつ信頼性高く反復アルゴリズムで解くにはどうすればよいか？
• RQ4 リーク境界条件の場合は、圧力場に加法的定数が存在するが、その影響は誤差解析にどのように現れるか？
• RQ5 適切な正則性および接線速度の符号条件のもとで、最適収束速度を達成できるか？

主な発見

• 適切な有限要素離散化と変分不等式定式化のもとで、スリップおよびリーク境界条件の両方において離散問題は一意解を有する。
• スリップ境界条件の誤差推定は $O(h^{\text{min}\u007b\epsilon,1\u007d})$ であり、接線速度に有利な符号条件があると $O(h^{\epsilon})$ に改善され、$\epsilon=2$ のとき最適な $O(h^2)$ 収束が達成される。
• リーク境界条件では、圧力定数の不確実性のため誤差推定が $O(h^{\text{min}\u007b\epsilon/2,1/4\u007d})$ となるが、圧力定数を適切に処理すれば $O(h^2)$ に改善される。
• 数値実験により理論的収束速度が確認され、$N=120$ の参照解と比較して、SBCFおよびLBCFの両ケースで $O(h^2)$ 収束が観察された。
• Uzawa型反復アルゴリズムは離散解に収束し、収束速度は問題パラメータおよび初期推定値に依存するが、数値実験による反復回数の結果からその有効性が示された。
• 境界上のラグランジュ乗数 $\lambda_h$ は、$u_h(M)=0$ のとき $|\lambda_h(M)| \leq 1$ を満たし、それ以外の場合は $\lambda_h(M) = \pm 1$ となる。これは理論的境界と整合的である。