[論文レビュー] An integral model based on slender body theory, with applications to curved rigid fibers

Helge I. Andersson, Elena Celledoni+3 more|arXiv (Cornell University)|Dec 21, 2020

Structural Analysis and Optimization参考文献 60被引用数 15

ひとこと要約

本稿では、非局所な細長い体理論に基づき、半径に依存する滑らかなカーネルを導出することで、粘性流れ内の剛体で曲がった細長い繊維の動的挙動をシミュレートする新しい積分モデルを提示する。この手法は、漸近的に整合性のある正則化を用いて、第一種フリードホルム積分方程式を第二種方程式に変換し、数値的安定性と可逆性を保証する。これにより、検証された収束性とスペクトル特性を有する、複雑な繊維幾何形状の高速かつ高精度なシミュレーションが可能となる。

ABSTRACT

We propose a novel integral model describing the motion of curved slender fibers in viscous flow, and develop a numerical method for simulating dynamics of rigid fibers. The model is derived from nonlocal slender body theory (SBT), which approximates flow near the fiber using singular solutions of the Stokes equations integrated along the fiber centerline. In contrast to other models based on (singular) SBT, our model yields a smooth integral kernel which incorporates the (possibly varying) fiber radius naturally. The integral operator is provably negative definite in a non-physical idealized geometry, as expected from PDE theory. This is numerically verified in physically relevant geometries. We propose a convergent numerical method for solving the integral equation and discuss its convergence and stability. The accuracy of the model and method is verified against known models for ellipsoids. Finally, a fast algorithm for computing dynamics of rigid fibers with complex geometries is developed.

研究の動機と目的

• ストークス流れ内における剛体で曲がった繊維を安定的かつ高精度かつ計算的に効率的なモデルでシミュレートすること。
• 細長い体理論における第一種積分方程式の不適切な定式化を、証明可能に安定化する正則化手法を導入することで解決すること。
• 追加の任意パラメータを導入することなく、繊維の断面半径を積分カーネルに自然に組み込むこと。
• 既知の回転楕円体の力学と比較してモデルを検証し、複雑な曲がった繊維幾何形状への応用を拡張すること。
• 繊維の剛性と行列・ベクトル演算を活用することで、時間ループ内での行列積に基づく高速な動的シミュレーションを可能にすること。

提案手法

• 繊維表面に沿った角項の漸近的キャンセレーションを用いて、非局所な細長い体理論から滑らかな積分カーネルを導出する。
• 力密度のための第一種フリードホルム方程式を第二種方程式に変換する正則化技術を導入し、数値的条件数の改善と可逆性の向上を図る。
• 第二種積分方程式を解くために、高次精度の数値積分法（Nyström法）を採用し、収束性と精度を保証する。
• 時間ループ内での行列・ベクトル積に基づく高速アルゴリズムを用い、繊維の剛性を活用して動的シミュレーションを加速する。
• 直線からの曲がりの小さな変動を有する繊維をシミュレートし、回転動的挙動を直線繊維の基準ケースと比較する。
• Jeffery [14] の標準的な回転楕円体モデルと比較することで、力、トルク、角運動量の観点からモデルを検証する。

実験結果

リサーチクエスチョン

• RQ1 特異なストークスレットとダブルレットから、滑らかで半径に依存する積分カーネルをどのように導出できるか？
• RQ2 得られた積分作用素のスペクトル的性質は何か？正則化は負定値性と安定性を保証するか？
• RQ3 第一種と第二種の定式化の間で、モデルの精度と条件数はどのように比較されるか？
• RQ4 直線繊維からの小さな曲がりのずれが回転動的挙動にどの程度影響を及ぼし、その影響を定量的に評価できるか？
• RQ5 複雑な形状を有する剛体繊維の動的挙動をシミュレートするための、高速で安定的かつ収束性を有する数値手法を開発できるか？

主な発見

• 正則化を施した第二種積分方程式の定式化により、細かいメッシュでもスペクトルがゼロから離れて保たれ、数値的可逆性と安定性が保証される。
• 数値的テストにより、空間方向で2次収束を示すことが確認され、第一種定式化に比べて条件数が改善されている。
• 直線的で周期的である繊維に対しては、積分作用素が証明可能に負定値であることが示され、PDE理論および期待される物理的挙動と整合的である。
• 曲がった繊維の角運動量および姿勢のずれは、曲がりの大きさδに比例して線形に増加し、100時間単位のシミュレーションでδ = 0.0003では1%のずれ、δ = 0.0015では7.5%のずれを示す。
• 1回転後の繊維の長軸とx軸のなす角度θは、δ = 0.0003では約3°、δ = 0.0015では約8°ずれることが確認され、測定可能な回転ずれが存在することが裏付けられた。
• 本手法はJefferyの楕円体の力学を高精度に再現でき、複雑な曲がった繊維への応用に移る前の標準的ケースにおけるモデルの正確性が検証された。