Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

本研究は、線形理論と直接数値シミュレーションを組み合わせることで、初期の鉛直シアーを持たない強層状流において、水平シアー不安定性が、鉛直シアーの直接的な発生、あるいは柱状渦への非線形進化のいずれかという2つの異なる経路を通じて、必然的にダイオピクニカル混合（密度成層に沿った混合）を駆動することを明らかにし、その両者が最終的には小規模なケルビン・ヘルムホルツ不安定性を誘発するものの、異なる鉛直スケールを励起することによって混合効率が異なることを示している。

要約

海洋や大気を、巨大な層状のケーキだと想像してみてください。その層は、異なる密度の流体（ケーキの異なるフレーバーのようなもの）でできており、互いに混ざり合うことを嫌います。通常、科学者たちは、これらの層を横方向と、同時に上下方向（垂直方向のシアー）の両方に押し出したときに何が起こるかを研究します（垂直シアー）。しかし、この論文は異なる問いを投げかけています。もし、非常に強く層状化された環境において、上下の動きを伴わずに、横方向だけに押し出したとしたら、何が起こるのでしょうか？

研究者たちは、たとえ最初は完全に水平な流れであったとしても、自然界には、混沌（乱流）を生み出し、層を混合させるための2つの明確な「レシピ」または経路が存在することを発見しました。自然がどちらのレシピを選ぶかは、実験にどのような「種」をまくか、つまり、最初に流体にどのような小さな刺激を与えるかに完全にかかっています。

以下に、比喩を用いた2つの経路の解説を記します。

設定：穏やかな川

幅の広い、穏やかな川が水平に流れている様子を想像してください。水はパンケーキを積み重ねたようになっています（強い成層）。最初、流れは滑らかで二次元的です（つまり、左右にしか動きません。上下には動きません）。

経路1：「混み合った部屋」効果（ダイレクト・ルート）

始まり方： 川のあらゆる場所に、一度に小さなランダムな刺激を与えます（空中に紙吹雪を投げ入れるようなものです）。
何が起こるのか：

1. さざ波： 層があるため、流体は単に左右に波打つだけでなく、即座にさまざまな大きさで上下にも波打ち始めます。これは、部屋の中にいる大勢の人々が皆同時に動こうとして、混沌とした多方向の揺れを生み出しているようなものです。
2. シアー： これらのさざ波が、強力な垂直方向の流れ（シアー）を作り出します。
3. 崩壊： これらの垂直方向の流れが非常に強くなると、それらは「スナップ」し、小さな激しい渦（小さな渦巻きのようなもの）を生み出します。これは「ケルビン・ヘルムホルツ不安定性」と呼ばれ、風が水面を吹き抜けるときにできる砕け波のように見えます。
   結果： 混合は効率的に行われます。エネルギーがさまざまなサイズのさざ波に分散されているため、「摩擦」（粘性散逸）が低くなり、混合プロセスが比較的効率的になります。

経路2：「シンクロナイズド・ダンス」効果（インダイレクト・ルート）

始まり方： 非常に特定の、組織化された刺激を与えます（指揮者がタクトを振って、全員に特定のパターンで動くよう指示するようなものです）。
何が起こるのか：

1. 渦： 混沌としたさざ波の代わりに、流体は長い垂直方向の柱状の渦へと自己組織化します（川の中に立つ巨大な竜巻のようなものです）。長い間、流れは完璧に二次元的なまま、これらの大きな渦の柱として存在し続けます。
2. 揺らぎ： やがて、これらの巨大な柱は不安定になります。これらは非常に特定の高周波の仕方で揺れ始めます。研究者たちはこれを「双曲型不安定性」と呼んでいます。回転している独楽（こま）が、倒れる直前に激しく揺れ始める様子を想像してください。
3. 崩壊： この激しい揺れが、非常に薄く鋭い垂直方向のシアーの層を作り出します。そして、これらの薄い層が、経路1と同様に、小さな激しい渦へとスナップします。
   結果： 混合は行われますが、効率は低くなります。なぜでしょうか？ この経路は、極めて薄く鋭い層を作り出すからです。これらの小さく鋭い層を作り出し、壊すためには、多くのエネルギー（摩擦）を必要とします。それは、鈍いナイフで厚いチーズの塊を切ろうとする場合（経路1）と、カミソリの刃を使う場合（経路2）の違いのようなものです。カミソリの刃は、より鋭く、よりエネルギーを消費する切り込みを作ります。

