[论文解读] Theoretical justification and error analysis for slender body theory

Yoichiro Mori, Laurel Ohm+1 more|arXiv (Cornell University)|Jun 30, 2018

Micro and Nano Robotics参考文献 40被引用 5

一句话总结

本文通过在细纤维上建立具有混合边界条件的新PDE问题，为斯托克斯流中的细长体理论（SBT）建立了严格的理论基础，证明了适定性，并推导出误差估计，表明SBT近似相对于真实解的收敛速率与纤维半径成正比（对数校正项内），在力密度和中心线几何的适度正则性假设下成立。

ABSTRACT

Slender body theory facilitates computational simulations of thin fibers immersed in a viscous fluid by approximating each fiber using only the geometry of the fiber centerline curve and the line force density along it. However, it has been unclear how well slender body theory actually approximates Stokes flow about a thin but truly three-dimensional fiber, in part due to the fact that simply prescribing data along a one-dimensional curve does not result in a well-posed boundary value problem for the Stokes equations in $\mathbb{R}^3$. Here, we introduce a PDE problem to which slender body theory (SBT) provides an approximation, thereby placing SBT on firm theoretical footing. The slender body PDE is a new type of boundary value problem for Stokes flow where partial Dirichlet and partial Neumann conditions are specified everywhere along the fiber surface. Given only a 1D force density along a closed fiber, we show that the flow field exterior to the thin fiber is uniquely determined by imposing a fiber integrity condition: the surface velocity field on the fiber must be constant along cross sections orthogonal to the fiber centerline. Furthermore, a careful estimation of the residual, together with stability estimates provided by the PDE well-posedness framework, allow us to establish error estimates between the slender body approximation and the exact solution to the above problem. The error is bounded by an expression proportional to the fiber radius (up to logarithmic corrections) under mild regularity assumptions on the 1D force density and fiber centerline geometry.

研究动机与目标

• 通过定义一个SBT所逼近的适定PDE问题，解决细长体理论中长期存在的理论模糊性。
• 解决仅在中心线上指定1D力密度时，无法在三维中定义适定斯托克斯问题的问题。
• 仅使用1D数据，为细纤维问题提供一个数学上严谨的“真实”解的定义。
• 推导SBT近似与所提出PDE公式的确切解之间的定量误差估计。
• 为计算模拟中在黏性流体中广泛使用的SBT提供严格的收敛保证。

提出的方法

• 提出一种新的斯托克斯流PDE公式，其中在纤维表面施加部分狄利克雷（无滑移）和部分诺伊曼（牵引力）边界条件。
• 施加纤维完整性条件：垂直于中心线的截面上表面速度必须恒定，以确保解的唯一性。
• 使用PDE适定性框架，建立细长体PDE解的存在性与唯一性，且常数与纤维半径ε无关。
• 通过将解分解为调和部分和一个满足边界条件的校正项，构造一个校正后的速度场。
• 利用迹不等式、科恩不等式和索伯列夫嵌入，以纤维半径ε表示SBT近似的残差。
• 通过构造一个对比场ũ并估计其与精确解的偏差，推导误差界，利用PDE框架的稳定性与残差估计。

实验结果

研究问题

• RQ1 给定仅1D力数据无法定义适定的3D斯托克斯问题，细长体理论实际上近似的是哪一个PDE问题？
• RQ2 如何仅使用1D力密度和中心线几何，为细纤维问题定义唯一解？
• RQ3 细长体近似相对于3D斯托克斯问题真实解的收敛速率是多少？
• RQ4 PDE框架中的适定性常数如何依赖于纤维半径ε？当ε → 0时是否保持稳定？
• RQ5 SBT近似与精确解之间的定量误差是多少？其随纤维半径如何变化？

主要发现

• 具有混合边界条件和纤维完整性条件的细长体PDE，对任意给定的1D力密度和中心线几何，均存在唯一解。
• PDE的适定性已建立，且常数与纤维半径ε无关，确保了当ε → 0时的稳定性。
• 细长体近似的残差被一个与ε成正比的项所控制（对数校正项内）。
• SBT近似与PDE精确解之间的误差被有界为Cε|f|_H^1，其中C与ε无关，f为力密度。
• 误差估计通过构造的对比场ũ推导得出，依赖于PDE框架的稳定性估计与残差分析。
• 在力密度和中心线的正则性条件假设下，收敛速率是最优的，证实了SBT的渐近准确性。