翼型・翼の空力解析 — トラブルシューティングガイド

代表的な問題を整理しよう。翼型のCFDは精度要求が厳しいぶん、落とし穴も多いんだ。

症状: 遷音速翼型で翼上面の衝撃波位置がCFDと実験で5--10%翼弦ずれる

考えられる原因:

• 風洞壁干渉の未考慮（特に閉塞比が大きい試験）
• 乱流モデルの予測精度不足
• メッシュ解像度の不足（衝撃波前後）

対策:

• 風洞壁を含めた解析、または壁干渉補正の適用
• SA/SSTの両方で比較し、感度を確認
• 衝撃波近傍にアダプティブメッシュ細分化を適用

症状: $C_D$が実験値と10%以上乖離する

考えられる原因:

• $y^+$が大きすぎて壁面せん断応力の精度が低い
• 後縁近傍のメッシュ不足
• 遠方境界が近すぎて人工的な閉塞効果が発生
• 風洞の支持装置干渉を未考慮

対策:

• $y^+ < 1$を厳守し、壁関数を使わない
• 後縁上下面にそれぞれ最低50点を配置
• 遠方境界を翼弦の50倍以上に拡大
• near-field法とfar-field法の両方で抗力を評価し比較

near-field法は翼面上の圧力と摩擦力を直接積分する方法。far-field法は翼の十分下流での運動量変化から抗力を算出する方法だ。両者の差がメッシュ品質の指標になる。差が大きいほどメッシュが不十分ということだよ。

症状: RANSで失速迎角が実験値より2--4度大きい

対策:

• SST k-omegaモデルに切り替え（SAより剥離予測が良好）
• 大規模剥離にはDDESまたはIDDESを使用
• 非定常解析に切り替え（定常RANSでは剥離後の振動挙動を捉えられない）
• 時間刻みを $\Delta t \cdot V_\infty / c \approx 0.01$ 程度に設定

症状: 残差が発散、または$C_L$が非物理的に振動

対策:

• CFL数を1--5に下げて再スタート
• 1次精度で初期解を作り、2次精度に切り替え
• メッシュ品質チェック: 非直交性 > 70度、スキューネス > 0.95のセルを修正
• 遷音速の場合、限定関数（Venkatakrishnan limiter）のパラメータを調整

• 密度ベースソルバーではCFL初期値を5以下に抑える
• Roe-FDSスキームでEntropy Correctionを有効にする（膨張衝撃波の防止）
• Reference Valuesを正しく設定しないと$C_L$/$C_D$が間違った値で出力される
• 2D解析でもAxial/Radialの座標方向設定に注意（デフォルトがZ軸上向きの場合あり）

翼弦長、翼面積、参照密度・速度が正しく設定されていないと、無次元係数の値が全くおかしくなる。CFDの結果が物理的に妥当か、まず概算値と比較する習慣をつけることが大事だよ。