大きな教訓

この論文は、垂直シアー（上下の動き）が、最初から存在している必要はないということを証明しています。流体が十分に厚い（高いレイノルズ数を持つ）場合、水平シアーから垂直シアーは避けられない副産物として生じるのです。

• ランダムなノイズから始めた場合： 経路1（ダイレクトで効率的な混合）になります。
• 特定のパターンから始めた場合： 経路2（インダイレクトで、効率の低い混合）になります。

研究者たちは強力なコンピュータ・シミュレーションを用いて、これら2つの経路が実在する明確に異なる経路であること、そして最初に選ぶ「レシピ」によって、どれだけのエネルギーが熱として浪費され、どれだけのエネルギーが実際に層を混合するために使われるかが決まることを示しました。

要約すると： 完全に静かで層状化された流体であっても、水平方向の押しは、最終的に垂直方向の混沌を生み出します。しかし、その押し方をどう始めるかによって、その混沌の見え方と、層を混合する効率は変わるのです。

技術概要

問題提起
大規模で強層状化された水平流の三次元（3D）摂動に対する安定性は、地球物理学的および天体物理学的システムにおけるエネルギーの下方への転送およびダイピクナル輸送の重要なメカニズムである。垂直シアーを持つ層状流はよく研究されているが、水平シアーを持つもの（初期の垂直シアーを持たないもの）については、これまであまり注意が払われてこなかった。中心的な問いは、当初二次元（2D）であり、かつ強層状化された（$Fr \ll 1$）水平シアー流が、どのようにして乱流へと遷移し、垂直シアーを生成するのかという点である。先行文献では、2つの異なる経路が示唆されている。一つは、一次モードの直接的な成長によってケルビン・ヘルムホルツ（KH）不安定性を引き起こす経路であり、もう一つは、柱状渦（columnar vortices）への非線形進化を経て二次不安定へと至る経路である。しかし、これらの経路が、初期条件の違い（例：ハイパボリックタンジェント型か正弦波型か）や問題設定の違い（初期値問題か強制駆動型か）による排他的なものなのか、あるいは同一の物理系において利用可能な代替ルートなのかについては、不明なままであった。さらに、選択された経路がダイピクナル混合効率に与える影響についても、系統的な定量化が行われていない。

手法
著者らは、線形安定性理論と直接数値シ法（DNS）を組み合わせ、三周期領域における垂直不変な水平平面平行コルモゴロフ流（$\bar{u}(y) = \sin(y)\mathbf{e}_x$）を対象として、この遷移を調査している。流体は定数浮力振動数 $N$ を持つ線形層状流である。シミュレーションは、高いレイノルズ数（$Re = 10,000$）およびペクレ数（$Pe = 10,000$）、低いフルード数（$Fr = 0.15$および$0.1$）で行われ、その結果として得られる浮力レイノルズ数（$Re_b = Fr^2 Re$）は 225 および 100 である。

初期条件が経路選択に与える影響を分離するため、以下の2組のDNSが行われた：

1. DNS1 (経路 P1): グリッドスケールのホワイトノイズで初期化され、一次水平シアー不安定性のすべての不安定モードの成長を許容する。
2. DNS2 (経路 P2): グリッドスケールのノイズに加え、一次不安定性の最も成長の速い2D ($k_z=0$) 固有モードに対する特定の「シード」を用いて初期化される。

さらに、著者らは、DNS2において出現する柱状流の「凍結された（frozen-in）」解析モデルに対して線形安定性解析を行っている。このモデルは、非線形に進化した2Dの流れを、テイラー・グリーン渦配列と正弦波流の重ね合わせとして近似しており、これにより固有モードと二次3D摂動の成長率の計算を可能にしている。

主要な貢献と結果

1. 単一の設定における2つの経路の確認：
   本研究は、初期摂動スペクトルのみによって区別される両方の経路が、全く同じモデル設定内で起こり得ることを示している：

