乱流モデル選択ガイド・y+計算機 戻る
流体解析

乱流モデル選択ガイド・y+計算機

k-ε・k-ω SST・Spalart-Allmaras・LESを比較し、最適なモデルと壁面y+を選択。入口乱流条件(k,ε,ω)もリアルタイム算出。

y+計算機
レイノルズ数 Re
基準速度 U [m/s]
m/s
密度 ρ [kg/m³]
kg/m³
動粘度 ν [×10⁻⁶ m²/s]
×10⁻⁶
目標 y+ 値
入口乱流条件
乱流強度 Tu [%]
%
長さスケール L [m]
m
k = — m²/s²   ε = — m²/s³   ω = — s⁻¹
計算結果
摩擦速度 u* [m/s]
セル高さ Δy [μm]
境界層厚さ δ [mm]
推奨モデル
壁面法則 u⁺ – y⁺ プロファイル
乱流モデル選択ガイド
モデル性能比較表
モデル自由せん断逆圧力勾配剥離熱移動回転計算コスト
k-ε 標準⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐
k-ε RNG⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐
k-ε Real.⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐
k-ω 標準⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐
k-ω SST⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐低〜中
SA⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐最低
LES⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐⭐非常に高
理論・主要公式

粘性底層($y^+<5$):

$$u^+ = y^+$$

対数則領域($y^+>30$):

$$u^+ = \frac{1}{\kappa}\ln(y^+) + B \quad (\kappa=0.41,\; B=5.2)$$

摩擦速度:$u^* = \sqrt{\tau_w/\rho}$, $C_f \approx 0.027 Re_L^{-1/7}$

入口乱流:$k = \tfrac{3}{2}(T_u U)^2$, $\varepsilon = C_\mu^{3/4}k^{3/2}/L$, $\omega = k^{1/2}/(C_\mu^{1/4} L)$

乱流モデル選択とy+計算とは

🙋
CFDで「y+」ってよく聞くけど、具体的に何ですか?シミュレーターの「初期セル高さΔy」とどう関係するんですか?
🎓
大まかに言うと、壁面近くのメッシュの細かさを決める無次元数だよ。y+が小さすぎると計算コストが増え、大きすぎると精度が落ちる。このツールの「レイノルズ数」や「動粘度」を変えると、推奨される「初期セル高さΔy」がリアルタイムで計算される。例えば、自動車の外気解析で流速を上げてみると、必要なΔyが小さくなることがわかるよ。
🙋
え、そうなんですか!じゃあ、下のグラフで見える「粘性底層」と「対数則領域」って何が違うんですか?k-εとk-ω SSTで推奨y+が違う理由も知りたいです。
🎓
壁面にべったりくっついて速度が直線的に増える層が粘性底層(y+<5)、その外側で速度が対数的に増えるのが対数則領域(y+>30)だ。k-εモデルは壁面関数を使うので、最初のセル中心を対数則領域(y+≈30〜300)に置く。一方、k-ω SSTは粘性底層も直接解く「低Re数モデル」だから、最初のセルをy+≈1にしないと精度が出ない。ツールで乱流モデルを切り替えると、推奨y+とグラフの関係がよくわかるよ。
🙋
実務では、どうやってこの情報を使うんですか?「乱流強度」とか「長さスケール」も入力する意味は?
🎓
現場では、このツールで計算したΔyを元に、壁面の最初のメッシュ層の厚さを決めるんだ。例えば、航空機の翼解析なら、SSTモデルを選び、y+≈1になるようにΔyを設定する。「乱流強度」と「長さスケール」は、kやεといった乱流変数の初期値を推定するために使う。風洞実験データがあればそれを使い、なければツールのデフォルト値から始めて、シミュレーション結果を見ながら調整するのが一般的だね。

よくある質問

ツールが出力するy+値は目標値です。第1セル高さΔyは、想定するレイノルズ数と摩擦速度から逆算します。例えばy+=1を目指す場合、Δy = y+ × ν / u* で計算し、メッシュ作成時の初期値としてください。
k-εは高レイノルズ数・非回転流れ(パイプ流など)に適し、壁面は標準壁関数(y+≈30~300)が必要です。k-ω SSTは剥離や逆圧力勾配を伴う流れ(翼周りなど)に強く、低y+(≈1)でも高y+(≈30)でも使えるため汎用性が高いです。
算出値は、設定した乱流強度(例:5%)と代表長さから推定した目安です。kは乱流運動エネルギー、εとωは散逸率を表します。実際の実験データや経験則と比較し、オーダーが合っているかを確認してから入力してください。
LESは非定常な大規模渦を直接解くため、渦放出や騒音予測に適しますが、壁面近傍で非常に細かいメッシュ(y+≈1)と小さいタイムステップが必要です。RANS(k-ε, k-ω SSTなど)は平均場の予測に優れ、定常解析や設計初期のパラメータスタディに実用的です。

実世界での応用

自動車の空力解析:車体周りの流れや抗力係数(Cd値)の予測にCFDが広く使われます。k-ω SSTモデルとy+≈1のメッシュが用いられ、ドアミラーやピラーまわりの複雑な剥離流れを精度良く捉えます。

航空機・翼設計:揚力と抗力の評価、失速特性の予測が主な目的です。逆圧力勾配下での境界層発達や剥離を扱えるSSTモデルが適しており、翼面上のメッシュはy+≈1に設定されます。

建物周りの風環境解析:都市計画やビル風の評価に利用されます。複雑な形状が多いため、比較的頑健なk-εモデルと標準壁面関数(y+≈30-300)の組み合わせがよく使われ、計算コストと精度のバランスが取られます。

ターボ機械(ポンプ・ファン)内部流れ:羽根車内の流れや効率予測にCFDが不可欠です。強い曲率と回転の影響を受けるため、SSTモデルが推奨され、壁面近傍のメッシュは非常に細かく(y+≈1)設定されます。

よくある誤解と注意点

まず、「y+さえ合っていれば全てOK」という誤解があります。確かにy+は壁面メッシュの第一歩ですが、それだけでは不十分。例えば、y+≈1に設定しても、壁面に垂直な方向のセル数が少なすぎると(例えば5層以下)、境界層全体の速度勾配を解像できず、剥離の予測を誤ります。目安として、境界層全体をカバーするのに10〜15層以上のメッシュを積むことが多いです。

次に、「推奨y+値は絶対的な神の数値」と思い込むこと。ツールが示す「k-εモデルならy+≈30」はあくまで出発点。実際の流れ場が複雑(強い圧力勾配、剥離・再付着がある)なら、より小さなy+(例えば10〜20)にした方が良い場合もあります。シミュレーション結果の壁面せん断応力や圧力分布を必ずチェックし、必要に応じてメッシュを再構成する柔軟さが重要です。

最後に、流入境界条件の「乱流強度」と「長さスケール」を適当に設定すること。これらは計算の初期値に直結し、特に流れの発達に影響します。例えば、風洞実験を模倣するなら、実験値を使うべきです。データがない場合、内部流路(配管など)では強度を数%、外部流れ(自動車周りなど)では0.5〜1%程度から始めることが多いですが、「とりあえずデフォルト」で済ませず、感度解析で影響を確認する姿勢がプロの証です。