[論文レビュー] 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

ひとこと要約

本稿は、高レイノルズ数における2次元チャネル内の対称的せん断流れのスペクトル不安定性について、漸近的マッチングを回避する新しい作用素に基づくアプローチを用いて、厳密な数学的証明を提供する。不安定な固有モードの成長率は $ e^{t/\sqrt{\alpha R}} $ であり、不安定性は3つの領域で発生する：$ \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.

研究の動機と目的

• 一般に非線形コウエットでない対称的せん断流れが高レイノルズ数で線形不安定であるという長年の物理的予想に対する完全な数学的証明を提供すること。
• ノースリップ境界条件を満たす2次元チャネルにおける非圧縮性ナビエ＝ストークス方程式のスペクトル不安定性問題を解明すること。
• $ R \to \infty $ の極限において、不安定固有モードの正確な成長率と不安定性領域を特定すること、特に臨界安定性領域を対象とすること。
• 従来の内層・外層漸近層のマッチングを回避する新しい作用素論的枠組みを構築すること。
• 臨界層における特異性を点ごとのグリーン関数のバインディングで分析することで、正確な不安定固有値および固有関数を導出すること。

提案手法

• ノースリップ境界を持つチャネルにおける線形化ナビエ＝ストークス方程式を定式化し、正規モードのアンザッツ $ e^{\lambda t} $ を用いてスペクトル問題を解析する。
• レイリー方程式と修正されたアイル方程式を用いて、解をゆっくりと変化するオール・モードと速く変化するオール・モードに分解することで、オール＝ゾンマーフェルト方程式を新しい作用素的アプローチで分析する。
• $ \alpha \neq 0 $ に対してレイリー方程式の正確な解を構成し、漸近的に境界条件を満たす2つの特別な解を導出する。
• 臨界層における特異性を制御するため、レイリー作用素およびアイル作用素のグリーン関数に対する点ごとのバインディングを導出する。
• ランガー変換を適用し、粘性補正を扱うために修正されたアイル方程式の近似グリーン関数を構成する。
• 畳み込み推定を確立し、修正されたアイル方程式を反復的に解くことで、臨界層付近における固有関数の挙動を制御する。

実験結果

リサーチクエスチョン

• RQ1 物理理論が予測するように、無限大のレイノルズ数の極限において、2次元チャネル内の対称的せん断流れのスペクトル不安定性は継続するか？
• 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}} $ を持つ不安定モードの存在が確認され、ハイゼンベルク、リン、トールマイエン、ドレイジン、レイトの物理的予測と整合的である。