• 経路 P1 (直接的3D成長): ホワイトノイスで初期化した場合、一次水平シアー不安定性の $k_z \neq 0$ 固有モードが直接成長する。これらのモードは広範なスケールで垂直シアーを生成し、その後、$Re_b$ が十分に大きい場合には二次的な小スケールのKH不安定性に不安定となる。この経路は、強制駆動型のコルモゴロフ流における観測結果と一致する。
• 経路 P2 (柱状進化): $k_z=0$ モードをシードした場合、それが流れを支配し、背景流を長寿命な時間依存的2D柱状（渦状）流へと駆動する。この2D状態は、その後に二次的な3Dモードに対して不安定となる。著者らは、これらの二次モードを、渦配列内の双曲停留点から生じる**双曲不安定性（hyperbolic instabilities）**であると特定した。これらの不安定性は高い垂直波数（$k_z \sim Fr^{-1}$）を持ち、垂直シアーの出現を駆動し、それが最終的に三次KH不安定性を引き起こす。

2. 経路 P2における二次不安定の特性評価：
   凍結された柱状流の線形安定性解析を通じて、著者らは、この二次不安定が高波数 $k_z$ の双曲不安定であることを示した。この不安定の成長率は、双曲点の局所的な歪み速度（$\Delta_h$）に比例し、大きな振幅においては層状化の強さにほとんど依存せず、$\lambda \approx 0.46 \Delta_h$ となる。この不安定性は、双曲点を中心とした垂直速度場の「リップル（さざ波）」として現れ、孤立した渦対に見られる伝統的な「ジグザグ」不安定とは明確に異なる。

3. 混合効率への影響：
   本研究は、ダイピクナル混合と混合効率（$\eta$）を両経路について定量化している。

• 浮力散逸 ($\chi$): 浮力分散のピーク散逸率は、両経路で同程度である。これは、両者の場合とも小スケールの運動（$k_z \gtrsim Fr^{-1}$）によって制御されるためである。
• 粘性散逸 ($\epsilon$): 経路 P2 は経路 P1 よりも有意に高い粘性散逸を示す。これは、P2 における双曲不安定性が狭い高 $k_z$ スペクトルを励起し、薄い垂直シアー層における激しい散逸を引き起こすためである。対照的に、P1 は広範な垂直スケールのスペクトルを励起し、散逸をより広く分散させる。
• 混合効率: その結果、ピーク混合効率は経路 P2（$\eta_{max} \approx 0.165$）よりも経路 P1（$\eta_{max} \approx 0.25$）の方が有意に高い。この差は、P2 が混合される浮力の量に対して、より多くの運動エネルギーを散逸させることに起因する。

意義と主張
本論文は、十分な $Re_b$ があれば、強層状化された流れにおいて、垂直シアーの出現が小スケールのKH不安定性を駆動することは、変曲点を持つ場合に水平シアー不安定性の不可避な副産物であると主張している。経路の選択は初期条件、具体的には $k_z=0$ モードの初期振幅と $k_z \neq 0$ モードの比率によって決定される。

著者らは、過去の文献（特にハイパボリックタンジェント型のせん断層の研究）における経路 P2 の普及は、計算を加速するために $k_z=0$ のシードをしばしば追加する初期化慣行によるアーティファクトである可能性を示唆している。これにより、3Dモードの直接的な成長（経路 P1）が抑制されているのである。逆に、一部の文脈で P1 が見られないのは、層状化がより弱い限界において、不安定な $k_z \neq 0$ モードの範囲が制限されているためである可能性がある。

最終的に、本研究は、水平シアー不安定性が、既存の垂直シアーが存在しない場合であっても、強層状化された流れにおける乱流生成とダイピクナル混合の鍵となるプレイヤーであることを確立した。これは、乱流への経路の選択（直接的な3D成長か、2Dから3Dへの遷移か）が、エネルギー散逸スペクトル、ひいては混合の効率を根本的に変えることを強調しており、気候・地球物理モデルにおける乱流のパラメータ化において考慮すべき要因である。著者らは、彼らの結果は理想化された流れにおいては堅牢であるが、現実の地球物理学的設定における回転や周囲の波動場との相互作用は、依然として未解決の課題であると述べている。