[論文レビュー] Theoretical justification and error analysis for slender body theory

Yoichiro Mori, Laurel Ohm+1 more|arXiv (Cornell University)|Jun 30, 2018

Micro and Nano Robotics参考文献 40被引用数 5

ひとこと要約

この論文は、ストークス流れにおける細長い体理論（SBT）の厳密な理論的基盤を確立する。細い繊維上での混合境界条件を有する新しいPDE問題を定式化し、適切な解の存在と一意性を証明した。さらに、力密度および中心線幾何の弱い正則性仮定のもとで、SBT近似が真の解へ、繊維半径に比例する速度で収束することを示す誤差推定を導出。収束速度には対数補正項が含まれる。

ABSTRACT

Slender body theory facilitates computational simulations of thin fibers immersed in a viscous fluid by approximating each fiber using only the geometry of the fiber centerline curve and the line force density along it. However, it has been unclear how well slender body theory actually approximates Stokes flow about a thin but truly three-dimensional fiber, in part due to the fact that simply prescribing data along a one-dimensional curve does not result in a well-posed boundary value problem for the Stokes equations in $\mathbb{R}^3$. Here, we introduce a PDE problem to which slender body theory (SBT) provides an approximation, thereby placing SBT on firm theoretical footing. The slender body PDE is a new type of boundary value problem for Stokes flow where partial Dirichlet and partial Neumann conditions are specified everywhere along the fiber surface. Given only a 1D force density along a closed fiber, we show that the flow field exterior to the thin fiber is uniquely determined by imposing a fiber integrity condition: the surface velocity field on the fiber must be constant along cross sections orthogonal to the fiber centerline. Furthermore, a careful estimation of the residual, together with stability estimates provided by the PDE well-posedness framework, allow us to establish error estimates between the slender body approximation and the exact solution to the above problem. The error is bounded by an expression proportional to the fiber radius (up to logarithmic corrections) under mild regularity assumptions on the 1D force density and fiber centerline geometry.

研究の動機と目的

• 1Dの力密度のみを指定しても3次元のストークス問題として適切に定式化できないことから、長年の理論的曖昧性を解消し、SBTが近似する適切なPDE問題を定義すること。
• 3次元ストークス問題として適切に定式化されない1D力密度指定の問題を解消すること。
• 1Dデータ（力密度および中心線幾何）のみを用いて、『真の』細い繊維問題の数学的に整合した定義を提供すること。
• 提案されたPDE定式化の正確な解とSBT近似との間の定量的誤差推定を導出すること。
• 粘性流体中での細い繊維の計算シミュレーションに広く用いられるSBTの正当性を、収束保証を伴って裏付けること。

提案手法

• 繊維表面に部分的ディリクレ（スライプなし）および部分的ノイマン（応力）境界条件を課すストークス流れの新しいPDE定式化を導入する。
• 繊維の整合性条件を課す：中心線に直交する断面において、表面速度が一定でなければならない。これにより解の一意性が保証される。
• PDEの適切な定式化フレームワークを用いて、繊維半径εに依存しない定数を用いて、細長い体PDEの解の存在と一意性を確立する。
• 解を調和部と補正項に分解することで、補正された速度場を構成する。
• トレース不等式、コルンの不等式、ソボレフ埋め込みを用いて、SBT近似の残差を繊維半径εの関数として推定する。
• 比較場ũを構成し、その真の解からのずれを安定性および残差推定を用いて評価することで、誤差境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 1Dの力データのみでは3次元ストークス問題として適切に定式化できないことから、細長い体理論が実際に近似しているPDE問題は何か？
• RQ2 1Dの力密度および中心線幾何のみを用いて、細い繊維問題の解を一意に定義することは可能か？
• RQ3 細長い体近似の真の3次元ストークス問題の解に対する収束速度は何か？
• RQ4 PDEフレームワークにおける適切な定式化の定数は繊維半径εに依存するか？ε → 0の極限でも安定か？
• RQ5 SBT近似と正確な解との間の定量的誤差は何か？また、これは繊維半径にどのように依存するか？

主な発見

• 混合境界条件および繊維の整合性条件を備えた細長い体PDEは、任意の与えられた1D力密度および中心線幾何に対して一意な解を有する。
• 繊維半径εに依存しない定数を用いてPDEの適切な定式化が確立され、ε → 0の極限でも安定性が保証される。
• 細長い体近似の残差は、対数補正項を除けば、εに比例する項で有界である。
• SBT近似とPDEの正確な解との間の誤差は、Cε|f|_H^1で有界であり、Cはεに依存せず、fは力密度を表す。
• 誤差推定は、構築された比較場ũを用い、PDEフレームワークからの安定性推定および残差解析に依拠している。
• 力密度および中心線の正則性仮定のもとで収束速度は最適であり、SBTの漸近的精度が裏付けられる。