Dynamic similarity of vortex shedding in a superfluid flowing past a penetrable obstacle

本文表明，超流体绕过可穿透障碍物时的动态相似性由基于马赫-1 轮廓定义的有效直径所确定的超流体雷诺数支配，而非障碍物的几何尺寸，该发现成功统一了不同障碍物参数下的尾流动力学、涡旋脱落转变及阻力特性。

要点

想象一种超流体，它是一条神奇且无摩擦的河流，其中的水由步调完全一致的原子组成。现在，想象你向这条河中投下一块石头。在普通河流中，水流会绕过石头，在其后方形成一片混乱的漩涡尾迹。而在这条神奇的超流体中，这些“漩涡”是微小的、量化的涡旋，称为涡旋。

长期以来，科学家们知道当“石头”是坚硬且不可穿透的（如巨石）时，这些涡旋如何表现。但如果这块“石头”实际上是一个幽灵般、半透明的屏障，流体可以部分地从中流过，会发生什么呢？这正是本文所解决的谜题。

以下是他们发现的故事，分解为简单的概念：

1. 问题：如何测量“幽灵石头”？

在常规物理学中，如果你想预测水流如何绕过物体，你需要知道其尺寸。如果物体是实心圆柱体，你只需测量其直径。

但在这个实验中，“障碍物”是一束激光。它不是固体墙壁，而是一座温和的能量山丘。超流体原子可以翻越它或流过它。由于流体穿透了障碍物，障碍物的“尺寸”并非固定不变。这就像试图测量云的尺寸；它会根据风吹的强度而变化。

研究人员意识到，仅仅测量激光束的宽度行不通。他们需要一种新的方式来定义对流体而言有意义的障碍物“尺寸”。

2. 解决方案：“速度极限”区域

团队发现，涡旋（即漩涡）仅在流体速度超过该特定位置的局部“声速”时才会出现。

这就像音爆。当喷气式飞机突破音障时，会产生冲击波。在这种超流体中，当流速快得足以突破局部“声障”时，流体会变得不稳定并喷出一个涡旋。

研究人员定义了一个障碍物的新“有效尺寸”。他们测量的不是激光束本身，而是障碍物周围流体速度足以突破声障的不可见区域的大小。

• 类比：想象一座灯塔。你很难测量光束的“大小”。但你可以测量水面上光线强烈到灼伤眼睛的区域大小。对于游过的鱼来说，这个“灼伤区域”才是关键。研究人员利用这个“灼伤区域”（即超音速区域）作为障碍物的真实尺寸。

3. 发现：通用规则手册

一旦他们使用了这个新的“有效尺寸”，神奇的事情发生了。他们能够将所有杂乱的数据整理成单一、清晰的规则手册，就像经典物理学对普通水所做的那样。

他们发现，尾迹的行为取决于一个单一的数字（“超流体雷诺数”）。

• 低数值（慢速流动）：障碍物以整齐、有节奏的行，成对地（一个正涡旋和一个负涡旋）吐出涡旋，如同行进乐队。
• 高数值（快速流动）：节奏被打乱。成对的涡旋变得拥挤，相互碰撞，并重组为同符号涡旋的混乱集群。

该论文表明，无论激光束有多大或多强，这种转变都发生在完全相同的“数值”上。无论障碍物是微小、微弱的幽灵，还是巨大、强烈的实体，一旦你考虑了“速度极限区域”，流体的行为方式都是一样的。

4. 阻力与节奏

研究人员还考察了另外两点：

• 阻力：障碍物使流体减速的程度。他们发现，如果将阻力与新的“超流体数”作图，所有不同尺寸的障碍物都会坍缩到一条平滑的曲线上。
• 节奏（斯特劳哈尔数）：涡旋脱落的频率。同样，当使用他们新的尺寸测量方法时，脱落频率遵循通用模式，就像普通流体（如烟囱后的烟圈）中看到的著名“冯·卡门涡街”一样。

核心结论

该论文声称，尽管超流体是奇怪的量子现象，但如果正确测量障碍物，它们仍然遵循“动力相似性”的古老规则（即小模型可以预测大流动）。

你不应该测量物理激光束。你应该测量流体变得太快而无法保持平静的区域。一旦你这样做了，超流体混乱的量子世界就会表现出与水流绕过岩石时一样的可预测秩序。

简而言之：他们找到了测量幽灵的正确“尺子”，证明了即使量子流体，只要观察流动的正确部分，也遵循与你浴缸中水相同的通用规则。

技术摘要

技术摘要：超流体流经可穿透障碍物时的涡旋脱落动力学相似性

