A Compressible Miles Instability for Wind-Wave Generation in a Neutral Atmosphere

技術要約：中立大気における風波生成のための圧縮性マイルズ不安定性

問題設定
本研究は、重力と表面張力を考慮した、大気中の圧縮性オイラー流と海洋中の非圧縮性オイラー流を隔てる平坦な界面の正規モード安定性を扱うものである。基底状態は、水の下に静止した水、その上に定数エントロピーの剪断流 $U(z)$ を持つ、高度 $U(0)=0$ の大気からなる。この構成は、風波の臨界レベルが存在する海面上の数十メートルの標準的な記述である中立成層大気を表している。

中心となる問いは、空気が非圧縮性流体ではなく圧縮性ガスとして扱われる場合、古典的なマイルズ（Miles）メカニズム（すなわち、臨界レベル $U(z_c) = c_r$ で平均風の剪断からエネルギーを抽出することで波が増幅するという仕組み）がどのように修正されるかである。非圧縮性理論（B"uhler, Shatah, Walsh, および Zeng によって厳密に確立されたもの）は、風速プロファイルの曲率 $U''(z_c)$ に基づいて不安定性を予測するが、圧縮性が導入する補正の大きさが、圧縮性ガス、非圧縮性液体、および剪断風を組み合わせた厳密な正規モードの枠組みの中で定量化されたことはなかった。

手法
著者は、二相自由境界問題に対して厳密な数学的手法を用いている：

1. 線形化と作用素の定式化： 非線形の二相オイラー系を、層流状態の周りで線形化する。移動する界面を扱うために、問題を固定された参照領域へと輸送し、Alinhacの「グッド・アンノウン（good unknown）」を用いた変換によって、Fr'echet微分を線形化系と特定する。これにより、閉じた、稠密に定義された線形作用素 $L_k$ を、$L^2$ 状態空間上で定義することが可能となる。この作用素の固有値は正規モードに対応する。
2. 正規モードへの還元： 解を $e^{ik(x-ct)}$ の形式で求める。空気側の方程式は、垂直方向の速度振幅 $w(z)$ に関する二階常微分方程式へと還元される。圧縮性かつ中立な場合、この方程式は発散形式 $(\bar{A} w')' = (\bar{\rho}_0 k^2 + G/V)w$ （ここで $V = U-c$）をとる。
3. 臨界層解析： 規則的な臨界レベル（$U(z_c) = c_r$）付近での解の挙動を扱う。著者は、実位相速度における解を、虚部 $\text{Im } c > 0$ を持つ解の極限として定義する「限定吸収原理」を確立する。これにより、臨界レベルにおける特異性が解消され、解の展開における対数項が明らかになる。
4. 留数公式と分散関係： 空気側の圧力応答の虚部 $Z_a^+(k, c_r)$ の境界値に関する留数公式を導出する。この量は、臨界層を横切る波運動量フラックスの跳躍（ジャンプ）に関連している。分散関係は、正則関数 $E(k, c, \epsilon) = 0$ （ここで $\epsilon$ は空気密度と水密度の比の小量）として定式化される。
5. 摂動論： ルーシェの定理を用いて、自由重力キャピラリー波速度 $c_k$ の近傍にある不安定根を特定する。

主要な貢献と結果

• 圧縮型マイルズ基準： 本論文は、圧縮性かつ中立成層の大気においてもマイルズ不安定性が存続することを確立している。しかし、不安定性の基準は変化する。非圧縮性の曲率項 $U''(z_c)$ は、以下の量に置き換わる：
  $$Q(z_c) = U''(z_c) + \frac{g}{a^2(z_c)} U'(z_c) = \bar{\rho}_0 \left( \frac{U'}{\bar{\rho}_0} \right)'$$
  ここで $g$ は重力加速度、$a$ は音速、$\bar{\rho}_0$ は正規化された基底状態の空気密度である。不安定条件は $Q(z_c) < 0$ （具体的には、すべての臨界レベルにおいて $Q \le 0$ であり、最も高いレベルにおいて $Q < 0$ であること）を要求する。この量 $Q$ は、単なる速度プロファイルの曲率ではなく、密度によって重み付けされた比渦度（specific vorticity）の垂直勾配（$U'/\bar{\rho}_0$）を表している。

