エグゼクティブサマリー
準一次元シミュレータで「粗大液滴の入口粒径分布を与え、キャリアガスが一様流として加速されて超音速に達する」問題を検証するなら、最も整合的な一次ベンチマーク群は、衝撃波管で単一液滴を衝撃後のほぼ一様流にさらす実験です。理由は、流れの主方向が明確で、Weber数・Ohnesorge数・Mach数を比較的独立に制御しやすく、液滴の変形・分裂開始・分裂様式・子液滴形成の時間発展を追えるためです。一方で、これらの古典的主力データの多くは単分散の単一液滴であり、ユーザーが最終的に必要とする入口PSDそのものは与えません。したがって、実務上は、単一液滴の物理カーネル検証 → 入口PSDのアンサンブル検証 → 実噴霧/二流体噴霧でのPSD発展検証という三段階が最も堅実です。
トップ候補の中では、Faeth系の一連の衝撃波管実験が、bag・shear・multimode の各分裂領域で、子液滴サイズ、速度、分裂時間、液滴形状変化を整理しており、粗大液滴入口からの二次微粒化モデルの核検証に最も有用です。特に、Chou 1997、Chou 1998、Dai 2001、Hsiang 1993 は、時間発展と子液滴特性を扱うため、単一液滴モデルの遷移判定・分裂時刻・代表粒径予測の較正に向きます。
超音速性そのものの影響を明示的に含めるには、Theofanous 系の Mach 3 希薄超音速流実験と、Wang 2020 の 一定 We = 約1100で M∞ を 0.3–1.19 に掃引した衝撃波管実験が特に重要です。前者は高We・高Mach域の分裂機構再解釈、後者はMach数だけを比較的分離して観察できるため、準一次元モデルの圧縮性補正の検証に最適です。
ただし、**「入口が粗大液滴PSD」「キャリアが超音速一様流」「下流でPSD進化を高分解能で計測」**を同時に満たす公開ベンチマークは非常に少ないです。したがって、最終的な推奨は次の通りです。
第1層：単一液滴衝撃波管データで局所分裂則を検証。
第2層：その局所分裂則で入口PSDをモンテカルロ/節点法で伝播し、代表統計量を検証。
第3層：Varga 2003、Menon 2022、Liu 2022 などの噴霧・二流体系で、SMD・D10/D50/D90・軸方向PSD進化を二次検証する。
選定方針と検証の考え方
本調査では、文献を次の優先順位で選びました。第一に、原著の実験論文であること。第二に、一様流に近い外力条件、またはそれに近い比較可能な条件を持つこと。第三に、検証に直接使える出力、つまり breakup time、変形比、子液滴径、子液滴速度、SMD、空間分布、モード遷移、あるいはその代理量が含まれることです。これに従うと、最上位は衝撃波管の単一液滴実験で、噴霧ノズルや二流体噴霧は補完扱いになります。これは、Varga 2003 でも一次破砕がほぼ一定ガス速度域で完了すると示されている一方、実噴霧系ではノズル近傍の三次元性・液膜形成・初期コア不安定が不可避だからです。
準一次元モデルにとって最も重要な比較量は、単一粒径に対する条件付き応答と、入口PSDに対する統計応答を分けることです。前者では、個々の初期粒径 (d_0) に対する分裂開始時刻、最大変形比、モード、子液滴代表径が主指標です。後者では、入口PSDから下流PSDへの写像、すなわち SMD、D10/D50/D90、数密度、質量分率、未分裂液滴の残存率が主指標になります。後者を直接与える理想ベンチマークは少ないため、実際には前者の高品質データを積み上げて後者を構成するのが現実的です。これは本報告の重要な推奨事項です。
推奨トップベンチマーク
| 順位 | 文献・一次ソース | 実験条件 | 入口液滴分布 | 検証に使える主出力 | 計測法・不確かさ | 準一次元粗大液滴入口への適合性と限界 | 推奨抽出データ・生データ状況 / 推奨比較指標 |
|---|---|---|---|---|---|---|---|
| 1 | Nicholls, J.A.; Ranger, A.A., Aerodynamic shattering of liquid drops, AIAA J. 7(2), 285–290 (1969), urlDOI:10.2514/3.5087https://doi.org/10.2514/3.5087. | 衝撃波により液滴を高速気流へ急加速。空気。Mach/速度範囲は検索断片では未特定。 | 単一液滴。平均径・分布型・初期径計測法は未特定。 | 分裂様式、変形/破砕の代表時間尺度。Ranger–Nicholls の特性時間の基礎。 | 古典的可視化。具体手法・不確かさは本検索では未特定。 | 古典基準として重要。準1Dの加速問題に近いが、PSDや子液滴統計は弱い。 | 抽出推奨: breakup onset、特性時間、最大変形比。生データ公開は未確認。比較指標: (t_b/t_c)、破砕モード、最大横径比。 |
| 2 | Hsiang, L.-P.; Faeth, G.M., Drop properties after secondary breakup, Int. J. Multiphase Flow 19, 721–735 (1993), urlDOI:10.1016/0301-9322(93)90039-Whttps://doi.org/10.1016/0301-9322(93)90039-W. | 衝撃波起因擾乱。bag / multimode / shear 領域。空気、常温常圧。 | 単一液滴。液種・初期径詳細は本検索では未特定。 | 分裂後の液滴特性。子液滴径・速度・残存母滴の情報が主用途。 | 検索断片では手法未記載だが、同系統研究では shadowgraphy / holography 系。直接の不確かさは未特定。 | 子液滴統計の一次ベンチマークとして非常に有用。Mach掃引は弱い。 | 抽出推奨: 子液滴代表径、速度分布、母滴残存率。機械可読データ公開は未確認で図 digitize 前提。比較指標: SMD、D32、D10/D50/D90、子液滴速度平均/分散。 |
| 3 | Hsiang, L.-P.; Faeth, G.M., Drop deformation and breakup due to shock wave and steady disturbances, Int. J. Multiphase Flow 21(4), 545–560 (1995), urlDOI:10.1016/0301-9322(94)00095-2https://doi.org/10.1016/0301-9322(94)00095-2. | 衝撃波擾乱と定常擾乱の比較。3試験装置を使用。詳細圧力・Machは本検索では未特定。 | 単一液滴。分布型なし。入口はモノドロップ。 | 変形特性、分裂境界、衝撃波条件と定常条件の比較。 | 複数設備比較。具体計測法・不確かさは本検索では未特定。 | 準1Dモデルに極めて重要。衝撃波管結果を定常一様流モデルへ写像する橋渡しになる。PSD出力は限定的。 | 抽出推奨: 変形履歴、モード遷移境界、steady/shock 等価性。生データ公開は未確認。比較指標: (d_\perp/d_0)、(d_\parallel/d_0)、breakup map in We–Oh。 |
| 4 | Chou, W.-H.; Hsiang, L.-P.; Faeth, G.M., Temporal properties of drop breakup in the shear breakup regime, Int. J. Multiphase Flow 23(4), 651–669 (1997), urlDOI:10.1016/S0301-9322(97)00006-2https://doi.org/10.1016/S0301-9322(97)00006-2. | 衝撃波擾乱、空気、常温常圧。We=125–375、Oh=0.003–0.040、液/気密度比 670–990、Re=3000–12000。 | 単一液滴。液体物性を変えて Oh を掃引。初期径平均・分布型は検索断片では未特定。 | shear領域の時間発展、分裂機構、断片径/速度情報。 | pulsed shadowgraphy と holography。定量不確かさは未特定。 | 高Weの周辺で最重要。粗大液滴入口から微粒化への主経路が shear になる条件で有力。 | 抽出推奨: breakup initiation time、断片サイズ/速度、strip 速度。生データ公開は未確認。比較指標: (t_b/t_c)、平均断片径、幅方向拡がり、未分裂率。 |
| 5 | Chou, W.-H.; Faeth, G.M., Temporal properties of secondary drop breakup in the bag breakup regime, Int. J. Multiphase Flow 24(6), 889–912 (1998), urlDOI:10.1016/S0301-9322(98)00015-9https://doi.org/10.1016/S0301-9322(98)00015-9. | 衝撃波擾乱、空気、常温常圧。水・エタノール・グリセリン混合。We=13–20、Oh=0.0043–0.0427、Re=1550–2150、密度比 633–893。 | 単一液滴。モノドロップ。 | bag 領域の時間発展、膜形成・破断過程。 | 検索断片では直接未記載だが同系統で pulsed optical 計測。不確かさ未特定。 | 入口PSDの大粒側でWeが低い場合の基準。超音速というより regime 基準。 | 抽出推奨: bag形成時刻、最大扁平化、bag破断時刻。生データ公開は未確認。比較指標: モード判定、(t_b/t_c)、最大投影面積、D50。 |
| 6 | Dai, Z.; Faeth, G.M., Temporal properties of secondary drop breakup in the multimode breakup regime, Int. J. Multiphase Flow 27, 217–236 (2001), urlDOI:10.1016/S0301-9322(00)00015-Xhttps://doi.org/10.1016/S0301-9322(00)00015-X. | 衝撃波管、空気、常温常圧。水・エタノール。液/気密度比 >500、Oh<0.1、We=15–150。 | 単一液滴。初期径分布はなし。 | multimode 領域の時間発展、分裂機構、分裂結果。 | pulsed shadowgraphy と holography。 | bag と shear の間の橋渡しデータ。入口PSDでWeが分布する場合に不可欠。 | 抽出推奨: multimode での breakup time、残存母滴、断片代表径。生データ公開は未確認。比較指標: モード別確率、(t_b/t_c)、SMD、質量分率。 |
| 7 | Theofanous, T.G.; Li, G.J.; Dinh, T.N., Aerobreakup in Rarefied Supersonic Gas Flows, J. Fluids Eng. 126(4), 516–527 (2004), urlDOI:10.1115/1.1777234https://doi.org/10.1115/1.1777234. | 希薄超音速空気 Mach 3。水またはtributyl phosphate。We=26–2600。 | 単一液滴。mm級。例として TBP 3.7–3.9 mm、We 28/57/109/2643 の図示あり。分布型なし。 | 界面不安定、混合、モード遷移、詳細可視化。 | 高解像可視化。具体的不確かさは未特定。 | 超音速性を強く含む最重要古典。ただし希薄流であり、常圧一様流への直接転用には注意。 | 抽出推奨: モード境界、波数、breakup onset、画像系列の digitize。生数値データは未確認。比較指標: breakup map in We, onset time, wave count/spacing, 最大広がり。 |
| 8 | Theofanous, T.G.; Li, G.J., On the physics of aerobreakup, Phys. Fluids 20, 052103 (2008), urlDOI:10.1063/1.2907989https://doi.org/10.1063/1.2907989. | 超音速ガス流中の液滴。LIF可視化。高We域では We>10^3、低We域では We<10^2 を強調。 | 単一液滴。水・TBP。例として 1 mm 水滴、TBP の複数ケース。 | 高Weでの SIE、低Weでの RTP/multiwave piercing の高解像画像。 | laser-induced fluorescence。従来 shadowgraph の誤認を修正。数値不確かさは未特定。 | 高We・超音速領域の機構検証に最重要。子液滴PSDは弱い。 | 抽出推奨: RT/KH 波長、sheet成長、モード識別。生の機械可読データは未確認、図/画像は利用可。比較指標: 波長、sheet厚さ、onset、mode classification。 |
| 9 | Wang, Z.; Hopfes, T.; Giglmaier, M.; Adams, N.A., Effect of Mach number on droplet aerobreakup in shear stripping regime, Exp. Fluids 61:193 (2020), urlDOI:10.1007/s00348-020-03026-1https://doi.org/10.1007/s00348-020-03026-1. | 24 m 衝撃波管。空気/CO2、水/エチレングリコール。We≈1100一定、M∞=0.3–1.19、Re=2.6×10^3–2.4×10^4。ポストショック一様流は約1.6 ms。 | 単一液滴。水 3.1 mm、EG 2.6 mm。分布型なし。生成は外径0.9 mm針。 | 横径変化、変形過程、breakup initiation、破片分布の形態変化、Mach依存性。 | shadowgraph/schlieren、HPV-X、100 kfps、0.087 mm/pixel。横径の 90%信頼区間を4反復から算出。 | 最有力の準1D圧縮性ベンチマーク。Mach効果をWe一定で分離。PSD全量は未公開。 | 抽出推奨: (d_c/d_0(t))、breakup onset、ケース別画像列、M依存。オープンアクセスPDFあり。比較指標: (t_b/t_c)、最大変形比、fragment spread、Mach補正誤差。 |
| 10 | Sharma, S.; Singh, A.P.; Rao, S.S.; Kumar, A.; Basu, S., Shock induced aerobreakup of a droplet, J. Fluid Mech. 929, A27 (2021), urlDOI:10.1017/jfm.2021.860https://doi.org/10.1017/jfm.2021.860. | 衝撃波 Mach 数 Ms≈1.1–1.8。水滴 Di=0.5, 2.5, 2.9 mm。We≈200–4133 を含む。 | 単一液滴。monodisperse single-drop。 | 波動段階と breakup 段階の連続観測、top-view のRT/KH様式、recurrent sheet breakup、補助表・動画。 | schlieren と shadowgraphy、20k–100k fps、局所高解像では 8 µm/pixel、補足 movie/table あり。 | 近年で最良の高時空間ベンチマークの一つ。機構検証に強いがフルPSDは弱い。 | 抽出推奨: 波長、sheet/ligament 発生周期、breakup onset、top-view画像系列。補足表と動画あり、完全生数値PSDは未確認。比較指標: RT/KH波長、onset、周期、モード遷移、最大広がり。 |
上の10件を、準一次元・粗大液滴入口という観点でさらにまとめると、最初に使うべき核セットは Wang 2020 + Chou 1997 + Chou 1998 + Dai 2001 + Hsiang 1993 です。ここで、Mach依存、shear/bag/multimode の領域別時系列、子液滴統計を一通り押さえられます。次に、Theofanous 2004/2008 と Sharma 2021 を加えると、高Mach・高We・高分解能の機構検証まで拡張できます。Nicholls 1969 と Hsiang 1995 は、基礎時間尺度と shock/steady の写像確認に有用です。
追加で有用な関連データセット
Hsiang & Faeth 1992, Near-limit drop deformation and secondary breakup, IJMF 18, 635–652, urlDOI:10.1016/0301-9322(92)90036-Ghttps://doi.org/10.1016/0301-9322(92)90036-G。衝撃波起因擾乱、空気、常温常圧で、**near-limit の変形と分裂開始条件を整理した高被引用古典です。粗大液滴PSDのうち、We がしきい値近傍にある成分の未分裂率**や分裂開始確率を校正するのに価値があります。
Varga, Lasheras & Hopfinger 2003, Initial breakup of a small-diameter liquid jet by a high-speed gas stream, JFM 497, 405–434, urlDOI:10.1017/S0022112003006724https://doi.org/10.1017/S0022112003006724。単一液滴ではなく**高速度同軸ガス流による液柱一次破砕**ですが、主要な atomization がガス速度がほぼ一定の領域で完了すると示し、液滴 sudden exposure の二次破砕と機構相似を議論しています。入口PSDを持つスプレー側ベンチマークへの橋渡しとして有効です。
Dorschner et al. 2020, On the formation and recurrent shedding of ligaments in droplet aerobreakup, JFM 904, A20, urlDOI:10.1017/jfm.2020.699https://doi.org/10.1017/jfm.2020.699。高倍率 shadowgraphy と3D数値を組み合わせ、ligament formation と recurrent shedding を解析した近年の高品質研究です。準一次元モデルで難しい sheet→ligament→droplet のサブモデル検証に適しますが、主力は機構画像であり、公開PSD表は限定的です。
Menon & Gurunadhan 2022, Droplet behavior in overexpanded supersonic two-phase jets, IJMF 152, 104076, urlDOI:10.1016/j.ijmultiphaseflow.2022.104076https://doi.org/10.1016/j.ijmultiphaseflow.2022.104076。**超音速共流二相噴流**で、shadowgraphy と PDPA により droplet size/velocity の軸方向変化を調べています。ショック列に沿った SMD の周期変動や、液滴速度の鈍感性が示されており、入口PSDを持つ噴霧の下流統計を見る二次ベンチマークとして有用です。
Liu et al. 2022, Experimental study of the spray characteristics of twin-fluid atomization: Focusing on the annular flow regime, Phys. Fluids 34, 123309, urlDOI:10.1063/5.0128231https://doi.org/10.1063/5.0128231。shadowgraph と PDPA により、軸方向 SMD、対数正規PSD、速度分布、GLR依存性、不確かさ表を与えます。超音速一様流そのものではありませんが、粗大液滴入口PSDを扱う統計モデルの下流PSD進化を確認する実用ベンチマークとして価値があります。
推奨比較指標と図表
モデル比較は、単一液滴条件付き指標とPSD統計指標の二層で行うのがよいです。前者では、(t_b/t_c)（breakup initiation の無次元時間）、最大横変形比 (d_\perp/d_0)、モード分類、sheet/ligament 発生位置と波長を主に見ます。後者では、SMD、D10/D50/D90、質量基準PSD・個数基準PSD、子液滴速度平均と分散、未分裂大滴の残存質量分率、必要なら breakup length または breakup completion time を見ます。Wang 2020 と Sharma 2021 は前者、Hsiang 1993 と Chou/Dai 系は両者、Menon 2022 と Liu 2022 は後者に向いています。
可視化は、少なくとも次の5種類を推奨します。
第一に、We–Oh–M の regime map。
第二に、軸方向または無次元時間に対する SMD・D50・未分裂率の推移。
第三に、PSD evolution heat map（横軸: (x/d_0) または (t/t_c)、縦軸: (d/d_0)、色: 質量または個数確率密度）。
第四に、joint PDF of droplet size and velocity。
第五に、mode atlas（同一縮尺の時系列画像を並べた比較図）です。レガシー論文の多くは機械可読データがなく、図の digitize が前提になりますが、Wang 2020 はオープンアクセス、Sharma 2021 は補足動画/表があり、追加抽出が比較的容易です。
生データが入手できるなら、最初に生成すべき図は PSD evolution plot と breakup onset / mean size / residual coarse fraction の3本です。具体的には、
- ( \mathrm{SMD}(x) ), (D50(x)), (D90(x))
- ( P(d|x) ) または ( P(d|t) ) のヒートマップ
- ( F_{\mathrm{coarse}}(x; d>d_{\mathrm{th}}) )
- ( \bar{u}d(d,x) ) と ( \sigma{u_d}(d,x) )
を推奨します。入口PSDを直接与えない単一液滴実験では、粒径クラスごとの条件付き結果を重み付き合成して、擬似的な (P(d|x)) を作るのが実務的です。これは本報告の推奨ワークフローです。
以下の mermaid は、採用候補を時系列に整理するための最小テンプレートです。年代の流れと「何を検証すべきか」を一目で見られます。
1969Nicholls & Rangershock-inducedshattering1993Hsiang & Faethchild-drop properties1995Hsiang & Faethshock vs steady1997Chou et al.shear breakuptemporal data1998Chou & Faethbag breakuptemporal data2001Dai & Faethmultimode temporaldata2004Theofanous et al.Mach 3 rarefiedaerobreakup2008Theofanous & LiLIF-based RTP/SIEphysics2020Wang et al.fixed-We Machsweep2021Sharma et al.high-resolutionshock-inducedaerobreakup液滴エアロブレークアップ主要実験の流れShow code
検索語は、英語では “shock tube droplet breakup”, “secondary atomization supersonic uniform flow”, “single droplet shock-induced breakup holography size velocity”, “Mach number effect droplet aerobreakup”, “rarefied Mach 3 droplet breakup”, “supersonic two-phase jet PDPA droplet size velocity”, “twin-fluid atomization annular flow PDPA SMD” を強く推奨します。日本語では 「衝撃波管 液滴 破砕 実験」, 「二次微粒化 液滴 衝撃波」, 「超音速 一様流 液滴 分裂 計測」, 「液滴 エアロブレークアップ マッハ数」, 「二流体噴霧 PDPA SMD 軸方向」, 「超音速 二相噴流 液滴径 速度」 が有効です。
制約と未解決点
最大の制約は、理想的な公開ベンチマークがほぼ存在しないことです。つまり、入口に粗大液滴PSDを与え、キャリアガスが一様流として超音速へ加速し、下流のPSD・速度場・数密度・不確かさが同時に公開されているデータは、今回確認した範囲では見当たりませんでした。そのため、単一液滴の原理検証と、噴霧系の統計検証を分ける必要があります。
また、Theofanous 2004 のような 希薄 Mach 3 データは、圧縮性の影響を強く含む点で貴重ですが、常圧・高密度の超音速一様流への直接転用には注意が必要です。逆に、Menon 2022 や Liu 2022 のような噴霧系は PSD 発展を見やすい一方、ノズル近傍の三次元性や coalescence を含むため、準一次元の純粋な breakup only モデルにはそのままは対応しません。したがって、両者を無理に同列比較するのでなく、役割分担して使うのがよいです。
最後に、レガシー主力論文の多くは、raw numeric data を公開していません。実務的には、Wang 2020、Sharma 2021、一部の近年 open-access 論文を優先し、古典論文は図の digitize を前提にするのが現実的です。もし著者データや補足表が入手できるなら、最初に数値化すべきは breakup onset, 最大変形比, 子液滴代表径, 速度–粒径相関 です。
コメント