编辑推荐：
研究人员研究了无限深流体自由表面上具有延伸到有限深度的常涡度剪切层的大振幅行进重力波。通过使用共形映射（conformal mapping）描述问题的两个自由边界，建立了边界积分（boundary-integral）公式，并使用伪谱方法（pseudospect
研究人员研究了无限深流体自由表面上具有延伸到有限深度的常涡度剪切层的大振幅行进重力波。通过使用共形映射（conformal mapping）描述问题的两个自由边界，建立了边界积分（boundary-integral）公式，并使用伪谱方法（pseudospectral method）进行了数值求解。线性波速 $c$ 满足三次色散关系（cubic dispersion relation），该关系本质上由两个重力模态和一个剪切模态组合而成。在某些参数范围内，这些模态的碰撞导致线性不稳定性。与 $c$ 的三个值之一相对应的非线性解曾被 Milinazzo & Saffman（1990 J. Fluid Mech. 216, 93–101）研究过。这些解由表面重力模态主导；然而，与其余两个 $c$ 值相关的解被忽视了。研究人员详细研究了这些新解，并发现了与 Milinazzo & Saffman 计算出的解不同的特性，例如涡度界面上大幅值波的出现以及两个自由边界之间相对相位的变化。此外，研究显示，在共振条件下，不同波数的解支之间会出现剪切模态与重力模态之间的分岔（bifurcation），这导致在非常小的参数空间区域内解剖面发生显著的定性变化。
**论文解读：剪切层上非线性表面重力波的数值研究与新解发现**
**一、 研究背景与问题**
当风吹过气-水界面时，水面受到的剪切应力会诱发剪切流，其不稳定性是海洋表面波产生的主要机制之一。在真实海洋中，这种剪切流的涡度通常集中于水面附近并随深度衰减。为了在数学上处理这一物理情境，研究人员构建了最简单的理论模型：一个由两层不同涡度组成的无粘等密度流体。下层流体无限深且无旋，上层流体则具有均匀的常涡度，即一个有限深度的剪切层。该模型在数学上构成一个具有两个自由边界的非线性问题：一个是涡度跃变处的内界面，另一个是剪切层顶部压强恒定的自由表面。
此前的研究已经建立了该模型的线性色散关系，它是一个关于波速 $c$ 的三次方程，物理上对应于两个表面重力模态与一个剪切模态的耦合。在某些参数范围内，这些模态的碰撞会导致线性不稳定性。在非线性层面，Milinazzo & Saffman（1990）曾利用边界积分方法计算了全非线性行进波解，但仅关注了三个波速解中对应于表面重力模态、且沿剪切流方向传播的其中一个分支，并将其用于与 Banner & Phillips（1974）关于最大波幅对表面漂移依赖关系的实验结果进行比较。然而，三次色散关系中剩余两个波速解对应的非线性解支被完全忽略了。本文研究的目的正是系统地研究这两个被忽视的解支所对应的非线性行波解，探索其全新的物理特性。该研究发表在《Journal of Fluid Mechanics》上。
**二、 关键技术与方法**
研究人员采用时变共形映射（conformal mapping）将物理平面上的两个自由边界映射到计算平面上的固定位置，并结合无旋下层流体的 Hilbert 变换（Hilbert transform）与有旋上层流体的 Cauchy-Riemann 方程，推导出稳态边界积分公式。数值上采用伪谱方法（pseudospectral method），以傅里叶级数展开未知函数，利用快速傅里叶变换（FFT）计算导数与积分变换，并通过牛顿迭代（Newton iteration）求解非线性代数方程组。为追踪解支随参数的变化，使用了弧长延拓（arclength continuation）技术。针对涡度界面主导的解中出现的严重混叠误差，引入了36阶傅里叶滤波器（$\exp {(-36 (1.2 m/N)^{36})}$）进行去混叠处理。
**三、 研究结果**
**3.1 有限振幅数值解**
在固定无量纲波长 $L=5$、振幅 $A=0.03$ 的条件下，研究人员通过数值延拓获得了与线性解支对应的三条 $n=1$ 非线性解支。其中波速 $c$ 最大的上支和右下支在定性上与线性解支一致，因为波速的一阶弱非线性修正为零。然而，左下支则发生了定量和定性上的显著变化：它不再是一条连续曲线，而是在多个位置分裂成若干与更大 $n$ 值（如 $n=2$、$n=3$）解支发生分岔的独立曲线，且可能存在无穷多条不相交的分支。这种行为源于不同波数线性解支的交叉导致的共振。在 $H \to 0$ 的极限下，该分支上的内波波峰显著变宽，形成狭窄而深的波谷，使解接近纯内孤立波（internal solitary wave）形态。研究还详细展示了 $(1,2)$ 和 $(1,3)$ 分岔附近的两个 $n=1$ 分支从同一 $n=2$ 或 $n=3$ 分支的不同位置或相同位置分岔的现象。这些分岔发生在极小的参数空间范围内，但会导致解剖面的显著变化，包括主导振幅界面的反转（由自由表面主导变为涡度界面主导）以及波峰极性的改变。
**3.2 极限振幅解**
为研究解的极限形态，研究人员选取了六个代表性解（对应于线性解中的 $(a,c,d,f,g,i)$），不断增加振幅 $A$ 直至傅里叶系数谱出现平台，以此判断解是否趋近于具有尖角的极限波。结果显示，不同解首先趋于极限的界面不同：解 $(d)$ 和 $(i)$ 的涡度界面 $\zeta_2(x)$ 首先形成尖角，而解 $(a)$、$(c)$、$(f)$、$(g)$ 的自由表面 $\zeta_1(x)$ 首先形成尖角。各解的极限振幅差异显著，自由表面和涡度界面的极限振幅分别在不同参数下表现出不同的量级。流线和驻点分析表明，在剪切层内部出现了闭合流线和驻点（stagnation point），这对于具有常涡度的流体而言，揭示了内部流动结构的复杂性。
**四、 结论与讨论**
研究人员通过共形映射的边界积分公式和伪谱方法，成功计算了剪切层上同时受常背景涡度和重力作用的三种行进波解支，重点研究了以往被忽视的两个新解支。研究发现，剪切模态的解中，涡度界面位移 $\zeta_2$ 远大于自由表面位移 $\zeta_1$；而重力模态的解中，$\zeta_2$ 与 $\zeta_1$ 的相对大小取决于无量纲剪切层厚度 $H$。在基态失稳前，两个界面具有相同极性（波峰对齐），而在基本流重新稳定后，两者极性相反。此外，当剪切模态与重力模态的波长成整数倍关系时会发生共振，导致非线性解支在极小的参数空间内发生分岔，例如 $(1,2)$ 和 $(1,3)$ 分岔。虽然这些分岔现象类似于重力-毛细波研究中的 Wilton 涟漪，但本文在解剖面和傅里叶谱中均未观察到指数小振荡成分。研究结论表明，该模型为理解风生剪切流作用下的非线性表面波及内部界面波提供了新的解族，其分岔特性对海洋表面波生成机制的研究具有重要意义。未来的工作将聚焦于推导完整的时变共形映射演化方程，以研究这些非线性行波解的线性稳定性及受扰剪切层的一般演化过程，并可能将模型扩展到密度分层海洋中薄涡度层对大振幅内波的影响研究。
