[論文レビュー] Accurate evaluation of integrals in slender-body formulations for fibers in viscous flow

Anna‐Karin Tornberg|arXiv (Cornell University)|Dec 23, 2020

Micro and Nano Robotics参考文献 9被引用数 4

ひとこと要約

本稿では、ストークス流れにおける柔軟繊維の非局所的細長い体理論に現れる有限部積分およびほぼ特異的積分を高精度で評価する数値法を提示する。積分法と事前に計算された修正された求積重みを組み合わせることで、追加コストなしに特異的積分に対して分光的精度を達成する。また、類縁の技法により、繊維表面付近での高精度なほぼ特異的ストークスレット積分を処理する。

ABSTRACT

A non-local slender body approximation for slender flexible fibers in Stokes flow can be derived, yielding an integral equation along the center lines of the fibers that involves a slenderness parameter. The formulation contains a so-called finite part singular integral, and can in the case of several fibers or evaluation of the flow field require the evaluation of nearly singular integrals. We introduce a numerical technique to accurately and efficiently evaluate the finite part integral. This technique can be applied combined with any panel based quadrature rule and will add no additional cost except for a small precomputation of modified quadrature weights. We also show how a related technique that was recently introduced can be applied for the evaluation of the nearly singular integrals.

研究の動機と目的

• 非局所的細長い体理論における細長い繊維のストークス流れにおける有限部積分を評価する、堅牢で効率的な数値法の開発。
• 繊維同士が近接している、または評価点が繊維中心線に近い場合に特異的またはほぼ特異的な積分を正確に計算する課題に対処する。
• 特異性を解析的に取り扱う求積技法を提供し、計算コストの増加を最小限に抑えつつ高精度を保証する。
• 特に繊維表面付近での流体速度評価に生じるほぼ特異的積分に対しても、同様の手法を拡張する。
• ヘリシスや円形ループなどのベンチマーク幾何形状に対して検証を行い、標準的な求積法と比較して優れた収束性と誤差制御を示す。

提案手法

• 被積分関数を滑らかな多項式因子と特異的カーネルに分解することで、有限部積分を再定式化し、各パネル上で解析的積分を可能にする。
• 滑らかな因子の多項式近似を用い、その係数は事前計算段階で転置 Vandermonde 方程式を一度だけ解くことで決定し、修正された求積重みを生成する。
• 生成された修正された求積重みはすべてのパネルで再利用可能であり、評価ごとに追加コストをかけずに、特異的積分を効率的かつ高精度に評価可能である。
• ほぼ特異的ストークスレット積分に対しては、同様の積分法を適用し、根の抽出から得られる複素共役ペアを用いて被積分関数を分割する。
• 得られた積分の解析的評価には再帰公式を用い、数値的実装に注意を払い、精度を維持する。
• 本手法はガウス＝ルジャンドル求積則を用いて実装され、高精度な基準解を用いた3次元繊維幾何形状での検証が行われた。

実験結果

リサーチクエスチョン

• RQ1 被積分関数が対角線上で特異的である場合、細長い体理論における有限部積分をどのように高精度に評価できるか。
• RQ2 特異性を解析的に取り扱いながらも、効率性を保ち、追加の計算コストを伴わない求積法を設計できるか。
• RQ3 評価点が繊維中心線に近い場合、本手法はほぼ特異的積分に対してどのように性能を発揮するか。
• RQ4 異なる繊維幾何形状および力密度関数に対して、本手法の収束特性はどのように変化するか。
• RQ5 自由端を持つ複数の相互作用する繊維が3次元ストークス流れに存在する場合、本手法は適用可能か。

主な発見

• 提案手法は、事前に計算された修正された求積重みを用いて特異性を解析的に取り扱うことで、有限部積分に対して分光的精度を達成する。
• ヘリカル繊維の周辺での流体速度評価において、同数量のパネルを用いた標準的なガウス＝ルジャンドル求積法と比較して、最大誤差が桁違いに低減される。
• $ M = 8 $ パネルの場合、特殊な求積法では繊維境界に近づくまで低誤差を維持するが、標準的求積法では特異性付近での誤差低減が失敗する。
• 本手法は急速な収束を示す：パネル数 $ M $ の増加に伴い最大誤差が著しく減少し、特に滑らかな力密度関数に対して顕著である。
• ほぼ特異的積分を扱う技術は、繊維に限りなく近い評価点に対しても効果的に機能し、距離が $ 2.2 \times 10^{-3} $ にまで小さくても誤差は小さく保たれる。
• 事前計算ステップでは、$ n = 16 $ の $ n $ 個の転置 Vandermonde 方程式（$ n \times n $）を解き、$ n^2 $ 個の修正重みを格納する。これにより、すべてのパネルで重みを再利用可能となり、実行時コストは一切かからない。