If you get the boundary layer mesh wrong, your simulation results can be way off. y+ is the main number to watch—it tells you if your mesh near the wall is in the right flow region. Use the Y-Plus Interactive Calculator below to get y+ itself, the wall distance you’ll need, friction velocity, and wall shear stress based on your input fluid properties, velocity, and characteristic length. This matters for aerodynamics, turbomachinery, and heat exchange—small missteps in y+ can mean huge errors in skin friction or local heat transfer. The page explains all the equations, gives a flat-plate example, runs through the relevant CFD wall theory, and includes an FAQ for practical reference.
What is y+ (y-plus)?
y+ is a way to measure how close your first mesh cell is to the wall, but it doesn’t use physical units like millimeters. Instead, it scales the distance by the thickness of the viscous layer in turbulent flow—a scale set by viscosity and wall shear. CFD solvers use this value to determine if they should actually resolve the near-wall flow, or use a so-called wall function to approximate it.
Simple Explanation
Think of y+ as a non-dimensional ruler set by the fluid’s wall friction, not by length units. If y+ is 1, your cell is inside the viscous sublayer—viscosity dominates. When y+ is 50, you’re deep in the turbulent layer, and wall functions are needed. The “messy” range is between about 5 and 30; neither direct resolution nor wall functions work properly there.
📐 Browse all 1000+ Interactive Calculators
Contents
Boundary Layer Diagram
Y-Plus Calculator
How to Use This Calculator
This calculator is intended for education, concept evaluation, and preliminary design. Results are based on the equations and assumptions described on this page, but cannot account for every real-world load case, tolerance, material property, environmental condition, installation detail, safety factor, code, or regulatory requirement. Verify all inputs, assumptions, units, and results independently before selecting components or using the result in a real application. Safety-critical, structural, medical, lifting, transportation, or regulated applications must be reviewed by a qualified engineer.
- Pick what you need to solve—y+, wall distance, friction velocity, or cell height.
- Enter the fluid properties: density (ρ), dynamic viscosity (μ), and characteristic length (L).
- Fill in the operating conditions for your chosen calculation mode—wall distance, velocity, target y+, etc.
- Hit Calculate to view your answer.
Y-Plus Interactive Visualizer
Visualize how wall distance and flow conditions affect y+ values in the boundary layer. Watch the viscous sublayer, buffer layer, and logarithmic regions change as you adjust parameters.
Y-PLUS
12.5
FRICTION VEL
1.85 m/s
WALL SHEAR
4.2 Pa
FIRGELLI Automations — Interactive Engineering Calculators
Equations & Variables
Here’s the core formula to get y+ from mesh spacing and wall friction velocity:
y+ = (y · uτ) / ν
uτ = √(τw / ρ)
τw = 0.5 · Cf · ρ · U∞2
Cf ≈ 0.058 · ReL-0.2 (turbulent flat plate)
y+ = Dimensionless wall distance (no units)
y = Physical distance from wall to first grid point (m)
uτ = Friction velocity (m/s)
ν = Kinematic viscosity = μ/ρ (m²/s)
τw = Wall shear stress (Pa or N/m²)
ρ = Fluid density (kg/m³)
μ = Dynamic viscosity (Pa·s)
U∞ = Free-stream velocity (m/s)
Cf = Local skin friction coefficient (dimensionless)
ReL = Reynolds number based on characteristic length (dimensionless)
L = Characteristic length such as plate length or chord (m)
Simple Example
Take air at 50 m/s across a 1 m flat plate. Density: 1.225 kg/m³, dynamic viscosity: 1.81 × 10-5 Pa·s. If your first cell is 0.05 mm (0.00005 m) from the wall:
Re = 1.225 × 50 × 1 / 1.81 × 10-5 ≈ 3.38 × 106 → Cf ≈ 0.00320 → τw ≈ 4.90 Pa → uτ ≈ 2.00 m/s
y+ = (0.00005 × 2.00) / (1.478 × 10-5) ≈ 6.8 — this is the buffer layer. Dropping to about 0.007 mm will get y+ ≈ 1.
Theory & Practical Applications
Boundary Layer Structure and the Origin of y+
y+ comes out of scaling the near-wall, turbulent boundary layer. Here, the right length scale isn’t the boundary layer thickness, but the viscous length ℓν = ν/uτ. Divide the mesh wall distance y by this to get y+. This coordinate shows why certain layers matter: y+ < 5 is the viscous sublayer (linear profile, dominated by viscosity), 5 < y+ < 30 is the buffer layer (both viscous and turbulent effects), and y+ > 30 is the log-law region.
Why care about y+ in practice? The structure of these wall layers is nearly the same by y+ across a wide range of Reynolds numbers. But in actual units, as wall shear rises, your physical sublayer gets thinner. For a commercial aircraft wing in cruise, y+ = 1 means Δy ≈ 0.003 mm—very hard to resolve. That's why wall functions are common: they sidestep the need for fine mesh by modeling this region instead of resolving it.
Wall Treatment Strategies in CFD
Your choice of wall treatment sets your mesh and model. If you want to fully resolve the wall (low-Re turbulence models, LES), you’ll need y+ < 1 for the first cell, and several points within y+ < 5. This gives the best friction and heat transfer result, but the mesh size quickly becomes a problem for high Reynolds numbers—a single industrial case could take trillions of cells. Wall functions, on the other hand, let you start with cells at y+ = 30–100 and use empirical laws for the near-wall behavior. The usual law of the wall: u+ = y+ (below y+ ≈ 11.6), u+ = (1/κ)ln(y+) + B for higher y+ (with typical κ ≈ 0.41, B ≈ 5.2).
The main failure is when the cell lands in the “buffer” region (5 < y+ < 30)—not quite viscous, not quite logarithmic. Here, both physical and empirical assumptions fall apart and errors are large. Some CFD codes blend models, but even those work poorly here. Trouble gets worse with pressure gradients or separation; in these settings, you often need wall-resolved meshes even though it hurts computing cost.
Estimating Required First Cell Height
If you’re planning mesh, you’ll need to estimate first cell height Δy to hit a target y+. That means estimating uτ (and that in turn needs τw, which you don’t always know before running the simulation). For simple cases—flat plates, streamlined bodies—a skin friction coefficient does the job: Cf ≈ 0.058Re-0.2 is a fast estimate in turbulent boundary layers. Once you’ve got Cf, it’s just τw = 0.5CfρU∞2, find uτ = √(τw/ρ), and then Δy = (y+·ν)/uτ.
For pipe flow, use a friction factor from Moody or Colebrook, then τw = (f/4)ρU2. If you don’t have a formula—complex shapes—sometimes you do a coarse run first to get rough wall shear, then remesh for proper y+. Adaptive mesh methods can help, but working with wall-normal cell refinement is tricky and not all software handles it well. Also, wall shear can change by 3–5× over a surface, so “matching” y+ everywhere with fixed cell height is rarely possible; most meshers don’t do this automatically.
Industry-Specific Applications and Non-Standard Cases
In aerospace, y+ control is central to accurate drag simulation—skin friction is a big piece of total drag (up to half). Wind tunnel data shows that going above y+ ≈ 2 at the wall grid can give a meter-scale error in drag, enough to knock new designs off target. Rough transition prediction (laminar-turbulent) may need y+ < 0.5 and tight streamwise nodes. For hypersonics, whole new physics sets in (rarely handled with conventional CFD) and y+ has to be replaced with other non-dimensional numbers.
For turbomachinery, y+ is still critical, but practical problems get worse: blades run at high Re and have difficult surfaces. Getting Δy < 0.01 mm on a finished blade can be unrealistic, so practical meshes use thin prism layers at the wall, coarser tetrahedral elements in the core, and often have to compromise. Combustion means changing density and new source terms. If density at the wall changes much (compressibility), you need alternate scaling (like y* = yuτρ/μ). Temperatures shift viscosity by large factors—important in supersonic and heat transfer problems. Keep these effects in mind for real engineering jobs.
Fully Worked Multi-Part Example: External Flow Over a Flat Plate
Problem Setup: Consider a flat plate 2.4 meters long in a wind tunnel at U∞ = 67.3 m/s, air at 288 K. The goal is to get y+ = 0.8 at the trailing edge (where wall shear is lowest), and check that you’re still below y+ = 1 at the leading edge (where wall shear is highest).
Given Data:
- Free-stream velocity: U∞ = 67.3 m/s
- Plate length: L = 2.4 m
- Temperature: T = 288 K (standard conditions)
- Air density: ρ = 1.225 kg/m³
- Dynamic viscosity: μ = 1.81 × 10-5 Pa·s (from Sutherland's law at 288 K)
- Target y+ at trailing edge: y+TE = 0.8
Step 1: Calculate Reynolds Number
ReL = (ρ U∞ L) / μ = (1.225 kg/m³)(67.3 m/s)(2.4 m) / (1.81 × 10-5 Pa·s) = 1.094 × 107
This tells us the flow is fully turbulent all the way down the plate.
Step 2: Estimate Skin Friction Coefficient at Trailing Edge
With the flat plate rule: Cf(x=L) ≈ 0.058 ReL-0.2 = 0.002744
Use this at the trailing edge (x = L).
Step 3: Calculate Wall Shear Stress at Trailing Edge
τw,TE = 0.5 × Cf × ρ × U∞2 = 7.617 Pa
Note: The bigger the boundary layer, the smaller τw at the trailing edge.
Step 4: Calculate Friction Velocity at Trailing Edge
uτ,TE = √(τw,TE / 1.225) = 2.495 m/s
Step 5: Determine Required First Cell Height
ν = μ / ρ = 1.478 × 10-5 m²/s
ΔyTE = (0.8 × 1.478 × 10-5 m²/s) / 2.495 m/s = 0.004738 mm
Step 6: Verify y+ at Leading Edge
At the leading edge (use x ≈ 0.1 m):
Rex = 4.558 × 105
Cf(x=0.1m) = 0.004209
τw,LE = 11.68 Pa
uτ,LE = 3.089 m/s
y+LE = (4.738 × 10-6 m × 3.089 m/s) / (1.478 × 10-5 m²/s) = 0.99
Conclusion: Using Δy = 0.00474 mm gives y+ = 0.8 (trailing edge) and y+ = 0.99 (leading edge). A 24% swing in y+ is quite normal. You’ll need 10–15 cell layers through the viscous sublayer—watch the normal spacing; growth ratios above ~1.2 are likely to miss needed resolution.
Advanced Considerations: Roughness, Compressibility, and Heat Transfer
Surface roughness adds its own scale. Once the non-dimensional roughness height ks+ gets above 5, wall friction goes up and the log-law profile shifts. Most CFD wall models assume smooth walls; using them on rough surfaces causes systematic errors unless you specify the equivalent sand-grain roughness. Real-world blades and pipes are rarely smooth, and roughness can hit performance hard if not modeled.
Compressible flows add further complexity: varying density through a boundary layer shifts everything. At Mach 3 and above, the wall temperature alters both density and viscosity, so the usual y+ and friction velocity don’t describe the scaling unless you use compressibility corrections (like Van Driest). For flows with strong wall heating/cooling, you’ll also need to resolve the thermal sublayer—sometimes thinner than the velocity sublayer, especially for Prandtl numbers far from 1. This matters a lot in high-speed heat exchangers, turbine blades, and electronics cooling, where missing the right y+ (or thermal y+) can throw off both local and integrated predictions.
With commercial CFD codes, understand how your wall models actually work. Some blend low-Re and wall function methods dynamically, but still break down in the 5–30 buffer region. Some codes (like OpenFOAM’s nutUSpaldingWallFunction) cover the whole range, but accuracy for y+ in the “gray” zone is still poor. For any serious CFD work, always report y+ distributions—spot-checking a few values isn’t enough. Run mesh refinements to be sure your results don’t change noticeably with even finer near-wall spacing.
Frequently Asked Questions
Free Engineering Calculators
Explore our complete library of free engineering and physics calculators.
Browse All Calculators →🔗 Explore More Free Engineering Calculators
- Reynolds Number Calculator — Laminar or Turbulent
- Hydraulic Pressure Drop Calculator — Pipe Hose
- Pneumatic Valve Cv Flow Coefficient Calculator
- Venturi Flow Meter Calculator
- Porosity And Permeability Calculator
- Poise Stokes Converter Calculator
- Oblique Shock Calculator
- Pump Head Calculator — Total Dynamic Head
- Lead Screw Torque and Force Calculator
- AGV/AMR Wheel Friction & Traction Calculator
About the Author
Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations
Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.
📹 Video Walkthrough — How to Use This Calculator
📹 Video Walkthrough — How to Use This Calculator
Need to implement these calculations?
Explore the precision-engineered motion control solutions used by top engineers.
