噴霧・微粒化
理論と物理
概要
先生、噴霧・微粒化の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以上の印刷を可能にしています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
噴霧シミュレーションの数値的なポイントを教えてください。
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をキャリブレーションすることが実務の標準手順です。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
噴霧シミュレーションの手順を教えてください。
ディーゼル噴射（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が必要です。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
噴霧・微粒化シミュレーションのツールを比較してください。
| ツール | 一次微粒化 | 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年代から実用化されています。
選定で最も重要な3つの問い
- 「何を解くか」：噴霧・微粒化に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
噴霧・微粒化の最新研究にはどんなものがありますか？
いくつかの方向性を見ていこう。
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マップを計算して問題箇所を特定することが実務的アプローチです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——噴霧・微粒化の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告