【物理】真空中へ噴出する超音速蒸気自由噴流への水スプレー微粒化・蒸発・衝突を予測するシミュレータ設計（１）

エグゼクティブサマリー
対象問題と前提条件
物理モデル設計
数値解法と計算設定
ソルバー候補と実装アーキテクチャ
入出力仕様とユーザーインターフェース
検証・妥当性確認と優先文献・データセット

エグゼクティブサマリー

本件は「飽和蒸気（1 ata ≒ 1 bar abs）を末広管（内径10 mm・長さ100 mm）から 1000 Pa の真空へ蒸気放出し、管出口でマッハ1（チョーク）とし、出口から 50 mm 下流の平板に衝突させる」条件で、上流から水を噴霧（噴射位置＝上流リザーバ／管入口／管壁途中）したときの 液滴粒径分布（SMD と全 PSD）、破砕レジーム、蒸発・凝縮、二相カップリング、平板衝突後の挙動を一体で予測する数値シミュレータの設計課題である。真空中への過膨張（超膨張）自由噴流では、衝撃波セルや膨張扇に加え、希薄化・凝縮・非理想性が効き得ることが古典的に指摘されており、平板は運動量束を直接測る“プローブ”にもなる一方、衝突により流れ場自体が大きく再編される。

設計としては、**多忠実度（multi‑fidelity）**が現実的である。
(1) 形状（A(x)）と境界条件の妥当性確認・パラメータスイープ用に、準一次元（quasi‑1D）＋自由噴流近似（MoC/等価噴流）＋簡易液滴追跡の“高速モデル”を用意し、(2) 主要成果物（PSD/蒸発/衝突）を担う中核として **圧縮性気相＋ラグランジュ液滴（Eulerian–Lagrangian）**を基本に据え、必要に応じて (3) 高忠実度（LES、希薄域のハイブリッド CFD‑DSMC、湿り蒸気の非平衡凝縮）へ段階的に拡張する。希薄気体には DSMC が有効で、並列実装のオープンソースとして SPARTA が利用できる。

蒸気物性は、産業用途で標準的な IAPWS‑IF97（密度・音速・比熱・エンタルピ等を高速計算）に基づく実装が堅牢である。
また、飽和蒸気の急膨張では過冷却・核生成・液滴成長を伴う 非平衡（homogeneous）凝縮が問題になり得る。湿り蒸気 CFD の国際比較（IWSMP）でも、凝縮モデルにより結果のばらつきが大きいことが報告されているため、凝縮を「無視できる」前提とせず、モデル選択を明示的に切り替えられる設計が重要である。

以下では、必要物理、数値手法、メッシュ・時間刻み、境界条件、カップリング戦略、入出力仕様、UI、ソルバー候補、V&V（検証・妥当性確認）計画、優先文献・データセットを、この装置条件に即して整理する。

対象問題と前提条件

幾何と運転条件の明示

蒸気：飽和蒸気、供給圧力 1 ata（特記ない場合 1 bar abs と解釈）。蒸気物性は IAPWS‑IF97 を基準に計算。
ノズル：直管（円管）、出口径 D = 10 mm、長さ L = 100 mm（L/D = 10）。管出口で Mach = 1（チョーク）を満たす。
下流：100 Pa の真空（背圧）。
平板：管出口から 50 mm 下流、軸に垂直（ノズル‐平板間隔 H=50 mm、H/D=5）。
水噴霧：噴射位置オプション＝(a) 上流リザーバ内、(b) 管入口、(c) 管壁途中（単孔/多孔、角度任意）。
追加自由度：管内で「Laval 的」な軸方向断面積変化 A(x) をユーザー指定可能（ただし出口直径 10 mm は固定とするのが自然）。

この条件が難しい理由（設計上の論点）