问题陈述
动力学相似性是经典流体动力学的基石，它指出当通过适当的无量纲参数（最 notably 雷诺数 $Re$）表达时，具有不同物理尺度的流动会表现出相同的行为。虽然这种普适性已扩展到流经不可穿透障碍物（即流体密度在边界处完全耗尽）的超流体，但针对可穿透障碍物的定量框架仍然缺失。在可穿透情形下（典型如原子玻色 - 爱因斯坦凝聚体 (BEC) 实验，其中障碍物势低于化学势），流体部分流过障碍物。这从根本上改变了涡旋成核过程，使得无法定义一个明确且动力学相关的特征几何长度尺度（如圆柱直径）。因此，尚不清楚动力学相似性是否适用于这些系统，以及如何定义控制尾流动力学的有效长度尺度。

方法论
作者通过数值求解二维 Gross-Pitaevskii 方程 (GPEs) 来研究 BEC 流经高斯势障碍物的问题。障碍物被建模为可穿透的 ($V_0 < \mu$)，确保不会形成完全耗尽的密度核心。

• 模拟设置：系统采用傅里叶伪谱法结合四阶龙格 - 库塔时间积分器进行模拟。采用边缘区域技术以吸收计算域边界处的密度波动，并采用涡旋解缠算法以防止涡旋在周期性边界条件下重新进入。
• 流动表征：研究分析了超流体涡度 ($\omega_s$) 的空间分布，以及阻力 ($F_D$) 和升力 ($F_L$) 的时间演化。
• 长度尺度的定义：为解决缺乏几何尺度的问题，作者定义了一个有效直径 ($D_{eff}$)。该参数源自障碍物周围的时间平均无旋流分量 ($v_{irr}$)。具体而言，$D_{eff}$ 是无旋流速度超过当地声速 ($|v_{irr}| > c_s$) 区域的横向宽度，由马赫 -1 等值线界定。该区域代表了由可压缩性驱动的失稳导致量子化涡旋成核的区域。
• 无量纲参数：利用 $D_{eff}$，作者构建了一个超流体雷诺数：
  $$Re_s \equiv \frac{(v_0 - v_c)D_{eff}}{\hbar/m}$$
  其中 $v_0$ 为流速，$v_c$ 为临界速度，$m$ 为粒子质量。他们还计算了斯特劳哈尔数 ($St = f_s D_{eff}/v_0$) 和超流体阻力系数 ($C_D$)。

主要结果

1. 偶极子到团簇的转变：随着流速增加，尾流动力学表现出明显的转变。在临界速度附近，系统发射周期性的涡旋偶极子。随着速度增加，转变为同符号涡旋团簇的脱落。这一转变以升力均方根值的突然增加为标志，且阻力标度从 $(v_0 - v_c)^{1/2}$ 变为 $(v_0 - v_c)^2$。
2. 与 $Re_s$ 的普适标度：当将转变速度 ($v_{th}$) 对提出的超流体雷诺数作图时，偶极子到团簇的转变始终发生在 $Re_s \approx 2$ 处，无论障碍物的几何尺寸 ($\sigma$) 或强度 ($V_0$) 如何。
3. 无量纲量的坍缩：
   • 斯特劳哈尔数：当斯特劳哈尔数 ($St$) 对$Re_s$作图时，数据坍缩到一条普适曲线上。对于$Re_s > 2$，$St$饱和于一个值 ($St_0 \approx 0.20\text{--}0.30$)，该值与经典湍流涡旋脱落的数值惊人地接近。
   • 阻力系数：阻力系数 ($C_D$) 对 $Re_s$ 作图时也坍缩到单条曲线上。在低 $Re_s$ 区域，$C_D$ 遵循幂律，指数 $\alpha \approx 0.6\text{--}0.7$（类似于二维斯托克斯定律），而在高 $Re_s$ 区域，指数降低至 $\alpha \approx 0.22\text{--}0.27$。
4. $D_{eff}$ 的有效性：结果表明，有效直径 $D_{eff}$ 与涡度分布的横向宽度 ($w_v$) 呈线性正比，证实了超音速区域（而非障碍物的几何尺寸）设定了涡旋脱落的相关长度尺度。

意义与主张
本文主张将超流体流动中的动力学相似性概念扩展到可穿透障碍物区域。通过基于无旋流的马赫 -1 等值线定义有效直径，作者证明了动力学相关的长度尺度由涡旋成核的局部超音速区域决定，而非由障碍物的几何尺寸决定。

其主要意义在于统一了不同障碍物尺寸和强度下的尾流动力学。结果表明，一旦将流动定义的长度尺度纳入雷诺数，复杂的脱落行为（转变点、斯特劳哈尔数和阻力系数）就表现出与经典流体类似的普适行为。该框架提供了一种量化工具，用于表征从可穿透到不可穿透区域的交叉，并表明马赫 -1 等值线可作为超流体中相关障碍物长度尺度的物理定义。作者指出，当流动变为全局超音速 ($v_0/c_s > 1$) 时，该框架失效，此时激波动力学取代涡旋脱落成为主导。