1. 問題設定
粗大液滴の噴霧・微細化シミュレータを検証する場合、理想的には次のようなデータが欲しくなる。
入口液滴粒径分布 PSD
↓
超音速加速流 / 衝撃波 / 膨張波
↓
下流液滴粒径分布 PSD
↓
SMD, D10, D50, D90, 速度分布, 粒径–速度相関しかし、実際にはこの条件を完全に満たす公開データは少ない。単一液滴の衝撃波管実験は多いが、入口から下流までポリディスパース噴霧のPSDを追ったデータは限られる。そこで本記事では、有料利用や著者への生データ依頼も含め、準一次元シミュレータ検証に最も使いやすい3つのデータ候補を整理する。
2. 推奨データセットの結論
| 優先 | 文献 | 主な特徴 | 使いどころ |
|---|---|---|---|
| 1 | Reuter & Tuttle, 2024 | 噴霧が衝撃波を通過する前後で、液滴径分布・速度分布・joint PDFを測定 | 入口PSD→衝撃波通過後PSDの最有力ベンチマーク |
| 2 | Menon & Gurunadhan, 2022 | 過膨張超音速二相噴流中でPDPAにより液滴径・速度を軸方向測定 | ショックトレイン中のPSD/SMD発展の検証 |
| 3 | Gobyzov et al., 2020 | 収束拡大ノズル内で液滴群が加速され、形状・サイズ分布・速度を計測 | 準一次元ノズル内加速流の検証 |
3. 第1候補：Reuter & Tuttle 2024
最も推薦したいのは、“Interactions between liquid sprays and shock waves in underexpanded flows” である。この研究は、underexpanded jet 中の衝撃波を通過する噴霧について、衝撃波上流・下流の液滴径分布と液滴速度分布をPDIとSchlierenで測定している。特に重要なのは、平均SMDだけでなく、液滴径分布、液滴速度分布、粒径–速度のjoint probability distributionまで扱っている点である。これは、準一次元モデルで「入口PSDを入れて、衝撃波・膨張波・圧縮波を通過した後のPSDを予測する」場合に非常に相性がよい。
この論文の価値は、単なる単一液滴のbreakup time検証ではなく、ポリディスパース噴霧の統計量変化を直接扱っている点にある。論文概要では、oblique shock、weak normal shock、strong normal shockで液滴径・速度分布の変化が異なることが示されており、強い垂直衝撃波では平均液滴径と平均液滴速度が一貫して低下する、と整理されている。
準一次元コードに落とし込むなら、このデータは次の検証に向く。
| 検証項目 | 比較対象 |
|---|---|
| 衝撃波前後のD32変化 | PDI測定の平均径 |
| D10/D50/D90の変化 | 液滴径分布 |
| 液滴速度分布 | PDI速度分布 |
| 粒径–速度相関 | joint PDF |
| 衝撃波種別ごとのモデル妥当性 | oblique / weak normal / strong normal shock |
このデータは、あなたのシミュレータにとって最も“入口PSD→下流PSD”に近い実験系である。ただし、論文PDFや図表だけでは完全なイベントデータが得られない可能性があるため、真に厳密な検証をするなら、著者にPDIの生イベントデータを依頼するのが望ましい。
4. 第2候補：Menon & Gurunadhan 2022
第2候補は、“Droplet behavior in overexpanded supersonic two-phase jets” である。この研究は、過膨張超音速ジェット中に噴霧を入れ、ShadowgraphyとPhase Doppler Particle Anemometry、さらに流れ場シミュレーションを組み合わせて、液滴サイズと速度の変化を評価している。論文概要では、shock train構造に対応して液滴サイズが周期的に増減し、その原因としてcoalescenceとbreakupが示唆されている。
この文献は、入口からノズル内で加速する問題というより、ノズル出口後の過膨張噴流・ショックセル・ショックトレイン中のPSD発展を見るのに向いている。準一次元モデルが、ノズル出口後の圧縮・膨張を擬似一次元的に扱うなら、かなり有用である。
使い方としては、次のように整理するとよい。
| 検証項目 | 使い方 |
|---|---|
| 軸方向SMD変化 | 計算SMD(x)と比較 |
| 液滴速度変化 | Stokes数が大きい粒子の慣性応答を検証 |
| shock trainの影響 | 圧縮・膨張領域ごとのPSD変化を比較 |
| breakup / coalescenceモデル | SMD増減の符号と周期を比較 |
この研究の特徴は、液滴速度がshock structureにあまり影響されず、中心線上の液滴がballistic particleのように振る舞うと説明されている点である。これは、準一次元モデルで液滴慣性を扱う際の重要な検証ポイントになる。
5. 第3候補：Gobyzov et al. 2020
第3候補は、“Study of Deformation and Breakup of Submillimeter Droplets’ Spray in a Supersonic Nozzle Flow” である。この論文はオープンアクセスで、収束拡大ノズル内の加速流中に導入された15–60 µm程度の液滴群について、PIVとshadow photographyで、ガス速度、液滴速度、液滴形状、サイズ分布を測定している。
この研究は、あなたの準一次元モデルに対して構造的にかなり近い。なぜなら、液滴が収束拡大ノズル内で亜音速から超音速側へ加速される流れを通過するためである。論文では、ノズルスロート寸法、チョーク流れ、Mach diskの位置、液滴導入位置なども記述されており、準一次元の流路計算と対応させやすい。
ただし、注意点もある。対象液滴は主に15–60 µmであり、「粗大液滴」というよりは小粒径側である。したがって、1 mm級の粗大液滴の二次破砕ではなく、すでに一次微粒化された液滴群が、超音速加速流中でさらに変形・破砕する問題として使うのがよい。
このデータで検証すべき項目は次の通り。
| 検証項目 | 比較対象 |
|---|---|
| ノズル軸方向の平均粒径 | 実験の平均粒径分布 |
| サイズ分散 | PSD幅の変化 |
| 液滴速度 | PTV/PIV由来の液滴速度 |
| Weber数履歴 | ガス速度・密度・液滴径から再計算 |
| breakup threshold | We ≈ 12付近の破砕有無 |
| breakup mode | bag / multi-bag / pulling breakup |
この論文は、各可視化領域で多数の液滴統計を集めており、サイズ分布と形状変化を検証しやすい。特に、平均粒径とサイズ分散が軸方向にどう変化するかを、準一次元モデルのPSD evolutionと比較するのに向いている。
6. 3つのデータの使い分け
3つのデータは、同じ目的で横並びに使うより、役割を分けるべきである。
| 目的 | 最適データ |
|---|---|
| 衝撃波前後のPSD変化 | Reuter & Tuttle 2024 |
| 過膨張ジェット・ショックトレイン中のPSD発展 | Menon & Gurunadhan 2022 |
| ノズル内加速流中の液滴PSD・速度・形状変化 | Gobyzov et al. 2020 |
| 粒径–速度joint PDF | Reuter & Tuttle 2024 |
| SMD(x)の検証 | Menon 2022、Gobyzov 2020 |
| 準一次元ノズルモデルとの対応 | Gobyzov 2020 |
| 衝撃波モデルとの対応 | Reuter 2024 |
| coalescenceを含む統計変化 | Reuter 2024、Menon 2022 |
7. 実装者向けの検証ワークフロー
準一次元シミュレータでは、いきなり全物理を合わせに行くと、誤差の原因が分からなくなる。したがって、次の順番がよい。
Step 1：ガス流れ場を合わせる
まず、ノズル内または噴流中のガス速度、密度、Mach数、衝撃波位置を合わせる。ここがずれると、Weber数履歴がすべてずれる。
x → M(x), U_g(x), ρ_g(x), p(x), T(x)Step 2：液滴速度を合わせる
次に、液滴慣性応答を確認する。特にMenon 2022では、液滴速度がshock trainの局所構造に強く追従しないことが示されているため、Stokes数が大きい領域では、液滴を過度にガス速度へ追従させないモデルが必要になる。
Step 3：breakup thresholdを合わせる
Gobyzov 2020では、Weber数が低い領域での小液滴破砕が扱われている。まずは、Weがしきい値未満のときに破砕しないこと、Weがしきい値を超えたときにbag / multi-bag / pullingなどのモードが発生することを確認する。
Step 4：PSD evolutionを合わせる
最後に、入口PSDを粒径ビンに分け、各ビンについて
drag → relative velocity → Weber number → breakup probability → child droplet distributionを計算し、下流のPSD、D32、D50、D90、速度分布と比較する。
8. 結論
有料でもよい前提で、“完全な入口PSD→下流PSD”に最も近いデータを選ぶなら、第一候補は Reuter & Tuttle 2024 である。理由は、衝撃波前後の液滴径分布・速度分布・joint PDFを扱っており、ポリディスパース噴霧の統計変化を直接検証できるためである。
一方、Menon & Gurunadhan 2022 は、過膨張超音速ジェット中のショックトレインに沿ったSMD変化・速度変化を見るのに向く。これは、ノズル出口後の超音速二相噴流モデルに強い。
Gobyzov et al. 2020 は、収束拡大ノズル内で液滴群が加速されるため、準一次元ノズル内モデルとの対応が取りやすい。粒径は粗大液滴より小さいが、ノズル内加速流中のPSD・速度・変形・破砕を検証するには非常に有用である。
実務的には、次の3本立てが最も堅い。
Reuter 2024 ：衝撃波前後の入口PSD→下流PSD
Menon 2022 ：過膨張ジェット中のSMD/PSD発展
Gobyzov 2020 ：ノズル内加速流中の液滴群挙動この3つを組み合わせれば、準一次元の粗大液滴微細化シミュレータに対して、流れ場、液滴速度、breakup threshold、PSD evolutionを段階的に検証できる。
コメント