翼型・翼の空力解析 — 理論と支配方程式
概要
先生、航空機の翼まわりの空力解析って、CFDでどうやるんですか？
翼型（エアフォイル）の揚力・抗力特性を予測し、失速挙動や高揚力装置の効果を評価する解析だ。航空機設計の根幹を担う技術だよ。
NACA翼型シリーズから始まり、現代ではスーパークリティカル翼型や自然層流翼型の設計にCFDが不可欠になっている。風洞試験だけでは探索しきれない設計空間をCFDで広げるわけだ。
風洞だけじゃダメなんですか？
風洞試験は1条件あたり数百万円規模のコストがかかる。CFDで設計候補を絞り込んでから風洞に持ち込むのが現代の標準的なワークフローだよ。
支配方程式
翼まわりの流れを記述する方程式を教えてください。
圧縮性Navier-Stokes方程式が基本だ。連続の式、運動量の式、エネルギーの式の3本セットで記述する。
揚力係数と抗力係数はこう定義される。
ここで $L$ は揚力、$D$ は抗力、$\rho_\infty$ は自由流密度、$V_\infty$ は自由流速度、$S$ は翼面積だ。
レイノルズ数はどのくらいの範囲になるんですか？
翼弦長ベースのレイノルズ数は旅客機巡航時で $Re_c \approx 2 \times 10^7$ 程度だ。マッハ数は遷音速域の $M \approx 0.78$--$0.85$ が主戦場になる。
遷音速域では翼上面に局所的な超音速領域が生じて衝撃波が発生する。この衝撃波と境界層の干渉がバフェット現象を引き起こし、失速の引き金になる。
なるほど。衝撃波の位置を正確に捉えることが重要なんですね。
乱流モデルの選択
翼まわりの解析ではどの乱流モデルが使われるんですか？
用途に応じて使い分ける。遷移予測が必要な場合と、完全乱流を仮定する場合で選択肢が変わるんだ。
| モデル | 特徴 | 翼解析への適性 |
|---|---|---|
| Spalart-Allmaras (SA) | 1方程式モデル。航空分野で広く使用 | 巡航条件で良好。失速近傍はやや苦手 |
| SST k-omega | 壁近傍k-omegaと遠方k-epsilonのブレンド | 逆圧力勾配・剥離に強い |
| gamma-Re_theta 遷移モデル | SA/SSTと組み合わせて自然遷移を予測 | 自然層流翼の設計に必須 |
| DDES/IDDES | RANS+LESのハイブリッド | 大規模剥離・バフェット解析 |
Spalart-Allmarasモデルって航空系ではすごくメジャーですよね。
そうだ。SA模型はもともとNASAで翼型解析のために開発されたモデルだからね。Boeing、Airbus両社とも広く使っている。ただし失速近傍の大規模剥離にはSST k-omegaやDDESが必要になる。
翼型の空力特性
具体的にどんな数値が出てくるものなんですか？
RAE 2822って風洞データが公開されてるからCFDの検証によく使われますよね。
その通り。Case 9（$M=0.73$, $\alpha=2.79°$, $Re=6.5 \times 10^6$）は業界標準のベンチマークだ。CFDで翼上面の圧力分布と衝撃波位置を風洞データと比較するんだよ。
実務上の注意点
翼の解析で特に気をつけるべきポイントは何ですか？
一番重要なのは境界層解像だ。壁面第一層の$y^+$値を1以下に設定して、境界層内に十分なプリズム層を確保する必要がある。
- $y^+ \approx 1$: SA/SSTモデルでの推奨値。壁関数なしで粘性底層を解像
- プリズム層数: 最低20層以上。成長率は1.2以下を推奨
- 後縁処理: 鋭い後縁は特異点になるため、有限厚さ（翼弦の0.1%程度）に丸める
- 遠方境界: 翼弦長の20--50倍離れた位置に配置
翼の解析は境界層の扱いが命なんですね。よく分かりました。
そうだ。抗力の予測精度は境界層メッシュの品質に直結する。$\Delta C_D = 0.0001$（1カウント）の精度を求められることもあるから、メッシュには細心の注意が必要だよ。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース（下向きの空力的な力）を発生します。つまり理論上、天井に貼り付けて走れる！ チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値例：円管内層流（d=10mm, L=1m, 水μ=0.001Pa·s, Q=0.1L/min）
Re = ρuD/μ = 998×0.021×0.01/0.001 ≈ 212（層流） 圧力損失 ΔP = 128μLQ/(πd⁴) ≈ 68.2 Pa
乱流モデル別の精度比較（後向きステップ、再付着長さ）：
k-ω SSTは精度とコストのバランスが良く、多くの実務で最初の選択肢になります。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「翼型・翼の空力解析をもっと効率的に解析できないか？」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →