Target y+ 1 or 30: Practical CFD Rules for Mesh and Heat
Y plus (y+) is a dimensionless number that tells you how close your first mesh cell sits to a wall, calculated as y+ = y uτ / ν. As a practical target, aim for y+ near the lower end of its range to resolve the viscous sublayer directly, and near or above the higher end of the range when you plan to rely on wall functions. A worked example and a setup checklist follow below.
TL;DR:
Target y+ around 1 to resolve the viscous sublayer directly, requiring a mesh that’s roughly 30 times thinner than a y+ of 30 target.
Local wall shear stress variations can cause significant y+ differences at the same mesh height, especially near leading edges or in accelerating flows.
Using thermal wall functions accurately depends on turbulent Prandtl number and flow conditions; buoyancy effects can invalidate standard assumptions.
Check the wall treatment flag in the solver matches your intended y+ range, as mesh defaults may not align with your setup goals.
Regularly verify the actual Yplus across the wall surfaces after each run, since flow changes can shift y+ into unintended regions.
Table of Contents
What y plus is: formula, variables, and solver outputs
The formula is straightforward once you unpack each term. According to ANSYS CFX documentation, y+ = y uτ / ν, where y is the distance from the wall to the centroid of the first mesh cell, uτ is the friction velocity defined as the square root of wall shear stress divided by fluid density, and ν is the kinematic viscosity of the fluid.
The tricky part is that uτ depends on wall shear stress, which varies across a model even when the geometric wall distance stays the same. For example, a flat plate near the leading edge and the same plate further downstream can show differing y+ values at the same first-cell height, simply because the local shear changes.
y: distance from the wall to the first cell centroid, set by your mesh.
uτ: friction velocity, derived from wall shear stress and fluid density.
ν: kinematic viscosity, a fluid property that depends on temperature.
ANSYS CFX draws a distinction worth remembering: the solver can report both a standard Yplus and an internal Solver Yplus. The documentation recommends checking the standard Yplus output when judging near-wall spacing, since the internal variant may reflect solver specifics rather than the physical quantity intended.
Why y plus regions matter for turbulence modeling
The value of y+ at a given point tells you which part of the near-wall flow that cell sits in, and each region behaves differently. The Sandia/INL Fuego wall-function documentation lays out the boundaries used across most CFD codes.
y+ ≤ 5: the viscous or sub-laminar sublayer, where molecular viscosity dominates and velocity varies almost linearly with distance from the wall.
5 < y+ < 30: the buffer region, a transition zone where neither the viscous nor the log-law assumptions hold cleanly.
y+ ≥ 30: the logarithmic or inertial region, where the log-law of the wall gives a reliable velocity profile.
Wall functions exist because resolving the viscous sublayer directly is expensive. Instead, they prescribe the near-wall shear and heat-transfer behavior based on empirical profiles, which works well when the flow is close to equilibrium. The same Sandia/INL documentation notes that this approximation becomes less reliable under separation or strong non-equilibrium conditions, where the assumed velocity and temperature profiles no longer match reality.
The NASA Turbulence Modeling Resource frames y+ selection as part of a broader verification and validation exercise rather than a fixed rule: the right target depends on the turbulence model, the flow case, and what you are validating against.
How y plus dictates solver setup: wall functions vs. low-y meshes
Once you know your target y+, the mesh and turbulence model follow from it, not the other way around.
Target y+ around 1 when you need to resolve the viscous sublayer directly, typically with low-Reynolds-number turbulence formulations or when wall-adjacent gradients (heat flux, skin friction) matter for your results.
Target y+ of 30 or higher when standard wall functions are acceptable, which assumes the flow near that wall is close to equilibrium and free of major separation.
Avoid the buffer zone (5 to 30) as a deliberate target. Neither the viscous-layer assumptions nor the log-law hold reliably there, so cells that land in this range give you the worst of both approaches.
Check your solver’s wall-function flags explicitly rather than trusting mesh defaults, since the same mesh can be interpreted differently depending on how the near-wall treatment is configured.
Choosing low-y meshing buys you more physical fidelity near walls at the cost of far more cells and smaller time steps. Wall functions trade some of that fidelity for a mesh you can actually run on a reasonable budget, which is why most production CFD work defaults to wall functions unless the near-wall physics is the point of the study, as explained in the practical insights on how near-wall behavior affects material processing from ASTRA CHEMICAL. For more on how wall treatment choices tie into boundary condition setup generally, see our practical guide to CFD boundary conditions.
Pro Tip: Check the standard Yplus field after every run, not just once during mesh setup. Local flow changes (separation, recirculation, heat load) can shift wall shear enough to move cells into a different y+ region than you planned for.
Calculating first-cell height for a target y plus
Working backward from a target y+ to a mesh is the most useful skill in this whole topic, and the arithmetic is simple once you have a shear estimate.
Estimate uτ from a reference velocity and a rough skin-friction correlation, or from a preliminary coarse run, as recommended by CFD-Wiki’s guidance on dimensionless wall distance.
Gather ν for your fluid at the relevant temperature; for air near room temperature this is roughly 1.5 × 10⁻⁵ square meters per second.
Solve for y by rearranging the formula to y = y+ ν / uτ.
Say your estimated friction velocity uτ is 0.5 meters per second and ν is 1.5 × 10⁻⁵ square meters per second. For a low target y+, the computed first-cell height will be a certain very thin value. For a higher target y+, this height will be roughly an order of magnitude larger. This serves as an illustrative calculation only; actual values depend on your specific uτ and ν.
One dimensionless check decides your mesh cost by an order of magnitude or more: the y+ = 1 target here demands a first cell roughly 30 times thinner than the y+ = 30 target, from the same dimensionless wall distance formula.
After solving, read the Yplus field across every wall surface, not just at a single probe point. A surface average can hide pockets where y+ spikes near separation or drops unexpectedly near stagnation points, both signs your mesh and wall treatment may not match in that region.
Thermal simulation caveats for y plus
Momentum y+ confirms your mesh matches your velocity wall treatment, but it says nothing directly about whether your wall heat flux is accurate. Thermal wall functions use their own blending across the sub-laminar, buffer, and log regions, and that blending depends on thermal diffusivity and the turbulent Prandtl number, according to the Sandia/INL wall-function documentation.
Turbulent Prandtl number governs how momentum and thermal boundary layers relate; getting this wrong skews predicted wall temperatures even with a well-resolved momentum y+. Our guide to computing the Prandtl number walks through that calculation.
Buoyancy-driven flows can break standard wall-function assumptions outright. A Purdue study on wall functions for indoor airflow found that standard thermal wall functions can underpredict heat transfer when buoyancy effects are significant, and that buoyancy-aware formulations improve results when calibrated against experimental data.
Validate against wall heat flux and local Nusselt number, not just velocity y+, since these metrics directly reflect how well your thermal model is capturing near-wall behavior.
Pro Tip: Treat a buoyant or high-heat-flux case as its own validation exercise. A mesh tuned purely for momentum y+ can still produce a wall temperature that is meaningfully wrong.
For broader context on matching solver setup to heat-transfer physics, see our heat transfer modeling guide.
Mesh and solver checklist for avoiding y plus problems
Most y±related headaches trace back to a handful of avoidable mistakes.
Pick your target y+ before meshing, not after, so the first-cell height is a design input rather than an afterthought.
Check cell aspect ratio and growth rate near the wall. Extremely skinny cells can cause numerical trouble even when y+ looks right.
Confirm the solver’s wall-function flag matches your intended treatment; a mismatch here silently invalidates the mesh you built.
Re-run and remesh locally where y+ swings wildly, rather than refining the entire domain.
Reconsider the turbulence model if separation or strong pressure gradients push y+ out of your intended range across a large area; a generalized wall function with pressure-gradient terms may fit better than the standard log-law assumption in those zones. Our notes on relaxation settings for conjugate heat transfer cover related solver tuning.
How Jewlztech thinks about y plus in practice
We treat y+ as a compatibility check between mesh and wall treatment, never as a standalone score of mesh quality. A model can hit the “right” y+ everywhere and still mispredict heat transfer if the thermal wall function or turbulence model doesn’t match the physics at hand. That framing shapes how we build calculators into our own tools: y+ is one input among several, not the finish line.
— Joel
Get hands-on with near-wall and thermal calculations
Working through a y+ calculation by hand is a good way to learn the physics, but running it across dozens of surfaces and design iterations gets tedious fast. We built the Jewlz Engineering Toolkit to take that repetition off your plate, with a property database covering fluid viscosity and thermal diffusivity across a wide temperature range so you are not hunting for ν values mid-project.
Near-wall and thermal calculators that sit alongside full conduction, convection, and radiation models.
A built-in property database so kinematic viscosity and other fluid properties are on hand without separate lookups.
CFD and thermal modules under one Engineering Toolkit, alongside our CFD simulation software for fluid flow analysis.
Try the free tools first, then explore the full toolkit when you need deeper thermal or CFD analysis.
FAQ
What is y plus in CFD, in plain terms?
Y plus is a dimensionless measure of how far your first mesh cell sits from a wall, relative to local flow conditions. It combines wall distance, friction velocity, and kinematic viscosity into one number so you can judge whether your mesh matches your intended near-wall turbulence treatment, as described in ANSYS CFX documentation.
What y plus value should I target?
Target y+ around 1 when resolving the viscous sublayer directly, and y+ of 30 or higher when using standard wall functions, based on the region boundaries in the Sandia/INL wall-function documentation. Avoid deliberately targeting the buffer zone between 5 and 30, since neither modeling approach works reliably there.
Why does the same mesh give different y plus values in different places?
Because y+ depends on local wall shear stress through the friction velocity term, not just on geometric distance. The CFD-Wiki entry on dimensionless wall distance notes that regions with higher local shear, like near a leading edge or in accelerating flow, produce higher y+ at the same first-cell height than regions with lower shear.
Does a good y plus guarantee accurate heat transfer results?
No. Momentum y+ only confirms your mesh suits your velocity wall treatment. Accurate wall heat flux also depends on the turbulent Prandtl number and the thermal wall function’s own assumptions, and a Purdue study on buoyant indoor airflow found standard thermal wall functions can underpredict heat transfer when buoyancy effects are strong.
How do I calculate the first-cell height for a target y plus?
Rearrange the formula to y = y+ ν / uτ, using your fluid’s kinematic viscosity and an estimated friction velocity from a reference velocity or a preliminary run. This approach is detailed in the CFD-Wiki wall distance guidance, and the result should be confirmed by checking the solved Yplus field afterward.
Sources
Recommended
Comments