Direct Lagrangian tracking simulation of droplet growth in vertically-developing turbulent cloud

この論文は、地上から雲頂までを縦長に設定した直接数値シミュレーションを用いた新しいラグランジュ粒子追跡モデルを開発し、乱流が雲滴の衝突・合体成長を促進して降水の発生を早め、雨滴のサイズを大きくすることを明らかにしたものである。

要約

この論文は、**「雲の中で雨粒がどうやって大きくなり、地面に降り注ぐのか」**という謎を、スーパーコンピュータを使って解き明かした研究です。

特に注目したのは、「風の乱れ（乱流）」が雨粒の成長にどんな影響を与えるかという点です。

難しい数式や専門用語を使わず、身近な例え話で解説しますね。

---

🌧️ 研究の舞台：巨大な「雲の縦断模型」

まず、研究者たちは雲をシミュレーションするための「実験室」を作りました。
これまでの研究では、箱の中をぐるぐる回るような「周期的な箱」を使っていましたが、これでは雲の「上から下までの高さ」や「地面からの距離」を正しく再現できませんでした。

そこで今回、**「地面から雲のてっぺんまで続く、細長い縦長のトンネル」**のような実験室を作りました。
これにより、雲の「下層（入り口付近）」「中層（真ん中）」「上層（天井付近）」で、それぞれ何が起きているかを詳しく観察できるようになりました。

🌪️ 2 つのシミュレーション：静かな風 vs 激しい乱れ

この実験室で、2 つのパターンを比較しました。

1. ラミナー（LAM）ケース：静かな風
   • 風はまっすぐ上に向かって吹くだけ。乱れはありません。
   • 雨粒はのんびりと成長します。
2. タービュランス（TURB）ケース：激しい乱れ
   • 風が激しく揺れ動き、渦（うず）が生まれます。
   • これが**「乱流」**です。

🏃‍♂️🏃‍♂️ 雨粒の成長ストーリー：3 つの段階

この 2 つのパターンで、雨粒（水滴）がどう成長するかを詳しく見てみましょう。

1. 赤ちゃん雨粒の誕生（雲の底）

まず、水蒸気が集まって小さな雨粒（雲の粒）が生まれます。この段階では、乱流があってもなくても、あまり違いはありません。

2. 中層での「お友達集め」競争（雲の真ん中）

ここが最大のポイントです！

• 静かな風の場合： 雨粒たちはまっすぐ上か下へゆっくり移動します。同じ大きさの雨粒同士が出会うチャンスは少ないです。
• 乱流の場合： 風が激しく揺れるので、同じくらいの大きさの雨粒同士が、まるでダンスのように激しくぶつかり合います。
  • これを**「オートコンバージョン（同サイズ同士の合体）」**と呼びます。
  • 乱流があるおかげで、雨粒たちは「お友達」を見つけやすくなり、早く大きな粒に成長し始めます。

3. 下層での「大物狩り」（雲の底へ降りる時）

中層で大きくなった雨粒は、重力に引かれて下へ落ちていきます。

• 乱流の場合： 中層で大きくなった「大物」が、下層へ降りてきます。そして、下層にまだ小さく残っている雨粒たちを**「集め（アクリション）」**て、さらに巨大化します。
• 乱流のおかげで、この「大物狩り」がスムーズに行われ、雨粒は驚くほど早く大きくなります。

🌧️ 結果：雨は早く、そして大きく降った！

このシミュレーションから、驚くべき結果が得られました。

• 雨の降り始めが早かった：
  乱流がある場合、地面に雨粒が到達する時間が、静かな風の場合より約 270 秒（4 分半）も早かったのです。
• 最初の雨粒がデカかった：
  最初に地面に届いた雨粒のサイズは、乱流がある場合の方が1.5 倍も大きかったのです。

🧩 隠された秘密：「成長の歴史」の違い

なぜ、同じ水蒸気から生まれたのに、乱流がある方が雨粒が大きくなるのでしょうか？

