噴霧・微粒化
噴霧・微粒化の理論基礎
概要
先生、噴霧・微粒化のCFDって何を計算するんですか？
液体がノズルから噴射されて微細な液滴に分裂する過程（atomization、微粒化）をシミュレーションする。ディーゼルエンジン燃料噴射、ガスタービン燃料噴射、スプレー塗装、農薬散布、消火器噴霧など、あらゆる噴霧プロセスの設計に使う。
液滴分裂モデル（secondary breakup）とは違うんですか？
噴霧の過程は2段階に分かれる。ノズル出口で液柱やシートが分裂して液滴になるのが一次微粒化（primary atomization）、生成した液滴がさらに細かくなるのが二次微粒化（secondary breakup）だ。この記事では一次微粒化から噴霧全体のモデリングを扱う。
噴霧の基本パラメータ
噴霧を特徴づけるパラメータは何ですか？
以下が主要なパラメータだ。
| パラメータ | 定義 | 意味 |
|---|---|---|
| Weber数 $We$ | $\rho_g U_{rel}^2 d / \sigma$ | 空気力 vs 表面張力 |
| Ohnesorge数 $Oh$ | $\mu_l / \sqrt{\rho_l \sigma d}$ | 粘性 vs 表面張力 |
| SMD $d_{32}$ | $\sum d_i^3 / \sum d_i^2$ | Sauter平均径 |
| 噴射圧力 $\Delta p$ | $p_{inj} - p_{amb}$ | 噴射エネルギー |
| 噴霧角 $\theta$ | コーン角度 | 噴霧の広がり |
SMD（Sauter Mean Diameter）は噴霧の代表径として最も使われる。蒸発や反応の速度は表面積に比例するため、体積/表面積比の代表径が有用だ。
一次微粒化モデル
CFDでは一次微粒化をどうモデル化するんですか？
直接的にノズル内部流れから液柱分裂を追うのは計算コストが膨大なので、工学モデルが使われる。
| モデル | 概要 | 適用 |
|---|---|---|
| Blob injection | ノズル径の液滴を投入 | ディーゼル噴射 |
| LISA (Linearized Instability Sheet Atomization) | シート状液膜の不安定性 | 圧力噴霧ノズル |
| ELSA (Eulerian-Lagrangian Spray Atomization) | Eulerian液相→Lagrangian液滴遷移 | 研究用 |
Blob injection法はReitz（1987）が提案した最も実用的な手法で、ノズル径と同じ大きさの液滴（Blob）をDPMで投入し、二次微粒化モデル（KHRT等）で細かくしていく。
Rayleigh崩壊——液柱はなぜ液滴になるのか
蛇口から落ちる水が連続流ではなく液滴に分れるのは、Rayleigh（1878年）が解明した「Rayleigh不安定」のためです。円柱状の液柱は表面張力によって波長が直径πD以上の撹乱が増幅する固有不安定性を持ち、最終的にほぼ均一な液滴列に分裂します。この不安定成長率は液柱のウェーバー数やOhnesorge数で決まり、スプレーノズルの設計に直接応用されます。単分散液滴を生成するインクジェットプリンタはRayleigh崩壊を意図的に制御している好例で、±1%以内の液滴径制御がA4用紙に600 DPI以上の印刷を可能にしています。
噴霧・微粒化の数値計算手法
数値解法の詳細
噴霧シミュレーションの数値的なポイントを教えてください。
Lagrangian噴霧計算の最大の課題はメッシュ依存性だ。ノズル近傍に大量のparcelが集中するため、CFDセルサイズがparcelの運動量ソースに影響する。
メッシュ依存性の問題
Abraham（1997）のガイドラインでは、各セルに占める液相体積分率が低いこと（理想的に1%以下）が推奨される。実務的にはノズル近傍のメッシュを0.5〜2 mmにし、AMRを併用する。
AMRが噴霧で特に重要な理由は何ですか？
噴霧先端が移動するので、細かいメッシュが必要な領域が時間変化する。AMR（Adaptive Mesh Refinement）で噴霧が存在する領域のみを自動的に細分化すれば、固定メッシュに比べて計算コストを大幅に削減できる。CONVERGEはこのAMRがネイティブに組み込まれているため、噴霧計算で強みを発揮する。
ノズル内部流れの影響
高圧噴射ノズルでは、ノズル内部でキャビテーションが発生し、これが噴霧の微粒化を促進する。ノズル内部流れを事前に計算し、出口での乱流プロファイルと液膜分布をLagrangian噴霧計算の入口条件にする手法が高精度化に有効だ。
ツール別の実装
| ツール | 一次微粒化 | 二次微粒化 | AMR | 特徴 |
|---|---|---|---|---|
| CONVERGE | Blob, KH-ACT | KHRT, TAB | ネイティブ | 噴霧計算のベンチマーク |
| Ansys Fluent | Blob, LISA, Flat Fan | TAB, KHRT, SSD | Gradient-based | VOF-to-DPM転換 |
| STAR-CCM+ | Blob, LISA | TAB, KHRT, Reitz-Diwakar | Table-based | Lagrangian/Eulerianハイブリッド |
| OpenFOAM (sprayFoam) | BlobInjection, ConeInjection | TAB, ETAB, ReitzKHRT | dynamicRefineFvMesh | 完全OSS |
CONVERGE は自動メッシュ生成とAMRの組み合わせにより、メッシュ設計の手間を大幅に削減できるため、エンジンメーカーでの噴霧・燃焼計算の標準ツールになりつつある。
Rosin-Rammler分布——スプレーCFDの液滴径設定の定番
スプレーCFDで液滴径分布を設定する際の最も一般的な手法がRosin-Rammler分布 F(d) = 1-exp(-(d/d_bar)^n) です。d_bar（特性径）とn（分布幅係数）の2パラメータで広範な噴霧特性を表現できます。このパラメータはレーザー回折（Malvern等）による実測から決定しますが、測定条件（液圧、ノズルからの距離）によって大きく変わるため「どの条件で測ったデータか」の記録が不可欠です。CFDとRosin-Rammlerを使う場合、D32（ザウター平均径）が実測と一致するようにd_barとnをキャリブレーションすることが実務の標準手順です。
噴霧・微粒化の実務適用
実践ガイド
噴霧シミュレーションの手順を教えてください。
ディーゼル噴射（ECN Spray A条件）を例にとろう。
1. 定圧容器の気相: 900 K、60 bar の窒素雰囲気
2. メッシュ: ノズル近傍0.25 mm（AMRで自動細分化）、遠方2 mm
3. 噴射条件: n-ドデカン、噴射圧1500 bar、ノズル径90 μm
4. 一次微粒化: Blob injection（Blob径 = ノズル径）
5. 二次微粒化: KHRTモデル
6. 蒸発: Abramzon-Sirignanoモデル
8. 検証: 液相到達距離、蒸気浸透長をECNデータと比較
ECN（Engine Combustion Network）
ECNって何ですか？
Sandia National Laboratoryが主導する国際的な燃焼研究ネットワークだ。実験条件と計測データが公開されていて、噴霧・燃焼CFDの検証ベンチマークとして世界中で使われている。
| ケース | 燃料 | 条件 | 計測量 |
|---|---|---|---|
| Spray A | n-ドデカン | 900 K, 60 bar | 液相到達距離、蒸気浸透長、着火遅れ |
| Spray H | n-ヘプタン | 各種温度 | 同上 |
| Spray G | iso-オクタン | GDI条件 | Flash boiling噴霧 |
KHRTモデルのパラメータチューニング
Breakupモデルの定数はどう調整するんですか？
| パラメータ | 標準値 | 感度 | 影響 |
|---|---|---|---|
| $B_0$ | 0.61 | 高 | KH子液滴径（大→径大） |
| $B_1$ | 10〜40 | 高 | KH分裂時間（大→分裂遅延） |
| $C_{RT}$ | 0.1〜1.0 | 中 | RT分裂速度 |
まず液相到達距離を $B_1$ で合わせ、次にSMD（Sauter平均径）を $B_0$ で微調整する。ECNのSpray Aデータとの比較がキャリブレーションの標準だ。
よくある失敗
| 症状 | 原因 | 対策 |
|---|---|---|
| 噴霧が短すぎる | 分裂が速すぎる | $B_1$を大きくする |
| SMDが大きすぎる | 分裂不足 | $B_0$を小さく or $C_{RT}$を大きく |
| 発散 | メッシュ内のparcel集中 | AMRを有効化、$\Delta t$を下げる |
| 噴霧角が合わない | ノズル内部流れ未考慮 | DRW調整 or ノズル内部計算 |
ガスタービン燃焼器のスプレー——寿命vs性能のトレードオフCFD
ガスタービン燃焼器のリキッドフュエルスプレーは、液滴の蒸発速度・混合・点火・安定燃焼に直接影響するため、燃焼効率と排出ガス（NOx、スモーク）の両方を決める最重要因子です。ロールスロイス・GEアビエーション・三菱重工のエンジン開発部門では、スプレーノズル設計をCFDで最適化してから実験確認するプロセスを標準としており、CFDなしでは認証取得まで2〜3年かかっていた設計検証が1年以下に短縮されています。特に高圧噴霧条件（100 bar以上）での液滴挙動は古典理論が当てはまらず、非定常VOF-LESが必要です。
噴霧・微粒化のソフトウェア比較
商用ツール比較
噴霧・微粒化シミュレーションのツールを比較してください。
| ツール | 一次微粒化 | AMR | 燃焼連携 | 特徴 |
|---|---|---|---|---|
| CONVERGE | Blob, KH-ACT | ネイティブAMR | SAGE詳細化学 | エンジン噴霧の業界標準 |
| Ansys Fluent | Blob, LISA, Flat Fan | Gradient AMR | Flamelet, EDC等 | VOF-to-DPM機能 |
| STAR-CCM+ | Blob, LISA | Table AMR | 詳細化学 | ポリヘドラルメッシュ |
| OpenFOAM (sprayFoam) | Blob, Cone | dynamicRefine | 基本モデル | 完全OSS |
| AVL FIRE | 独自モデル | 対応 | 詳細化学 | エンジン特化 |
用途別推奨
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| ディーゼル噴射+燃焼 | CONVERGE | AMR + 詳細化学の統合 |
| GDI（ガソリン直噴） | Fluent, STAR-CCM+ | Flash boilingモデル |
| ガスタービン燃料噴射 | Fluent, STAR-CCM+ | 旋回流+噴霧 |
| スプレー塗装 | Fluent | DPM + 液膜 + 蒸発 |
| 農薬・消火噴霧 | Fluent, OpenFOAM | 低圧噴霧モデル |
| 学術・モデル開発 | OpenFOAM | ソースコードアクセス |
CONVERGEがエンジン噴霧で強い理由は何ですか？
3つの理由がある。(1) 自動メッシュ生成でCADからメッシュ設計なしで計算開始できる。(2) AMRが噴霧・燃焼領域を自動追従するのでメッシュチューニングが不要。(3) SAGE化学反応ソルバーが大規模反応メカニズム（数千反応）を効率的に解ける。これにより、エンジンの噴射-微粒化-蒸発-着火-燃焼-排ガスの一連を1つのシミュレーションで回せる。
CONVERGE CFD vs GT-SUITE——エンジン噴霧シミュレーションの分業
内燃機関の燃料噴霧解析では、三次元噴霧CFD（CONVERGE、Fluent）と一次元エンジンサイクルシミュレータ（GT-SUITE、AVL BOOST）の使い分けが確立しています。GT-SUITEはエンジン全体の吸排気系・燃焼サイクルを高速に計算できますが、噴霧の空間分布は扱えません。CONVERGEは自動AMRで噴霧の詳細を三次元で解けますが、エンジン全サイクルを計算するには数日の計算時間が必要です。実務では「GT-SUITEで全稼働域マップを計算→問題点をCONVERGEで詳細解析」という協調ワークフローが主流で、両ツールのAPI連携機能が2020年代から実用化されています。
噴霧・微粒化の先端研究
先端技術と研究動向
噴霧・微粒化の最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
VOF-to-DPM転換（Primary Breakup直接計算）
ノズル近傍の液柱分裂をVOF法で直接解像し、液滴がメッシュ解像限界以下になったらDPM粒子に自動変換する手法だ。Fluent 2020以降で搭載されており、一次微粒化を経験モデルに頼らず物理的に正確に計算できる。
一次微粒化の直接計算は計算コストが膨大ではないですか？
AMRと組み合わせれば現実的になりつつある。ノズル近傍のVOF領域のみ超高解像度（数μm）にし、下流でDPMに転換すれば計算量を抑えられる。Shinjo & Umemura（2010, 2011）のDNS研究が先駆的で、液柱表面のKH不安定性からリガメント形成、液滴分離までを直接捕捉している。
超臨界噴射
ロケットエンジンの液体酸素/水素噴射や超臨界CO2の噴霧では、臨界点を超えると液-気界面が消失する。Real-Fluid EOS（Peng-Robinson, SRK等）を用いた超臨界噴射のCFDが活発に研究されている。
Flash Boiling噴霧
Flash boilingは先の蒸発モデルの記事でも出てきましたね。
GDIエンジンの高温条件でノズル内が沸騰する現象だ。通常の微粒化とは異なるメカニズムで超微粒化が起こる。ECN Spray Gのデータベースが検証に使われている。Fluent、STAR-CCM+ともにFlash Boiling Spray専用モデルを搭載している。
機械学習による噴霧特性予測
噴射条件（圧力、温度、ノズル形状）からSMDや噴霧角を瞬時に予測するサロゲートモデルが研究されている。CFDで数百ケースのパラメトリックスタディを行い、Deep Learningで写像を構築する。エンジン制御マップの最適化への応用が期待されている。
電気流体力学（EHD）アトマイゼーション——電場で液滴を制御する
電場を印加した液体ノズル先端で「Taylorコーン」と呼ばれる円錐形の液体が形成され、先端から数nmサイズのイオン化液滴が放出されるエレクトロスプレー（ESI）は、生体分子の質量分析（Fenn, 2002年ノーベル化学賞）の根幹技術です。このEHDアトマイゼーションは印刷エレクトロニクスや薬物搭載ナノ粒子製造にも応用されており、印加電圧・液体流量・ノズル径で液滴径を100 nm〜10 µmの範囲でチューニングできます。CFDでのEHD解析は流体・電磁場・荷電粒子輸送の三連成が必要で、OpenFOAMのEHDFoamソルバーが研究グループの基盤ツールとして使われています。
噴霧・微粒化のトラブル対応
トラブルシューティング
噴霧シミュレーションでよくあるトラブルを教えてください。
順番に見ていこう。
1. 噴霧到達距離が合わない
短すぎる場合: 分裂が速すぎる。$B_1$（KH分裂時間定数）を大きくする。
長すぎる場合: 分裂が遅すぎる。$B_1$を小さくする、メッシュを細かくして気相の運動量フィードバックを改善する。
2. 液滴径（SMD）が実験と合わない
対策:
- KH子液滴径パラメータ $B_0$ を調整（小→SMD減少）
- RT分裂定数 $C_{RT}$ を調整（大→SMD減少）
- Blob injection のBlob径がノズル径と一致しているか確認
- 蒸発モデルが有効な場合、蒸発によるSMD変化も考慮
3. 噴射開始直後に発散する
噴射初期で計算が壊れます…
対策:
- タイムステップを十分小さくする（噴射初期は $\Delta t = 10^{-7}$ s以下）
- DPMのunder-relaxationを下げる（0.1〜0.3）
- まず1-way couplingで安定後に2-way couplingへ
- ノズル近傍のメッシュを細かくする（parcel集中の緩和）
4. 噴霧角が合わない
対策:
- ノズル内部のキャビテーションを考慮（噴霧角が増大する原因）
- DRW（乱流分散）のパラメータを調整
- LISA噴射モデルの場合、シート厚さとコーン角の設定を確認
- Blob injectionの場合、速度プロファイルの設定を確認
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| CONVERGE | AMR Embed levelが噴霧領域の解像度を決定。Level 3〜4推奨 |
| Fluent | Parcel数を十分に多くする（最低10,000以上）。少ないと統計誤差大 |
| STAR-CCM+ | Injector Rate Profileの時間解像度が噴射初期の精度に影響 |
| OpenFOAM | sprayFoamのatomizationModelとbreakupModelの互換性を確認 |
スプレーが壁に当たって飛散する——インパクト角の見落とし
スプレーCFDで実機と最も乖離しやすい現象が「液滴の壁面衝突後の挙動」です。液滴は衝突角・衝突速度・壁温・液滴径に応じて「付着」「広がり」「跳ね返り」「分裂（splashing）」の4つのモードを取ります。Splashing（We_crit > 3000程度）が発生すると1次液滴から多数の小径二次液滴が生成され、追跡粒子数が爆発して計算が重くなります。K = We^0.5 × Re^0.25 でK > 57.7のとき飛散が発生するという経験則（Mundo criterion）を使い、事前にスプレー全体のKマップを計算して問題箇所を特定することが実務的アプローチです。
関連トピック
なった
詳しく
報告