パラメータ設定
完全気体・断熱・1次元・定常を仮定。M_1n = M_1·sin(β) \gt 1 のときのみ衝撃波として意味があり、それ以外は `—` 表示になります。
楔まわりの斜め衝撃波
下=水平流入の上流(寒色・矢印長=高速)/斜線=衝撃波(角度 β)/楔表面に沿って θ 偏向した下流(暖色・矢印短)
θ-β 曲線(複数 M_1)
横軸=流れ偏向角 θ/縦軸=衝撃波角 β/曲線下半分=弱解、上半分=強解/黄丸=現在の (θ, β)
理論・主要公式
斜め衝撃波は、上流マッハ数 M_1 の流れが衝撃波角 β で立つ薄い不連続面を横切るときに、角度 θ だけ偏向される現象です。3 量の関係を θ-β-M 関係式と呼びます。
θ-β-M 関係(M_1, β から θ を陽に求める):
$$\tan\theta = \frac{2\cot\beta\,(M_1^2\sin^2\beta - 1)}{M_1^2(\gamma + \cos 2\beta) + 2}$$
衝撃波の法線成分 M_1n と、その下流側 M_2n(垂直衝撃波関係式の流用):
$$M_{1n} = M_1\sin\beta, \qquad M_{2n}^2 = \frac{1 + \tfrac{\gamma-1}{2}M_{1n}^2}{\gamma M_{1n}^2 - \tfrac{\gamma-1}{2}}$$
下流マッハ数と圧力比(M_1n に対して垂直衝撃波と同じ形):
$$M_2 = \frac{M_{2n}}{\sin(\beta - \theta)}, \qquad \frac{P_2}{P_1} = 1 + \frac{2\gamma}{\gamma+1}\bigl(M_{1n}^2 - 1\bigr)$$
M_1n \gt 1 でなければ衝撃波として物理的に意味がありません。同じ θ に対して β は 2 つあり、小さい方が弱解、大きい方が強解です。
斜め衝撃波 シミュレーターとは
🙋
超音速の機体の主翼の写真でよく見る、ピンと斜めに伸びた線、あれが斜め衝撃波ですか?
🎓
そう、それが斜め衝撃波だよ。垂直衝撃波が「流れに対して垂直に立つ不連続面」だったのに対して、斜め衝撃波は楔やコーンに当たって角度 β で斜めに立つ。流れは衝撃波を通ったあと角度 θ だけ偏向して楔の表面に沿って流れる。シミュレーターのデフォルト(M_1=3, β=32°, γ=1.4)で、流れ偏向角 θ は 14.77°、下流マッハ M_2 は 2.27 と表示されるはずだよ。
🙋
下流もまだ超音速なんですね?垂直衝撃波だと必ず亜音速になると習いました。
🎓
するどい。斜め衝撃波の場合、衝撃波の法線成分 M_1n = M_1·sin(β) だけに垂直衝撃波の関係式が効く。接線成分は素通りするから、全体の M_2 は超音速のままになりやすい。今のデフォルトで M_1n は 1.59 だから、これが垂直衝撃波として効いて、その後ろの法線成分 M_2n は 0.67 まで落ちるけど、接線成分が大きいので M_2 はまだ 2.27 と超音速だ。
🙋
右のグラフ、M_1 ごとに曲線が膨らんでいて、上下に分岐していますね。これは何ですか?
🎓
それが弱解と強解の分岐だよ。同じ θ と M_1 に対して、β は 2 つ存在する。下半分が弱解(β が小さい)、上半分が強解(β が大きい)。実在の楔まわりの流れではほぼ常に弱解が選ばれる。曲線の頂点が θ_max で、楔角がそれを超えると衝撃波が楔から離れる「離脱衝撃波」になる。シミュレーターで β を大きくしていくと、同じ θ なのに点が上の枝(強解)に乗るのが見えるはずだ。
🙋
β を 35° より小さくしたら下流マッハがすごく増えて、もっと大きくしたら下流が遅くなりますね。
🎓
そう、それが斜め衝撃波の本質だ。β が小さい(衝撃波が水平に近い)ほど、法線成分 M_1n が小さくなって衝撃が弱くなり、下流のマッハ数も圧力も上昇が小さい。逆に β が 90° に近づくと垂直衝撃波そのものに近づくから、M_1n = M_1 で最大強度になって M_2 はぐっと落ち込む。超音速インテークの設計では、なるべく弱い斜め衝撃波で段階的に減速して総圧損失を抑える、というのが鉄則だ。
よくある質問
β が arcsin(1/M_1) より小さくなると M_1n = M_1·sin(β) が 1 を下回り、衝撃波の法線成分が亜音速になるため、衝撃波として物理的に成り立ちません。この極限である β = arcsin(1/M_1) は「マッハ角」と呼ばれ、その方向に伝播する弱い圧力擾乱(マッハ波)に相当します。シミュレーターでは M_1n ≤ 1 になる β に対しては全 stat 値が `—` 表示になり、衝撃波の代わりに灰色の点線が描画されます。
θ-β-M 関係式から d(tan θ)/dβ = 0 となる β を解析的に求めると、各 M_1 に対して θ_max が決まります。γ = 1.4 では M_1 = 2 で θ_max ≈ 22.97°、M_1 = 3 で約 34.07°、M_1 = 5 で約 41.12°、M_1 → ∞ で約 45.58° に漸近します。シミュレーター右側の θ-β 曲線で、各 M_1 の曲線の最右端の頂点がこの θ_max に相当し、ここを境に弱解(下枝)と強解(上枝)が分岐します。
θ_max を超える楔角では、楔の先端に張り付く(attached)斜め衝撃波としての解が存在せず、衝撃波は楔の前方に湾曲した形で離脱します。これを離脱衝撃波(detached shock)またはボウショックと呼びます。楔の中心線上では局所的に垂直衝撃波として振る舞い、その背後ではよどみ点圧力・温度が極大になります。再突入カプセル・極超音速の鈍頭物体では基本的にこの離脱衝撃波が支配的です。
外部流れ(楔・コーンなど周囲が開いた流れ)では、上流境界条件のみで決まるため、ほぼ必ず弱解が選ばれます。強解はディフューザのように下流の背圧条件が高く、流れを亜音速まで落とす必要がある場合に現れます。境界層の影響や 3 次元効果でわずかに強解が混じることもありますが、設計上は弱解を仮定するのが標準です。シミュレーターでは β を直接スライダーで指定するため、両方の解を行き来して挙動を比較できます。
実世界での応用
超音速インテーク(多段ショック):SR-71 ブラックバードやコンコルド、F-15 などの超音速機のエンジンインテークでは、複数の斜め衝撃波で段階的に減速し、最後だけ弱い垂直衝撃波で亜音速に落とす多段ショック・インテークが使われます。1 枚の垂直衝撃波で全減速するピトー型に比べて総圧損失が大幅に小さく、エンジン効率を大きく改善します。θ-β-M 関係はランプ角(楔角)と各段の β を決める設計の核心ツールです。
超音速ノズル外部の衝撃ダイヤモンド:ロケットや戦闘機のアフターバーナー出口で見られる「ショックダイヤモンド(マッハディスク)」は、出口圧力と外気圧のミスマッチで生じる斜め衝撃波と膨張波の繰り返し模様です。出口圧が外気圧より高い「過膨張」ではプラントル・マイヤー膨張で広がり、外気圧より低い「不足膨張」では斜め衝撃波で収束し、これが交互に発生します。各斜め衝撃波の角度・強度は θ-β-M 関係で予測できます。
極超音速機の機体周り流れ:X-15、SR-71、最近の HTV-2 や中国の DF-17 などの極超音速ビークルでは、機首・主翼前縁・コントロールフィン周りに複雑な斜め衝撃波系が形成されます。各衝撃波の交差や相互作用(ショック-ショック干渉)は局所的に極めて高い熱・圧力を生じ、構造設計の難所になります。CFD で詳細解析する前段階のクイック評価として θ-β-M 関係が今も活用されます。
超音速翼型・ダイヤモンド翼:超音速旅客機やミサイルの翼型として用いられるダイヤモンド翼は、前縁・後縁ともに鋭い楔状の断面を持ち、上下面で 2 枚ずつの斜め衝撃波と膨張波が形成されます。各面の流れ偏向角と衝撃波角を θ-β-M 関係で求め、表面圧力分布から揚力・抗力を解析的に計算する「ショック-膨張理論」は、超音速空力の基礎です。
よくある誤解と注意点
最も多い誤解は、「斜め衝撃波の前後で M_2 は必ず亜音速になる」と思い込むことです。垂直衝撃波の場合は確かに常に M_2 \lt 1 ですが、斜め衝撃波では衝撃波の接線成分が保存されるため、下流もしばしば超音速のままです。デフォルト条件(M_1=3, β=32°)で M_2 = 2.27 と超音速を保つのが典型例です。M_2 が亜音速になるのは、強解側か、β が 90° に近い垂直衝撃波に近づいた場合に限られます。
次に多いのが、θ_max を超える楔角でも斜め衝撃波の解があると思い込むことです。θ-β-M 関係式は各 M_1 で許容される最大の流れ偏向角 θ_max を持っており、これを超えると数学的に解が存在しません。実流れでは衝撃波が楔から離脱して湾曲したボウショックに遷移します。シミュレーター右の θ-β 曲線で各 M_1 の曲線の最右端が θ_max で、ここから先は曲線が無いことが視覚的に分かります。
もう一つ注意したいのが、β を直接指定するモードと θ を指定するモードの違いです。実用設計では「楔角 θ が決まっていて、それに対応する β を逆算する」のが一般的ですが、本シミュレーターは数値逆解法を避けるために β を直接入力する設計になっています。同じ θ に対して弱解と強解の両方の β を確認できる利点があり、教育用途には適しています。実際の楔まわり流れでは弱解(下枝)が選ばれることを念頭に、結果を読み取ってください。