[论文解读] Time dependent simulation of the Driven Lid Cavity at High Reynolds Number

Nuno Vieira Cardoso, Pedro Bicudo|ArXiv.org|Sep 18, 2008

Computational Fluid Dynamics and Aerodynamics参考文献 11被引用 3

一句话总结

本研究采用1024×1024的网格和SSPRK(5,4)格式，在涡量-流函数形式下对雷诺数高达30,000的二维驱动方腔流进行了时间依赖性模拟。通过双精度CUDA加速计算，结果收敛至稳态解，解决了以往争议：此前研究中报告的非定常解实为粗网格导致的数值伪影，而非物理行为。

ABSTRACT

In this work, numerical solutions of the two dimensional time dependent incompressible flow, in a driven cavity at high Reynolds number Re, are presented. At high Re, there is a controversy. Some studies predicted that the flow is steady, others found time dependent non-steady flow, either periodic or aperiodic. In this study, the driven lid cavity is successfully solved using a very fine grid mesh, for Re up to $30 000$. We discretize the Vorticity-Stream formulation of the Navier-Stokes equation with the SSPRK(5,4) scheme in a $1024 imes1024$ grid. Using this very fine grid, the results obtained converge to a stationary solution. Detailed results for Re between $5 000$ and $30 000$ are presented. The driven lid cavity problem is solved with a NVIDIA GPU using the CUDA programming environment with double precision.

研究动机与目标

• 解决长期存在的争议：即驱动方腔流在高雷诺数下是否变得非定常。
• 确定在极细网格下进行的时间依赖性模拟是否收敛至稳态解，而非振荡或周期性行为。
• 研究文献中先前报告的非定常解是否由于空间分辨率不足而成为数值伪影。
• 验证足够细的网格是捕捉真实物理解的必要条件，可避免粗网格引起的虚假振荡。
• 展示使用GPU加速双精度计算的高分辨率直接数值模拟的有效性。

提出的方法

• 采用涡量-流函数形式求解二维不可压Navier-Stokes方程，自动满足连续性方程。
• 使用SSPRK(5,4)龙格-库塔格式对时间依赖的涡量输运方程进行时间离散，以保证精度与稳定性。
• 采用1024×1024结构化笛卡尔网格，确保高空间分辨率，尤其在小涡旋形成的角落区域。
• 在NVIDIA GPU上使用CUDA实现，采用双精度算术以最小化截断误差并提高解的保真度。
• 应用一致的边界条件：除顶部滑动壁面以恒定速度驱动流场外，其余壁面均为无滑移边界。
• 时间推进直至解达到静止状态，通过残差范数和涡量场演化进行监控。

实验结果

研究问题

• RQ1 在高雷诺数（高达30,000）下，驱动方腔流是否表现出时间依赖性行为，还是收敛至稳态？
• RQ2 以往直接数值模拟（DNS）中报告的非定常解是物理真实现象，还是由于网格粗化导致的数值伪影？
• RQ3 极细网格（1024×1024）是否能够解析所有相关涡旋结构（包括小角落涡），并稳定解？
• RQ4 与低精度方法相比，使用GPU上的高精度双精度算术是否能提高模拟的精度与收敛性？
• RQ5 在稳态Navier-Stokes方程中存在稳态解是否足以断定其具有物理现实性，还是需要时间依赖性模拟来验证？

主要发现

• 在雷诺数5,000至30,000范围内，时间依赖性模拟在1024×1024网格上收敛至稳态解，表明真实物理解为稳态。
• 在Re = 5,000时，仅能清晰识别主涡和部分角落涡（BL1、BR1、TL1、TR1）；BL3和TL2涡旋因分辨率不足而缺失。
• 在Re = 15,000时，所有主要涡旋（包括BL3和TL2）均被完全解析且解稳定。
• 在Re = 20,000、25,000和30,000时，解保持静止，证实高雷诺数下收敛至稳态。
• 本研究将先前报告的非定常或周期性流动归因于粗网格引起的数值振荡，而非物理不稳定性。
• 结果支持Erturk的观点：极细网格对于解析小尺度涡旋并获得正确稳态解至关重要。