エグゼクティブサマリー
液柱を一様気流に直交噴射して微粒化させる問題で、CFD/VOF/Level-Set/CLSVOF/VOF-LPT の検証に本当に使いやすい一次実験データは、単一論文ではなく、役割の異なる複数系列を組み合わせるのが最も実務的です。近接場の液柱形状と波動・衝撃波相互作用には、Mach 1.94 水噴流の古典系列と、それを拡張した高速度 X 線計測が有効です。速度場には Mach 2.1 水噴流の PIV、粒径・粒子速度には Mach 1.86 ケロシン PDPA、Mach 2.0 キャビティ内 PDA、Mach 2.85 水噴流の PDA 系列が特に有用です。さらに、公開 PDF で再現しやすいケースとして、Purdue University の Mach 1.94 エタノール噴流実験は、噴霧侵入深さ・幅・不確かさの整理が非常に良く、いま一番 “手を動かして回せる” ベンチマーク群です。
一方で、完全に公開された機械可読の生データは、このニッチではまだ少なく、公開入口の多くは論文本文、図、補助資料、あるいは受理原稿への導線です。したがって、実務上は、まず公開 PDF・要旨・受理原稿で境界条件と観測量を固定し、そのうえで図のデジタイズと著者問い合わせを組み合わせる方針が現実的です。今回確認できた範囲でも、オープン PDF が明確に取れる代表例は Purdue の ICLASS 2021 論文で、受理原稿入口は OSTI.GOV にあり、X 線論文には Argonne National Laboratory の publication page があります。
1960年代1966初期の超音速二次噴射研究1970年代1971超音速中の液膜・液柱breakup1980年代1980波動現象と平均滴径1985 plume 横断分布2000年代2002 高速 crossflowpenetration2004 Mach 1.94水噴流構造2004/2007亜音速の正準breakup 基準2010年代2014 Ma 2.1 PIV速度場2015 Ma 1.86 PDPA滴径分布2019 Ma 2.85 PDA系列2020年代2020 X線診断2021 断面滴分布2021 オープンなMach 1.94エタノール像計測2024 キャビティ内PDA液柱噴流 in 高速・超音速気流 主要研究の流れShow code
上の流れで見ると、古典研究は penetration と平均滴径、2010 年代の研究は速度場と滴径空間分布、2020 年代の研究は X 線・キャビティ内計測・時系列統計へと重心が移っています。検証用の観測量も、単なる軌跡から、SMD 分布、断面形状、速度 PDF、液相 path length へと広がっています。
優先度付き文献とデータ源
優先順位は、原著性、定量量の豊富さ、再現に必要な境界条件の明瞭さ、公開入手性の順で付けています。
- 最優先: “Experimental investigation on structures and velocity of liquid jets in a supersonic crossflow” (2014) — 水噴流を Ma = 2.1 交差流へ噴射し、PIV で構造と速度を計測した原著です。対称面速度場、近傍表面波、遠方渦構造を同時に押さえており、液柱初期加速、波動発生、近接場速度回復の検証に最適です。要旨では、x/d < 15 で噴霧速度が主流速度の 66% まで急増すると報告されています。テスト部は長さ 200 mm、断面 60 mm × 40 mm、ノズル径 0.5 mmです。url公式ページhttps://doi.org/10.1063/1.4893008
- 最優先: “Experimental study of spray characteristics of liquid jets in supersonic crossflow” (2019) — Ma = 2.85 の水噴流に対して、PDA/PDPA 系で滴径と 2 速度成分を測った代表的な原著です。中心対称面の滴径・速度分布に加え、平板壁とキャビティ前方噴射の比較があり、VOF-LPT や Eulerian-Lagrangian の far-field 検証に非常に使いやすいです。公開要約では、さまざまなノズル径とキャビティ効果を扱うことが読めますが、各ノズル径の正確な値は要約では未記載です。url公式ページhttps://doi.org/10.1016/j.ast.2019.105426
- 最優先: “Cross-sectional droplets distribution of a liquid jet in supersonic crossflow” (2021) — 2019 系列をさらに進め、断面形状と断面内粒径分布に焦点を当てた原著です。断面がΩ形を示し、高さと幅がオリフィス径とともに線形増加し、断面滴径は 0–30 μm の範囲で、主として横断方向位置に依存すると要旨が述べています。断面 spray-foot 形状の比較には、現時点で最も便利な supersonic baseline 群です。url公式ページhttps://doi.org/10.1016/j.actaastro.2021.05.024
- 最優先: “Experiments investigation on atomization characteristics of a liquid jet in a supersonic combustor” (2024) — Ma = 2.0、全圧 0.55 MPa、全温 300 K の supersonic combustor/cavity 条件で、PDA によりキャビティ内部の滴径と速度を測った原著です。SMD は 30–55 μm、平均 streamwise 速度は −20〜150 m/s、平均 vertical 速度は −20〜30 m/s と明示され、キャビティ内再循環と ignition-friendly zone の議論まで含みます。キャビティ付きスクラムジェット模擬器の検証には第一候補です。url公式ページhttps://doi.org/10.1063/5.0204890
- 高優先: “Characterization of the steady and unsteady spray structures of a liquid jet in supersonic crossflow” (2021) — Misentropic = 1.94、エタノール、d = 1 mm、L/D = 12、J = 3–13、20 kHz の backlit imaging、Mie scattering、shadowgraphy を用いたオープン PDF の会議論文です。x = 10 cm, 15 cm における penetration 深さの線形回帰やspray widthの図が公開されており、図表からそのままベンチマーク化できます。再現性重視なら、今すぐ動かしやすいケースです。
- 高優先: “Exploration of Water Jets in Supersonic Crossflow Using X-Ray Diagnostics” (2020) — 近接場の高速度 X 線 imaging と path-length-integrated X-ray fluorescence を使った原著で、dense near-field を光学法より頑健に捉えられるのが最大の利点です。後続の高忠実度シミュレーション論文でも、この実験がX 線で妥当性確認された基準として使われています。近傍液柱形状や projected liquid mass/path length を検証したい場合の最有力候補です。url公式ページhttps://doi.org/10.1615/AtomizSpr.2020034448
- 高優先: “Structures of Water Jets in a Mach 1.94 Supersonic Crossflow” (2004) — Mach 1.94、水、L/d₀ = 20、d₀ = 0.5 mm と 1.0 mm の古典系列で、normal injection、広い J、freestream velocity、aeration level を振った実験です。以後の Mach 1.94 研究の“祖形”で、形状・近接場 breakup 機構の基準として外せません。url公式ページhttps://doi.org/10.2514/6.2004-971
- 高優先: “Investigations on the droplet distributions in the atomization of kerosene jets in supersonic crossflows” (2015) — Ma = 1.86 のケロシン噴流に対する PDPA の原著で、SMD 分布が C 型から I 型へ変化すること、大滴 breakup と小滴 coalescence が atomization 全域で共存することを示しています。燃料液での粒径分布を検証するには重要です。url公式ページhttps://doi.org/10.1063/1.4930817
- 補助的だが重要: “Breakup of Round Nonturbulent Liquid Jets in Gaseous Crossflow” (2004) と “Primary Breakup of Turbulent Round Liquid Jets in Uniform Crossflows” (2007) — どちらも亜音速ですが、液物性、ノズル長さ、Re/We/Oh の整理が非常に丁寧で、コード検証の checkout case として圧倒的に有用です。超音速本番へ行く前に、breakup length、ligament/drop size、drop velocity、regime transition をここで合わせる価値が大きいです。urlnonturbulent 公式ページhttps://doi.org/10.2514/1.3749 ／ urlturbulent 公式ページhttps://doi.org/10.2514/1.19397
- 歴史的核: 1966–1985 の超音速原著群 — “Secondary liquid injection into a supersonic airstream” (1966)、“Breakup of liquid sheets and jets in a supersonic gas stream” (1971)、“Wave phenomena in liquid jet breakup in a supersonic crossflow” (1980)、“Distributions across the plume of transverse liquid and slurry jets in supersonic airflow” (1985) は、penetration、平均滴径、波動、plume 横断分布の出発点です。現代の再現ケースとしては情報不足がちですが、どの指標が古くから本質的かを知るには重要です。
実験条件の比較
下表は、シミュレータ検証に使いやすい順に近い supersonic 群と、事前コード検証に強い subsonic 正準群を並べたものです。
※「公開要約では未記載」は、今回確認できた公式アブストラクト・オープン PDF・公開スニペットの範囲で値を特定できなかった、という意味です。
| 事例 | 液体・ノズル | 気流条件 | 無次元数 | 計測法 | 公開される主な定量量 | 再現メモ |
|---|---|---|---|---|---|---|
| Mach 2.1 水噴流 PIV (2014) | 水、d = 0.5 mm、低乱れノズル、テスト部 200 mm, 60 × 40 mm | Ma = 2.1 | J は公開要約では未記載 | PIV | 対称面の速度場、近接場波、遠方渦、x/d < 15 で 0.66 U∞ まで加速 | 速度検証に最適。液物性は水で扱いやすい。 |
| Mach 2.85 水噴流 PDA (2019) | 水、複数ノズル径（正確値は公開要約では未記載）、平板壁・キャビティ前方 | Ma = 2.85 | 公開要約では未記載 | PDA/PDPA, 高速度撮影 | 滴径、2 速度成分、壁・キャビティ効果 | Far-field 粒子統計に強い。 |
| Mach 2.85 断面滴分布 (2021) | 液体噴流、オリフィス径可変（正確値は公開要約では未記載） | supersonic crossflow | 公開要約では未記載 | PDA | 断面が Ω形、高さ・幅の径依存、滴径 0–30 μm | spray-foot / cross-section の benchmark family。 |
| キャビティ内 PDA (2024) | 燃料液体、ノズル詳細は公開要約では未記載 | Ma = 2.0, p₀ = 0.55 MPa, T₀ = 300 K | 公開要約では未記載 | PDA | SMD 30–55 μm, uₓ = −20〜150 m/s, uᵧ = −20〜30 m/s, cavity 内再循環 | combustor/cavity 検証に非常に有用。 |
| Purdue open case (2021) | エタノール 100%, d = 1 mm, L/D = 12, T_liq ≈ 290 K | M_isen = 1.94, p_stat = 0.319–0.434 bar, air mass ≈ 0.3–0.4 kg/s | J = 3.95–13.17 | 20 kHz backlit, Mie scattering, shadowgraph | penetration depth, spray width, downstream growth, 2σ 不確かさ、回帰式 | PDF がオープン。今すぐ再現しやすい。 |
| Mach 1.94 水噴流古典系列 (2004) | 水、d = 0.5, 1.0 mm, L/d₀ = 20, flush-mounted normal injection | M = 1.94, total temperature 最大 533 K | J, freestream velocity, aeration level を広く掃引 | 画像計測、PDPA 系の後続研究と連結 | 近接場構造、penetration、古典的 breakup 機構 | Mach 1.94 系列の基準。 |
| X 線近接場 (2020) | 水、flush injection | blowdown supersonic wind tunnel、正確な Mach は今回の公開スニペットでは未特定 | 公開スニペットでは未記載 | high-speed X-ray imaging, X-ray fluorescence | 近接場形状、投影 liquid path length / 分布 | dense near-field に最適。 |
| 亜音速 nonturbulent baseline (2004) | 水・エタノールほか、広範囲の液物性 | 空気、NTP | We_g 0–200, We_l 0–1.6×10⁶, Re_l 0–8×10⁵, ρ_l/ρ_g 700–1100, Oh 6×10⁻⁵–0.3 | pulsed shadowgraph | regime, onset time, column/surface wavelength, breakup | 超音速前の code checkout に最適。 |
| 亜音速 turbulent baseline (2007) | 水・エタノール、d = 1, 2 mm, L/D > 100 | 空気、98.8 kPa, 298 K, u_g = 11–142 m/s | 水では Re_l = 3800–59000, We_g = 0–159, q = 3–200 | pulsed shadowgraph, holograph | breakup length, ligament/drop size, drop velocity, breakup rate | 最も物性と無次元数が揃った正準 case。 |
推奨ベンチマークケース
まず、典型的な supersonic crossflow 実験構成を図にすると、だいたい次のようになります。これは Mach 1.94–2.85 級の代表系列を抽象化したもので、CD ノズルまたは blowdown facility で一様高速流を作り、壁面 flush オリフィスから液柱を横断噴射し、近傍ではボウショックと表面波、下流では液滴群と recirculation/counter-rotating structures を測ります。
貯気槽 / Settling chamberCDノズル or blowdown section一様高速気流 M≈2壁面オリフィス液柱ボウショック表面波・ligament液滴群 / spray plume計測面 x/d, y/d, z/dPIV / PDA / PDPA / backlight / Mie / X-rayShow code
実際にベンチマークとして推すケースは、次の 5 つです。
- SJ-PIV-21
水、d = 0.5 mm、Ma = 2.1、テスト部 200 mm・60 × 40 mm。
まずこれを使って、近接場の液柱加速、x/d に沿った平均速度回復、表面波と far-field vortices を合わせるのがよいです。評価量は、中心面での (u_l/U_\infty)、主流方向・横方向の液相速度、噴霧輪郭、x/d<15 での加速率です。J、p₀、T₀ は今回確認した公開要約では未記載なので、厳密な single-shot benchmark というより、速度場 benchmark として使うのが適しています。 - SJ-OPEN-194-J6
エタノール、d = 1 mm、L/D = 12、(M_{\mathrm{isen}} = 1.94)、(J \approx 6.05)。
x = 100 mm 面での penetration を図から直接比較でき、さらに x = 60–200 mm の downstream growth も追えます。Purdue の open PDF では、10 cm 面での penetration 深さの回帰が (y \approx 1.086J + 8.755) mm、15 cm 面では (y \approx 0.973J + 14.098) mm と読めます。case-by-case の静圧は図上の色分けで示される一方、各 run の完全な CSV は公開 PDF には付属していません。侵入深さの一次検証に非常に良いです。 - SJ-OPEN-194-J12.8
上と同一装置で (J \approx 12.8) を選ぶケースです。
J を二段階に絞るだけで、penetration slope の J 依存、spray depth の不確かさ、spray width の下流増大が比較できます。VOF/CLSVOF では liquid core の侵入と spread、VOF-LPT では handoff 後の plume width を見るのに向いています。 - SJ-PDA-20-Cavity
Ma = 2.0、(p_0 = 0.55) MPa、(T_0 = 300) K の cavity-based combustor。
粒径・速度の定量比較を主目的にするなら、この case family は強いです。SMD 30–55 μm、(u_x=-20 \sim 150) m/s、(u_y=-20 \sim 30) m/s が要旨で明示されます。逆に、ノズル形状、液物性、cavity geometric detail の exact value は、今回確認した要約だけでは不足しているので、strict single-case benchmark ではなく family benchmark とみなすのが安全です。 - SB-WATER-1MM-NTP
水、d = 1 mm、L/D > 100、空気 98.8 kPa・298 K、(u_g = 11–142) m/s。
超音速そのものではありませんが、Re/We/Oh が揃い、breakup length・ligament/drop size・drop velocity が揃うので、本番 supersonic case に入る前の code checkout として不可欠です。特に、界面捕捉法の breakup onset、droplet identification、velocity statistics の妥当性確認に使えます。水の物性は (\rho_l=995) kg/m³、(\mu_l=8.94\times 10^{-4}) kg/m·s、(\sigma=70.8) mN/m、Re_l = 3800–59000、We_g = 0–159、q = 3–200 です。
結論だけ言うと、「近接場形状」= Mach 1.94 水/X-ray、 「速度場」= Ma 2.1 PIV、 「粒径・速度統計」= Ma 2.0 / 2.85 PDA、 「誰でも回せる公開ケース」= Purdue 2021、 「コードの基礎検証」= 亜音速 Sallam/Lee という役割分担が最も実務的です。
検証指標と後処理
検証指標は、実験が見ているものに合わせて定義しないと意味がありません。supersonic LIJCF/LJSC では、penetration、breakup length、SMD、断面形状、速度 PDF、液相投影量が主役です。PDA 系は粒径と速度、PIV は液相追跡可能な速度場、backlit/Mie は閾値化された輪郭、X 線は投影された液体 path length に対応します。したがって、シミュレーション後処理でも、計測器と同じ “見え方” を模擬するのが大事です。
| 指標 | 実験対応 | シミュレーションでの計算方法 | 実務上の注意 |
|---|---|---|---|
| penetration depth / trajectory | backlit, shadowgraph, Mie | 計測面・投影方向を実験に合わせ、液相体積率 (\alpha_l) を line-of-sight に投影して閾値化し、外縁座標 (y(x)) を取る | 3D 体積の“真の外縁”ではなく、画像と同じ閾値輪郭を作ること。 Purdue 2021 はこの比較に最適。 |
| breakup length (L_b) | column breakup location, end of continuous core | injector と連結した液相連結成分が途切れる最遠軸方向位置を定義し、時間平均する | threshold と connected-component 条件に敏感。 subsonic Sallam/Lee 系で先に定義を固めるのが安全。 |
| SMD / (D_{32}) | PDA/PDPA | 離散滴では (D_{32}=\sum d_i^3/\sum d_i^2)。VOF では detached liquid component を球等価直径へ変換して同式を使う | number-weighted か volume-weighted かを明示する。 2024 cavity case と 2019/2021 Mach 2.85 系列で直接比較しやすい。 |
| 粒径ヒストグラム / PDF / CDF | PDA/PDPA, laser particle analyzer | 実験と同じ bin 幅で粒径分布を作る。必要なら number-weighted と mass-weighted を両方出す | instrument sample volume と小粒径 cut-off を模擬する。 kerosene Ma 1.86 では C→I 型分布遷移が重要。 |
| 液滴速度平均・PDF | PIV, PDA/PDPA | 計測面内の液滴または detached component の重心速度を集計し、(u_x,u_y,u_z) の平均と PDF を出す | far-field では Eulerian 相平均と Lagrangian 粒子平均を混同しない。 Ma 2.1 PIV と 2024 cavity PDA が基準。 |
| 断面幅・高さ・Ω/C/I 形状 | Mie scattering, PDA 断面計測 | 特定 x 断面で (\alpha_l) あるいは液滴数密度を投影し、輪郭幅・高さ・形状指標を取る | 2021 Acta の Ω 形 benchmark と、Purdue 2021 の width-vs-J が使いやすい。 |
| 波長・ligament 統計 | 高速度可視化、X-ray | 界面抽出後、windward surface の波峰間距離、ligament 長さ・太さ・向きを時系列統計化する | 近接場は X 線計測が強い。解像度不足だともっとも壊れやすい指標。 |
| 投影 liquid path length / extinction | X-ray radiography, extinction diagnostics | (\int \alpha_l,dl) や Beer–Lambert 型透過率を模擬して実験画像と直接比較する | dense near-field では binary edge 比較より robust。 X-ray と extinction 論文が直接の対応先。 |
| mass flux profile | isokinetic sampling, PDA/LPA の面内分布からの再構成 | 断面セルごとに (\dot m”=\frac{1}{A\Delta t}\sum \rho_l V_i u_{n,i}) を構成し、面内 map を出す | 直接 mass flux がない論文では、粒径・速度・数密度からの再構成になるので、仮定を明示する。 historical papers と PDA 系で有効。 |
実装上の要点は、**“raw CFD output をそのまま比べない”**ことです。たとえば、Purdue の輪郭比較なら synthetic backlight / synthetic Mie を作り、X 線論文と比べるなら projected path length に落とし、PDA 論文と比べるなら measurement volume を想定した plane sampling で binning を合わせるべきです。これをやるだけで、モデル誤差と計測定義のズレをかなり切り分けられます。
数値上の難所と解像度・不確かさ
supersonic liquid-jet atomization の数値計算が難しい理由は、圧縮性衝撃波、界面張力、一次 breakup、二次 breakup、壁面近傍輸送、場合によってはキャビティ再循環が同時に出るからです。Ma 2.1 PIV 論文は、近傍表面波と遠方渦が同時に効くことを示しており、Purdue 2021 は injection point でunsteady bow shockを観察しています。近年のレビューでも、liquid column breakup を主に Rayleigh–Taylor 型の windward wave 発達、surface/droplet breakup を Kelvin–Helmholtz 型の周方向せん断不安定が支配すると整理されています。
さらに、解像度の影響はかなり大きいです。Kuhn–Desjardins の 2022 論文要旨は、jet shape と surface-instability wavelength の収束性を resolution study で確認したことを明示しています。つまり、少なくとも supersonic LJSC では、軌跡が合うだけでは不十分で、波長や breakup 由来の統計量が収束しているかまで見ないと危険です。
そのうえで、実務向けの推奨指針を整理すると、次のようになります。
- 界面解像度
penetration と平均輪郭だけを目的にするなら、ノズル径あたり 32–64 cells を最低線に置き、近接場の波長・ligament・droplet production まで狙うなら 64–128 cells/d を推奨します。これは文献の resolution sensitivity を踏まえた実務上の推奨値です。 - AMR の当て先
少なくとも、界面 (|\nabla \alpha_l|)、衝撃波 (|\nabla p|)、高せん断 / 高渦度領域の三つに refinement criterion を置くのが有効です。supersonic case では bow shock と shear layer の両方を同時に追わないと、liquid column の deformation と droplet formation が壊れやすいです。 - 時間刻み
compressible solver では、少なくともshock/interface 近傍の gas-phase CFL を 0.2–0.5 程度以下に抑え、さらに capillary time-step 制約も同時に満たす設定が安全です。もし VOF-LPT handoff を使うなら、handoff 前後で mass conservation の監視を必須にしてください。 - VOF-LPT handoff
detached liquid structure が少なくとも数セル以上で直径表現でき、形状がほぼ球状に移った段階で handoff し、その閾値を変えて統計量が不変か確認するのがよいです。粒径 PDF を本気で合わせるなら、この handoff criterion が実はかなり効きます。 - 統計サンプリング
Purdue 2021 では、penetration 抽出に平均 836 枚で収束し、図では 2,000 枚規模の統計と 2σ error bar を出しています。したがって、シミュレーション側も、数百〜数千スナップショットの block average を前提にした方が安全です。瞬時場 1 枚で比較すると、実験ノイズより大きいサンプリング誤差を自分で作ってしまいます。
不確かさ評価は、離散化誤差, モデル誤差, measurement-matching 誤差の 3 層に分けるのが扱いやすいです。離散化誤差には少なくとも 3 段階格子の GCI、モデル誤差には breakup/handoff/subgrid-capillary の感度解析、measurement-matching には閾値・bin 幅・plane 位置・sampling window の摂動を使います。粒径ヒストグラムや速度 PDF には bootstrap、penetration/width には block-average と信頼区間をつけると、実験側の 2σ 表示と比較しやすくなります。
公開資料と入手法
入手しやすさの観点では、まず Purdue 2021 のオープン PDF を起点にするのが最も実用的です。装置、撮影条件、ノズル径、L/D、J 範囲、静圧範囲、回帰式、spray width の点群が無料で取れます。
次に、数値検証側の実務資料としては、Kuhn–Desjardins 2022 の OSTI 記録が有用です。ここでは、accepted manuscript への入口が明示されており、X 線計測に対する高忠実度計算の位置づけ、解像度感度、評価対象の整理が分かります。
近接場 X 線データを追うなら、Argonne National Laboratory の publication page も入口になります。そこから原著と著者情報に辿れるので、raw radiography / fluorescence frames の有無を問い合わせる際の起点として使えます。
その他の主要原著への公式導線は次の通りです。Mach 1.94 古典系列は url2004 Mach 1.94 水噴流https://doi.org/10.2514/6.2004-971、**Ma 2.1 PIV** は url2014 Ma 2.1 PIVhttps://doi.org/10.1063/1.4893008、**Ma 1.86 ケロシン PDPA** は url2015 Ma 1.86 PDPAhttps://doi.org/10.1063/1.4930817、**Ma 2.85 PDA** は url2019 Ma 2.85 PDAhttps://doi.org/10.1016/j.ast.2019.105426、**断面滴分布** は url2021 断面滴分布https://doi.org/10.1016/j.actaastro.2021.05.024、**キャビティ内 PDA** は url2024 cavity PDAhttps://doi.org/10.1063/5.0204890、**X 線原著** は url2020 X線診断https://doi.org/10.1615/AtomizSpr.2020034448 です。
今回確認した公開入口の多くは、論文本文、受理原稿、図表、institutional publication page であり、Zenodo / Mendeley Data / Dryad 型の独立した機械可読 benchmark repository は、少なくとも主要 supersonic 正準ケースでは目立ちませんでした。したがって、生データが必要なら、著者へ依頼する際に、CSV 化した penetration curve、PIV ベクトル場、PDA の bin 定義と sample-volume 情報、X 線校正情報、ノズル断面図、境界層厚さ、各 run の (p_0, T_0, \dot m_g, \dot m_l) 時系列をセットで求めるのが効率的です。公開 PDF だけで始めるなら、Purdue 2021 と 2024 cavity PDA が最も着手しやすく、近接場の厳密化には X 線原著を足す、という順番がおすすめです。
コメント