🙋
飛行機が亜音速で飛ぶときって、低速のときと圧力分布って違うんですか？マッハ 0.3 までは非圧縮と思っていいって聞いたんですけど。
🎓
マッハ 0.3 程度までは確かに非圧縮で良い近似なんだけど、それ以上、特に M = 0.5〜0.7 だと密度変化が無視できなくなる。そこで使うのがプラントル・グラウェルトの補正則だ。シミュレーターで C_p,inc = −0.50、M = 0.60、γ = 1.40 にすると、補正係数 1/β = 1.250、C_p,compr = −0.625 と表示されるはずだ。約 25% 増しになるんだよ。
🙋
えっ、25% も変わるんですか？じゃあ揚力係数も増える方向なんですね？
🎓
そう、揚力係数も C_L,compr = C_L,inc / β で増える。デフォルト値で C_L,inc = 0.50 → C_L,compr = 0.625 となる。同じ迎角でも圧縮性があるほうが揚力が増えるから、亜音速ジェット旅客機の巡航設計では「圧縮性は味方」として扱う場面と「逆風」として扱う場面が両方ある。たとえば設計揚力が増えれば抗力・モーメントも変わるから、トリム計算で必須なんだ。
🙋
補正係数のグラフを見ると、M = 0.9 あたりでめちゃくちゃ大きくなりますね。これはどこまで信用できるんですか？
🎓
いいところに気づいたね。1/β は M → 1 で発散するけれど、現実にはそんなに増えない。プラントル・グラウェルトは線形化の仮定があるから、M ≤ 0.7 程度までしか実験と合わない。それ以上では局所超音速領域が出てきて衝撃波が立ち、衝撃失速や造波抗力の影響が出てくる。そこからは Karman-Tsien 補正やフル CFD の領域になる。シミュレーターのカーブはあくまで補正則そのもので、物理的に意味があるのは M ≤ 0.7 の左半分なんだ。
🙋
全圧/静圧比 P_0/P もスライダーで表示されてますが、これは何に使うんですか？
🎓
これは等エントロピー完全気体の流れの状態比で、ピトー管の指示値からマッハ数を逆算するときの基礎式だよ。M = 0.6、γ = 1.4 なら P_0/P ≈ 1.276。風洞試験では翼面の静圧と上流の全圧を測って局所マッハ数を求めるから、補正則と並んで実験データ処理の必需品なんだ。圧縮性補正の C_p とは別物の量だから混同しないようにね。
プラントル・グラウェルトは線形化された亜音速ポテンシャル流の最も単純な補正則で、C_p,compr = C_p,inc / √(1−M²) という形をとります。Karman-Tsien 則は等エントロピー関係を一部取り入れて遷音速域への適用範囲を広げた補正で、Laitone 則はさらに局所マッハ数の効果まで考慮します。実務では M ≤ 0.6 までは P-G、0.6〜0.85 では Karman-Tsien を選び、それ以上の遷音速・超音速近傍はフルポテンシャル解や Euler 解析が必要です。シミュレーターは最も基本的な P-G 則を扱っており、結果の物理的限界を理解する出発点になります。
プラントル・グラウェルト補正の原則として、亜音速の線形化空気力学では、局所 C_p と積分量である C_L、ピッチングモーメント係数 C_m すべてが 1/β だけ増加します。一方、誘導抗力係数 C_Di は揚力の自乗に比例するため、近似的に 1/β² 倍になります。粘性抗力（摩擦抗力・形状抗力）はこの線形化に乗らないので、別途レイノルズ数依存の経験式や境界層解析が必要です。造波抗力（M ≥ 0.7 で顕著）は P-G では一切扱えず、衝撃波解析が必須になります。
プラントル・グラウェルトは薄翼・小擾乱・小迎角を仮定しています。実用的には、最大翼厚比 12% 以下、迎角 ±5° 以下、平均キャンバー 4% 以下程度なら、M ≤ 0.6 の範囲で実験と 5% 以内の精度で一致します。それを超えると非線形効果が無視できなくなります。厚翼や高迎角では局所マッハ数が翼面で 1 を超えてしまい（超臨界状態）、衝撃波が発生して P-G の前提が破綻します。設計初期の翼型比較やパラメータスタディには有効ですが、最終的な空力特性は風洞試験や CFD で検証する必要があります。
はい、これは「非圧縮極限で補正は無効」という物理的に正しい挙動です。M → 0 で 1/β → 1 となり C_p,compr = C_p,inc・C_L,compr = C_L,inc となります。実用的にはこれは M ≤ 0.3 の低速空気力学で補正不要という意味で、ハングライダー・ドローン・自動車などの空力解析は非圧縮ナビエ・ストークス（または非圧縮ポテンシャル流）で十分です。逆に旅客機の巡航 M = 0.78〜0.85 では補正係数が 1.6〜1.9 に達し、無視すれば設計揚力を 40〜90% も過小評価してしまいます。
プラントル・グラウェルト補正則は、亜音速の小擾乱ポテンシャル流方程式 $(1-M_\infty^2)\phi_{xx} + \phi_{yy} = 0$ の解析から導かれます。座標変換 $x = X$、$y = Y\sqrt{1-M_\infty^2}$ により、圧縮性流れの方程式が非圧縮ラプラス方程式 $\phi_{XX} + \phi_{YY} = 0$ に帰着できることを利用します。境界条件と圧力係数の関係を線形化することで、最終的に次の補正式が得られます：
$$C_{p,\text{compr}}(x, y) = \frac{C_{p,\text{inc}}(x, y/\beta)}{\sqrt{1 - M_\infty^2}}$$
ここで $\beta = \sqrt{1 - M_\infty^2}$ はプラントル・グラウェルト因子です。揚力係数・モーメント係数も同じ補正係数で線形に増加します：
$$C_L^{\text{compr}} = \frac{C_L^{\text{inc}}}{\sqrt{1 - M_\infty^2}}, \qquad C_m^{\text{compr}} = \frac{C_m^{\text{inc}}}{\sqrt{1 - M_\infty^2}}$$
等エントロピー完全気体の関係から、全圧 $P_0$ と静圧 $P$ の比は：
$$\frac{P_0}{P} = \left(1 + \frac{\gamma - 1}{2}M^2\right)^{\gamma/(\gamma - 1)}$$
これは補正則とは別の独立した式で、ピトー管測定からマッハ数を求める場合や、全温・静温の比 $T_0/T = 1 + (\gamma-1)/2 \cdot M^2$ とともに使われます。シミュレーターでは両者を同時に表示することで、亜音速圧縮性流れの全体像が直感的に把握できるようにしています。
亜音速ジェット機の予備設計：巡航マッハ数 0.7〜0.85 の旅客機・ビジネスジェットの翼設計初期段階では、非圧縮の翼型データベース（NACA・SC シリーズ等）にプラントル・グラウェルト補正を適用して、圧縮性込みの揚力・抗力係数を見積もります。最終的な空力特性は風洞や CFD で検証しますが、トリム・配置・主翼面積の初期推算には P-G 補正のシンプルさが今でも重宝されています。設計巡航 M = 0.78 では補正係数が 1.6 程度になるため、無視すると設計揚力が大きくずれます。
風洞試験データの圧縮性補正：低速風洞（M ≤ 0.3）で得た翼型データを、目的の巡航マッハ数（M = 0.5〜0.7）の値に換算する場合に P-G 則を使います。逆に、亜音速風洞のデータを地上対地速度（M ≈ 0）の感覚で議論する場合には除算で「非圧縮等価」値に戻します。風洞補正の標準作業として、JAXA や AIAA のテストレポートにも明示的に P-G 補正が登場します。
プロペラ翼端速度の評価：大型プロペラの翼端では局所マッハ数が 0.7〜0.9 に達することがあり、翼根と翼端で補正係数が大きく異なります。各半径位置で P-G 補正を行い、局所揚力分布を補正することで、回転翼の効率・推力・騒音を予測する標準手法です。ヘリコプターのメインローター翼端や、無人機のプロペラ設計でも使われます。
UAV・小型機の高高度運用：高高度で巡航する偵察 UAV や成層圏プラットフォームでは、低密度のため対気速度が高く、結果としてマッハ数も 0.5〜0.7 に達することがあります。気圧 P の低下による静圧計算と、P-G 補正を組み合わせて翼面の局所圧力分布を予測し、構造設計に反映させます。地上低速試験のデータを直接使うと過小評価になるため、補正が不可欠です。
最も多い誤解は、「プラントル・グラウェルトは遷音速にも使える」と考えてしまうことです。1/β は M → 1 で発散するため、補正式は数学的には任意の M で値を返しますが、物理的には M ≥ 0.7 程度で線形化の前提が破綻し、局所超音速領域・衝撃波の発生によって補正値は実測と大きく外れます。シミュレーターで M = 0.9 にすると 1/β ≈ 2.29 となり、表面上は「2 倍以上に補正される」と読めますが、実際の翼面では衝撃波抵抗が支配的になっていて、薄翼線形理論で計算した C_p,compr は意味を失います。
次に多いのが、補正係数 1/β と全圧/静圧比 P_0/P の混同です。両者は M に依存する量ですが、物理的な意味が全く異なります。1/β は「翼面の局所 C_p をどれだけ強化するか」という線形化補正の係数で、P_0/P は「気流自体の全圧と静圧の比」という等エントロピー流の状態量です。シミュレーターのデフォルト値（M = 0.6、γ = 1.4）で、1/β = 1.250、P_0/P ≈ 1.276 と偶然似た値になりますが、M が変わると挙動が違います（1/β は M → 1 で発散、P_0/P は M = 0.6 → 1.0 で 1.28 → 1.89 と緩やかに増加）。
もう一つの誤解は、「補正は風洞データだけに使う」というものです。実際には、非圧縮 CFD（簡易パネル法・渦格子法）の出力を圧縮性込みの値に変換する場合にも頻繁に使われます。MIT 開発の XFOIL や Drela 系のパネル法は基本的に非圧縮ですが、Karman-Tsien 補正を組み合わせて圧縮性込みの C_p 分布を出力します。設計初期の高速計算ループで、フル CFD の代わりにこうした補正則ベースのツールを使うのは依然として一般的です。シミュレーターで補正係数のマッハ数依存性を視覚的に理解しておくと、これらのツールを使いこなすときの直感が養われます。