透過 MCP 連接 →

輸入計算

數學公式

數學公式: y+ 計算器
Show calculation steps (1)
  1. Friction velocity

    Friction velocity: y+ 計算器

    Friction (shear) velocity from wall shear stress and density

廣告

結果

無因次壁面距離
0.9669
y+
摩擦速度 u_τ 0.285714 m/s

什麼是 y+?

無因次壁面距離(寫作 \(y^+\),讀作 y-plus)是計算流體力學(CFD)中極為關鍵的參數,用來描述紊流邊界層中第一層網格相對於黏性長度尺度的高度。透過 \(y^+\),你可以判斷近壁網格究竟解析到了黏性次層、緩衝層,還是落在對數律區、必須仰賴壁面函數(wall function)。

近壁網格層示意圖,顯示從固體壁面到流場中第一個單元的距離 y
\(y^+\) 表示第一個網格單元到壁面的無因次距離。

如何使用這個計算器

請輸入流體密度 \(\rho\)、動力黏度 \(\mu\)、壁面剪應力 \(\tau_w\),以及壁面到第一層網格中心的距離 \(y\)。本工具會先算出摩擦速度 \(u_\tau\),再求得 \(y^+\)。一般建議:採用全解析(低雷諾數)模型時,目標 \(y^+ \approx 1\);若使用壁面函數,則約落在 30~300 之間較為理想。

公式說明

首先由壁面剪應力求得摩擦速度:

$$u_\tau = \sqrt{\frac{\tau_w}{\rho}}$$

這個速度尺度可用來描述近壁紊流的特性。接著

$$y^+ = \frac{\rho \, u_\tau \, y}{\mu}$$

等於把實際距離 \(y\) 以黏性長度 \(\nu/u_\tau\) 加以無因次化,其中 \(\nu = \mu/\rho\) 為運動黏度。換句話說,也可寫成 \(y^+ = u_\tau \, y / \nu\)。

近壁速度分佈圖,依 y+ 分為黏性底層、緩衝層與對數律區
邊界層的各區域對應不同的 \(y^+\) 範圍。

範例計算

以空氣為例,\(\rho = 1.225 \ \text{kg/m}^3\)、\(\mu = 1.81 \times 10^{-5} \ \text{Pa}\cdot\text{s}\)、\(\tau_w = 0.1 \ \text{Pa}\)、\(y = 5 \times 10^{-5} \ \text{m}\)。則

$$u_\tau = \sqrt{\frac{0.1}{1.225}} = 0.28571 \ \text{m/s}$$$$y^+ = \frac{1.225 \times 0.28571 \times 5 \times 10^{-5}}{1.81 \times 10^{-5}} \approx 0.967$$

由於目標是 \(y^+ \approx 1\),在此情況下,你的第一層網格位置已經相當理想。

常見問題

\(y^+\) 該設定多少才合適? 全解析模型(例如 k-ω SST 低雷諾數模型)約為 1;若採用標準壁面函數,則建議落在 30~300。

為什麼 \(y^+\) 是無因次的? 因為它是壁面距離除以黏性長度尺度 \(\nu/u_\tau\),單位互相抵消,自然成為無因次量。

如果不知道 \(\tau_w\) 怎麼辦? 可由表面摩擦係數估算:\(\tau_w = \tfrac{1}{2} \rho U^2 C_f\),其中 \(C_f\) 可採用經驗關係式,例如 \(0.058 \cdot Re^{-0.2}\)。

最後更新: