[论文解读] Simple and efficient representations for the fundamental solutions of Stokes flow in a half-space

Zydrunas Gimbutas, Leslie Greengard+1 more|arXiv (Cornell University)|Apr 1, 2015

Lattice Boltzmann Simulation Studies参考文献 30被引用 39

一句话总结

本文通过结合一个镜像 Stokeslet 与基于标量调和函数的单个 Papkovich-Neuber 位势，提出了一种简化且高效的半空间 Stokes 流基本解表示方法。该方法仅使用一个镜像力和一个调和校正项，即可精确满足无滑移边界条件，相较于 Blake 的复杂公式，提供了更直观且计算上更易处理的替代方案。

ABSTRACT

We derive new formulas for the fundamental solutions of slow, viscous flow, governed by the Stokes equations, in a half-space. They are simpler than the classical representations obtained by Blake and collaborators, and can be efficiently implemented using existing fast solvers libraries. We show, for example, that the velocity field induced by a Stokeslet can be annihilated on the boundary (to establish a zero slip condition) using a single reflected Stokeslet combined with a single Papkovich-Neuber potential that involves only a scalar harmonic function. The new representation has a physically intuitive interpretation.

研究动机与目标

• 与经典的 Blake 型公式相比，开发一种更简单、更高效的半空间 Stokes 流基本解表示方法。
• 解决现有镜像方法中涉及多个偶极子和四极子项所带来的计算复杂性问题。
• 提供一种物理解释直观的框架，仅通过最少数量的组件即可实现无滑移边界条件的强制满足。
• 利用现有的自由空间基本解和调和位势的快速求解器，实现高效计算。
• 将该方法扩展至应力偶极子、旋转偶极子和 Stokes 双偶极子，并给出统一且紧凑的公式表达。

提出的方法

• 利用 Papkovich-Neuber 表示法构造速度场与压力场，该方法通过 u(x) = x₃∇φ(x) − [0,0,φ(x)]ᵀ 和 p(x) = 2∂φ/∂x₃ 将一个调和函数 φ(x) 映射为 Stokes 解。
• 在镜像点 yI = (y₁,y₂,−y₃) 处使用一个反射 Stokeslet，以抵消边界 x₃ = 0 上的切向速度分量。
• 应用一个基于调和位势 φ(x) 的单个 Papkovich-Neuber 校正项，该位势由位于 yI 处的源与偶极子共同导出，用于抵消法向速度分量。
• 对于应力偶极子与双偶极子，校正位势 φ 由位于 yI 处的偶极子与四极子调和源构成，其系数由力矢量与取向矢量导出。
• 利用现有的自由空间格林函数（Stokeslet、偶极子、四极子）与调和位势计算的快速算法，使该方法可无缝集成至基于 FMM 的求解器中。
• 推导了二维 Stokes 流与线弹性力学的类似公式，表明该方法具有广泛的应用潜力。

实验结果

研究问题

• RQ1 能否为半空间中的 Stokeslet 推导出一种更简单的镜像表示方法，以避免 Blake 解中复杂的多偶极子结构？
• RQ2 能否仅通过一个调和位势校正项而非多个高阶奇异性，来实现无滑移边界条件的强制满足？
• RQ3 Papkovich-Neuber 表示法是否能为半空间中所有基本解（Stokeslet、应力偶极子、旋转偶极子、双偶极子）提供物理解释直观且计算高效的公式？
• RQ4 新公式能否利用现有的自由空间基本解快速求解器实现高效计算？
• RQ5 新公式在复杂度与计算成本方面与 Blake 的原始公式相比如何？

主要发现

• 所提出的方法将 Blake 复杂的镜像结构（需多个偶极子与四极子项）替换为基于标量调和函数的单个 Papkovich-Neuber 位势。
• 对于 Stokeslet，校正位势 φ(x) 是位于镜像点 yI 处的源与偶极子的线性组合，使得在 x₃ = 0 处有 u₃ = −φ(x)，从而抵消法向速度。
• 该方法仅使用一个反射 Stokeslet 与一个 Papkovich-Neuber 校正项，即可精确满足无滑移条件（x₃ = 0 处 u = 0）。
• 应力偶极子的校正使用调和位势 φT(x)，其由位于 yI 处的偶极子与四极子构成，系数取决于 νI·gI 与 y₃。
• 旋转偶极子的校正使用两个偶极子，其系数分别为 −2ν₃^I 与 2g₃^I，而 Stokes 双偶极子的校正将对称与反对称部分合并为一个 φD(x)，包含三项。
• 该公式可自然推广至二维 Stokes 流与线弹性力学，其中使用二维调和位势的类似 Papkovich-Neuber 校正项。