液滴分裂モデル
理論と物理
概要
先生、液滴分裂モデルって何に使うんですか？
燃料噴射、スプレー塗装、消火器噴霧など、液滴が空気力により分裂（二次分裂, secondary breakup）する過程のモデルだ。噴射直後の大きな液滴がより細かい液滴に分裂していく過程を記述する。
液滴の分裂にはどんなメカニズムがあるんですか？
Weber数 $We$ によって分裂モード（レジーム）が変わる。
| レジーム | Weber数範囲 | 特徴 |
|---|---|---|
| Vibrational | $We < 12$ | 振動のみ、分裂なし |
| Bag breakup | $12 < We < 50$ | 薄膜状に膨らんで破裂 |
| Multimode | $50 < We < 100$ | Bag + Stripping |
| Sheet stripping | $100 < We < 350$ | 表面から薄膜が剥離 |
| Catastrophic | $We > 350$ | 爆発的分裂 |
Ohnesorge数 $Oh = \mu_d / \sqrt{\rho_d \sigma d}$ も重要で、粘性が高いと分裂が遅延する。
代表的なBreakupモデル
CFDで使われるモデルを教えてください。
| モデル | 概要 | 適用範囲 |
|---|---|---|
| TAB (Taylor Analogy Breakup) | 液滴をばね-質量-ダンパ系に類推 | $We < 100$、低速噴霧 |
| KHRT (Kelvin-Helmholtz / Rayleigh-Taylor) | KH不安定性 + RT不安定性の競合 | 高速ディーゼル噴射 |
| SSD (Stochastic Secondary Droplet) | 確率論的にサイズ分布を生成 | 汎用 |
| ETAB (Enhanced TAB) | TABの改良、分裂後の子液滴分布を改善 | 中速噴霧 |
TABモデルは液滴の変形を強制振動方程式で記述する。
$y$ は液滴変形量の無次元パラメータで、$y = 1$ で分裂が発生する。$C_F$, $C_b$, $C_k$, $C_d$ はO'Rourke & Amsden（1987）の定数だ。
KHRTモデルはどういう考え方ですか？
Kelvin-Helmholtz不安定性（液滴表面の波の成長）とRayleigh-Taylor不安定性（加速度による界面不安定性）を競合させる。高速噴射（ディーゼルエンジン）ではKH不安定性が支配的で、減速領域ではRT不安定性が重要になる。
Weber数が支配する——液滴は「いつ」割れるのか
液滴の分裂を支配する無次元数がWeber数 We = ρ_g u_rel^2 d / σ です。Weが12以下なら表面張力が復元力として働き液滴は球形を保ちますが、Weが100を超えると「壊滅的分裂（Catastrophic Breakup）」が起き、液滴は瞬時に微細なミストに分散します。この臨界Weは1940年代にHinzeが実験的に求めた値とほぼ一致しており、75年後の今もCFD液滴分裂モデルのベンチマーク基準として使われています。エンジン燃料噴射の設計では噴霧液滴径の予測精度が燃費・排出ガスに直接影響するため、分裂モデルの選択は経営的インパクトを持つ技術判断です。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
液滴分裂モデルはCFDにどう組み込まれるんですか？
Lagrangian粒子追跡法（DPM）の中で、各計算粒子（parcel）について毎タイムステップ、分裂条件を評価する。分裂が発生したら子液滴のサイズ・速度・数を計算して新しいparcelを生成する。
TABモデルの実装
TABモデルでは液滴ごとに変形量 $y$ と変形速度 $\dot{y}$ を追跡する。$y \geq 1$ になった時点で分裂が発生し、子液滴径はエネルギー保存から決まる。
子液滴は何個できるんですか？
質量保存から子液滴の数が決まる。実際にはparcel概念を使うので、parcel内の液滴数が更新され、各parcelの代表液滴径が変わる形で実装される。
KHRTモデルの実装
KHRTモデルではKH不安定性による表面波の成長率 $\Omega$ と波長 $\Lambda$ をDispersion relationから求める。
ここで $T = Oh \sqrt{We}$ はTaylor数だ。KH分裂で生成される子液滴径は $r_{child} = B_0 \Lambda$ で、$B_0 = 0.61$ が標準値だ。
RT不安定性は液滴の減速度 $a_{decel}$ に依存し、最速成長波数から子液滴径を決める。KH分裂とRT分裂は競合し、先に条件を満たした方が適用される。
Fluent・OpenFOAMでの設定
実際のソフトウェアではどう設定しますか？
FluentのWaveモデルはKH部分のみで、KHRT（KH + RT）が推奨だ。ディーゼル噴射ではKHRTが最もよく使われる。OpenFOAMのsprayFoamソルバーはLagrangian噴霧計算に特化しており、breakupModelとatomizationModelを別々に選択できる。
TABモデルとKH-RTモデル——噴霧CFDを支える2大分裂モデル
エンジン噴霧シミュレーションで支配的な分裂モデルはTAB（Taylor Analogy Breakup）とKH-RT（Kelvin-Helmholtz / Rayleigh-Taylor）の二系統です。TABは液滴を弾性球として振動方程式で解き、振動振幅が臨界値を超えたとき分裂と判定します。計算が軽い反面、大径液滴の「stripping breakup」を再現しにくい弱点があります。KH-RTは流体力学的不安定理論から導出され、高We数域での精度が高い一方、モデル定数のキャリブレーションが必要です。商用エンジンCFDでは両モデルをハイブリッドさせる「KH-RT hybrid」が主流となっています。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
液滴分裂モデルを使った噴霧解析の手順を教えてください。
ディーゼル噴射を例にとろう。
1. メッシュ作成: 噴射ノズル近傍は0.5 mm以下、噴霧領域は1〜2 mm
2. 噴射条件設定: 噴射圧力、ノズル径、噴射期間
3. 一次分裂モデル: Blob法（ノズル径の液滴を投入）またはLISAモデル
4. 二次分裂モデル: KHRTを選択（ディーゼルの場合）
5. 蒸発モデル: 必要に応じてRanz-Marshall蒸発モデルを有効化
6. 乱流分散: Stochastic tracking（DRW）で液滴の乱流拡散を考慮
メッシュの影響
メッシュサイズは噴霧結果に影響しますか？
大きく影響する。Lagrangian噴霧計算の最大の課題はメッシュ依存性だ。1つのCFDセルに大量のparcelが存在すると、気相への運動量フィードバックが過剰になる。
Abraham（1997）のガイドラインでは、各セルに占める液相体積分率が1%以下になるようメッシュサイズを設定するのが理想的だ。実務的にはAMRを使って噴霧先端近傍を自動的に細かくする方法が有効。
パラメータチューニング
Breakupモデルの定数はどう調整するんですか？
KHRTモデルの主要パラメータとその感度を示そう。
| パラメータ | 標準値 | 影響 |
|---|---|---|
| $B_0$ (KH子液滴径) | 0.61 | 大きくすると子液滴が大きくなる |
| $B_1$ (KH分裂時間) | 10〜40 | 大きくすると分裂が遅くなる |
| $C_{RT}$ (RT分裂定数) | 0.1〜1.0 | 大きくすると分裂が速くなる |
| $C_\tau$ (RT分裂遅延) | 1.0 | 大きくすると分裂開始が遅れる |
実験データ（液滴径分布、噴霧到達距離、噴霧角）と比較してチューニングする。SMD（Sauter Mean Diameter）の実測値との一致が最もよく使われる検証指標だ。
実験データとの比較
検証に使える実験データはありますか？
代表的なベンチマーク実験はこちらだ。
| 実験 | 条件 | 計測量 |
|---|---|---|
| ECN Spray A | n-ドデカン噴射、900 K雰囲気 | 液相到達距離、蒸気浸透長 |
| Hiroyasu & Kadota | ディーゼル条件 | 噴霧到達距離 |
| PDPA/PDA計測 | 位相ドップラー法 | 液滴径・速度分布 |
農薬散布ドローンのスプレーCFD——漂流と付着のバランス
農業用ドローンの農薬散布では、液滴が作物葉面に到達するか漂流（ドリフト）して周辺に飛散するかが、農薬効果と環境規制の両面で重要です。液滴径50 µm以下では風速2 m/sでも数メートル漂流し、作物への到達率が30%を下回ることが報告されています。CFDで液滴サイズ分布と気流（ローターダウンウォッシュ）を組み合わせた解析により、最適な飛行高度（地面から1.5〜2 m）と散布速度（5〜7 m/s）が導出されており、大手農機メーカーが製品認証の根拠として採用しています。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
液滴分裂モデルを実装しているツールを教えてください。
| ツール | Breakupモデル | 噴射モデル | 特徴 |
|---|---|---|---|
| Ansys Fluent | TAB, KHRT, SSD, Madabhushi | Solid Cone, Hollow Cone, Flat Fan | VOF-to-DPM転換機能あり |
| STAR-CCM+ | TAB, KHRT, Reitz-Diwakar | Cone, Blob, LISA | Lagrangian-Eulerian切り替え |
| OpenFOAM (sprayFoam) | TAB, ETAB, ReitzKHRT, PilchErdman | BlobInjection, ConeInjection | 完全オープンソース |
| CONVERGE | KH-ACT, KHRT, TAB | Blob-injection | AMRネイティブ対応 |
用途別推奨
用途によってツールは変わりますか？
CONVERGEって何ですか？
Convergent Science社の内燃機関専用CFDソルバーだ。自動メッシュ生成（AMR）が最大の特徴で、噴霧・燃焼計算でメッシュ設計の手間を大幅に削減できる。エンジンメーカーでの採用が急増している。
CONVERGE vs Fluent——エンジン噴霧解析ツールの選択
エンジン燃料噴霧CFDの分野では、CONVERGE（Convergent Science）が圧倒的な存在感を示しています。自動適応メッシュ細分化（AMR）によって燃料噴霧とシリンダー壁を別々のメッシュ解像度で扱える独自機能が、エンジン開発サイクルの短縮に直結するためです。Toyota、BMW、Cumminsを含む主要エンジンメーカーのほとんどがCONVERGEを噴霧解析のプライマリツールとして採用しています。一方、一般的な工業スプレー（農薬、塗装、製薬）ではFluentの汎用性とサポート体制が選ばれる傾向があります。
選定で最も重要な3つの問い
- 「何を解くか」：液滴分裂モデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
液滴分裂モデルの最新研究にはどんなものがありますか？
重要な方向性をいくつか紹介しよう。
VOF-to-DPM転換
Fluent 2020以降で搭載された機能で、VOF法で一次分裂を直接解像し、液滴がメッシュ解像限界以下になったらDPM粒子に自動変換する。
一次分裂と二次分裂を別の手法で扱うんですね。
その通り。ノズル近傍の液柱分裂（一次分裂）はVOF法の方が物理的に正確で、下流の液滴分裂（二次分裂）はLagrangianモデルの方が効率的だ。この切り替えを自動化するのがVOF-to-DPM転換だ。
DNSによる液滴分裂の直接計算
個々の液滴の分裂過程をVOF法やLevel Set法で直接計算するDNSが進んでいる。分裂モードの遷移メカニズムの解明や、既存モデルのクロージャ改良に活用される。
Jalaal & Mehravaran（2014）、Shinjo & Umemura（2010）らの研究が先駆的で、液柱分裂からの子液滴サイズ分布を統計的に抽出している。
超臨界噴射
超臨界条件では液滴分裂は起きないんですか？
超臨界条件（$p > p_c, T > T_c$）では液相-気相の界面が消失し、通常の液滴分裂は起きない。代わりに密度の大きく異なるジェットの混合過程になる。Real-fluid EOS（Peng-Robinson、SRK等）を使った超臨界噴射のCFDが活発に研究されている。
ロケットエンジンの液体酸素/水素噴射がこの条件に該当し、ECN（Engine Combustion Network）のSpray AやSpray Hのベンチマークデータが利用可能だ。
機械学習によるBreakupモデル
DNSデータを教師データとして、液滴分裂の確率やサイズ分布をニューラルネットワークで予測する研究が増えている。Weber数・Ohnesorge数・Mach数の多次元パラメータ空間を効率的にカバーできる利点がある。
二次分裂のLES直接解析——マルチスケール界面の挑戦
液滴分裂の「二次分裂」は、サブミリスケールの界面不安定が支配するため、Reynolds Averaged手法では本質的に再現が困難です。2020年代のフロンティアはLESと高解像度VOF/Level Setの組み合わせによる「直接界面解析」で、1024^3の計算格子でKelvin-Helmholtz不安定から二次分裂までを直接計算し、TABモデルが過小評価していた小径液滴生成量を2.3倍に修正する結果が報告されています。この規模の計算にはスーパーコンピュータ「富岳」相当の1万コア以上が必要で、産業応用にはまだ距離があります。
トラブルシューティング
トラブルシューティング
液滴分裂モデルでよくあるトラブルを教えてください。
順番に見ていこう。
1. 噴霧到達距離が実験と合わない
短すぎる場合: 分裂が速すぎる。$B_1$ を大きくしてKH分裂を遅延させる。
長すぎる場合: 分裂が遅すぎる。$B_1$ を小さくするか、メッシュを細かくして気相の運動量フィードバックを改善する。
2. SMD（Sauter Mean Diameter）が大きすぎる
対策:
- KH分裂定数 $B_0$ を小さくして子液滴を細かくする
- RT分裂定数 $C_{RT}$ を大きくしてRT分裂を促進
- 一次分裂の初期液滴径が適切か確認（Blob法のBlob径 = ノズル径）
3. 計算が発散する
噴霧計算は発散しやすいですか？
特に噴射開始直後に発散しやすい。高速液滴が1つのCFDセルに大きな運動量ソースを与えるためだ。
対策:
- タイムステップを十分に小さくする（噴射初期は $10^{-7}$ s以下）
- DPMのunder-relaxationを下げる（0.3〜0.5）
- メッシュを十分に細かくしてparcel密度を下げる
- 2-way couplingを最初は1-wayに設定し、安定後に切り替え
4. 噴霧角が合わない
対策:
- ノズル内部流れの影響を考慮（キャビテーションによる噴射角増大）
- 乱流分散モデル（DRW: Discrete Random Walk）のパラメータを調整
- 噴射条件（コーン角、速度プロファイル）の入力を再確認
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Parcel数を十分に多くする（最低10,000以上）。少ないと統計誤差が大きい |
| STAR-CCM+ | Injector Rate Profileの時間解像度に注意 |
| OpenFOAM sprayFoam | atomizationModelとbreakupModelの組み合わせ互換性を確認 |
| CONVERGE | AMRのEmbed levelが噴霧解像度を決定。Level 3以上を推奨 |
液滴が消える——数値蒸発バグの診断法
スプレーCFDで「液滴粒子数が計算途中で急減する」現象は、蒸発モデルのバグではなくメッシュサイズの問題であることが多いです。1つのセルに液滴体積分率が高くなりすぎると（α > 0.1程度）、ガスへの質量輸送量が過大評価されて液滴が瞬時に蒸発してしまいます。対策はスプレー領域のメッシュをより細かくすることで「injection cell size < 0.3 x spray cone diameter」が一つの目安です。液滴が消える直前のセル位置と体積分率を粒子トラッキングで記録することが診断の第一歩です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——液滴分裂モデルの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告