超過膨張・真空膨張＋平板衝突：真空への自由噴流は早期に“ソース状（source‑like）”に近づき得る一方、凝縮、粘性、希薄化、内部自由度の凍結など、理想気体 MoC からの逸脱要因が列挙されている。
圧力比の極端さ：背圧 100 Pa に対し供給 1 bar は圧力比 (p_0/p_b \approx 1000)。自由噴流のマッハディスク位置は圧力比平方根スケーリング（係数 ≈0.67）がよく使われる。
これを当てはめると、自由噴流ならマッハディスクは 数十 D 以上下流になり得るため、H/D=5 の平板は shock‑cell/膨張領域の途中で流れを遮り、スタンドオフ衝撃波・再循環・壁面ジェットなど、典型的な超音速衝突噴流の構造が現れる可能性が高い。
希薄化の判定：自由噴流では、マッハディスク周りの現象が希薄化で“消える”目安として、マッハディスク径と平均自由行程に基づく Knudsen 数的な指標が議論されている。
→ シミュレータは 連続体（NS）だけでよいか／滑り境界が必要か／CFD‑DSMC が要るかを、ユーザーが判断できる診断（Kn 推算）を必須にする。
単一成分（H₂O）の相変化：周囲ガスが水蒸気そのものなので、一般的な「水滴が乾燥空気へ蒸発する」モデル（拡散支配）をそのまま使うと物理がずれる。低圧では界面の分子運動論（Hertz‑Knudsen‑Schrage）に基づく蒸発・凝縮フラックスや、Knudsen 層の熱抵抗が支配的になり得る点が整理されている。

装置概略と流れ領域（模式図）

（模式図は、後段のメッシュ分割・境界条件・後処理の“共通言語”として必須。）

mermaidCopyflowchart LR
  R[上流リザーバ\n飽和蒸気 p0=6 bar abs] -->|管入口| P[直管ノズル\nL=100 mm, Dexit=10 mm\nA(x)を任意指定]
  W1((水噴射\nOption A: リザーバ)) --> R
  W2((水噴射\nOption B: 管入口)) --> P
  W3((水噴射\nOption C: 管壁途中)) --> P

  P -->|出口 Mach=1| J[自由噴流（真空100 Pa）\n膨張・加速・衝撃波セル\n希薄化/凝縮の可能性]
  J --> I[衝突領域\nスタンドオフ衝撃波\n再循環/混合]
  I --> PL[平板（軸直交）\n出口から50 mm]
  PL --> WJ[壁面ジェット（放射状）\n液膜形成/再飛散/蒸発]

物理モデル設計

ここでは「何を解くか（方程式）」と「現象をどう閉じるか（サブモデル）」を、優先順位と交換可能性が分かる形で定義する。

気相（蒸気）モデル：圧縮性＋（必要なら）湿り蒸気

基礎方程式：3D 圧縮性 Navier–Stokes（質量・運動量・全エネルギ）を基本とし、乱流モデル（RANS/LES）を組み込む。外部境界での反射防止には特性（Characteristic）ベースの境界条件（NSCBC）が古典的手法として確立している。

物性：水・蒸気の物性は IAPWS の工業式 IF97 を基準にし、状態方程式（EOS）と派生量（音速、(c_p)、粘度・熱伝導等）を一貫させる。

湿り蒸気（非平衡凝縮）：飽和蒸気が急膨張すると過冷却状態となり、核生成と液滴成長により非平衡凝縮が起き得る。タービン翼列の湿り蒸気解析では、核生成率モデルと液滴成長モデルを組み込み、IAPWS‑IF97 を熱力学量に用いる枠組みが示されている。
実装の要点は以下。

支配変数：湿り度（液相質量分率）、核数密度、平均半径（または粒径分布近似）
ソース項：核生成による液相生成、成長による質量移動、潜熱のエネルギーソース
タービン分野の知見：凝縮開始位置・過冷却度・圧力分布がモデルに敏感で、国際比較でもばらつきが大きいので、モデル係数・感度解析を前提にする。

液滴モデル：Eulerian–Lagrangian（基本）＋ PBM（必要時）

本課題は「SMD と全 PSD」を要求しているため、**液滴の個別履歴を持てるラグランジュ粒子（parcels）**が最も直接的である。オープンソースでは OpenFOAM が SprayParcel （噴霧液滴の破砕・蒸発を含む）を提供している。

ただし、液滴負荷が高い（体積分率↑、二相音速↓、群としての気液連成が強い）場合は、(a) 2‑way coupling を強化する、または (b) Eulerian‑Eulerian＋人口バランス（PBM）へ段階移行する設計が必要になる。PBM の代表的なモーメント法として QMOM と DQMOM が確立している。

