Analysis of inviscid shear instability of axisymmetric flows

本文通过建立更精确的稳定性充分条件和一种新颖的不稳定性条件，推导出了管道和环形空间中轴对称流动无粘稳定性的改进解析判据，这两个条件均能有效预测经特征值计算验证的中性参数。

要点

想象你是一位流体力学家，试图预测管道或环形通道内平滑旋转的水流何时会突然变得混乱并转为湍流。通常，这需要运行耗时数小时甚至数天的大型复杂计算机模拟。

本文介绍了一组新的“经验法则”，使科学家能够利用简单的数学甚至纸上的快速草图，更快速地预测稳定性。以下是他们所做工作的分解，并辅以日常类比。

问题：“临界点”

想象流体流经管道或环形区域（如同甜甜圈）。有时，流动完全平滑（稳定）；有时，微小的涟漪会发展成巨大的波浪，导致湍流（不稳定）。

科学家们早已知道如何判断流动是否“绝对安全”（稳定），但很难找到一个简单的规则来判断流动何时“绝对不安全”（不稳定）。这就像确切知道桥梁何时不会坍塌，却没有任何简单的方法能预测它何时“会”坍塌，除非测试每一辆驶过的卡车。

新工具：两道“安全网”

作者开发了两项新的分析工具（定理）作为安全网。

1. “安全天花板”（稳定性条件）

• 旧方法： 科学家曾使用 1962 年巴切勒（Batchelor）与吉尔（Gill）提出的规则，它像一道低矮的天花板。如果流动保持在这道天花板之下，就是安全的。但这道天花板往往太低，导致许多实际上安全的流动被遗漏。
• 新方法： 作者构建了一道更高、更智能的天花板（基于“凯尔文 - 阿诺尔德第二定理”）。想象一位空中飞人演员。旧规则说：“如果你保持在这道低矮的横杆之下，你就不会坠落。”新规则则说：“实际上，你可以荡得高得多才处于危险之中。”
• 工作原理： 他们观察代表流动的具体数学曲线。如果该曲线保持在某条“安全线”之下（该线随管道形状而变化），则该流动被保证是稳定的。

2. “跨栏跳跃”（不稳定性条件）

• 概念： 这是本文最令人兴奋的新想法。想象一名运动员试图跨过一个栏架。
  • 在直管（平行流）中，栏架是一根平直的横杆。
  • 在环形或带中心的管道中，栏架呈山丘状或曲线状。
• 规则： 如果流动的数学曲线跃过了这个栏架，则该流动被保证会变得不稳定（混乱）。
• 特殊性： 在此之前，寻找“不稳定”流动需要复杂的计算。现在，你只需绘制曲线并观察它是否越过了栏架。如果越过了，你立刻就能知道湍流即将来临。

“栏架”形状至关重要

作者意识到，“栏架”的形状取决于几何结构：

• 在环形区域（Annulus）中： 栏架是平坦的，高度恒定。就像田径比赛中的标准栏架。
• 在管道中： 栏架很棘手。靠近管道中心处，规则发生变化。栏架不是平坦的，而像是一个斜坡，越靠近中心越陡峭。如果流动曲线试图跃过这个斜坡，它就会失败（变得不稳定）。

测试规则

为了证明这些规则有效，作者在两个特定的“模型”流动上进行了测试：

1. 环形流动： 水在两个圆柱体之间流动，外圆柱加热，内圆柱冷却，且两个圆柱体相互滑动。
2. 管道流动： 水在从内部加热的管道中流动。

他们将简单的“跨栏”预测与大规模的计算机模拟（“金标准”）进行了比较。

• 结果： 他们的简单规则出奇地准确。他们正确识别了流动从平滑转为混沌的“临界点”（中性稳定性）。
• 益处： 科学家不再需要为每种可能的情况运行计算机模拟，现在可以使用这些简单的图表来缩小搜索范围。这就像在开始挖掘之前，先用金属探测器寻找埋藏的宝藏。

他们未声称的内容

作者谨慎地说明了他们的规则不能做到什么：

• 粘度（粘性）： 这些规则假设流体没有“粘性”（无粘）。在现实世界中，流体是有粘性的。虽然这些规则在粘性影响较小的快速流动中效果良好，但它们并未考虑仅由粘性引起的特定类型的不稳定性（如著名的托尔明 - 施利希廷波）。
• 射流： 这些规则非常适用于管道和环形区域，但尚未完全解决“射流”问题（流体射入开放空间的射流，如花园水管）。开放空间的数学要困难得多，因为那里的“栏架”没有清晰的边界。

总结

本文为流体动力学家提供了一种新的、简单的方法，用于预测管道和环形区域中的旋转流动何时会失控。通过将复杂的计算机模拟替换为简单的“跨栏跳跃”检查，他们可以迅速识别哪些流动是安全的，哪些注定会变为湍流。

技术摘要

技术摘要：轴对称流动无粘剪切稳定性分析

问题陈述
本文探讨了环形区域和管道中轴对称基流的无粘稳定性问题，该问题在从圆射流到内部管道流动的各类应用中具有重要意义。尽管高雷诺数流动通常允许忽略粘性，从而简化为欧拉方程，但在轴对称几何构型中推导精确的解析稳定性判据一直颇具挑战。Batchelor 与 Gill（1962）的先前工作将 Rayleigh-Fjørtoft 条件（Kelvin-Arnol'd 第一定理，KA-I）推广至轴对称流动，但针对这些几何构型，第二 Kelvin-Arnol'd 稳定性定理（KA-II）以及不稳定的简单充分条件尚未被明确推导出来。此外，现有的识别中性稳定性点的方法通常依赖于复杂的积分方程或完整的特征值计算，限制了其在快速参数估算中的实用性。

