“没有黏性”究竟少掉了哪条路
搅动一杯普通液体，停止搅拌后，流动会逐渐衰减。大团运动可以把能量交给更小的运动结构，最后由黏性把机械运动能转成内能。超流体令人困惑的地方是：既然超流组分没有这种普通剪切黏性，为什么低温实验里的湍动仍然会减弱？
答案的第一步，是把“没有某一种损失项”与“所有能量出口都关闭”分开。在超流氦中，运动可以改变自身的组织形式：大尺度流动、量子涡旋、涡线上传播的波、声音以及其他准粒子并不是同一种能量库存。有限温度下还存在正常组分，容器壁也可能交换能量和角动量。1
因此，本篇要解释的不是能量如何凭空消失，而是我们观测的流动能如何减少，减少的那一部分去了哪里。只看流速变小，不能立即知道能量是否已经变成热；它也可能暂时进入了更细的结构，或离开了观测区域。
这里还有一个容易混淆的对象差别：液态氦-4在约2.17 K以下进入超流态；氦-3的超流转变在毫开尔文量级，而且其凝聚对象是成对的费米原子。两者都有量子涡旋，但涡核结构和低能激发不同。“氦的实验”不能不带同位素、温度和流动形态就合并成一条证据。1
两个组分可以互相传能，摩擦的总损失却不一定大
有限温度的两流体描述，把同一份液体的响应写成超流组分与正常组分。它不是把液体想成两种互不混合的化学物质。正常组分代表热激发的集体响应；当量子涡旋存在时，热激发与涡旋的相互作用，使两个组分之间产生相互摩擦。13
先看一个足以把能量账讲清的小模型。把某个流动模式的两个速度振幅记为$v_s,v_n$，正常组分质量份额为$r_n$，超流份额为$r_s=1-r_n$。令$g$为两者间的耦合率，$d$为正常组分的阻尼率：
这是本库自选参数的单模式模型，保留相反方向的内部力与正常组分的损失项，不生成三维湍流，也不模拟涡核。相应的单位总质量流动能是：
这行式子回答了两个问题。第一，超流速度方程里没有直接黏性项，超流运动仍可以经耦合被拖慢。第二，两个组分之间的能量交换不等于整体能量损失：分别记账时，一个组分可能大量失去能量，另一个得到其中大部分；把两边相加，只剩由相对运动决定的那部分摩擦损失，以及正常组分的阻尼损失。
本库令$r_n=0.1,g=10,d=1$，两个速度初始都为1。积分到无量纲时间100时，初始0.5单位流动能中，约8.18%经相互摩擦变成热，约91.82%经正常组分阻尼变成热，余下约一亿分之一。把耦合$g$改为0，正常组分仍会减速，超流组分却保持初始速度，最终留下90%的流动能。15
如果连正常组分的阻尼也去掉，相互摩擦只能消除两者的相对运动，不能把共同运动一并消除。例如初始$v_s=1,v_n=0$，仍取$r_n=0.1$，最后两个速度都趋于0.9。流动能从0.45降至0.405，剩下的共同运动继续存在。这说明“有相互摩擦”也不等于“所有运动必然完全停下”。
公开模拟能把这笔账算到哪里
2025年Polanco等研究了约1.44—2.16 K的氦-4湍流。他们使用HVBK粗粒化两流体方程，并把各尺度的能量传递、组分间交换和损失分别列账。粗粒化的意思，是把许多离散涡线平均成连续场；适用尺度不能任意缩到涡核内部。3
本库按照作者公开绘图代码的定义，从三组模拟的公开谱统计重新计算图中的能量通量。超流组分的摩擦项表现为损失，正常组分对应项可以为负——这里负值表示正常组分通过耦合得到能量，不是“负发热”。按两个组分的质量加权后，整体的摩擦损失仍为正。4
但是，这个模型有一条必须带着走的限定：作者给超流组分加入了一个人为设定的有效阻尼项，用来代表粗粒化以下没有显式解析的损失。 它不是测到了超流组分的普通分子黏性。所选三组模拟中，这个系数分别是正常组分运动黏性系数的0.2、0.1和0.2倍。34
把每个算例的整体注入功率定为1，重新算出的末端累计损失如下。表中保留两种损失项的区别：
| 模拟温度 | 正常组分黏性项 | 超流组分有效阻尼项 | 整体相互摩擦项 |
|---|---|---|---|
| 约2.16 K | 0.9617 | 0.0238 | 0.0064 |
| 约1.96 K | 0.7896 | 0.1484 | 0.0998 |
| 约1.44 K | 0.2403 | 0.7561 | 0.0333 |
它支持“组分间交换可以很强，整体摩擦损失却较小”，但不能被统一改写成“所有能量最终主要通过正常组分的普通黏性耗散”。尤其在所选1.44 K算例中，超流有效阻尼项承担了很大一部分损失。这个结果说明的是这套粗粒化模型如何分配未解析过程，不能替代对真实低温微观出口的验证。
三行的和也没有被强行归一到1。作者绘图脚本选取的是有限时间段，其间流动能量仍有增减。本库从同一批谱数据计算窗口两端的能量变化，并对损失功率作时间积分，得到“储存变化率＋损失率”相对注入的比值分别为0.9983、0.9984、0.9982。主要偏差可以由有限窗口中的储存变化解释，剩余约0.2%仍未完全闭合，本库未进一步确定其误差来源。15
这次重算使用的是350份已处理谱统计，没有重跑三维模拟，也没有重建作者的原始速度场。它核对的是记账方式和已发表数值的关系。
涡线密度告诉你有多少涡线，却不告诉你它们怎样组织
超流速度与凝聚态相位梯度有关。沿着避开涡核的闭合路径走一圈，相位要回到自身，只能累计整数倍的$2\pi$，因而环量取离散值。单量子涡旋的环量为：
氦-4的$\kappa$约为$9.97\times10^{-8}\,\mathrm{m^2/s}$。远离一条孤立直涡线的涡核，绕线速度大致按$1/r$减小。把周围流动的动能积分起来，每单位涡线长度的能量带有对数因子$\ln(\ell/a_0)$：内侧截止是涡核半径$a_0$，外侧截止与相邻涡线间距等尺度有关。12
通常把单位体积里的涡线总长度叫作$L$，单位是$\mathrm{m^{-2}}$；平均间距估计为$\ell\sim L^{-1/2}$。两组涡线可以有相同$L$，却有很不一样的方向关系：一组局部方向杂乱，另一组形成同向涡束，能够支持更大尺度的旋转运动。
这就像只知道一群人的总数，还不知道他们是四散走动，还是整队向同一方向行进。对涡旋而言，差别具体落实为速度场的空间相关和大尺度能量。反向涡线的粗粒化涡量还可能彼此抵消，而涡线长度不会抵消；因此，连续场的涡量大小也不能无条件代替真实的$\kappa L$。这正是粗粒化模型会遗漏部分无序涡线的原因之一。3
为什么会看到两种不同的衰减律
常用的现象学记法是把单位质量能量损失率写成：
这里的$\nu'$叫作有效运动黏性，是根据流动损失和涡线密度定义、拟合出的系数。它的量纲与普通运动黏性相同，却不是一项与流动组织无关的材料常数。2
如果没有大量大尺度运动，能量主要存在线附近，可以近似写成$E\simeq B\kappa^2L$，其中$B=\ln(\ell/a_0)/(4\pi)$。再把缓慢变化的$B$视为常数，就得到：
这类线能主导、缺少显著大尺度有序运动的衰减，通常称为Vinen或超量子湍流。这里$L$与能量近似成正比，后期二者都约按$1/t$减小。
另一类流动把很多能量存放在比涡线间距更大的结构里。如果形成Kolmogorov型谱、最大含能尺度已经受容器限制且近似固定，并且这部分大尺度能量占主导，那么大涡能量约按$(t+t_0)^{-2}$衰减。结合上一条损失率定义，就得到涡线密度的另一个指数：
这里$C$是谱常数，$d$代表固定最大尺度的估计。它不是说“量子世界突然变成经典世界”，而是说涡旋集体运动在大尺度上可以呈现近经典的能量分布，而小尺度仍受环量量子化约束。2
2008年Walmsley和Golov在低温氦-4中，通过带电涡环生成涡团，并用探测脉冲穿过涡团时的衰减估计$L$。代表性曲线测于0.15 K，低温比较涉及约0.5 K以下。在作者认定已充分扩展、进入晚期的时间段，短时注入形成的涡团显示接近$1/t$的衰减；与大尺度驱动形成的近经典湍流相比，所推得的低温有效黏性明显不同：前者较可靠的水平探测标定约为$0.08\kappa$，后者约为$0.003\kappa$。2
这支持组织形式会影响有效传能效率。但它没有让一条幂律成为万能分类器：论文自己指出，过渡阶段也可能暂时像$1/t$；探测方向和捕获截面的标定会影响绝对数值；近经典公式还依赖最大尺度和谱常数。根据上式反推$\nu'$时，假定的$d$若翻倍，推断值会改变四倍。这是计算条件，不是液体的材料黏性突然变了。215
涡线上确实有波，但“看见波”与“看见级联”有距离
Kelvin波可以想成一条涡线自身的螺旋状弯曲。波数越大，弯曲尺度越短。在细涡核、较小振幅等条件下，其色散关系大致为：
对数因子变化较慢时，常在有限范围内写成$\omega\approx C_K\kappa k^2$。不同波模之间的非线性相互作用，才有可能把输入某些波数的能量转移到更高波数。“波沿着涡线传播”与“能量持续跨越一系列尺度”是两个需要分别验证的问题。
Minowa等2025年发表的实验，在氦-4中利用附着于涡线的带电硅纳米颗粒激发波动，并从两个观察方向重建螺旋形状。实际取得的作者全文给出的温度是1.4 K，属于有限温度实验；正式公开图1及配套数据包含0.8—3 Hz的六个色散点。57
本库重新绘出六点及作者的99%波数区间。采用“频率方向残差最小”的未加权拟合，得到$C_K=0.526$；改成按公开波数区间宽度加权的波数方向残差，得到约0.569，接近作者使用的0.56。两种计算的误差模型不同，后者只作拟合敏感性检查，没有把99%非对称区间擅自解释成独立高斯标准差，也没有宣称重现了作者未公开说明的全部拟合步骤。15
不能因为曲线大致像$k^2$，就把绝对系数的差别略掉。在$k=10^4\,\mathrm{m^{-1}}$、涡核半径取0.1 nm的对数近似下，$C_K$约为1.10。论文补充材料明确讨论了差距：斜向观察会把波长投影缩短；示踪颗粒也可能影响有效涡核或运动。作者把“如果全部差距都解释成有效涡核变化”的反算写成约88 nm，但同时指出还需要进一步研究。6
所以，这项研究的直接贡献是可控激发、三维形状和传播色散的证据。88 nm不是独立测得的裸涡核半径，这六个受迫频率点也不是一条自发能量级联谱。
旋转超流体给出了更接近级联的证据
2023年Mäkinen等选择了另一种配置：让氦-3 B相中的涡线在旋转容器里大体平行，再对角速度作周期调制；在这类近乎平行的阵列中，间距近似满足$\ell\simeq\sqrt{\kappa/(2\Omega)}$。强旋转使涡线重联受到抑制，长尺度的集体惯性波仍可以被激发，继而与单条涡线上的短尺度波动联系起来。8
这是一项关键对照：重联可以产生Kelvin波，但并不是所有Kelvin波级联都必须不断依靠重联来供能。 若只是把“重联—Kelvin波—发热”画成唯一链条，就会漏掉这种近乎平行涡束中的波动路线。
实验利用NMR读出涡线取向相关的频移。停止调制后，先出现一段倾斜保持的阶段，然后才转入指数恢复。作者把前段与整体旋转中剩余能量继续经底部边界供给波动联系起来；后一段的衰减率记为$\tau^{-1}$。这里的NMR频移还要经过涡线取向模型解释，并非直接测得一条涡线每个波数的振幅。
本库取得的公开文件有142条已处理的衰减时间。在其中与作者斜率表对应的23个温度—压力组、共129条数据上，重新拟合$\tau^{-1}=A\Omega_f$。所有23个$A$都能算回，最大相对差约$4.34\times10^{-6}$。文件名虽有“Raw”，说明文件明确说这里已经是提取后的衰减时间；本次没有重新拟合原始NMR信号。9
为什么这些数据与级联有关？对一条固定波数的Kelvin波，相互摩擦给出的能量衰减率近似为$2\alpha\nu_\ell k^2$，其中$\nu_\ell$是带有环量和涡核对数的线张力参数，不是普通黏性，也不是前文由整体衰减定义的$\nu'$。如果波数相对涡线间距固定，降温使$\alpha$减小时，$A$应大致随$\alpha$一起减小。
实际最低温数据却趋于平台。作者的解释是：摩擦变弱后，波动能量可以抵达更大的波数、更短的尺度，$k^2$增大部分补偿了$\alpha$减小。因此，在有限温度实验中，总衰减对温度不敏感，并不等于最终损失项已经不包含相互摩擦。8
把这个推论写得更清楚，设波能谱在$k_0$至$qk_0$之间为$E_k\propto k^{-\eta}$，并暂时固定这段谱的形状。把各波数的摩擦损失加起来，瞬时的能量衰减率为：
这里$1<\eta<3$、$q>1$，$q\to1$时$F\to1$。它是对指定有限波数区间的加权平均，不是本库另行解出了一条自洽级联。若把它用于指数衰减时间，还要假定相关谱形和端点在拟合阶段足够稳定。
以9.6 bar、约$0.13T_c$的最低温组为例，取作者给出的$k_0=2.3/\ell$和近似线张力参数，其实测衰减率约为固定单尺度预测的9.03倍。采用$\eta=5/3$的有限谱积分，反推$q\approx7.34$；换成$\eta=7/5$，同一读数可对应$q\approx6.51$。采用论文的大带宽渐近式时，前一种谱会给出约8.76。三个数的差别来自明示的谱形和积分近似，均不是独立测得的截止波数。15
这就是证据的力量和边界：最低温衰减、涡线倾角、旋转速度依赖及模拟共同支持能量向较短尺度波动传递；但不能仅凭这一个衰减率，在不同Kelvin谱理论之间作唯一选择。作者也明确承认没有独立测量截止尺度、不能据此区分各谱模型。8
本库统一取$\nu_\ell=4\times10^{-8}\,\mathrm{m^2/s}$、$k_0=2.3/\ell$的近似时，23组中只有11组能反推出$q\geq1$，其余12组仍完整保留在中图。它说明这组统一参数和谱假设不能覆盖全部温压范围。论文讨论了最高温偏离可能涉及耗散尺度越出Kelvin波区间，但本库没有证据把这12组无解全部归于同一物理原因；参数近似和校准条件也需要分别检查。
此外，温度转成$\alpha$还存在系统误差：论文估计部分范围可达约30%，最低温又受温度计本底宽度影响。本篇沿用其校准，不将斜率复算的一致性当作温度校准也已独立重做。2023年的勘误补回了温度计压力外推公式，本篇已核对更正说明。14
2026年另有直接测量Kelvin波湍流、并辨别多波相互作用的工作，研究对象是水中的宏观涡核。它为波动理论提供了有价值的经典实验平台，但水涡环量不量子化，涡核也不是超流氦的涡核，不能把其证据移植成“氦的微观耗散链已经测全”。13
更短尺度的能量，最后怎样离开涡旋运动
Kelvin波级联解决的是传能问题，最后还要有能量出口。对于可压缩的玻色超流模型，涡旋运动和重联能激发密度波。Leadbeater等2001年的Gross–Pitaevskii计算显示，两个涡环发生重联后，出现向外传播的声脉冲；研究通过声脉冲离开一个指定测量体积时的能量下降，估算声辐射。11
这不是给无黏性方程偷偷补一个普通摩擦项。涡旋运动与密度波本来就是这个模型里可以互相转换的自由度。测量体积里的能量下降，是因为声波把能量带出了边界；对封闭、无外加损失的理想模型，把所有自由度一起记账，总能量仍守恒。 声能随后是否热化、返回或被容器吸收，取决于系统和边界，不能从一张涡线减少的图直接读出。
这项计算也有物理范围：GP模型能研究重联与声耦合，但不能完整代表强相互作用液氦的全部微观性质。它证明了特定模型中的传能通道，不等于给出了所有氦-4实验里各通道的贡献比例。
氦-3 B相还有不同的涡核内部结构。Silaev提出的低温理论中，加速运动的涡核可以驱动核内束缚准粒子偏离平衡，并向外释放非平衡准粒子携带的能量。即使远处热浴温度很低，运动本身也能产生这条非平衡出口；并不要求一开始就有大量正常组分。12
这是针对费米超流的理论机制，不能直接拿来解释氦-4，也不能因标题使用“普适”就把它视为每种低温实验都已独立确认的最终答案。本文未重解该微观动力学模型，更没有据此宣称实际实验已经达到严格的零温极限。
容器壁传递的角动量，与流体损失的能量不是同一本分账
边界条件在这里并不只是仪器细节。2013年Hosio等在光滑壁面的旋转氦-3 B相容器中研究传播的涡旋前沿，比较前沿速度和其后流体的旋转速度。要同时解释两种读数，他们分别引入与能量收支、角动量耦合有关的有效参数。10
模型拟合的低温残余参数约为$\widetilde\alpha_{\rm en}=0.20$与$\widetilde\alpha_{\rm am}=0.0019$，相差约两个数量级。它们是现象学参数，不是直接把能量损失率和角动量损失率相除；后两者的量纲本来就不同。
其物理含义是：流体内湍动可以有效地把有序运动能交给小尺度自由度，而流体与光滑容器之间的转矩耦合仍较弱。于是，流体不必紧紧跟着容器一起转。作者并未唯一确定残余角动量交换到底来自壁面附着涡旋、湍流贡献还是小涡环到达边界等过程。
这与2023年实验强调粗糙底部向波动供能并不冲突：接触条件、涡线组织和测量阶段不同，边界承担的作用也不同。一项研究中的“壁面耦合很弱”，不能扩展成“边界在超流湍流中都不重要”。
把条件换掉，哪些判断还站得住
目前较稳的认识是，超流组分缺少普通剪切黏性，并不阻止宏观流动经其他自由度衰减；量子涡旋及其组织形式，是把大尺度运动与小尺度过程联系起来的关键。直接波观测、旋转衰减实验、涡环实验与模型计算，分别约束这条链的不同环节。
但要判断某一次衰减究竟由谁控制，至少还要核对四件事：
- 温度和物质是什么。 有限温度氦-4的两组分交换、超低温氦-3的涡核准粒子，不能不带条件互换。
- 能量最初放在哪里。 大尺度涡束、杂乱涡团、受迫单波和旋转波动的库存不同；同一个$L$或相同外加功率不足以让它们可比。
- 测到的是哪一笔账。 涡线长度、速度谱、NMR频移、量热信号和容器转矩各有读出模型。流动衰减、传出观测区和转成热要分别说明。
- 最小尺度是否真的被解析。 涡丝模型中的重联处理、删点或有限网格，可能引入数值能量损失；HVBK有效阻尼也只代表未解析过程。图上能量下降必须和实际损失项、尺度收敛及边界对应，不能直接命名为真实微观耗散。310
仍未解决的部分，包括不同涡旋组织中各小尺度通道的定量份额、级联谱和截止尺度的独立测量，以及低温边界耦合的具体来源。这些缺口不妨碍理解已观察到的衰减，却限制了把某一条路线推广到所有超流体的程度。