一次・二次破砕とレジーム判定

水噴射の入口条件（初期 PSD）が不明な場合、実務 CFD では「blob 注入＋一次破砕モデル＋二次破砕モデル」で PSD を生成する。二次破砕は Weber 数・Ohnesorge 数で破砕モード（bag / multimode / shear など）を整理するレジームマップが古くから使われ、衝撃波擾乱に対する二次破砕の整理がレビュー形式でまとめられている。
また、加速度誘起破砕の破砕時間・安定破片最大径を与える相関（いわゆる Pilch–Erdman）が広く引用される。

実装上は、ユーザーの用途に応じて複数モデルを切替可能にするのが安全である（“正解が一つ”になりにくい）。

TAB（振動アナロジ）モデル：エンジンスプレー向けの実装・文献が整備されている。
KH‑RT（Kelvin‑Helmholtz + Rayleigh‑Taylor）系：商用でも実装が体系化され、KH と RT の双方で支配不安定を評価する説明が公開されている。
Pilch‑Erdman 相関：衝撃波・高相対速度の加速度誘起破砕で破砕時間を与える。

蒸発・凝縮・熱伝達：連続体から低圧まで

重要な注意：本系は「水滴 in 水蒸気」であり、拡散支配の蒸発相関（乾燥空気中の水蒸発）を無条件に適用すると、駆動力設定（蒸気分圧差）や輸送機構が崩れる。
設計としては、圧力・Kn・相平衡から (A) 連続体・対流伝熱支配 と (B) 界面分子運動論支配を切り替える。

(A) 連続体・対流支配：工程シミュレーションでは Ranz–Marshall 型の Nu/Sh 相関が“参照相関”として扱われることが多い（引用例・利用例が多数）。
(B) 低圧・界面支配：液体/蒸気界面の蒸発・凝縮質量流束を Hertz‑Knudsen‑Schrage（Schrage）式で与える枠組み、Knudsen 層が数平均自由行程厚であること、熱伝導が律速になり得ることが NASA 技術報告で整理されている。
真空への蒸発そのものの補正式（Hertz‑Knudsenとの差）も議論されているため、将来的に高度化する余地を残す。

平板衝突・付着・再飛散（impingement outcomes）

噴霧衝突は、(i) 付着（deposition）、(ii) 跳ね返り（rebound）、(iii) スプレッド、(iv) スプラッシュ（再飛散）、(v) 液膜形成・剥離、(vi) 蒸発促進、を同時に扱う必要がある。
実装戦略としては「壁面相互作用モデル（splash 判定＋運動量/エネルギ保存）＋壁面液膜（2D薄膜）モデル」を組み合わせるのが定石で、保存則に基づくパーセル衝突モデルの一例が公開されている。

さらに、超音速衝突噴流では、マッハディスク背後の亜音速域と外側の超音速域、スタンドオフ衝撃波、再循環域、壁面ジェットが現れ、これらの非定常が圧力・速度変動と結びつくことが実験で示されている。
→ 液滴衝突モデルは、時間平均場だけでなく、非定常圧力・せん断の影響を受ける（特に splash 境界付近）前提で選ぶ。

モデル選択肢の比較表

