翼型周りの数値風洞実験をANSYS Fluentで回す際、Cl／Cdを精度よく取り出すには前処理段階で三つの壁を突破する必要がある：メッシュ品質、境界層解像度、そして参照値の一致である。以下では各ステップで詰まりやすいポイントを図解しながら解説する。

1. メッシュ品質を数値化する

まずTUIで品質チェックを実行し、skewnessとorthogonal qualityを数値で確認する。

/mesh/quality/print-quality-summary

出力例：

Minimum Orthogonal Quality = 0.15
Maximum Skewness           = 0.82

翼型表面に沿って伸ばす境界層は、y+が1～5になるよう層厚を自動計算させる。Inflation設定でFirst Aspect Ratioを5に固定し、Growth Rateを1.2に抑えることで急激なセル拡大を回避できる。

2. 境界条件の落とし穴

流入境界では速度値だけでなく、Turbulent Viscosity Ratioを明示的に与えると収束が安定する。推奨値は10である。流出境界は圧力固定で問題ないが、Backflow Turbulent Intensityを5 %にしておくとリサーキュレーションが抑制される。

3. リアルタイムで係数を取得するスキーム

毎イテレーションでCl／Cdを出力するには、Schemeスクリプトを利用する手軽さがある。以下をcalc-lift-drag.scmとして保存し、File > Read > Scheme...で読み込む。

(define (report-cl-cd)
  (let* ((rho 1.225) (vel 30) (area 0.12) (chord 0.3)
         (forces (force-report 'wall-1 '(0 1 0))) ; wall-1 は翼型壁面
         (lift (list-ref forces 1))
         (drag (list-ref forces 0))
         (cl (/ lift (* 0.5 rho vel vel area)))
         (cd (/ drag (* 0.5 rho vel vel area))))
    (ti-menu-load-string
     (format #f "display set-window 1\nplot residuals\n~%"))
    (format #t "Iter ~a  Cl=~5f  Cd=~5f~%"
            (%iterate) cl cd)))

(cx-add-item "Solve" "Execute Commands"
  '(lambda () (report-cl-cd)) 10)

これにより10ステップごとにコンソールに値が表示される。

4. 可視化テンプレートの自動化

複数の迎角を回す際、手作業で図を保存すると抜け漏れが発生しやすい。以下のPythonスニペットはPyFluentを使い、pressure、velocity-magnitude、turbulent-kinetic-energyの三種を連続出力する。

import os, ansys.fluent.core as pyfluent

session = pyfluent.launch_fluent(precision="double", processor_count=8)
session.read_case("wing_aoa_6deg.cas.h5")
session.solution.initialization.hybrid_initialize()
session.solution.run_calculation.iterate(400)

views = ["pressure", "velocity", "tke"]
for v in views:
    session.results.graphics.contour[v].display()
    session.results.graphics.picture.save(
        os.path.join("output", f"{v}_aoa6.png"),
        width=1920, height=1080, format="png")
session.exit()

出力ディレクトリに迎角ごとの画像が整列されるため、比較が容易になる。

5. トラブルシューティングの指針

• 残差が振動する：密度ベースソルバーで圧縮性効果が顕著な場合、Courant Numberを1→0.5へ段階的に下げる。
• 翼端でk値が異常に高い：メッシュが急激に粗くなる領域を特定し、Adapt > Gradientで局所的に細分化する。
• 実験値と5 %以上ズレる：Reference Valuesで入力した面積・長さが実機スケールと一致しているか再確認。特に3D翼では翼端効果を含めた実効面積を用いる。

最後に、Solve > Monitors > Create > Force MonitorでWindow番号を固定しておくと、複数の迎角ケースを連続計算してもグラフが上書きされず、履歴比較が可能になる。