方法论
作者在圆柱坐标系中构建了轴对称基流 $U(r)$ 的线性化无粘稳定性问题。该问题简化为关于类流函数特征函数 $G(r)$ 的微分方程，由复相速度 $c$ 和波数比 $N = n/k$ 控制。

核心方法论涉及基于以下框架构建两个新的解析定理：

1. 稳定性（KA-II 的推广）： 遵循 Arnol'd（1966）基于伪能量的分析方法，作者推导出了稳定性充分条件（KA-II），该条件改进了经典的 Batchelor 与 Gill（1962）结果。这涉及定义量 $W_{\alpha,N}(r) = rQ'/(U-\alpha)$，并建立源自适用于环形和管道几何构型的 Poincaré 型不等式的上界 $H(r)$。
2. 不稳定性（障碍定理的推广）： 通过将 Deguchi 等人（2024）针对平行流动提出的“障碍定理”进行推广，作者推导出了不稳定性充分条件。该条件指出，如果在临界层（即 $U = \alpha$ 处）$W_{\alpha,N}$ 超过了特定的“障碍”函数 $h(r)$，则存在不稳定模态。推导采用两步法：(i) 利用试验函数最小化 Rayleigh 商来证明中性解的存在性；(ii) 证明波数的微小扰动会导致具有正增长率的模态。

理论预测通过与无粘特征值问题的数值解（采用 Chebyshev 配置法）进行验证，并与高雷诺数下的线性化 Navier-Stokes 计算结果进行了比较。

主要贡献
本文提出了两个主要定理，提供了评估稳定性的简单解析判据：

• 定理 1（稳定性）： 确立了若对所有波数比 $N$ 均满足 KA-I 或 KA-II 条件，则基流是稳定的。
  • KA-I： 要求存在实数 $\alpha$，使得 $W_{\alpha,N} \leq 0$ 处处成立。
  • KA-II： 要求存在实数 $\beta$，使得 $W_{\beta,N}$ 位于特定范围 $[0, H(r)]$ 内。阈值函数 $H(r)$ 已针对环形流动（依赖于半径比 $\eta$）和管道流动（依赖于贝塞尔函数零点）显式推导得出。这将 KA-II 框架扩展到了轴对称几何构型。
• 定理 2（不稳定性）： 确立了若对于某些 $N$，$W_{\alpha,N}$（其中 $\alpha = U(r_c)$ 为临界点处取值）超过了障碍函数 $h(r)$，则基流是不稳定的。
  • 对于环形流动，障碍 $h$ 是一个源自窄间隙极限的常数。
  • 对于管道流动，障碍被修正为幂律形式（$Cr^p$），以考虑管道轴线（$r=0$）附近的奇异性行为，其中 $W_{\alpha,N}$ 表现为 $r^2$。

结果
作者将这些定理应用于两个模型流动：

1. 环形流动： 由壁面滑动和热对流（垂直环形腔内的自然对流）驱动的同心圆柱间流动。在浮力参数 $\chi$ 和半径比 $\eta$ 的参数空间稳定性图中，蓝色区域（由定理 1 保证稳定）和红色区域（由定理 2 保证不稳定）包围了数值计算得到的中性曲线。这些定理有效地预测了中性点，特别是在窄间隙极限下，结果收敛于平行流动障碍定理。
2. 管道流动： 带有均匀内部加热的垂直 Hagen-Poiseuille 流动。分析表明，虽然 KA-I 条件保证了 $\chi < 4$ 时的稳定性，但定理 2 识别出了 $\chi \in [5.72, 7.42]$ 范围内的不稳定区域。数值确定的中性点出现在 $\chi \approx 7.86$ 处。作者指出，对于 $\chi \in [4, 6]$，不稳定性可能源于奇异中性模态，这超出了基于正则模态的定理 2 的适用范围，从而解释了理论不稳定区域与精确中性曲线之间的差距。

在这两种情况下，解析判据成功预测了从完整特征值计算中获得的中性参数位置，并与高雷诺数下的粘性线性稳定性结果一致。

意义与主张
本文主张，这些定理提供了“简单的解析工具”，用于估算特征值 $c$ 的行为及参数空间中中性点的位置。其主要意义在于能够：

• 推广 KA-II： 明确推导轴对称流动的第二 Kelvin-Arnol'd 稳定性条件，这是此前文献中缺失的结果。
• 推广障碍定理： 将针对平行流动的障碍定理适应于轴对称几何构型，提供了一种无需求解完整特征值问题即可识别不稳定性的实用图解方法。
• 降低计算成本： 通过提供稳定性边界的粗略估计，这些判据可以限制数值特征值问题中需要探索的参数范围。

作者对适用范围保持谦逊，指出结果严格适用于无粘近似。他们承认，即使在理论稳定区域内，粘性也可能诱发不稳定性（例如 Tollmien-Schlichting 波）。此外，他们明确指出，由于缺乏适用的 Poincaré 型不等式以及流动渐近行为的差异，这些定理不直接适用于圆射流等无界域，尽管他们建议未来工作可为此类情况设计合适的试验函数。