サブシステム	モデル候補	強み	弱み/注意	推奨の使い分け
気相乱流	RANS（例：SST） / URANS / LES	RANSは安価、LESは衝撃波セル・衝突噴流の非定常を捉えやすい	LESは格子・時間刻み要求が厳しい。衝撃波＋乱流で数値散逸設計が重要	まず軸対称RANSで設計空間を絞り、重要ケースのみLESへ
気相物性	理想気体 / IF97ベース実在物性	IF97は水蒸気/水の産業標準に近い	実在物性は実装コスト増。状態領域（飽和線近傍）で頑健性が要る	基本はIF97、診断用に理想気体モードも併設
湿り蒸気	なし（乾き蒸気） / 非平衡凝縮（核生成＋成長）	飽和蒸気急膨張の逸脱要因を扱える	モデル感度が大きいことが報告されている	凝縮ON/OFFを切替可能にし、V&Vで係数同定
液滴輸送	Eulerian–Lagrangian（パーセル） / Eulerian–Eulerian＋PBM	パーセルはPSDが出しやすい。PBMは高負荷に強い	PBMは実装が重く、粒径分布の表現自由度に制約	要求が「全PSD」なので基本はパーセル。高負荷時はPBMを追加
二次破砕	TAB / KHRT / Pilch‑Erdman 等	文献・実装が多い（商用・OSS両方）	係数調整が必要になりがち	複数選択可にし、実験データで校正する
蒸発・凝縮	Ranz–Marshall型 / 界面運動論（HKS/Schrage）	連続体⇄低圧を扱える	単一成分蒸気なのでモデル適用条件の判定が必須	Kn/圧力で切替、真空域ではHKSを優先
壁面相互作用	splash判定＋保存則パーセルモデル / 液膜モデル	付着～再飛散の最小機能を満たす	表面粗さ・濡れ性・温度の入力が必要	まず保存則型＋液膜。必要に応じ粗さ/濡れ性拡張

数値解法と計算設定

解析領域分割と境界条件

領域は少なくとも 2 つに分けるのが実装上有利：
(1) 管内（A(x) を反映）、(2) 外部真空領域（自由噴流＋平板＋周辺）。

入口条件（蒸気）：上流リザーバを「全圧 (p_0=6) bar abs、全温＝飽和温度（乾き度1）」の貯気槽として扱い、圧縮性入口境界（Total pressure/temperature）を指定する。熱力学量は IF97 から与える。

外部境界（真空 100 Pa）：外部計算領域の遠方境界に (p = 100) Pa を与える。ただし超音速域では単純圧力境界が反射源になり得るため、NSCBC（特性境界）やスポンジ層（sponge）を使うのが定石。

平板：ノースリップ壁、熱境界はユーザー入力（断熱／一定温度／熱伝達境界）にする。衝突後の液膜・跳ね返りは壁面相互作用モデルへ受け渡す。

圧縮性数値スキーム（衝撃波・接触面の扱い）

基本は有限体積（FVM）＋ Riemann solver 系フラックス（HLLC/Roe 等）＋ MUSCL/TVD 制限で衝撃波を捕捉する（LESでは数値散逸管理が要点）。
音速・密度・温度が強く変化するため、エネルギ式の一貫性（熱力学整合）を最優先し、IF97 を使う場合は物性計算の反復・テーブル化・範囲外対策を用意する。

希薄化の診断とハイブリッド化（必要時）

真空自由噴流では、連続体→遷移→自由分子へ移行し得ること、また連続体から自由分子へ移っても密度場が MoC と大きく外れない場合があることが述べられている。
加えて、自由噴流のマッハディスク近傍では希薄化によりマッハディスクが“見えなくなる”指標（マッハディスク径と平均自由行程で作った指標）が議論されている。

設計要件：

シミュレータ起動時に、外部真空圧・温度スケール・代表長さ（D、H）から Kn の概算を出し、
- 連続体 NS（Kn≪1）
- 滑り境界（Kn ~ 0.01–0.1 程度）
- ハイブリッド CFD‑DSMC（遷移域以上）
  の候補を提示する。
DSMC を採用する場合、SPARTA（並列DSMC）のような既存コードを外部モジュールとして結合するのが現実的。
液滴運動も希薄域では滑り補正（drag/熱伝達）を要することがあり、Kn 依存の slip correction 形が実験論文で扱われている。

メッシュ・時間刻みの要求（実装ガイドライン）

ここは厳密な“標準値”が存在しにくいため、**推奨値は「概算指針」**として記す（実際は格子収束・時間刻み収束で裏付ける。V&V 節参照）。

外部噴流：衝撃波セルとせん断層を解像するため、少なくとも噴流直径方向に 20–40 セル相当から開始し、LESではさらに増やす（AMR が有効）。
平板近傍：壁面ジェット・再循環・熱/質量移動を扱うなら、最初は壁関数 RANS（y+ 設計）で開始し、必要に応じて低Re解像または WMLES。
時間刻み：LES/非定常衝突噴流は音速スケールの CFL 制約が支配し得る。まず URANS で周期性と統計収束時間を掴み、次に LES の Δt を設定する（時間刻み収束もV&V対象）。