• 留数公式とフラックスの跳躍： 圧力応答の虚部は、各臨界レベルからの寄与の和として示される：
  $$\text{Im } Z_a^+(k, c_r) = -\frac{\pi}{k^2} \sum_{j} \frac{\bar{\rho}_0(z_j) Q(z_j)}{|U'(z_j)|} |w(z_j)|^2$$
  この結果は、成長率が臨界層を横切る波誘起フォームストレス（エリアセン・パルム・フラックス）の跳躍と厳密に結びついていることを示しており、圧縮性の場合において運動量沈降メカニズムが比渦度の勾配によって制御されていることを裏付けている。

• 不安定根の存在と一意性： 単調な風速プロファイルと一つの規則的な臨界レベルを持ち、$Q(z_c) < 0$ である場合、著者は、密度比 $\epsilon$ が十分に小さいとき、 $c_k$ の近傍にただ一つの不安定根 $c_\epsilon$ が存在することを証明する。この根は、線形作用素 $L_k$ の孤立した代数的に単純な固有値である。成長率は $\gamma_\epsilon = k \text{Im } c_\epsilon > 0$ である。

• 低マッハ数極限と非圧縮理論の回復： 本論文は、非圧縮近似の精度を定量化している。マッハ数 $M \to 0$ となるにつれ、圧縮性圧力応答 $Z_a^+$ は非圧縮性の対応物へと $O(M^2)$ の速度で収束する。したがって、不安定根は $|c_\epsilon(M) - c_\epsilon(0)| \le C \epsilon M^2$ を満たす。これにより、風波が発生する海面のような低マッハ数において、非圧縮理論が正確であることが確認される。

• 強い圧縮性下での不安定性の消失： 著者は、$Q$ の符号が $U''$ の符号と異なる可能性があることを示している。強い圧縮性を想定した特定のモデル領域（空気の音速よりもはるかに低い音速を必要とする設定）では、$U'' < 0$（これは非圧縮極限では不安定性を意味する）でありながら $Q > 0$ となることが起こり得る。このような場合、マイルズ分岐は発生せず、$\epsilon$ が増加しても自由水面波から成長モードは現れない。これは、非圧縮性の符号基準が圧縮性大気において普遍的に適用できるわけではないことを示している。

意義と主張
本論文は、圧縮性ガス、非圧縮性液体、剪断風、およびマイルズ臨界層を組み合わせた、初の厳密な正規モード理論を提供すると主張している。これは、B"uhler, Shatah, Walsh, および Zeng による厳密な非圧縮理論を、圧縮性領域へと拡張するものである。

主な意義は以下の通りである：

1. 数学的厳密性： 線形作用素の構成を含む、圧縮性二相問題の完全なスペクトル解析を提供したこと。
2. 物理的洞察： 単純な速度曲率 $U''$ に代わり、特定の渦度勾配 $Q$ が、中立かつ圧縮性の大気における臨界層不安定性の支配パラメータであることを特定したこと。
3. 定量的検証： 低マッハ数における風波生成に対する非圧縮近似の妥当性を、$M^2$ 次の誤差項を伴う明示的な誤差評価とともに証明したこと。
4. 古典的理論の限界： 強力な圧縮性条件下では、古典的なマイルズ基準（$U'' < 0$）が不安定性を予測できない可能性があること（静水圧的な密度勾配が、速度プロファイルが凹状であっても流れを安定化させ得るため）を浮き彫りにしたこと。

本研究は、乱流、粘性、または非線形発展をモデル化することを目的としたものではなく、あくまで非粘性、線形な界面のスペクトル不安定性に焦点を当てている。内部重力波を排除するために、大気を中立成層として扱っており、臨界層メカニズムに対する圧縮性の影響を孤立させている。