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

研究动机与目标

• 开发一种稳健且高效的数值方法，用于评估斯托克斯流中细长纤维非局部细长体公式中的有限部分积分。
• 解决当纤维彼此靠近或求值点靠近纤维中心线时，积分具有奇异或近似奇异特性时的精确计算挑战。
• 提供一种通过解析处理奇异性来实现高精度的积分技术，且计算成本仅在预计算阶段略有增加，运行时无额外开销。
• 将该方法扩展至流速计算中出现的近似奇异积分，特别是在纤维表面附近的情形。
• 在螺旋线和圆形环等基准几何结构上验证该方法，结果表明其收敛性与误差控制优于标准积分方法。

提出的方法

• 通过将被积函数分解为光滑多项式因子与奇异核函数，对有限部分积分进行重构，从而实现在每个面板上的解析积分。
• 采用光滑因子的多项式逼近，其系数通过一次求解转置范德蒙德线性系统（在预计算阶段完成）获得，以生成修正积分权重。
• 将生成的修正积分权重在所有面板间重复使用，从而实现高效且精确的奇异积分计算，且每次求值无需额外计算成本。
• 对于近似奇异的斯托克斯勒特积分，采用类似的乘积积分方法，利用根求解得到的共轭复数对来分解被积函数。
• 对所得积分的解析求值依赖于递推公式，并通过细致的数值实现以保持精度。
• 该方法采用高斯-勒让德积分规则实现，并在具有高精度参考解的三维纤维几何结构上进行了验证。

实验结果

研究问题

• RQ1 当被积函数在对角线上奇异时，如何在细长体理论中高精度地评估有限部分积分？
• RQ2 能否设计一种积分方法，通过解析处理奇异性来保持效率，且不增加额外计算成本？
• RQ3 当求值点靠近纤维中心线时，该方法在近似奇异积分上的表现如何？
• RQ4 该方法在不同纤维几何形状和密度函数下的收敛行为如何？
• RQ5 该方法能否应用于三维斯托克斯流中具有自由端的多个相互作用纤维？

主要发现

• 通过预先计算的修正积分权重对奇异性进行解析处理，所提方法在有限部分积分上实现了谱精度。
• 与相同面板数的高斯-勒让德积分相比，该方法在螺旋线纤维附近计算流速时，最大误差降低了多个数量级。
• 当 $ M = 8 $ 个面板时，该特殊积分方法在整个接近纤维边界区域均保持低误差，而标准积分方法在奇异点附近无法进一步减小误差。
• 该方法表现出快速收敛：随着面板数 $ M $ 增加，最大误差迅速衰减，尤其在光滑力密度分布下表现显著。
• 针对近似奇异积分的技术可成功处理距离纤维任意接近的求值点，即使在距离小至 $ 2.2 \times 10^{-3} $ 的测试案例中，误差仍保持较小。
• 预计算步骤包括求解 $ n $ 个大小为 $ n \times n $ 的转置范德蒙德线性系统（$ n = 16 $），并存储 $ n^2 $ 个修正权重，使得所有面板均可无运行时成本地复用。