[论文解读] The vortex method for 2D ideal flows in exterior domains

Diogo Arsénio, Emmanuel Dormy+1 more|arXiv (Cornell University)|Jul 5, 2017

Navier-Stokes equation solutions参考文献 24被引用 4

一句话总结

本文严格证明了在具有光滑边界的外部区域中，使用点涡模拟二维理想流体流动的涡方法的合理性，其中通过沿边界布置点涡来近似涡度，以确保流体切向流动和涡量守恒。本文建立了涡近似的收敛性，并提出了流体电荷方法，该方法在数值条件和模拟稳定性方面表现更优。

ABSTRACT

The vortex method is a common numerical and theoretical approach used to implement the motion of an ideal flow, in which the vorticity is approximated by a sum of point vortices, so that the Euler equations read as a system of ordinary differential equations. Such a method is well justified in the full plane, thanks to the explicit representation formulas of Biot and Savart. In an exterior domain, we also replace the impermeable boundary by a collection of point vortices generating the circulation around the obstacle. The density of these point vortices is chosen in order that the flow remains tangent at midpoints between adjacent vortices and that the total vorticity around the obstacle is conserved. In this work, we provide a rigorous justification of this method for any smooth exterior domain, one of the main mathematical difficulties being that the Biot--Savart kernel defines a singular integral operator when restricted to a curve (here, the boundary of the domain). We also introduce an alternative method---the fluid charge method---which, as we argue, is better conditioned and therefore leads to significant numerical improvements.

研究动机与目标

• 为在具有光滑边界的外部区域中，通过边界离散点涡表示来严格数学证明二维理想流的涡方法。
• 解决由限制在边界曲线上Biot–Savart核所引起的奇异积分算子的挑战。
• 建立涡近似对欧拉方程真实解的动态与静态收敛性。
• 提出并分析一种替代方法——流体电荷方法，以改善数值条件性和稳定性。
• 证明涡近似与流体电荷近似的收敛速率，特别是在Riemann和离散化的情境下。

提出的方法

• 通过在边界曲线上布置有限个点涡来近似外部区域中的涡度，以确保流体切向流动并保持总涡量守恒。
• 利用Biot–Savart定律从涡度分布重构速度场，将边界视为具有特定密度的涡片，以满足无通量条件。
• 应用Plemelj和Poincaré–Bertrand公式分析由速度边界积分表示所引出的奇异积分算子。
• 使用均匀或拟均匀网格对边界积分方程进行离散化，并证明离散算子的有界性与谱收敛性。
• 提出流体电荷方法作为替代方案，其中边界由流体电荷（双层势）表示，而非点涡，从而获得更良态的线性系统。
• 利用谱分析和Riemann和的收敛估计，量化边界积分离散化中的误差，尤其针对边界上光滑周期函数。

实验结果

研究问题

• RQ1 尽管Biot–Savart核在边界上具有奇异性，是否可以严格证明二维理想流在光滑外部区域中使用涡方法的合理性？
• RQ2 在这些区域中，涡近似对欧拉方程真实解的收敛速率是多少？
• RQ3 在数值条件性和稳定性方面，流体电荷方法与经典涡方法相比有何差异？
• RQ4 离散化边界积分算子的谱性质是什么？它们如何影响数值格式的收敛性？
• RQ5 在边界上光滑周期函数的假设下，能否对离散化中使用的Riemann和的收敛性进行量化，并获得最优收敛速率？

主要发现

• 涡方法在外部区域中对二维欧拉方程真实解的收敛是均匀的，误差界为 $ \mathcal{O}(N^{-\beta}) $，其中 $ \beta > 0 $，具体取决于涡度的正则性与网格质量。
• 与涡方法相比，流体电荷方法生成的线性系统具有更好的条件性，显著提升了模拟中的数值稳定性和性能。
• 在被积函数光滑的假设下，边界积分离散化中使用的Riemann和的收敛速率被量化为 $ \mathcal{O}(N^{-\beta}) $，其中 $ \beta \geqslant 2 $。
• 由涡方法产生的离散系统满足离散版本的Poincaré–Bertrand恒等式，这对证明一致性和稳定性至关重要。
• 在适当的网格条件下，与涡方法相关的离散算子的谱保持远离零点，从而确保线性系统的适定性。
• 在紧致时间区间上，流体电荷近似的动态收敛性已建立，收敛速率为 $ \mathcal{O}(N^{-\beta}) $，在 $ L^∞ $ 范数下，其中 $ \beta > 0 $。