Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

高分解能大渦シミュレーションは、成層均一乱流における受動スカラー混合が横方向には非成層の場合と類似した挙動を示す一方で、成層が大規模な攪拌を抑制することで鉛直混合を著しく阻害し、スカラーフラックスに対する特定のモデル化手法を導く明確な成長限界と変動強度をもたらすことを明らかにしている。

要約

インクの一雫がコップの水の中へ広がる様子を想像してください。水をそっとかき混ぜると、インクは均一に広がります。これが、科学者が流体の混合について通常考えていることです：乱流は巨大なスプーンのように働き、すべてを均一になるまでかき混ぜます。

しかし、その水が層状になっている場合はどうなるでしょうか？底の水は重くて塩分を含み、上の水は軽くて真水だと想像してください。これを成層と呼びます。現実世界では、海洋（深い水ほど密度が高い）や大気（高度が上がると空気が薄く軽くなる）でこの現象が起こります。

この論文は、高度なコンピュータ実験であり、**「この層状構造は、乱流中の流体（汚染物質や煙のような）の『染み』の広がり方をどのように変えるか？」**という問いに答えるものです。

彼らが発見したことを、簡単な概念に分解して物語りましょう。

設定：2 種類の「染み」

研究者たちは、最初はカオス的に渦巻いている（乱流状態の）仮想流体を作成しました。そして、2 種類の異なる「染み」（受動スカラー）を導入して、その挙動を観察しました。

1. 水平な染み：水に浮かぶ平らな紙のように、横方向に広がったインクの層。
2. 垂直な染み：色の垂直な壁のように、上下方向に広がったインクの層。

彼らは 2 つのシミュレーションを実行しました。1 つは通常の水中（非成層）、もう 1 つは「層状」の水（成層）中です。

大きな発見：「上下」対「左右」の違い

1. 左右の染み（横断層）

何が起こったか：染みが横方向に広がっていた場合、層状構造は広がりを妨げませんでした。実際、成層水の中では、通常の水よりもわずかに速く広がりました。
比喩：廊下を走る人々の群れを想像してください。床が完全に平ら（非成層）であれば、彼らはあらゆる方向に走ります。床に目に見えないわずかな傾斜（成層）があっても、彼らは横方向には問題なく、むしろ少しエネルギーを込めて走り続けます。「インク」はどちらの場合も広範囲に広がります。
注意点：全体的な広がりは似ていましたが、成層水の中の「インク」はより「棘状」でした。滑らかな勾配の代わりに、より鋭く、ギザギザした縁を持っていました。より「断続的」であり、純粋なインクのポケットと純粋な水のポケットがあり、その中間の滑らかな領域は少なかったのです。

2. 上下の染み（垂直層）

何が起こったか：ここで魔法（そして制限）が起きました。通常の水中では、垂直な染みは横方向のものと同じように、上下に容易に広がりました。しかし、成層水中では、広がりがほぼ完全に停止しました。
比喩：スプーンで濃厚なミルクシェイクをかき混ぜようとするのを想像してください。スプーンを上下に動かそうとすると、シェイクの層が抵抗します。重いものは底に残り、軽いものは上に残ろうとするからです。「かき混ぜる」動きは押しつぶされてしまいます。
結果：垂直な染みはごく初期にわずかに成長しましたが、その後「天井」にぶつかりました。流体の安定した層が蓋のように作用して、乱流が上下に混合するのを防いだため、これ以上広がることができませんでした。流体は横方向に渦を巻くことはできましたが、上下に混合することはできませんでした。

なぜこれが重要なのか？（「何」の背後にある「なぜ」）

研究者たちは、垂直方向において流体がバネのように振る舞うことを発見しました。乱流が重い層を押し上げたり、軽い層を押し下げたりしようとすると、重力がそれを引き戻します。これが「かき混ぜる」動きを止めてしまいます。

