Wall y+ Calculator
This y+ calculator will assist you in determining the necessary wall distance for your CFD analysis. The fluid flow near the wall is a complicated phenomenon. Several models and approximations are suggested to estimate the fluid flow parameters in the vicinity of a surface or wall. A CFD package utilizes such approximations and solutions to give us an overall impression of fluid flow around an object or along a surface.
Parameters like skin friction coefficient, y+ wall thickness, and shear stress are crucial to any fluid flow computation. Establishing a reasonably small mesh sizing strategy is necessary, especially in finite element-based tools like ANSYS FLUENT, OpenFOAM, Simscale, or XFlow. The mesh should be small enough in the proximity of walls for the resolution of the boundary layer phenomenon but coarse enough to be solved on the available computational resources.
A structured meshing is used to maintain this delicate balance between the boundary layer resolution and computation time. It is refined near the crucial areas like walls and surfaces, but coarse as we move away from it. But despite localized mesh strategy, often the computational time is very high. The y+ wall thickness method is often used in the case of turbulent flows to solve the near-wall flow problems. This method does not require a very fine mesh and is more effective in turbulent fluid flows that are not fully developed.
Read on to understand how this y+ wall distance is estimated.
What is boundary layer?
Before we get into the modeling of the boundary layer, let's understand what the boundary layer is. Consider a case of fluid flowing over an object, say, a flat plate. The edges of the plate in contact with the fluid are known as boundaries. When a fluid, let's say air, is flowing over the plate, certain particles will stick to the boundary. In that case, the velocity of those particles in the vicinity of the boundary would be the same as the boundary — in some cases even zero, if the boundary is at rest.
This very narrow region of fluid particles near the boundary is called the boundary layer. The boundary layer consists primarily of 3 subregions or layers:
- Laminar boundary layer;
- Transition region; and
- Turbulent boundary layer.
As one moves away from the boundary, the velocity of the particle is governed by the gradient du/dy, which ultimately increases up to the fluid velocity, otherwise known as the free-stream velocity. This velocity gradient causes shear resistance and forms a very thin layer of fluid near the boundary. The thickness of the boundary layer begins from zero at an edge and increases as one moves away from it.
The region where the thickness is very small is known as the laminar boundary layer. It is only considered up to the point where the flow is laminar. The type of flow can be identified using Reynolds number (refer to the Reynolds number calculator for more information). Past this point of the laminar boundary layer, the retarded fluid flow increases the thickness of the boundary layer to such a point that it becomes unstable and leads to the irregular motion of fluid particles. This layer is known as the turbulent boundary layer.
The zone between the two layers where the transition occurs is known as the transition zone. The y+ wall thickness, or wall functions approach, or modeling strategy is primarily used to find the approximate solution to the behavior of fluid particles in the near-wall region of the turbulent boundary layer by computing the shear stress along the wall.
Y+ wall functions approach
The wall function approach is introduced because one of the most common turbulence models, k−ε, is only valid when the turbulent flow is fully developed. The wall functions are introduced to solve for the fluid flow in the near-wall or boundary region. They are used to connect the fluid flow in the near-wall region with fully developed turbulent flow. This method also helps reduce elements in the wall region, assisting in reducing computation time.
The boundary layer is usually meshed with prism-shaped cells to capture the boundary layer and its viscous sublayer effectively. In CFD solvers, this method is often known as inflation layer meshing. The viscous sublayer is often so thin that it is not possible to resolve it using a mesh. Therefore, the y+ wall function method is used.
Now, what is y+? The y+ parameter is the dimensionless wall thickness parameter used to implement the wall functions approach. It is defined as:
y⁺ = (ρ · u* · Δy) / μ = (u* · Δy) / ν
where ρ is the fluid density, μ is the dynamic viscosity, ν = μ/ρ is the kinematic viscosity, Δy is the wall distance (the height of the first cell), and u* is the friction velocity at the nearest wall.
How to calculate wall distance?
To find the first-cell wall distance Δy for a target y+, this calculator works through the following chain of equations for turbulent flow over a flat plate:
- Reynolds number — based on the reference length L:
Re = (ρ · U · L) / μ
- Skin friction coefficient — from the 1/7-power-law correlation:
Cf = 0.026 / Re1/7
- Wall shear stress:
τw = ½ · Cf · ρ · U²
- Friction velocity:
u* = √(τw / ρ)
- Wall distance — rearranging the y+ definition:
Δy = (y⁺ · μ) / (u* · ρ)
As a rule of thumb, choose y+ ≈ 1 when you want to fully resolve the viscous sublayer (wall-resolved mesh, low-Reynolds turbulence models), and y+ ≈ 30–300 when you intend to use standard wall functions.
Example: Using the wall y+ calculator
Suppose air at 20 °C (ρ = 1.204 kg/m³, μ = 1.825 × 10⁻⁵ Pa·s) flows over a 1 m long flat plate at a free-stream velocity of 10 m/s, and you want the first cell to give y+ = 1.
- Re = (1.204 × 10 × 1) / (1.825 × 10⁻⁵) ≈ 6.6 × 10⁵
- Cf = 0.026 / Re1/7 ≈ 0.0049
- τw = ½ × 0.0049 × 1.204 × 10² ≈ 0.295 Pa
- u* = √(0.295 / 1.204) ≈ 0.495 m/s
- Δy = (1 × 1.825 × 10⁻⁵) / (0.495 × 1.204) ≈ 3.06 × 10⁻⁵ m ≈ 0.031 mm
So you would set the first prism-layer height to roughly 0.03 mm. Just pick your unit system, enter the values (or use one of the built-in fluid presets), and the calculator does all the steps for you — showing the result in both metric and American/imperial units.
Units
This calculator supports both the metric (SI) system (m/s, kg/m³, Pa·s, m) and the American / imperial system (ft/s, lb/ft³, lb/(ft·s), ft). Whichever unit system you choose for the inputs, the results are presented in both systems so you can compare instantly.
FAQs
What is a good y+ value for my mesh?
It depends on the turbulence treatment. For wall-resolved simulations (e.g. SST k−ω with low-Reynolds correction), aim for y+ ≈ 1. For wall-function-based modeling (e.g. standard k−ε), the first cell should fall in the log-law region, typically 30 ≤ y+ ≤ 300.
Why is my wall distance so small?
The viscous sublayer is extremely thin, so resolving it to y+ ≈ 1 often requires first-cell heights on the order of micrometers, especially at high Reynolds numbers. This is normal — that is exactly why wall functions were developed, to avoid such fine near-wall meshes.
Which length should I use as the reference length L?
Use a characteristic length representative of the boundary-layer development, such as the plate/chord length, the pipe diameter, or the distance from the leading edge to the location of interest. The Reynolds number — and therefore the wall distance — scales with this choice.
Does this account for laminar flow?
The skin-friction correlation used here (Cf = 0.026 / Re1/7) is intended for turbulent boundary layers, which is where the y+ wall-function approach is relevant. For predominantly laminar flow, a dedicated laminar boundary-layer analysis is more appropriate.