CFD near-wall y⁺計算機 戻る
流体解析

CFD near-wall y⁺計算機

CFD壁面境界層解像度の指標 y⁺ と必要な第1セル高さ Δy を Reynolds 数・流体・代表長さから即計算。壁関数法と低レイノルズ数モデルの選定基準も一覧で確認できます。

流れ条件

10⁴10⁸
0.01 m100 m
0.1300
必要第1セル厚さ Δy vs Reynolds数(対数-対数)
第1セル厚さ Δy
摩擦速度 u_τ (m/s)
壁面せん断力 τ_w (Pa)
推奨 y⁺ 範囲
理論・主要公式
$C_f \approx 0.058 \cdot Re^{-0.2}$
$\tau_w = C_f \cdot \rho U^2 / 2$
$u_\tau = \sqrt{\tau_w / \rho}$
$\Delta y = y^+ \cdot \nu / u_\tau$

CFD壁面y⁺計算機とは

🙋
CFDで「y⁺」ってよく聞くけど、何ですか?シミュレーションの設定で非常に気にするみたいだけど。
🎓
大まかに言うと、壁面から一番近い計算セルの中心が、壁面からどれだけの「無次元距離」にいるかを表す値だよ。これが適切でないと、壁面近くの摩擦抵抗や熱伝達の計算が大きく大きく狂ってしまうんだ。上のシミュレーターで、レイノルズ数を変えて確認してみて。y⁺が1の時と100の時で、必要な第1セル厚さが大きく異なるのがわかるよね。
🙋
え、そうなんですか!でも、なんでy⁺の「目標値」を自分で決めなきゃいけないんですか?乱流モデルをk-εからk-ω SSTに変えたら、推奨範囲が変わった!
🎓
良いところに気が付いたね。実務では、使う乱流モデルによって壁面の流れをどう扱うかが決まってるんだ。例えばk-εモデルは壁関数を使うのでy⁺を30〜300くらいに設定することが多い。逆にk-ω SSTモデルは壁面まで直接解くのでy⁺≒1にする必要がある。シミュレーターでモデルを切り替えると、推奨範囲が変わるでしょ?これがメッシュ設計の第一歩なんだ。
🙋
なるほど!でも、この計算で出てくる「第1セル厚さ」が例えば0.001mmとか、すごく小さくなってしまうことありますよね。現実的にそんな薄いメッシュ作れるんですか?
🎓
そこがCFDの難しいところで、特に航空機の翼やF1マシンみたいにレイノルズ数が大きいと、セル厚さがミクロン単位になることもある。現場では、全体のメッシュ数を抑えつつ、壁面近くだけ極端に細かくする「境界層メッシュ」を生成するよ。このツールで流体を「水」から「空気」に変えてみ。同じ条件でも必要な厚さが変わるから、物性の影響も体感できるね。

よくある質問

壁関数法を使用する場合は30〜300、低レイノルズ数モデルを使用する場合は1以下(理想的には0.2〜1)を目安にしてください。乱流境界層の粘性底層(y⁺<5)を解像するか、対数則領域(y⁺>30)をモデル化するかで選択します。
Δyが小さすぎると格子数が膨大になります。代表速度や代表長さを適切に見直すか、壁関数法(y⁺=30〜50)を採用してΔyを大きくしてください。また、流体の動粘度が小さい場合は、低レイノルズ数モデルを諦めて壁関数法に切り替える現実的な判断も重要です。
本ツールでは平板乱流境界層の代表的な経験式(Cf≈0.058Re⁻⁰·²)を使用しています。これはRe=10⁵〜10⁷程度の乱流平板流れで有効です。複雑な形状や強い圧力勾配がある流れでは誤差が大きくなるため、あくまで初期設計の目安としてご利用ください。
代表速度には主流の流速、代表長さには平板の長さや翼弦長などの流れ方向の代表寸法を使用します。内部流れの場合は水力直径を用いることが一般的です。レイノルズ数は乱流遷移が生じる目安(Re>5×10⁵)以上で計算してください。

実世界での応用

自動車・航空機の空力設計:車体や機体表面の摩擦抵抗を正確に予測するために、境界層メッシュの設計は極めて重要です。特にモータスポーツでは、数%の抵抗低減が勝敗を分けるため、y⁺に基づいた高品質なメッシュ生成が必須です。

ターボ機械の設計:タービンブレードやポンプインペラーの表面は複雑な流れと高いせん断応力にさらされます。適切なy⁺設定により、剥離や失速現象を精度良く捉え、効率と耐久性を両立させます。

建築・環境風工学:高層ビル周辺の風環境や構造物に働く風荷重を評価する際、地面や壁面近くの乱流を正しくモデル化する必要があります。大規模計算では計算コストと精度のバランスを取りつつy⁺を決定します。

電子機器の熱設計:基板やヒートシンクの強制空冷シミュレーションでは、固体壁面近くの熱伝達率の計算精度が全体の温度予測に直結します。y⁺を適切に設定することで、現実に近い冷却性能を評価できます。

よくある誤解と注意点

まず、「y⁺の目標値は絶対守らなければならない神の数値」と思い込むことです。確かにガイドラインは重要ですが、例えばk-ω SSTモデルでy⁺≒1を目指しても、複雑な幾何形状では全ての壁面でこれを達成するのは不可能に近い。実務では、流れの物理的に重要な領域(剥離点、再付着点、高せん断応力部)で優先的に目標y⁺を満たすようにメッシュを設計します。それ以外の場所では多少外れても許容されることが多いです。

次に、計算した第1セル厚さΔyをそのままメッシャーに入力すればOK、と考える落とし穴。このツールで出てくるΔyは「セル中心」までの距離です。多くのメッシュ生成ソフトでは、「第1層の厚さ」はセルそのものの厚さを指すため、計算値の約2倍を入力する必要があります(セル中心がΔyなら、セル厚さは約2Δy)。ここを間違えると、意図したy⁺の倍の値になってしまいます。

最後に、物性値や代表速度・長さの設定を軽視すること。例えば、自動車の外気解析で代表長さを「全長」にするか「ホイールベース」にするかで、レイノルズ数が大きく変わり、必要なΔyが一桁違うことも。また、温度変化で空気の動粘性係数νが変わることも忘れがちです。夏と冬の条件で、同じ速度でも必要なメッシュサイズが異なってきます。

使い方ガイド

  1. レイノルズ数(Re)と代表長さ(m)を入力し、動粘性係数ν(m²/s)を設定します
  2. 流体密度ρ(kg/m³)と動粘性係数μ(Pa·s)から摩擦速度u*と壁面せん断応力τwを計算します
  3. 目標y⁺値(壁関数法は30~300、低レイノルズ数モデルは0.1~1)に応じて壁面第1セル高さΔy(m)を算出します

具体的な計算例

空気流れ(ρ=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%)。初期メッシュ設計の目安であり、最終解像度はソルバーのマニュアル推奨で確認のこと。

実務での注意点

  1. 層流領域(Re<5×10⁵)では低レイノルズ数モデル採用時にy⁺<1管理が必須です
  2. 乱流境界層(Re>5×10⁵)では壁関数法でy⁺=30~100、粗さを考慮する場合は相当粗さks値との比較が重要です
  3. ANSYS Fluent、OpenFOAMなど解析ツール別に推奨y⁺範囲が異なるため、マニュアルで確認してください