Navier-Stokes hydrodynamics near the critical point with out-of-equilibrium modes

技术摘要：临界点附近的纳维-斯托克斯流体动力学与非平衡模态

问题陈述
标准的流体动力学描述（包括纳维-斯托克斯方程）依赖于系统处于局部热平衡以及相关长度（$\xi$）远小于涨落特征波长（$\lambda$）（即 $\xi \ll \lambda$）的假设。在流体的临界点（CP）附近，相关长度会发散，违反了这一条件。因此，流体脱离了平衡态，使得标准的流体动力学预测失效。虽然已有研究提出通过梯度展开方案（如二阶 Müller-Israel-Stewart 理论）来扩展流体动力学，但这些方案往往难以收敛，或无法充分描述具有大规模涨落的流体。 “Hydro+” 框架提供了一种替代方案，通过引入一个参数化慢模（$\phi$）来解释这些大规模的非平衡涨落。然而，该慢模对动态结构因子的具体影响，以及在非相对论纳维-斯托克斯框架内由热力学一致性所施加的约束，仍需进一步研究。

方法论
作者通过将一个缓慢的、非流体动力学的标量模 $\phi$ 引入熵的定义中，扩展了非相对论纳维-斯托克斯流体动力学框架。其方法如下：

1. 热力学公式化： 熵密度被推广为 $s_+(\epsilon, \rho, \phi)$，其中 $\epsilon$ 是能量密度，$\rho$ 是质量密度。吉布斯-杜亥姆（Gibbs-Duhem）关系被修改，以包含共轭力 $\pi = (\partial s_+ / \partial \phi)_{\epsilon, \rho}$。
2. 控制方程： 作者基于质量、动量和能量守恒定律，并辅以一个关于慢模 $\phi$ 的演化方程，推导出一组闭合的线性化流体动力学方程。该演化方程为弛豫类型：
   $$\frac{\partial \phi}{\partial t} + \nabla \cdot (\phi \mathbf{v}) = -\Gamma_\phi \pi - A_\phi (\nabla \cdot \mathbf{v})$$
   其中 $\Gamma_\phi$ 是弛豫速率，$A_\phi$ 代表该模态对压缩/膨胀的敏感度。
3. 热力学约束： 推导过程中的关键步骤是强制执行热力学第二定律（$\partial_\mu S^\mu_+ \geq 0$）。作者证明，外延性条件（extensivity condition）对慢模与速度散度之间的耦合施加了严格的约束。具体而言，为了确保在任意流场下熵产生均为正定，与 $(\nabla \cdot \mathbf{v})$ 项相关的系数 $S_n$ 必须消失（即 $S_n = 0$）。这导致了热力学变量与敏感度 $A_\phi$ 之间存在特定的关系。
4. 线性模态分析： 在均匀平衡态附近对控制方程进行线性化。作者在傅里叶空间中求解密度涨落 $\delta \rho(k, \omega)$，以导出表征密度涨落功率谱的动态结构因子 $S_{\rho\rho}(k, \omega)$。
5. 数值评估： 使用来自文献 [22] 的物态方程和来自文献 [23] 的输运系数（粘度、热导率）及响应函数的标度律，对氦-4 ($^4$He) 进行数值评估。计算分别在远离临界点（$T = 20$ K）和靠近临界点（$T = 6$ K，$T_c = 5.2$ K）的情况下进行。

主要贡献与结果

• 对耦合的约束： 本研究确立了热力学外延性要求对慢模 $\phi$ 与速度四散度之间耦合的约束。这一约束对于维持熵产生的正定性是必要的，并将 Hydro+ 构架与任意的纳维-斯托克斯理论扩展区分开来。
• 动态结构因子的修正：
  • 远离临界点时： 在不存在慢模（$\phi=0$）的情况下，$S_{\rho\rho}(k, \omega)$ 呈现标准的三峰结构：一个中心瑞利峰（热扩散）和两个对称的布里渊峰（声传播）。引入非零的 $\phi$ 后，分布的对称性得以保持，但表现为：
    • 布里渊峰向零频率方向移动，表明由于压力响应变慢，有效声速降低。
    • 所有峰的宽度均减小，意味着涨落的衰减速率降低。
    • 布里渊峰的高度出现不对称性。
  • 靠近临界点时： 随着系统接近临界点，由于压缩率的发散和声速的消失（$c_s \to 0$），布里渊峰完全消失。瑞利峰依然存在，但在存在 $\phi$ 的情况下，其宽度显著减小。这种变窄现象表明了“临界减慢”（critical slowing down），即涨落的衰减被显著延迟。
• 慢模 ($\phi$) 与耦合 ($b$) 的作用：
  • 增加平衡态慢模的值 ($\bar{\phi}$) 会导致瑞利峰和布里渊峰的高度单调增加。作者将其归因于慢模充当了涨落能量的储库，在释放能量之前暂时捕获能量，从而增强了谱权重。
  • $\phi$ 的效应在临界点附近比在远离临界点时更为显著，这与临界减慢现象一致。
  • 耦合参数 $b$（代表慢模与热通量之间的相互作用）会抑制峰值高度。较强的耦合引入了额外的弛豫通道，增加了热涨落和声学涨落的阻尼。

意义与主张
本文声称，引入非平衡模（OEM）从根本上重塑了临界点附近动力学密度涨落的功率谱。其主要意义在于展示了受热力学第二定律约束的 Hydro+ 框架，能够提供一种一致的流体动力学描述，捕捉到涨落的减慢以及谱权重的重新分布。

具体而言，作者主张：

1. 外延性条件是一个非平凡的约束，它决定了慢模与流体动力学梯度的耦合方式。
2. 慢模 $\phi$ 的存在降低了涨落的衰减速率（使峰宽变窄），为临界减慢现象提供了流体动力学层面的解释。
3. 慢模作为一个源，重新分配了涨落能量，增强了动态结构因子中各峰（特别是临界点附近的瑞利峰）的幅度。

研究结论认为，虽然慢模不会产生新的流体动力学模态，但它显著地改变了现有模态的动力学特性，为在扩展的流体动力学框架内理解临界现象提供了理论基础。作者建议未来的工作可以探索非线性涨落、多重慢模以及与实验光散射或中子散射数据的直接对比。