しかし、流体は横方向にはまだ渦を巻くことができます。したがって、乱流の「垂直長さ」は特定のサイズに固定され（重力と流体の層の強さによって決定されます）、染みはそのサイズを超えて成長することはできません。

予測のための「レシピ」

この論文は、スーパーコンピュータを使わずにこれらの染みの広がりを予測するための簡単な数学的「レシピ」の作成も試みました。

• 染みの形状が分かっている場合：横方向の広がりの速さを予測するために、単純な 1 数式の公式を使用できます。これは非常にうまく機能します。
• 染みの形状が分からない場合：形状を推測する必要があります（滑らかな曲線のように見えると仮定します）。この場合、2 数式の公式が必要になります。これは、染みが渦巻く流体とリズムを合わせて落ち着く「後」には非常にうまく機能します。

結論

• 横方向の混合：安定した層（深海や成層圏など）は横方向の混合を止めません。むしろ、それを少し激しく、ギザギザしたものにするかもしれません。
• 垂直方向の混合：安定した層はブレーキのように働きます。流体が上下に混合するのをほぼ完全に防ぎます。
• 「かき混ぜる」対「混合」の区別：流体は横方向にまだ「かき混ぜる」（動き回る）ことができますが、層が入れ替わるのを抵抗するため、垂直方向に「混合」（混ぜ合わせる）ことはできません。

著者は、彼らの実験が特定の流体特性（プラントル数 0.7）を使用していたことに言及しています。もし流体が「厚い」場合や異なる特性（プラントル数 > 1）を持っていた場合、混合が自身の浮力を生み出す「逆」効果のために結果が変わる可能性があると警告しています。しかし、彼らがテストした条件においては、「横方向は自由、上下は遮断」というルールが真実です。

技術概要

技術的サマリー：成層および非成層均等乱流における受動スカラー混合層の進化

問題提起
安定成層環境（対流圏上部、成層圏、海洋混合層など）における粒子およびスカラーの拡散は、大気および海洋科学における重要な課題である。安定成層が密度場および乱流の減衰に及ぼす影響はよく文書化されているが、受動スカラー（エアロゾル、汚染物質、マイクロプラスチックなど）の拡散については依然として理解が浅い。具体的には、安定密度勾配の影響下で減衰する均等等方乱流（HIT）に導入された受動スカラー混合層（PSML）がどのように進化するかは不明である。中心的な問いは、浮力スケールと異方性の導入が、非成層 HIT において観測されるスカラー統計の自己相似進化をどのように変化させるかである。

手法
著者らは、非成層および安定成層均等乱流の減衰を調査するために、高解像度大渦シミュレーション（LES）を採用した。シミュレーションには、非静水圧ボウシネスク近似を用いた非発散ナビエ・ストークス方程式が使用された。

• 流れの構成: 2 つのシミュレーションが実施された。最初のものは非成層 HIT である。2 つ目は初期条件は同一であるが、一様な安定密度勾配（$g \neq 0$）を含んでおり、これにより流れは約 1 浮力周期後に強成層乱流へ遷移する。
• 受動スカラー: 2 つの受動スカラー混合層が初期化された。1 つは鉛直方向に平均勾配（$\xi_z$）を持ち、もう 1 つは横方向に平均勾配（$\xi_y$）を持つ。これらは混合分率（$0 \le \xi \le 1$）としてモデル化された。
• 数値的詳細: 自己相似からの逸脱を引き起こす拡散効果を最小化するため、シミュレーションでは等しい係数を持つ超粘性および超拡散が使用された。分子プラントル数は $Pr = 0.7$に設定された。レイノルズ数は高く（非成層の場合$Re_h > 15,000$および$Re_\lambda > 450$）、慣性対流領域および慣性拡散領域の存在が保証された。解像度は浮力長さスケール（$L_b/\Delta \gg 1$）を解像するのに十分であった。
• 解析: 本研究では、スカラープロファイル幅、分散、高次モーメント（歪度、尖度）、フラックス、混合効率の進化を追跡した。また、2 つの条件下での横方向スカラーフラックスに対する渦拡散モデルを評価した。1 つは平均スカラープロファイルが既知の場合、もう 1 つは仮定しなければならない場合である。

