Fully Turbulent Wakes at Low Reynolds Numbers: the Case of the Thin Flat Plate

本文通过直接数值模拟与实验对比表明，薄二维平板背后的尾迹流在相对较低的雷诺数400时即进入完全湍流状态，并表现出与高雷诺数湍流尾迹无法区分的统计与谱特性，这一转变路径与典型圆形或方形圆柱的转变路径显著不同。

要点

想象一下，你正拿着一块薄而平的硬纸板（就像一张扑克牌）在强风中。当风吹向这张卡片时，它会在身后形成一条混乱、翻腾的气流轨迹，称为“尾迹”。长期以来，科学家们一直认为，要让这种尾迹真正变得混乱和“湍流化”，风速必须非常快，或者物体必须是像圆柱或方块这样的特定形状。

这篇论文讲述了一个不同的故事。研究人员发现，如果使用一个薄而平的平板，其背后的空气会在比预期低得多的风速下变得完全混乱且具有湍流性。事实上，在其他形状仍处于相对平静和有序的状态时，平板背后的空气已经进入了湍流状态。

以下是他们研究结果的简单类比分解：

1. “湍流阈值”带来的惊喜

把湍流想象成一个拥挤的舞池。

• 旧有的观念（圆柱体）： 如果风吹过一根圆柱形杆子，杆子背后的空气最初是平静、有节奏的舞蹈（左右摇摆）。需要大量的能量（高风速）才能让舞者们开始互相碰撞、疯狂旋转并创造出一片混乱的景象（湍流）。这种转变在一个很宽的速度范围内是缓慢发生的。
• 新的发现（薄平板）： 研究人员发现，对于薄平板而言，舞池从平静到狂欢的过程几乎是瞬间完成的。即使在相对较低的风速下（雷诺数 400），平板背后的空气已经完全陷入了混乱。它不会经历圆柱体那种缓慢的、有节奏的阶段，而是直接跳到了“派对”模式。

2. 他们是如何证明的

为了确保自己不是在凭空想象，团队像侦探对比犯罪现场一样进行了对比。

• 模拟实验（虚拟实验室）： 他们利用超级计算机模拟了低速风撞击平板的情景（Re 150 和 Re 400）。
• 现实世界测试（风洞实验）： 他们还观察了风速要快得多的真实实验情况（Re 12,500 和 Re 19,700）。
• 完美匹配： 当他们将低速计算机模拟（Re 400）与高速现实实验进行对比时，发现两者的模式完美契合。湍流的“指纹”——即空气如何移动、拥有多少能量以及如何旋转——是完全一致的。
• 对照组： 当他们观察更低速度下的模拟实验（Re 150）时，发现模式完全不同。此时仍处于“平静”阶段，尚未进入混沌状态。这证明了从有序到混沌的转变发生在 150 到 400 之间。

3. 湍流的“指纹”

如何判断一种流动是否真正具有湍流性？论文通过寻找数据中的特定“生命迹象”来进行判断：

• 能量谱（风的声音）： 在平静的流动中，能量集中在少数几个特定的音符上（就像长笛吹奏单一的音调）。而在湍流中，听起来则像是白噪音或静电噪声，能量分布在极广的频率范围内。研究人员发现，在 Re 400 时，平板背后的风“声音”已经充满了这种混乱的静电噪声，就像在高速实验中观察到的一样。
• “间歇性”（偶尔的尖叫）： 在真正的湍流中，空气不仅仅是轻轻地旋转，还会出现突然且剧烈的速度和旋转爆发。研究人员在 Re 400 的数据中发现了这些“尖叫”，但在 Re 150 时却并未发现。

4. 为什么会有所不同？

论文指出，这种突然跳跃的原因在于物体的形状。

• 圆形/方形物体： 当风撞击圆形或方形物体时，物体的后部起到了类似“屏蔽罩”的作用，稳定了物体背后的气流。打破这种稳定性需要大量的能量。
• 薄平板： 由于平板非常薄，没有任何“后部”来屏蔽空气。压力波动（空气的推力和拉力）从一开始就直接与旋转的涡流相连。这就像尝试用笔尖平衡铅笔与用保龄球平衡物体一样；铅笔（薄平板）本质上是不稳定的，会更快地倾斜并陷入混沌。

核心结论

这篇论文改变了我们对空气绕过平面物体流动方式的理解。它证明了薄平板在极低的速度下就能产生完全湍流的尾迹，其速度远低于圆形或方形物体。对于这些形状而言，从有序到混沌的转变不是一个缓慢、渐进的过程，而是一个发生在极早期速度范围内的突然且根本性的转变。

研究人员并未讨论这如何应用于建造桥梁、设计汽车或医疗器械。他们严格专注于证明这一现象确实存在，以及空气流动的物理机制如何与我们之前的认知有所不同。

技术摘要

技术摘要：低雷诺数下的完全湍流尾迹：以薄平板为例

