マッハ円錐 シミュレーター — 超音速飛行のソニックブーム

マッハ円錐は、超音速で運動する点源が引きずる球面波の包絡面（ホイヘンス的構成）として定義されます。Ernst Mach（1838-1916）と Peter Salcher が 1887 年にウィーンで発射体（弾丸）の影絵写真を撮影し、衝撃波面を初めて可視化しました。

飛行体が時刻 $-t$ に位置 $(-Vt, h)$ にいたとすると、その瞬間に発生した球面波は時刻 $0$ までに半径 $ct$ に広がります。一方、飛行体は時刻 $0$ に原点 $(0, h)$ にいます。波面の包絡面が円錐となり、その半頂角 $\mu$（マッハ角）は次式で与えられます。

$$\sin\mu = \frac{ct}{Vt} = \frac{c}{V} = \frac{1}{M}$$

ここで $M = V/c$ がマッハ数です。$M < 1$ では右辺が 1 を超え、解が存在しません（円錐は形成されない）。$M = 1$ で $\mu = 90°$（平面波）、$M > 1$ で $\mu < 90°$ の鋭い円錐となり、$M \to \infty$ で $\mu \to 0$ に漸近します。

飛行体が高度 $h$ で水平飛行しているとき、円錐の側面が地表 $y=0$ と交わる点の、飛行体直下からの水平距離 $\ell$ は幾何学的に

$$\ell = \frac{h}{\tan\mu} = h\sqrt{M^2 - 1}$$

第二式は $\sin\mu = 1/M$, $\cos\mu = \sqrt{1-1/M^2}$ から $\tan\mu = 1/\sqrt{M^2-1}$ を代入して得られます。観測者が飛行体直下にいる場合、ソニックブームは飛行体が真上を通過してから $t_{\text{delay}} = \ell/V = h/(M\,c\,\tan\mu)$ 秒遅れて到達します。本ツールの既定値 $M=2$, $c=343$ m/s, $h=1000$ m では $\ell=1732$ m, $t_{\text{delay}}=2.52$ 秒となります。

実際のソニックブームの圧力波形は典型的に「N 字波形」と呼ばれ、機首衝撃波（前方衝撃）で急激な圧力上昇 $+\Delta p$、機体後方の膨張で滑らかに $-\Delta p$ まで降下、機尾衝撃波で再び急上昇して大気圧に戻る形になります。$\Delta p$ は機体長・揚力・飛行高度に依存し、F-18（M=1.4, h=10km）で約 50 Pa、コンコルド（M=2, h=18km）で約 100 Pa が典型値です。

超音速航空機の設計と認証：FAA・EASA・JCAB は陸上での超音速飛行を Δp 制限（米国は事実上禁止）で規制しています。新世代 SST（NASA X-59、Boom Overture）は機体全長を 30m 以上に延ばし、機首・主翼・機尾の衝撃波を分散させる「低ブーム形状」で Δp を 25 Pa 以下に抑えます。これらの設計は CFD（OVERFLOW、Cart3D）と圧力波伝播コード（PCBoom、ZEPHYRUS）を組み合わせて最適化されています。マッハ円錐の幾何学はこの低ブーム解析の出発点です。

大気圏再突入：スペースシャトルや HTV（こうのとり）の再突入時、機体は M=25 程度で大気に突入し、極めて鋭いマッハ円錐に弓形衝撃波（bow shock）が重なります。CFD では DSMC（直接シミュレーションモンテカルロ）と化学非平衡 Navier-Stokes を併用し、表面熱流束（ピーク 1MW/m²）を予測します。マッハ円錐の理論はこの解析の境界条件を決める基礎です。

弾丸・ライフルの内弾道 / 外弾道：軍事・狩猟用ライフル弾は M=2〜4 の超音速で飛翔し、地表に接近すると円錐状の衝撃波が「クラッキング音」として聞こえます。Ernst Mach の最初の実験はまさにこの弾丸の影絵写真でした。現代では狙撃手の音響探知システム（Boomerang、PILAR）がこのマッハ円錐の到達時間差から発射点を逆算し、市街戦での対策に使われています。

キャビテーション・水中発射：水中での超音速はほぼ不可能ですが、超キャビテーション魚雷（露 VA-111 シュクヴァル）は機首から発生するガスでキャビティを形成し、水中で M=0.07 相当（200 km/h）の高速航行を実現します。気液境界面の挙動はマッハ円錐の理論とも類似しており、ホモジニアス・モデルによる解析対象です。

最も多い誤解は 「ソニックブームは音速突破の瞬間だけ起きる」というものです。実際には超音速飛行中ずっと円錐が引きずられているので、地表のブーム軌跡（footprint）は飛行体が通過した経路全体に沿って残ります。一直線飛行では幅 80km 程度の細長いベルト状になります。1947 年のチャック・イェーガー初の音速突破飛行で「音速の壁」という言葉が生まれましたが、壁を破って一度大きな音がするわけではなく、音速を超えている間ずっと地上で観測されます。

次に多いのが 「マッハ円錐は機体形状によらない」という思い込みです。点源の理論ではそうですが、実際の機体は機首・主翼前縁・キャノピー・主翼後縁・機尾の各所から個別の衝撃波を発生させ、地表ではこれらが N 字波形に合成されます。NASA の SonicBAT 飛行試験（2017 年、F/A-18B）では、地表 25 点での Δp 測定値と CFD 予測が ±5% 以内で一致することが確認されています。マッハ円錐は「主衝撃波の位置」を与えますが、波形の詳細には CFD が必須です。

最後に、「マッハ角は飛行高度に依存する」という誤解です。式 sin μ = 1/M は機体周りの局所マッハ数のみで決まり、高度は直接は関係しません。ただし高度によって音速 c が変わる（標準大気で地表 c=340 m/s、高度 11km で c=295 m/s）ため、同じ対地速度 V でもマッハ数 M = V/c が変わり、結果としてマッハ角も変わります。本ツールでは音速 c をスライダーで独立に変えられるので、高度による音速変化の影響を疑似的に体験できます。