[論文レビュー] A Locally Corrected Multiblob Method with Hydrodynamically Matched Grids for the Stokes Mobility Problem

Anna Broms, Mattias Sandberg+1 more|arXiv (Cornell University)|Jul 22, 2022

Lattice Boltzmann Simulation Studies被引用数 4

ひとこと要約

本稿では、3次元ストークス移動度問題における剛体粒子の精度を向上させるために、流体力学的に整合したグリッドを備えたローカル補正付きマルチブロブ法を提示する。ブロブグリッドの配置を最適化し、理想粒子の自己相互作用を再現するとともに、ストークシアンダイナミクスにインspiredされた局所的ペア補正を適用することで、自己相互作用誤差および近接領域誤差を低減し、球体および軸対称ロッドに対しても粗い解像度でも高精度なシミュレーションを可能にする。

ABSTRACT

Inexpensive numerical methods are key to enable simulations of systems of a large number of particles of different shapes in Stokes flow. Several approximate methods have been introduced for this purpose. We study the accuracy of the multiblob method for solving the Stokes mobility problem in free space, where the 3D geometry of a particle surface is discretised with spherical blobs and the pair-wise interaction between blobs is described by the RPY-tensor. The paper aims to investigate and improve on the magnitude of the error in the solution velocities of the Stokes mobility problem using a combination of two different techniques: an optimally chosen grid of blobs and a pair-correction inspired by Stokesian dynamics. Optimisation strategies to determine a grid with a certain number of blobs are presented with the aim of matching the hydrodynamic response of a single accurately described ideal particle, alone in the fluid. Small errors in this self-interaction are essential as they determine the basic error level in a system of well-separated particles. With a good match, reasonable accuracy can be obtained even with coarse blob-resolutions of the particle surfaces. The error in the self-interaction is however sensitive to the exact choice of grid parameters and simply hand-picking a suitable blob geometry can lead to errors several orders of magnitude larger in size. The pair-correction is local and cheap to apply, and reduces on the error for more closely interacting particles. Two different types of geometries are considered: spheres and axisymmetric rods with smooth caps. The error in solutions to mobility problems is quantified for particles of varying inter-particle distances for systems containing a few particles, comparing to an accurate solution based on a second kind BIE-formulation where the quadrature error is controlled by employing quadrature by expansion (QBX).

研究の動機と目的

• ストークス流れにおけるマルチブロブ法の主な自己相互作用誤差を、理想粒子の流体力学的応答に一致するようにブロブグリッド配置を最適化することで低減すること。
• ストークシアンダイナミクスにインスパイアされた低コストで局所的なペア補正技術を用いて、中程度に離れた粒子および近接した粒子の精度を向上させること。
• 自己相互作用に起因する基本誤差を最小限に抑えることで、粗いブロブ離散化を用いた剛体粒子の高精度なシミュレーションを可能にすること。
• 不均一な懸濁系や複雑な粒子形状に適用可能な、制御可能な精度を実現するフレームワークを確立すること。
• 将来の大規模シミュレーションを支援するため、補正済みの移動度行列が正定値であり、確率的ダイナミクスに適していることを保証すること。

提案手法

• 理想粒子の並進および回転移動度に一致するように、ブロブグリッドの幾何配置（例：rtグリッド、r-およびtグリッド）を最適化し、自己相互作用誤差を最小化すること。
• 並進と回転のためのグリッドを別々に最適化し、その解を統合することで自己相互作用誤差を低減する、統合的ソルブ戦略を用いること。
• ストークシアンダイナミクスにインスパイアされた局所的かつ低コストのペア補正を用い、粒子間の近接領域流体力学的相互作用を補正すること。
• 誤差評価のための高精度基準として、四次元展開（QBX）を用いた第二種境界積分方程式（BIE）定式化を採用すること。
• 球体およびロッドを用いた少数粒子系の数値実験を通じて、粒子間距離や配置の変化に伴う誤差を評価すること。
• 補正済み手法を基準解と比較して、遠方および近接領域の両方において速度精度の向上を定量すること。

実験結果

リサーチクエスチョン

• RQ1 ストークス移動度のマルチブロブ法において、自己相互作用誤差を最小化するためのブロブグリッド配置をどのように最適化できるか。
• RQ2 ストークシアンダイナミクスにインスパイアされたペア補正は、中程度から近接した距離に置かれた粒子の速度誤差をどの程度低減できるか。
• RQ3 並進と回転のための別々のグリッドを用いる統合的ソルブ戦略は、標準的なマルチブロブグリッドと比較して、自己相互作用誤差を顕著に低減できるか。
• RQ4 自己相互作用誤差とペア補正の有効性との間にはどのような関係があるか。
• RQ5 特に滑らかなキャップを持つ軸対称ロッドに対して、粒子形状がマルチブロブ法の精度にどのように影響するか。

主な発見

• 特にロッドに対してrtグリッドを最適化することで、自己相互作用誤差を、遠方シミュレーションにおける主な誤差床レベルまで低減できる。
• 統合的ソルブ戦略により、自己相互作用誤差が顕著に低減され、特に粒子間隔が大きい場合に顕著である。これは並進と回転のためのグリッドを別々に最適化することで達成される。
• ペア補正は、近接領域相互作用誤差が自己相互作用誤差床を上回っている場合にのみ有効である。つまり、自己相互作用がすでに良好に一致している場合に最も効果的である。
• 粗い離散化において、統合的ソルブにペア補正を適用することで、顕著な精度向上が得られる。これは、低い自己相互作用誤差により補正の影響が明確に現れるからである。
• すべての数値実験において、補正済みの移動度行列が正定値であることが、明示的な固有値計算により確認された。これは、確率的ダイナミクスへの応用が可能であることを支持する。
• グリッドが流体力学的に整合されており、適切にペア補正が適用されていれば、粗いブロブ解像度でもストークス移動度問題の高精度な解が得られる。