[论文解读] A Simple Redistribution Vortex Method (with Accurate Body Forces)

O.R. Tutty|arXiv (Cornell University)|Sep 1, 2010

Fluid Dynamics Simulations and Interactions参考文献 13被引用 5

一句话总结

本文提出了一种简单、显式的重新分布涡方法，将涡量从涡量元转移到固定节点，而非在涡量元之间重新分布，从而在避免核心重叠的同时实现高分辨率的黏性流模拟。该方法利用贴体网格实现了精确的点面对力计算，并与基准数据表现出极好的一致性，尤其在雷诺数高达9500时的阻力和涡量分布图方面表现优异。

ABSTRACT

A circulation redistribution scheme for viscous flow is presented. Unlike other redistribution methods, it operates by transferring the circulation to a set of fixed nodes rather than neighbouring vortex elements. A new distribution of vortex elements can then be constructed from the circulation on the nodes. The solution to the redistribution problem can be written explicitly as a set of algebraic formulae, producing a method which is simple and efficient. The scheme works with the circulation contained in the vortex elements only, and does not require overlap of vortex cores. With a body fitted redistribution mesh, smooth and accurate estimates of the pointwise surface forces can be obtained. The ability of the scheme to produce high resolution solution is demonstrated using a series of test problems.

研究动机与目标

• 开发一种简单、高效的黏性流涡方法，避免核心重叠并保持高分辨率。
• 通过使用贴体重新分布网格，实现精确的点面对力计算。
• 提供涡量重新分布的显式代数解，简化实现过程并提高稳定性。
• 通过基准测试案例（包括雷诺数为150和9500的圆柱绕流）展示该方法的精度和分辨率能力。
• 表明该方法在较粗网格下仍能与有限差分法在阻力和涡量方面保持良好一致性，即使存在噪声。

提出的方法

• 每个涡量元的涡量被转移到一组固定、结构化的节点，而非在邻近涡量元之间重新分布。
• 通过代数方程显式求解重新分布过程，确保计算效率和稳定性。
• 从固定节点上的涡量值重构新的涡量元分布，保持涡量和涡量守恒。
• 该方法无需涡量核心重叠，消除了粒子强度交换（PSE）方法的关键限制。
• 使用贴体重新分布网格，以实现精确的点面对力计算。
• 该方法可扩展至三维，并支持局部网格加密而无需全局重网格化。

实验结果

研究问题

• RQ1 能否开发一种简单、显式的重新分布方案，避免涡方法中的核心重叠并保持高分辨率？
• RQ2 在拉格朗日涡方法中，使用贴体重新分布网格能否实现高精度的点面对力计算？
• RQ3 该方法在高雷诺数流动中，对阻力和涡量是否能与基准解保持良好一致性？
• RQ4 通过避免随机游走和PSE方法中常见的噪声，该方法能否在较粗网格下保持精度？
• RQ5 与迭代或基于网格的替代方法相比，涡量重新分布的显式代数解是否更具计算效率和稳定性？

主要发现

• 使用贴体重新分布网格时，该方法能产生精确的点面对力，噪声极小且分辨率高。
• 在 Re = 150 时，采用贴体网格的 DVM 解在涡量分布图和阻力历史方面与 FDS 解高度一致。
• 在 Re = 9500 时，高分辨率 DVM 模拟（3000 个面板，h = 1/1000）在阻力和涡量分布方面与 FDS 解表现出极佳的一致性，仅在极早期瞬态阶段存在微小偏差。
• 即使在较粗网格（1600 个面板，h = 1/400）下，DVM 的阻力历史在 t ≈ 1.5 之前与 FDS 保持良好一致，但后期因分辨率不足而出现偏差。
• 该方法在无核心重叠的情况下实现高分辨率结果，与 PSE 方法不同，且只要黏性扩散长度小于网格间距，就能保持稳定性。
• 涡量重新分布的显式代数解计算效率高，便于实现和扩展至三维。