问题陈述
针对垂直于均匀流场的薄平板绕流进行准二维（2D）数值模拟，在历史上一直具有挑战性。计算流体力学（CFD）研究在低雷诺数（$Re < 1000$）下的尾迹动力学方面得出了相互矛盾的结果。虽然实验数据表明，在三维（3D）机制下，积分参数（如阻力系数和斯特劳哈尔数）对$Re$的依赖性较弱，但模拟研究却报告了截然不同的尾迹行为。一个核心问题仍然存在：薄平板后的尾迹是在三维不稳定性刚出现后就立即进入完全湍流状态（类似于典型的圆柱尾迹），还是遵循一种独特的、多阶段的转变路径？作者研究了在$Re = 400$这一相对较低的雷诺数下，薄平板后的尾迹是否已经处于完全湍流状态；在该机制下，典型的圆形和方形柱尾迹通常表现出强烈的$Re$ 依赖性和非湍流波动特征。

研究方法
本研究采用了一种结合直接数值模拟（DNS）与实验验证的多维度方法：

1. 数值模拟： 为了确保稳健性，使用了两种独立的 DNS 代码：

   • 谱/hp 元方法 (SEM)： 在 Nektar++ 中实现，用于进行 $Re = 400$ (DNS400) 和 $Re = 150$ (DNS150) 的主要高精度 DNS。
   • 有限体积法 (FVM)： 在 OpenFOAM 中实现，用于网格收敛性研究和交叉验证 (DNS400V)。
   • 计算域经过截断，并在平板上施加无滑移边界条件，设置自由流入口条件，以及由 Dong 等人提出的特定出口条件。在展向方向实施了周期性边界条件。
   • 进行了严格的网格收敛性测试，确保分辨率达到或优于现有文献中针对类似几何形状所要求的 Kolmogorov 长度尺度 ($\eta$) 的 $\Delta x/\eta$ 值。

2. 实验数据： 使用了来自加州大学风洞的随时间解析立体粒子图像测速技术（PIV）获取的两组实验数据集：

   • EX12K： $Re = 12,500$
   • EX20K： $Re = 19,700$
     这两个数据集代表了完全湍流机制，并作为低 $Re$ 模拟的基准。

3. 分析： 研究对比了不同 $Re$ 情况下的全局积分量、平均速度场、雷诺应力、湍流动能（TKE）收支、能量谱以及速度梯度（特别是非高斯性）的统计特性。

主要贡献与结果

• 湍流的早期发生： 主要发现是，薄平板的尾迹在 $Re = 400$时即进入完全湍流状态。这显著低于典型圆形（$Re \approx 10^5$）或方形（$Re \approx 6000$）柱尾迹的转变过程。
• **与高 $Re$实验的定量一致性：**$Re = 400$ 的 DNS 结果 (DNS400) 在以下方面与 $Re = 12,500$和$19,700$ 的实验数据 (EX12K, EX20K) 表现出极高的吻合度：
  • 全局量： 平均阻力系数 ($C_D$) 和斯特劳哈尔数 ($St$) 是$Re \ge 400$时关于$Re$ 的弱变化函数。
  • 流动结构： 平均流线、速度剖面以及雷诺应力分布（$u'^2, v'^2, w'^2$）在测量误差范围内，与高 $Re$ 实验在统计上是不可区分的。
  • 能量收支： $k$-输运方程中 TKE 产生项和平流项的空间分布在 $Re=400$与高$Re$ 情况之间保持一致。
• **低 $Re$下的湍流特征：**$Re=400$的情况展现了$Re=150$ 情况中所缺失的确定性湍流特征：
  • 能量谱： 脉动速度分量显示出宽频谱，其惯性子区遵循近似为 $f^{-5/3}$ 的标度（Kolmogorov 定律），表明存在能量级联。相比之下，$Re=150$ 的能谱仅在少数几个模态之间进行能量传递。
  • 非高斯性： 脉动应变率 ($\Sigma$) 和旋转率 ($\Omega$) 的概率密度函数 (PDF) 呈现重尾特征，标志着湍流特有的间歇性和罕见的高强度事件。而 $Re=150$ 的情况则表现为高斯分布。
• 独特的转变路径： $Re=150$的转变涉及三维不稳定性，但缺乏$Re=400$中观察到的相干、完全湍流结构。$Re=150$的尾迹特征为较弱的背流、较慢的尾迹恢复以及非单调演化的尾迹半宽，这与典型圆柱的行为有本质区别——在圆柱中，三维不稳定性（模态 A, B, C）会在较宽的$Re$ 范围内演化，随后才达到完全湍流。

意义与结论
本文认为，薄平板向湍流转变的路径与典型圆柱有着本质的不同。虽然圆柱尾迹在广泛的雷诺数范围内经历了一系列长时间的不稳定性序列（模态 A, B, C）后才达到完全湍流状态，但薄平板的尾迹在三维不稳定性出现后（在 $Re=150$到$Re=400$ 之间）几乎立即转变为完全湍流状态。

作者认为，这种差异可能归因于该几何形状缺乏流向后体（afterbody）。对于圆柱而言，后体稳定了流动，并将压力波动与分离点解耦。相比之下，对于薄平板，压力波动直接与涡旋形成区耦合，从而可能驱动更早且更快速的向湍流的转变。这些发现有助于调和先前文献中关于薄平板尾迹对长径比和雷诺数敏感性的争议，确立了对于该类几何形状而言，“完全湍流”机制的起始 $Re$ 比此前假设的典型钝体要低得多。