ポテンシャル流れ理論 — トラブルシューティングガイド
パネル法のよくある問題
パネル法を使っていてハマりやすいポイントを教えてください。
代表的なトラブルと対策を解説しよう。
1. 圧力分布に不自然なスパイクが出る
翼型のパネル法解析で、前縁付近に圧力のスパイクが出るんですが。
原因と対策は以下だ。
| 原因 | 対策 |
|---|---|
| パネル分布が不均一 | 前縁・後縁にパネルを集中させるコサイン分布を使用: $x_i = \frac{c}{2}(1 - \cos(\pi i/N))$ |
| パネル数不足 | $N \geq 100$ に増やす（前縁に20パネル以上） |
| 後縁の処理が不適切 | 後縁を開いた翼型（有限厚み後縁）の場合、後縁パネルの特別処理が必要 |
| 座標データの精度 | 翼型座標の小数点以下6桁以上を確保 |
コサイン分布は前縁に自動的にパネルが集中するんですね。
前縁はよどみ点で圧力勾配が急峻だから、均等配分では解像度が全く足りない。コサイン分布にするだけで精度が劇的に改善されるよ。XFOILでは内部的にこの分布を使っているんだ。
2. 揚力が理論値と合わない
薄翼理論では $C_l = 2\pi\alpha$ ですが、パネル法の結果が大きくずれます。
以下をチェックしよう。
- Kutta条件の実装: 後縁で上下面の速度が等しくなる条件を正しく実装しているか。これが揚力を決定する最も重要な条件
- 迎角の定義: 翼弦線の定義と迎角のゼロ点が一致しているか
- 翼型の厚み: 薄翼理論は厚み0の仮定。厚い翼型（t/c > 15%）では $C_l$ が $2\pi\alpha$ から乖離して当然
- パネルの向き: 法線ベクトルが全て外向き（反時計回り）に統一されているか。一部が逆だと解が狂う
3. 影響係数行列が特異になる
連立方程式を解こうとすると「singular matrix」エラーが出ます。
影響係数行列が特異になる原因は
- パネルが重複している: 翼型座標の始点と終点が重複。後縁のパネルが長さゼロになっている可能性
- パネルが極端に小さい: 長さ比が $10^4$ 以上になるとDouble精度でも条件数が悪化
- 閉じた物体で全Neumann条件: 湧き出しのみのパネル法では解が不定（一意性のための条件が必要）。内部点で $\phi = 0$ を追加するか、二重極定式化を使う
XFOILのトラブルシューティング
XFOILでよくあるトラブルも教えてください。
XFOIL特有の問題をまとめよう。
| 問題 | 原因 | 対策 |
|---|---|---|
| 「VISCAL: Convergence failed」 | 粘性解が収束しない | VACC で収束閾値を緩和（0.01→0.05）、迎角を小さいところから徐々に増やす |
| 失速角を超えると計算不能 | 大規模剥離はXFOILの適用外 | 失速後はCFD（Fluentなど）を使用 |
| 低Reynolds数で不正確 | $Re < 10^5$ で遷移モデルの限界 | Ncrit を適切に設定（風洞: 9、飛行: 11〜13） |
| 翼型データ読み込みエラー | 座標フォーマットの問題 | Selig形式（上面→下面の順、後縁始まり）に統一 |
MACH 設定で不安定 | 高亜音速で圧縮性補正が不安定 | $M < 0.5$ 程度に留める。遷音速はCFDを使用 |
Ncrit の設定って結果に大きく影響しますか？
Ncrit は $e^N$ 遷移モデルのN因子の臨界値で、乱流遷移位置を左右する。低乱風洞では $N_{crit} = 9$（デフォルト）、通常風洞では $N_{crit} = 5$〜$7$、飛行環境では $N_{crit} = 11$〜$13$ が目安だ。これを間違えると抗力が数十%変わるので要注意だよ。
CFDとポテンシャル解の不一致
Fluentの結果とXFOILの結果が合わないとき、どちらを信じればいいですか？
不一致の原因を切り分けることが重要だ。
- $C_l$ が合わない場合: まず迎角の定義を確認。FluentのReference Valuesの方向ベクトルとXFOILの迎角定義が一致しているか
- $C_d$ が合わない場合: XFOILは粘性抗力のみ。FluentのRANS解は圧力抗力+粘性抗力の合計。剥離がある場合はXFOILの$C_d$が過小
- 圧力分布の差: FluentでY+が不適切だと壁面近傍の速度場が不正確になり圧力分布がずれる
- 両方とも間違い: メッシュ収束していないFluentと粗いパネルのXFOILの比較には意味がない。それぞれ独立に収束確認してから比較すること
結局、どちらも検証が必要ってことですね。
実験データがゴールドスタンダードだ。Abbott & von Doenhoffの「Theory of Wing Sections」にNACA系翼型の詳細な風洞データが掲載されている。まずそれと合わせるのが第一歩だよ。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース（下向きの空力的な力）を発生します。つまり理論上、天井に貼り付けて走れる！ チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」（どの残差が下がらないか）を特定し、次に「何が詰まっているか」（メッシュ品質？境界条件？乱流モデル？）を調べ、最後に「どう直すか」（メッシュ修正？緩和係数？）を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ポテンシャル流れ理論の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
CAEの未来を、実務者と共に考える
Project NovaSolverは、ポテンシャル流れ理論における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →