流れ条件
$\tau_w = C_f \cdot \rho U^2 / 2$
$u_\tau = \sqrt{\tau_w / \rho}$
$\Delta y = y^+ \cdot \nu / u_\tau$
CFD壁面境界層解像度の指標 y⁺ と必要な第1セル高さ Δy を Reynolds 数・流体・代表長さから即計算。壁関数法と低レイノルズ数モデルの選定基準も一覧で確認できます。
自動車・航空機の空力設計:車体や機体表面の摩擦抵抗を正確に予測するために、境界層メッシュの設計は極めて重要です。特にモータスポーツでは、数%の抵抗低減が勝敗を分けるため、y⁺に基づいた高品質なメッシュ生成が必須です。
ターボ機械の設計:タービンブレードやポンプインペラーの表面は複雑な流れと高いせん断応力にさらされます。適切なy⁺設定により、剥離や失速現象を精度良く捉え、効率と耐久性を両立させます。
建築・環境風工学:高層ビル周辺の風環境や構造物に働く風荷重を評価する際、地面や壁面近くの乱流を正しくモデル化する必要があります。大規模計算では計算コストと精度のバランスを取りつつy⁺を決定します。
電子機器の熱設計:基板やヒートシンクの強制空冷シミュレーションでは、固体壁面近くの熱伝達率の計算精度が全体の温度予測に直結します。y⁺を適切に設定することで、現実に近い冷却性能を評価できます。
まず、「y⁺の目標値は絶対守らなければならない神の数値」と思い込むことです。確かにガイドラインは重要ですが、例えばk-ω SSTモデルでy⁺≒1を目指しても、複雑な幾何形状では全ての壁面でこれを達成するのは不可能に近い。実務では、流れの物理的に重要な領域(剥離点、再付着点、高せん断応力部)で優先的に目標y⁺を満たすようにメッシュを設計します。それ以外の場所では多少外れても許容されることが多いです。
次に、計算した第1セル厚さΔyをそのままメッシャーに入力すればOK、と考える落とし穴。このツールで出てくるΔyは「セル中心」までの距離です。多くのメッシュ生成ソフトでは、「第1層の厚さ」はセルそのものの厚さを指すため、計算値の約2倍を入力する必要があります(セル中心がΔyなら、セル厚さは約2Δy)。ここを間違えると、意図したy⁺の倍の値になってしまいます。
最後に、物性値や代表速度・長さの設定を軽視すること。例えば、自動車の外気解析で代表長さを「全長」にするか「ホイールベース」にするかで、レイノルズ数が大きく変わり、必要なΔyが一桁違うことも。また、温度変化で空気の動粘性係数νが変わることも忘れがちです。夏と冬の条件で、同じ速度でも必要なメッシュサイズが異なってきます。
空気流れ(ρ=1.225 kg/m³、μ=1.81×10⁻⁵ Pa·s)が長さ0.5mの平板を流速10 m/sで通過する場合、Re≈338,000となります。摩擦係数 Cf≈0.058×Re⁻⁰·² から摩擦速度u*≈0.48 m/sが得られ、y⁺=1を目指すなら壁面第1セル高さΔy≈3.1×10⁻⁵mが必要です。一方、同じ条件(L=0.5 m、U=10 m/s)の水流れ(ρ=998 kg/m³、μ=0.001 Pa·s)ではy⁺=100の壁関数法で約2.8×10⁻⁴mとなります
準拠/参考: 壁法則 \(y^+ = u_\tau\, y/\nu\)、\(u_\tau=\sqrt{\tau_w/\rho}\)。壁面せん断は平板乱流の摩擦係数経験式 \(C_f \approx 0.058\,Re^{-0.2}\)(Schlichting/Prandtl 1/7乗則系)、\(\tau_w=C_f\tfrac12\rho U^2\)、第1セル高さ \(\Delta y = y^+\nu/u_\tau\)。
モデルの前提: 非圧縮・平板境界層・滑面・ゼロ圧力勾配を仮定。物性は一定。算出 Δy は「セル中心」までの距離(多くのメッシャーの第1層厚さは約2Δy)。乱流モデル別の推奨 y⁺ は参考値の表示のみ。
適用範囲・限界: \(Re=10^5\)〜\(10^7\) の乱流平板で有効。曲面・強い圧力勾配・剥離・粗面では \(C_f\) 式が外れ誤差増大。係数 0.058 は丸め値(Schlichting局所値 0.0592、差 <3%)。初期メッシュ設計の目安であり、最終解像度はソルバーのマニュアル推奨で確認のこと。