[论文解读] Spectral instability of symmetric shear flows in a two-dimensional channel

Emmanuel Grenier, Yan Guo+1 more|arXiv (Cornell University)|Feb 6, 2014

Navier-Stokes equation solutions参考文献 13被引用 23

一句话总结

本文通过一种新颖的基于算子的方法，避免渐近匹配，对高雷诺数下二维通道内对称剪切流的谱不稳定性提供了严格的数学证明。研究发现不稳定特征模态的增长率为 $ e^{t/\sqrt{\alpha R}} $，不稳定性存在于三个区域：$ \alpha \approx R^{-1/7} $，$ R^{-1/7} \ll \alpha \ll R^{-1/11} $，以及 $ \alpha \approx R^{-1/11} $，其中增长率分别按 $ \sim A^{-1/3}R^{-2/7} $ 和 $ \sim A^{-3/2}R^{(3\beta-1)/2} $ 规律变化。

ABSTRACT

This paper concerns spectral instability of shear flows in the incompressible Navier-Stokes equations with sufficiently large Reynolds number: $R o \infty$. It is well-documented in the physical literature, going back to Heisenberg, C.C. Lin, Tollmien, Drazin and Reid, that generic plane shear profiles other than the linear Couette flow are linearly unstable for sufficiently large Reynolds number. In this work, we provide a complete mathematical proof of these physical results. In the case of a symmetric channel flow, our analysis gives exact unstable eigenvalues and eigenfunctions, showing that the solution could grow slowly at the rate of $e^{t/\sqrt {αR}}$, where $α$ is the small spatial frequency that remains between lower and upper marginal stability curves: $α_\mathrm{low}(R) \approx R^{-1/7}$ and $α_\mathrm{up}(R) \approx R^{-1/11}$. We introduce a new, operator-based approach, which avoids to deal with matching inner and outer asymptotic expansions, but instead involves a careful study of singularity in the critical layers by deriving pointwise bounds on the Green function of the corresponding Rayleigh and Airy operators.

研究动机与目标

• 为长期存在的物理猜想提供完整的数学证明，即除线性库埃特流外，一般对称剪切流在高雷诺数下是线性不稳定的。
• 解决无滑移边界条件下二维通道中不可压缩纳维-斯托克斯方程的谱不稳定性问题。
• 精确刻画当 $ R \to \infty $ 时不稳定特征模态的增长率与不稳定性区域，特别是临界稳定区域。
• 构建一种新的算子理论框架，避免传统内层与外层渐近展开的匹配。
• 通过分析临界层中瑞利与艾里算子格林函数的奇点，推导出精确的不稳定特征值与特征函数。

提出的方法

• 在无滑移边界条件下建立通道内的线性化纳维-斯托克斯方程，并通过正规模假设 $ e^{\lambda t} $ 研究谱问题。
• 采用新颖的基于算子的方法分析奥尔-索默菲尔德方程，通过瑞利方程与修正艾里方程将解分解为慢速与快速奥尔模态。
• 为 $ \alpha \neq 0 $ 构造瑞利方程的精确解，并推导出两个满足边界条件渐近行为的特解。
• 推导瑞利与艾里算子格林函数的逐点界，以控制临界层中的奇点。
• 应用兰格变换，并为修正艾里方程构造近似格林函数，以处理粘性修正项。
• 建立卷积估计，并通过迭代方法求解修正艾里方程，以控制特征函数在临界层附近的行为。

实验结果

研究问题

• RQ1 对称剪切流在二维通道中的谱不稳定性是否在无限雷诺数极限下依然存在，如物理理论所预测？
• RQ2 当 $ R \to \infty $ 时，不稳定特征值虚部（增长率）在不同 $ \alpha $ 频率区域的精确标度是什么？
• RQ3 临界层如何影响不稳定性？其奇异性行为能否在不匹配内层与外层渐近展开的前提下得到控制？
• RQ4 能否为非库埃特对称剪切剖面（包括中间与上临界稳定分支）构建完整的数学证明？
• RQ5 瑞利与艾里算子格林函数在确定不稳定性阈值与增长率方面起什么作用？

主要发现

• 当 $ \alpha \approx R^{-1/7} $ 时，特征值虚部按 $ \sim A^{-1/3}R^{-2/7} $ 标度，当常数 $ A $ 超过临界阈值 $ A_c $ 时出现不稳定性。
• 在中间区域 $ R^{-1/7} \ll \alpha \ll R^{-1/11} $，增长率按 $ \sim A^{-3/2}R^{(3\beta-1)/2} $ 标度，且对所有足够大的 $ A $，不稳定性必然成立，其中 $ \beta \in (1/11, 1/7) $。
• 当 $ \alpha \approx R^{-1/11} $ 时，色散关系中的 $ \alpha^4 \log \alpha $ 项变得显著，并导致随着 $ A $ 超过临界值 $ A_{2c} $，不稳定性向稳定性转变。
• 临界层位置 $ z_c $ 在 $ R^{-1/7} $ 区域满足 $ z_c / \delta \approx A^{4/3} $，且当 $ R \to \infty $ 时保持有界，表明其靠近边界。
• 该方法通过瑞利与艾里算子格林函数的逐点界，成功避免了内层与外层渐近展开的匹配。
• 分析证实了存在增长率为 $ e^{t/\sqrt{\alpha R}} $ 的不稳定模态，与海森堡、林、托尔米恩、德雷泽因和雷德的物理预测一致。