[論文レビュー] Regularized Stokes Immersed Boundary Problems in Two Dimensions: Well-posedness, Singular Limit, and Error Estimates

Jiajun Tong|arXiv (Cornell University)|Apr 21, 2019

Lattice Boltzmann Simulation Studies参考文献 39被引用数 4

ひとこと要約

本稿では、特異な力および流れの結合を滑らかにするために、滑らかなデルタ関数を用いた2次元ストークス埋め込み界面問題の正則化版を導入する。グローバルな適切性を確立し、正則化パラメータが0に近づくにつれて正則化されていない問題への正則化されたストリング運動の収束を証明し、適切に選ばれた滑らかさ関数を用いることで精度が向上するという誤差推定を導出する。

ABSTRACT

Inspired by the numerical immersed boundary method, we introduce regularized Stokes immersed boundary problems in two dimensions to describe regularized motion of a 1-D closed elastic string in a 2-D Stokes flow, in which a regularized $δ$-function is used to mollify the flow field and singular forcing. We establish global well-posedness of the regularized problems, and prove that as the regularization parameter diminishes, string dynamics in the regularized problems converge to that in the Stokes immersed boundary problem with no regularization. Viewing the un-regularized problem as a benchmark, we derive error estimates under various norms for the string dynamics. Our rigorous analysis shows that the regularized problems achieve improved accuracy if the regularized $δ$-function is suitably chosen. This may imply potential improvement in the numerical method, which is worth further investigation.

研究の動機と目的

• 2次元ストークス埋め込み界面問題の連続的正則化定式化を厳密に分析すること。
• 2次元ストークス流れにおける1次元弾性ストリングに対して、正則化系のグローバルな適切性を確立すること。
• 正則化パラメータが0に近づくにつれて、正則化されたストリング運動が正則化されていない問題に収束することを正当化すること。
• 正則化された解と正則化されていない解の間の定量的誤差推定を導出し、正則化パラメータおよび滑らかさ関数の選択に依存することを示すこと。
• 解析的誤差境界に基づいて、数値的埋め込み界面法の潜在的な改善策を提示すること。

提案手法

• ディラックのデルタ関数の代わりに、滑らかで compact な台を持つ滑らかさ関数を用いて、特異な力および速度場を正則化する2次元ストークスIB問題の定式化を行う。
• 流体変数を消去することでストリング運動の輪郭ダイナミクス方程式を導出し、2次元ストークス方程式の基本解を含む特異積分としてストリング速度を表現する。
• 周期的ソボレフ空間におけるホルダーおよびソボレフ型推定を用いて、ストリング配置およびその導関数の差を制御する。
• 特徴線法および積分推定を用いて、時間経過に伴う正則化および正則化されていないストリングダイナミクスの差を評価する。
• 速度カーネルの差に対する鋭い点での推定および$L^2$型推定を適用し、ストリング発展における誤差を制御する。
• 正則化パラメータ$ρ$、滑らかさ関数の性質、および解の正則性に明示的に依存する、さまざまなノルム（例：$L^2$、$;tilde{H}^1$）における誤差境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 初期ストリング配置が$H^{5/2}$に属する場合、正則化された2次元ストークス埋め込み界面問題はグローバルな解をもつか？
• RQ2 正則化パラメータ$ρ \to 0$のとき、正則化された問題の解は正則化されていない問題の解に収束するか？
• RQ3 正則化された解と正則化されていない解の間の誤差推定は、滑らかさ関数の選択にどのように依存するか？
• RQ4 正則化された定式化は、標準的な数値的埋め込み界面法よりも高い精度を達成できるか？
• RQ5 誤差を最小化する最適な正則化されたデルタ関数は何か？

主な発見

• 初期ストリング配置が$H^{5/2}(\bbT)$に属する場合、正則化問題はグローバルに適切であり、すべての時間において解の存在と一意性が保証される。
• 正則化パラメータ$ρ \to 0$のとき、正則化された問題におけるストリングダイナミクスは、$L^2$および$;tilde{H}^1$ノルムにおいて、正則化されていない問題の解に収束する。
• 誤差推定が導出され、ストリングダイナミクスの差が$C\lambda^{-2}\|Y_{1}^{\prime\prime}\|_{\dot{H}^\theta}^2$および高次項によって上界で抑えられることを示している。ここで$λ$は滑らかさ関数の幅を制御する。
• 誤差境界は滑らかさ関数の選択に強く依存する：より高い滑らかさを持つ滑らかさ関数（例：高次の滑らかさ）は、収束速度の向上をもたらす。
• 解析により、正則化された定式化が標準的な数値法よりも高い精度を達成できることを示しており、手法の改善への道筋を示唆する。
• より鋭い$̅dot{H}^1$型ノルムにおける誤差推定が導出され、$\|g_{Y_1} - g_{Y_2}\|_{\dot{H}^1} \leq C\lambda^{-3}(\|Y_1\|_{\dot{H}^{2+\theta}} + \|Y_2\|_{\dot{H}^{2+\theta}})^3\|\delta Y\|_{\dot{H}^{3/2}}$となる。これは弱い代替推定よりも鋭い。