研究者は、**「雨粒の成長履歴」を詳しく追跡しました。
乱流があるおかげで、雨粒同士が早くぶつかり合い、数が減って「大物」になりました。
すると、「残った雨粒一人一人が、より多くの水蒸気を独占して成長できる」**という状態になりました。

まるで、**「お菓子を分け合う子供たち」**のような話です。

• 静かな風： 子供たちがバラバラで、お菓子のかけら（水蒸気）を少ししかもらえない。
• 乱流： 子供たちがすぐに集まって大きなお菓子（大きな雨粒）を作ってしまう。すると、その大きなお菓子は、周りの小さな子供たちからさらに多くのお菓子を奪い取って、さらに巨大化する。

つまり、**「乱流は、雨粒同士の衝突を促し、それが結果として『水蒸気の取り合い』を変え、雨粒をより大きく成長させる」**という、意外な連鎖反応が起きていることがわかりました。

💡 まとめ

この研究は、**「風の揺らぎ（乱流）は、雲の中で雨粒を『急成長』させる魔法のスイッチ」**であることを証明しました。

• 乱流があると、雨粒同士が早く出会う。
• それで大きな雨粒が早くできる。
• 結果として、雨は早く降り始め、最初の雨粒はより大きくなる。

これは、気象予報の精度を上げたり、人工降雨の仕組みを理解したりする上で、非常に重要な発見です。雲という複雑な世界を、スーパーコンピュータという「時間と空間を操る魔法の鏡」を使って、鮮明に映し出した素晴らしい研究だと言えます。

技術概要

以下は、Masaya Iwashima および Ryo Onishi による論文「Direct Lagrangian tracking simulation of droplet growth in vertically-developing turbulent cloud（垂直に発達する乱流雲における液滴成長の直接ラグランジュ追跡シミュレーション）」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

雲の微物理過程（雲水から雨滴への成長）を数値シミュレーションする際、従来の手法には以下のような課題がありました。

• 周期的境界条件の限界: 従来の直接数値シミュレーション（DNS）とラグランジュ粒子追跡を組み合わせた研究の多くは、周期的境界条件を持つ立方体ドメインを使用していました。これは計算上便利ですが、実際の雲が持つ「高度による鉛直構造（過飽和度の勾配や液滴粒径分布の変化など）」を表現できず、現実的な雲の発達過程を再現するのには不十分でした。
• 乱流の影響の解明不足: 乱流が雲微物理過程（特に液滴の衝突・合体）に与える影響は理論的に指摘されてきましたが、鉛直構造を有する現実的な雲の発達過程において、乱流がどのように液滴成長や降水開始に影響するかを、すべての微物理過程（核生成、凝結、衝突合体、沈降）を含めて包括的に解析した研究は不足していました。

2. 手法 (Methodology)

本研究は、直接数値シミュレーション（DNS）とラグランジュ粒子追跡法を組み合わせた新しい明示的雲微物理モデルを開発し、以下の構成でシミュレーションを行いました。

• 計算ドメイン: 地面から雲頂までをカバーする「垂直に伸びた準 1 次元ドメイン」を採用。これにより、実際の雲の鉛直構造を明示的に捉えることを可能にしました。
• 乱流の導入: 均一等方性乱流（HIT）の瞬間的なスナップショットを、この垂直ドメインに周期的に積み重ねる形で導入しました。これにより、乱流風変動を明示的に解像しました。
• 微物理過程のモデル化:
  • CCN 活性化: 確率モデル（Twomey モデル）を使用。
  • 凝結/蒸発成長: 周囲の過飽和度と液滴サイズに依存する半径成長方程式。
  • 衝突・合体成長: 幾何学的な軌道計算とストークス流れの解を用いた流体力学的相互作用の考慮。合体効率を 1 と仮定（分裂過程は考慮せず、本研究の液滴サイズ範囲では無視可能）。
  • 沈降: 重力による沈降を考慮。