ソルバー候補と実装アーキテクチャ

既存ソルバー活用の現実解

ゼロから完全実装するより、圧縮性CFD＋噴霧（DPM）＋壁面相互作用が揃った既存ソルバーを“計算エンジン”にし、周辺をシミュレータ（入力生成・実行管理・後処理・同定）として構築する方が、研究開発速度・信頼性・V&V で有利になりやすい。

オープンソース系

OpenFOAM：SprayParcel が「破砕・噴霧機能を持つ反応性スプレーパ―セル」であること、KHRT や Pilch‑Erdman 等の breakup 実装が API ドキュメントとして参照可能。
長所：ソース改変が可能、モデル切替・追加が容易。
注意：湿り蒸気（非平衡凝縮）や低圧界面相変化（HKS）などは追加実装が要る可能性が高い。湿り蒸気モデルは商用CFDにUDFで実装された事例がある。
DSMC：SPARTA（Sandia開発の並列DSMC）。低密度ガスの衝突・表面衝突を扱う設計が明記されている。

商用系（代表）

ANSYS：Fluent の理論/ヘルプ内で KHRT breakup の説明が公開されており、KH と RT の不安定に基づく破砕を扱う。
長所：DPM、壁面相互作用、収束支援、前後処理。
注意：湿り蒸気・真空相変化はUDF/カスタムが必要になりやすい（実装例あり）。
CONVERGE：噴霧（liquid parcel breakup）に KH など複数モデルがあることを公式ページが述べる。
Siemens：Simcenter STAR‑CCM+ は CFD/マルチフィジックスとして公式に位置づけられ、スプレー等の離散相（Lagrangian multiphase）を含むマルチフェーズ機能の資料が流通している。

ソルバー比較（機能観点）

候補	圧縮性衝撃波	噴霧（破砕/蒸発）	壁面衝突/液膜	湿り蒸気（非平衡凝縮）	低圧（DSMC/滑り）	拡張性
OpenFOAM	可（実装/設定次第）	KHRT/Pilch‑Erdman等の実装参照可	追加実装/既存モデル組合せ	追加実装（研究例は文献に存在）	別エンジン結合が現実的	高
ANSYS Fluent	可	KHRT説明あり	標準機能あり（範囲は設定依存）	UDF実装例あり	滑り/希薄は限定的、DSMCは別系	中
CONVERGE	可（用途次第）	KH等スプレーモデルを公式に言及	あり（用途次第）	追加実装が必要な場合	低圧は別法が要る場合	中
STAR‑CCM+	可	Lagrangian multiphaseを含む	あり（用途次第）	追加実装が必要な場合	低圧は別法が要る場合	中

望ましい内部アーキテクチャ（データ構造と結合）

「シミュレータ」としての完成度は、CFD ソルバー単体よりも 周辺のデータ管理で決まる。

コアデータ構造（例）

GeometryProfile：A(x) もしくは r(x) を B‑spline/区分線形で保持。制約（x=0…L、出口 r=5 mm）を付ける。
InjectionSpec：位置（reservoir/inlet/wall）、方向、温度、質量流量、初期 PSD（分布種別：lognormal/rosin‑rammler/ユーザーCSV）、初速度、噴霧角、連続/パルス。
SteamStateModel：IF97 呼び出し（(p,T \to \rho,h,a,c_p) 等）＋範囲外処理。
DropletParcel：多数粒子をパーセルとして束ねる（代表径、数密度、温度、速度、位置、分裂履歴、質量変化）。OpenFOAM/商用DPMと互換を意識する。
WallFilmState：平板上の液膜厚、速度、温度、蒸発量、再飛散フラックス。保存則型衝突モデルと整合。

結合方式（推奨）

2‑way coupling：粒子→気相へ運動量・エネルギ・質量ソースを反映（PIMPLE/PISO などの反復ごとに source を更新）。
強結合が必要な領域（衝突域・高液滴負荷）は、アンダーリラクゼーションや Aitken 加速で安定化。
湿り蒸気（凝縮）を入れる場合は、気相方程式のソースとして統一（液滴噴霧による相変化と“同じ会計”になるように）。