主要な結果

1. 横方向混合（$\xi_y$）: 重力に垂直な方向において、成層ケースにおける受動スカラーの進化は、非成層ケースと概似しているが、わずかに速く拡がる。

   • スカラー変動は成層ケースでより激しく、ピーク分散は約 15% 高い。
   • 乱流/非乱流界面は、成層ケースでより断続的である。
   • スカラー場が速度場と準平衡に達すると、成層および非成層の両方のケースにおいて正規化されたフラックスプロファイルが収束し、成層にもかかわらず同様の輸送メカニズムが働いていることを示している。

2. 鉛直混合（$\xi_z$）: 成層は鉛直混合を劇的に抑制する。

   • 鉛直層幅（$\delta_z$）は初期に成長するが、約 1 浮力周期（$t^* \approx 3$）後に急速に拡大を停止する。
   • この成長の停止は、浮力力による大規模な鉛直攪拌の抑制に起因する。
   • 層の最終幅は、鉛直フルード数がオーダー 1 に制限される（$Fr_v \sim O(1)$）ように制約された水平速度の鉛直積分長さスケール（$L_v$）に比例する。
   • 平均プロファイルは分子拡散により単調に成長するが、フラックス解析は、内部波のような挙動と一致し、攪拌が時折スカラーのひも状構造を元の側へ掃き戻す「逆フラックス」現象を明らかにする。

3. 混合効率: 成層ケースにおける累積混合（スカラー分散の散逸）は、後期において最終的に非成層ケースを上回る。これは、混合が攪拌に追随し、成層流れにおける水平攪拌率が水平長さスケールに比例するためであり、そのスケールは非成層ケースよりも大きくなることに起因する。

4. フラックスのモデル化:

   • 既知のプロファイル: 平均スカラープロファイルが既知の場合、単一定数渦拡散モデル（$F_y = -c_1 \delta_y v' \partial \xi / \partial y$）は、$c_1 \approx 0.4$ として、成層および非成層の両方のケースに対してフラックスを効果的に予測する。
   • 仮定されたプロファイル: プロファイル形状（誤差関数と仮定）およびスカラー長さスケールを速度統計から推定する必要がある場合、2 定数モデルが必要となる。このモデルは、スカラーが速度場と準平衡にある後期において良好に機能し、スカラー長さスケールを運動エネルギー散逸率からスケーリングすることを可能にする。

意義と主張
本論文は、減衰する成層乱流における受動スカラー混合の詳細な特徴付けを提供し、異方性拡散の理論的予測と数値的観測の間のギャップを埋めると主張している。

• 異方性: 本研究は、成層が混合において明確な二極化を生み出すことを確認している。すなわち、水平混合は活発であり、むしろ強化される一方、流れが強成層となると鉛直混合は実質的に抑制される。
• モデル化の実用性: 著者らは、スカラーが速度場と準平衡にある場合、標準的な渦拡散アプローチが成層流れに対しても有効であることを実証している。彼らは、プロファイル形状を仮定するよりも既知のプロファイルを仮定することが、モデル化の要件を大幅に単純化することを強調している。
• 限界と将来の作業: 著者らは謙虚に、彼らの結果はプラントル数 0.7 に特化したものであると指摘している。彼らは明示的に、最近の文献はより高いプラントル数（$Pr > 1$）が小スケールで逆浮力フラックスを誘発し、それが流れの減衰率、ひいては受動スカラー混合率を変化させる可能性があると示唆していると述べている。したがって、現在の知見は、さらなる調査なしには高プラントル数シナリオに直接適用されない。

本研究は、成層が混合層の鉛直ダイナミクスを根本的に変化させる一方で、横方向混合は断続性と変動強度の増加を伴うものの、非成層乱流の多くの特性を保持することを結論付けている。