• シミュレーションケース:
  1. LAM-case（層流ケース）: 乱流変動なし、上昇気流のみを prescribed（規定）した流れ場。
  2. TURB-case（乱流ケース）: 上昇気流に加え、HIT の乱流速度場を合成した流れ場。
• 初期・境界条件: 海洋性条件を想定した CCN 分布、KiD warm-1 ケースに基づく温度・水蒸気プロファイル、および上昇気流の時間変化（最初の 600 秒で上昇後、停止）。

3. 主な貢献 (Key Contributions)

• 新しい計算枠組みの確立: 周期的境界条件の立方体ドメインではなく、地面から雲頂までの垂直構造を再現する「準 1 次元ドメイン」に乱流を導入した、世界初の DNS 結合ラグランジュ追跡シミュレーション手法を提案しました。
• すべての暖かい雲の微物理過程の統合: CCN 活性化から降水開始までの全過程を、乱流の影響下で直接追跡可能にしました。
• 乱流と微物理過程の相互作用の定量的評価: 従来の DNS 研究では得られなかった、「高度と時間依存性を持つ微物理統計量」を初めて取得し、乱流が降水形成に与える具体的なメカニズムを解明しました。

4. 結果 (Results)

乱流ケース（TURB-case）と層流ケース（LAM-case）の比較から、以下の重要な知見が得られました。

• 衝突頻度の増大: 乱流の存在により、液滴間の衝突頻度が大幅に増加しました。特に、上昇気流が存在する初期段階（$t < 600$秒）において、雲の中層で衝突頻度が 5000 $m^{-3}s^{-1}$ を超え、層流ケース（1000 $m^{-3}s^{-1}$未満）と比較して顕著に高まりました。
• 成長メカニズムの時間的・空間的変化:
  • 初期段階（自動変換）: 乱流は、同程度のサイズの液滴同士の衝突（自動変換：autoconversion）を促進しました。これにより、中層で粒径分布の広がり（100 $\mu m$ 以上の液滴の出現）が早期（約 400-450 秒）に発生しました。
  • 後期段階（吸着成長）: 初期に中層で成長した大きな液滴が重力で落下し、下層でより小さな液滴を捕獲する吸着（accretion）が活発化しました。これにより、下層での平均粒径が急激に増大しました。
• 降水開始の早期化と雨滴サイズの増大:
  • 地表への降水開始時刻は、乱流ケースの方が層流ケースより約 270 秒早まりました。
  • 地表に到達する最初の雨滴の直径は、乱流ケースの方が層流ケースより約 1.5 倍（最初の 10 滴の平均）大きく、全体的にも 1.2 倍程度大きくなりました。
• 凝結成長への間接的影響: 乱流による衝突合体の促進が液滴数密度を減少させ、結果として単位液滴あたりの凝結水量が増加しました。これにより、衝突合体だけでなく、凝結成長そのものも乱流の影響を受けたことが示唆されました。

5. 意義 (Significance)

本研究は、雲微物理と乱流の相互作用を理解する上で重要な進展をもたらしました。

• 降水メカニズムの解明: 乱流が「自動変換（同サイズ液滴の衝突）」を促進し、それが「吸着成長」を加速させることで、降水開始を早期化し、初期雨滴を大型化させるという具体的なメカニズムを、鉛直構造を考慮した現実的なシナリオで実証しました。
• 数値モデルの高度化: 従来の周期的ドメインの限界を克服し、実際の雲の鉛直構造を考慮した高精度な微物理モデルの構築可能性を示しました。
• 気象・気候モデルへの応用: 乱流効果をより正確に表現できる微物理パラメータ化手法の開発や、降水予測精度の向上に寄与する基礎データを提供しました。

総じて、本研究は「垂直に発達する乱流雲」における液滴成長のダイナミクスを、ラグランジュ追跡法を用いて初めて詳細に可視化・定量化した画期的な研究です。