シミュレーション手順（タイムライン/フローチャート）

mermaidCopyflowchart TD
  A[入力読込\n幾何A(x), 蒸気p0, 真空pB, 水噴射仕様, 平板条件] --> B[整合性チェック\n単位/範囲/出口Mach=1成立性\nKn診断: 連続体/滑り/CFD-DSMC]
  B --> C[前処理\nメッシュ生成(管内+外部)\n初期場(準1D/等エントロピー推定)]
  C --> D[気相CFD\n圧縮性RANS/URANS/LES\n非反射境界(NSCBC等)]
  D --> E[噴霧・相変化\n注入(パーセル)\n一次/二次破砕\n蒸発/凝縮(連続体↔HKS切替)]
  E --> F[壁面相互作用\n付着/跳ね返り/スプラッシュ\n液膜更新+再飛散]
  F --> G[収束/統計評価\n定常収束 or 時間平均の統計収束]
  G --> H[後処理\nSMD・PSD(空間/時間)\n蒸発率/凝縮率\n平板付着率・再飛散率\n熱流束/圧力分布]
  H --> I[校正/同定(任意)\n実験データと比較\nモデル係数推定\n感度解析]

入出力仕様とユーザーインターフェース

未指定パラメータの扱い（ユーザー入力化＋推奨デフォルト）

ご提示条件だけでは、噴霧・蒸発・衝突の結果が桁違いに変わり得る“未指定”が残る。シミュレータは次を 必須入力として要求し、未入力時は「合理的デフォルト＋根拠」を出す。

蒸気：乾き度（x=1 を既定＝乾き飽和蒸気）、供給配管の圧力損失（無視/指定）、壁面温度（断熱/指定）。IF97に基づく状態量計算。
真空：背圧 100 Pa は指定済みだが、チャンバー寸法と排気能力が未指定（実際は噴出で圧力上昇し得る）。Cassanova らが述べるように、真空噴流では連続体→自由分子への遷移が関与し得るため、チャンバー条件も入力にするのが望ましい。
水噴霧：
- 水温（例：室温～飽和近傍）、水供給圧、質量流量（蒸気に対する比が重要）
- 注入方式（単孔ジェット、スワール、エアアシスト相当など）
- 初期液滴径分布（未知なら「ノズル径に基づく blob＋一次破砕」）
平板：温度、熱伝導（材質）、表面粗さ/濡れ性（接触角）、液膜の有無（乾き/濡れ）。保存則型壁面モデルは dry/wet でモデルが変わることが明記されている。

出力（要求仕様に直結するもの）

粒径分布：ASTM の液滴径分布データ処理手順（サンプル数、階級幅、代表径、分散の扱い）を参照し、“全PSD”の定義（数基準/質量基準/表面積基準）を UI 上で明示する。
計測互換として ISO 13320（レーザ回折）ベースの出力（体積分布 etc.）も提供すると実験比較が楽になる。

推奨出力セット（最小）

噴流：中心線 (p,T,M,\rho)、衝撃波セル推定、壁面圧力分布（平板）
液滴：
- 空間分布：SMD（D32）、D10/D50/D90、分布ヒストグラム（標準化）
- レジーム：We, Oh, Re から破砕モード分類（bag/multimode/shear 等）
- 蒸発/凝縮：質量変化率、潜熱寄与（エネルギ収支）
衝突：付着率、再飛散率、液膜厚分布、二次飛散 PSD（もしモデル化するなら）

UI 機能（研究用に効く“差分”）

A(x) 入力：GUI での編集（区分線形・B‑spline）＋CSV/JSON import。制約（出口径固定、最小半径、曲率制限）をチェック。
噴射位置ウィザード：reservoir / inlet / wall‑ports をテンプレで作れる。
パラメータスイープ：
- サンプリング（Latin Hypercube）を標準装備し、感度解析・UQ に繋ぐ。
- 外部ツール連携：Sandia National Laboratories の Dakota（最適化/UQ）や NASA 系 OpenMDAO（MDAO）で、CFD実行をブラックボックスとして回す構成が取りやすい。
- Python では SALib（Sobol/Morris/FAST）が使える。

