🙋
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%程度から始めることが多いですが、「とりあえずデフォルト」で済ませず、感度解析で影響を確認する姿勢がプロの証です。