[论文解读] 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 能否为模拟复杂形状的刚性纤维动力学，开发一种快速、稳定且收敛的数值方法？

主要发现

• 正则化后的第二类积分方程形式确保了数值可逆性与稳定性，即使在精细离散化下，也能将谱系远离零点。
• 数值测试证实，该模型在空间上具有二阶收敛性，且相比第一类形式具有更优的条件数。
• 对于具有周期性的直纤维，积分算子被严格证明为负定，与PDE理论及预期的物理行为一致。
• 弯曲纤维的角动量与取向偏差随弯曲幅度δ线性增加，δ = 0.0003时偏差为1%，δ = 0.0015时偏差达7.5%，持续100个时间单位后。
• 在完成一次旋转后，纤维长轴与x轴之间的夹角θ在δ = 0.0003时偏离约3°，在δ = 0.0015时偏离约8°，证实了可测量的旋转偏差。
• 该方法准确再现了Jeffery的椭球体动力学，验证了模型在典型情况下的准确性，随后成功扩展至复杂弯曲纤维。