検証・妥当性確認と優先文献・データセット

V&V の枠組み（何を“正しい”とするか）

CFD の信頼性は、(i) **Verification（方程式を正しく解いているか）**と、(ii) **Validation（現実を正しく再現できるか）**に分けて管理するのが標準的である。ASME V&V 20 は CFD/伝熱のV&V枠組みを規定し、AIAA のガイドは格子・時間刻みの外挿（Richardson）等の実務指針を与える。

Verification（コード・解の検証）

コード検証：製造解（MMS）で、圧縮性NS＋ソース項（粒子カップリング、相変化）の実装が収束次数通りに落ちるか確認（単体テスト）。
解検証：格子独立性（少なくとも3段階）、時間刻み独立性（非定常なら）を実施し、主要 QoI（平板力、付着率、SMD、蒸発量）の数値不確かさを推定する。

Validation（妥当性確認）：優先すべき実験・文献

本件は「全現象を一度に」検証するのが難しいため、分解して段階的に妥当化する。

真空自由噴流・平板（気相のみ）

Cassanova & Stephenson（1967）：“jet into vacuum” の噴流‐平板実験（平板に働く力・流線ベーン）と MoC との比較、逸脱要因（凝縮・粘性・希薄化・内部自由度凍結）の列挙が本件に直撃。
Ashkenas & Sherman（NASA, 1965）：自由噴流の構造、マッハディスク位置のデータ整理、希薄化でマッハディスクが消える指標の議論。

超音速衝突噴流（気相・流れ構造）

Henderson ら（NASA NTRS, 2002）：衝突領域がマッハディスク背後の亜音速域と外側超音速域に分かれること、再循環域や壁面ジェットの非定常が示されている。
Orescanin ら（2014）：自由噴流と衝突噴流でマッハディスクの振る舞いが変わり得る点を扱う。

液滴二次破砕（単一液滴 or 希薄スプレー）

Pilch & Erdman（1987）：加速度誘起破砕の破砕時間・安定破片最大径推定相関。
二次破砕レジーム整理：We/Oh に基づく破砕モード整理がレビューにまとめられている。
TAB 法：エンジンスプレーでの破砕モデルとして提示。
KHRT：商用の理論説明が公開されており、実装比較の基準にできる。

湿り蒸気・非平衡凝縮（蒸気ノズル）

湿り蒸気 CFD の方法と方程式形：J‑STAGE公開の蒸気タービン翼列解析で、核生成・液滴成長・IAPWS‑IF97 の利用が明示されている。
IWSMP（国際比較）：凝縮ノズルケースで計算結果のばらつきが報告され、モデル感度の大きさを示す。
近年のノズル実験データ（例：Barschdorff形状の測定）：凝縮開始位置を高分解能で測る試験が報告されている。

壁面衝突（スプレー/液膜）

保存則ベースの壁面衝突モデル（dry/wet の splash 判定を含む）実装例が公開されている。

スプレーモデルの校正用に便利な公開データ（参考）

本件は蒸気＋真空で特殊だが、噴霧の“モデル係数”を同定するには、豊富な公開データセットで基礎校正するのが現実的。Engine Combustion Network は Spray A 等の実験データ検索・推奨診断指標を公開している。

計算コスト見積り（小／中／大：概算）

以下はオーダー見積り（ハードウェア・並列度・収束条件で大きく変動）。V&Vに基づき、QoIごとに必要解像度を決めて更新する前提。

忠実度	目的	代表モデル	メッシュ規模（目安）	時間積分	計算資源（目安）
小	スイープ・傾向把握	2D軸対称RANS＋パーセル（1‑way）＋簡易破砕/蒸発	0.3–2 M cells	定常 or URANS短時間	WS～小型サーバ（数時間～）
中	PSD/衝突の実用予測	3D URANS/簡易LES＋2‑way＋壁面衝突＋液膜	5–50 M cells	非定常（統計平均）	計算クラスタ（数日～）
大	非定常衝突噴流＋希薄化/相変化の高忠実度	3D LES（必要ならCFD‑DSMC分割）＋HKS相変化＋高度壁面モデル	100 M cells 以上	非定常（長時間統計）	HPC（週オーダー～）