壁関数 — トラブルシューティングガイド
より充実した内容を wall-function.html でご覧いただけます。
よくあるトラブルと対策
壁関数がらみでよくあるトラブルって何ですか？
現場で最も多い問題をまとめよう。
1. y+ がバッファ層に落ちる
症状: 壁面 $y^+$ が5〜30の範囲になっている
原因: メッシュの第1層が薄すぎる（壁関数用なのにLow-Re用の厚さにしてしまった）
対策:
- 第1層厚さを増やして $y^+ > 30$ にする
- あるいは逆に $y^+ < 5$ にしてEnhanced Wall Treatment に切り替える
- Scalable Wall Function（Fluent）やAll y+ Treatment（STAR-CCM+）を使えば、$y^+$ がバッファ層でもある程度補正される
2. 壁面摩擦係数が実験値と合わない
$C_f$ が実験と20%以上ずれることがあるんですが。
考えられる原因:
- $y^+$ が対数層の範囲外
- 強い逆圧力勾配で対数則が崩れている
- 粗面モデルの粗さパラメータ $k_s$ が不適切
- メッシュが壁に沿って粗すぎ（流れ方向の解像度不足）
対策:
- $y^+$ のコンター図を確認し、全壁面で30〜100に収まっているか検証
- 圧力勾配が強い領域ではNon-Equilibrium WFを検討
- 実験の粗さデータを $k_s$ に正しく変換する（等価砂粒粗さ）
3. 壁面熱伝達係数が過大/過小
Nu数が実験の半分くらいしか出ないことがあるんですけど。
温度の壁関数は速度の壁関数以上にメッシュ感度が高い。
原因: 高Pr数流体（オイル等）で $y^+$ が大きすぎると、温度境界層（速度境界層より薄い）が第1セル内に完全に埋もれてしまう
対策:
- Pr > 1 の流体では $y^+$ を小さめに設定する（目安: $y^+ < 50 / \text{Pr}^{0.5}$）
- Enhanced Wall Treatmentに切り替える
- STAR-CCM+ではTwo-Layer All y+ Wall Treatment が有効
4. 局所的な逆流で発散する
症状: 壁隣接セルに逆流があり、壁関数の計算が不安定化
対策:
- 剥離領域の $y^+$ を確認。剥離点近傍では対数則が成立しないため壁関数が適さない
- Enhanced WTやAll y+ Treatmentに変更
- Under-relaxation factor を下げる（Fluentでは Momentum を 0.5 → 0.3）
デバッグの手順
壁関数の問題を効率よく診断する手順を教えてください。
1. まず壁面 $y^+$ のコンター図を出す。全壁面が適切な範囲にあるか確認
2. 壁面 $C_f$ 分布を理論値や実験と比較。ずれが大きい箇所を特定
3. 疑わしい壁面近傍の速度プロファイルを $u^+$ vs $y^+$ でプロットし、対数則と比較
4. メッシュ感度チェック: 第1層を1/2と2倍にして結果の変化を確認
$u^+$ vs $y^+$ のプロットって、ソルバーの標準機能で出せますか？
Fluentでは壁面に垂直なラインプローブで $U$ と $y$ を取得し、$u_\tau$（壁面摩擦速度）で無次元化する。ParaViewでも同様に計算可能だ。対数則 $u^+ = 2.44 \ln(y^+) + 5.5$ と比較して、解が対数層を正しく再現しているかが一目瞭然になる。
y+=11.25という「魔法の数字」が生む悲劇
壁関数の適用限界として有名な「y+=11.25」は、粘性底層と対数層の交差点に対応する値です。この値付近のメッシュは「粘性底層を解像するには粗すぎ、対数層に合わせるには細かすぎ」という最悪のゾーンに入ります。大型モデルを一括でメッシュ生成すると、一部の局所的な高速域でy+がちょうどこの値に入ってしまうことがあり、その部分だけ壁面せん断応力がおかしな値になります。「圧力損失が特定の流量条件だけ大きくズレる」という症状はy+=11付近の壁面セルが原因である場合があり、y+コンターの確認が重要です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——壁関数の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告