翼型・翼の空力解析 — 実践ガイドとベストプラクティス

典型的なワークフローを示そう。

1. 翼型座標の取得: UIUC Airfoil Coordinates DatabaseやAirfoilToolsから翼型座標をダウンロード

2. 後縁処理: 尖った後縁を0.1%翼弦厚に丸める

3. 計算領域の設定: C型またはO型の外部境界を翼弦の30--50倍に設定

4. メッシュ生成: Pointwiseやsnappy HexMeshで境界層メッシュ込みの格子を作成

5. 境界条件設定: 遠方はfarfield（特性線ベース）、翼面はno-slip壁

6. 乱流モデル選択: 巡航条件ならSA、失速解析ならSST k-omega

7. 求解と収束確認: $C_L$/$C_D$のヒストリを監視

8. メッシュ収束性確認: 粗・中・密の3水準で$C_L$/$C_D$の収束を検証

Richardson外挿法を使うんだ。3水準のメッシュで解を求めて、格子収束指数GCI（Grid Convergence Index）を算出する。

ここで $F_s=1.25$（安全係数）、$\epsilon$ は粗密メッシュ間の解の差、$r$ はメッシュ細分化率、$p$ は収束次数だ。

遷音速翼型（RAE 2822 Case 9）を例に説明しよう。

遠方境界の流れ方向ベクトルで指定する方法が推奨だ。メッシュを回転させると各迎角でメッシュ再生成が必要になるが、ベクトル指定ならパラメトリックスイープが容易になる。Fluentでは flow-direction-vector で指定できるよ。

多要素翼型の解析は難易度が跳ね上がる。各要素間のギャップ（隙間）とオーバーラップが空力性能を決めるんだ。

• スラット/主翼/フラップ間のギャップ領域に高密度メッシュが必要
• 各要素面にそれぞれプリズム層を配置
• 剥離・再付着を伴うため、SST k-omega以上の乱流モデルを推奨
• 非定常解析（URANS/DDES）が必要になることも多い
• NASA Trap Wingが3要素翼型のベンチマークとして広く使われている

実際そうだ。STAR-CCM+のオーバーセットメッシュ機能を使うと、フラップ角度のパラメトリック変更が容易になるよ。各要素を個別メッシュで作成して重ね合わせるアプローチだ。

特にRANSでは失速迎角を2--3度過大に予測する傾向がある。失速近傍の精度が必要ならDDESかIDDESへの移行を検討すべきだ。