[论文解读] 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

一句话总结

本论文提出一种具有流体动力学匹配网格的局部校正多blob方法，以提高求解刚性粒子三维斯托克斯迁移率问题的精度。通过优化blob网格构型以匹配理想粒子的自相互作用，并应用受斯托克斯动力学启发的局部成对校正，该方法减少了自相互作用和近场误差，即使在球体和轴对称杆的粗分辨率下也能实现高精度模拟。

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).

研究动机与目标

• 通过优化blob网格构型以匹配理想粒子的流体动力学响应，减少多blob方法在斯托克斯流中主导的自相互作用误差。
• 通过一种低成本、局部的成对校正技术（受斯托克斯动力学启发），提高中等距离和近距离粒子的精度。
• 通过最小化自相互作用带来的基准误差，实现在粗略blob离散化下对刚性粒子的高精度模拟。
• 建立一个适用于异质悬浮液和复杂粒子形状的可控精度多blob方法框架。
• 通过确保校正后的迁移率矩阵保持正定且适用于随机动力学，为未来大规模模拟提供支持。

提出的方法

• 优化blob网格几何构型（如rt网格、r-和t-网格），以匹配理想粒子的平动和转动迁移率，从而最小化自相互作用误差。
• 采用联合求解策略，分别优化平动和转动的网格，然后组合解以减少自相互作用误差。
• 应用基于斯托克斯动力学的局部、低成本成对校正，以校正粒子之间的近场流体动力学相互作用。
• 采用第二类边界积分方程（BIE）公式，结合展开求积法（QBX），作为误差量化的高精度参考。
• 通过少量粒子系统（球体和杆体）的数值实验，评估在不同粒子间距和构型下的误差。
• 通过与参考解对比验证校正方法，量化在远场和近场区域速度精度的提升。

实验结果

研究问题

• RQ1 如何优化blob网格构型，以最小化多blob方法在斯托克斯迁移率问题中的自相互作用误差？
• RQ2 基于斯托克斯动力学的成对校正能在多大程度上减少中等至近距离粒子的速度误差？
• RQ3 与标准多blob网格相比，采用分别针对平动和转动优化的网格的联合求解策略，是否能显著降低自相互作用误差？
• RQ4 自相互作用误差与成对校正有效性的关系如何？
• RQ5 多blob方法的精度如何依赖于粒子形状，特别是具有光滑端帽的轴对称杆体？

主要发现

• 优化blob网格，特别是针对杆体的rt网格，可将自相互作用误差降低至远场模拟中的主导误差下限水平。
• 联合求解策略显著降低了自相互作用误差，尤其在粒子间距较大时，通过分别优化平动和转动的网格实现。
• 成对校正仅在近场相互作用误差超过自相互作用误差下限时才有效；因此，当自相互作用已良好匹配时，其效益最大。
• 在粗略离散化下，对联合求解应用成对校正可带来显著的精度提升，因为较低的自相互作用误差使校正产生可测量的影响。
• 所有数值测试中，校正后的迁移率矩阵均保持正定，经显式特征值计算验证，支持其在随机模拟中的应用。
• 只要网格实现流体动力学匹配且适当地应用成对校正，该方法即使在粗略blob分辨率下也能实现斯托克斯迁移率问题的高精度求解。