ポテンシャル流れ理論 — 実践ガイドとベストプラクティス

NACA 4412 翼型の揚力特性を求める例で説明しよう。

1. XFOILを起動し、NACA 4412 と入力して翼型を読み込む

2. OPER コマンドで解析モードに入る

3. VISC 500000 でReynolds数を設定（Re = 500,000）

4. ITER 200 で反復回数の上限を設定

5. ALFA 5 で迎角5度の解析を実行

6. CPLT で圧力分布をプロット

XFOILはポテンシャル流れ（パネル法）と境界層の積分法を連成させている。外部流れはパネル法で解き、壁面近傍の粘性効果は境界層方程式で別途計算し、両者を反復的に整合させるんだ。この粘性-非粘性連成（viscous-inviscid interaction）により、摩擦抗力や剥離の予測も可能になる。

主要な確認項目はこうだ。

段階的に検証することをおすすめする。

Step 1: 円柱周りの流れ

• 解析解が既知: $C_p = 1 - 4\sin^2\theta$
• パネル数を $N = 20, 40, 80, 160$ と増やして収束確認
• $N=80$ で $C_p$ の誤差が 0.1% 以下になるはず

Step 2: 楕円柱周りの流れ

• 厚み比を変えてテスト。よどみ点圧力 $C_p = 1.0$ の精度確認

Step 3: Joukowski翼型

• 解析解が既知の翼型。揚力がKutta-Joukowski定理と一致するか確認
• $C_l = 2\pi(\alpha + \alpha_0)$ との比較（薄翼理論の範囲で）

Step 4: NACA翼型

• XFOILの結果と比較。座標データはUIUC Airfoil Data Siteから入手可能

最初から複雑な形状で試すと、バグがあっても原因特定が困難になる。円柱から始めて徐々に複雑にしていくのが鉄則だ。

XFLR5やOpenVSPを使えばGUIベースで3D翼の解析ができる。XFLR5の典型的なワークフローはこうだ。

1. Direct Foil Design で翼型を定義（XFOILデータを読み込める）

2. Wing and Plane Design で翼平面形を設定（スパン、テーパー比、後退角、ねじり）

3. 解析手法を選択: LLT（揚力線理論）、VLM（渦格子法）、3D Panel

4. 迎角を変えて揚力・誘導抗力を計算

初期設計ではVLMで広いパラメータスペースを探索し、有望な候補をXFOILの2D解析やCFDで詳細検討するのが効率的だよ。

大いにある。具体的な活用場面を挙げよう。

• 概念設計段階: 数百の翼型・翼形状を数分で評価（N-Sだと1ケース数時間）
• 最適化ループ: XFOILをPythonから呼び出して自動最適化。数千ケースの評価が現実的
• CFDの境界条件: ポテンシャル解をCFDの初期条件や遠方境界条件に使用
• 結果の妥当性検証: CFD結果が薄翼理論の $C_l = 2\pi\alpha$ と大きく乖離していたら何かがおかしい

XFOILで1翼型・1迎角の計算は1秒未満だ。Fluentで同じ問題を解くと数分〜数十分かかる。この速度差がパラメトリックスタディや最適化で決定的な違いを生むんだ。