Congratulations! You’ve made it through the first step of the CFD analysis process and your CAD geometry is ready for meshing! One of the first questions CFD engineers ask themselves at this phase is: “What are my grid requirements going to be?” The focus of this blog will be on the application of the wellknown y^{+} value and how it may be used to estimate your first cell height normal to the wall.
To begin to understand the importance and application of the y^{+} calculation, we need to consider a generalized view of the structure of the boundary layer, and understand a few simple concepts. The boundary layer itself is either classified as laminar, turbulent or transitional as illustrated in the Figure 1. If we view the boundary layer on a semilog scale like the one found in Figure 1, we see that the flow in this region may typically be decomposed into a viscous sub layer, buffer layer, log layer and outer layer based on the distance from the wall and u^{+} value. The nondimensional velocity parallel to the wall is given by u^{+} and is divided by the friction velocity u_{f}.
For simplicity’s sake, let’s assume our flow is fully turbulent. We need to decide if we want to resolve the viscous sublayer, which is common for flow exhibiting adverse pressure gradients, external aerodynamics and any flows were the wall bounded effects are of importance. If our problem of interest lies in this arena, we must consider y^{+ }values on the order of unity.
Figure 1: Boundary Layer Details (Tutor Help Desk/Wikipedia)

If we decide wall bounded effects are of secondary importance, we typically march forward with wall treatments  commonly referred to as wallfunctions. For these applications, it is desirable to place the first cell no closer to the wall than a y^{+} of approximately 10 to ensure we are adequately outside the viscous sublayer. The wallfunction algorithm will model the shape of the relevant boundary layer profiles using an approach such as lawof thewall.
For example, let’s consider a highspeed Mach=1.5 external flow over a cylinder of diameter 1 meter with air as the working fluid at operating at sealevel conditions. To estimate the first cell height, we begin by calculating the Reynold’s number:
where is the fluid density, U the freestream velocity, L is the characteristic length and the absolute viscosity. The y^{+} value is then given by
where the first cell height is given by . The frictional velocity and wall shear stress will be found from
The skin friction coefficient (C_{f}) is determined empirically and substituted accordingly. The resulting Reynold’s number indicates that the flow is fully turbulent so I will assume C_{f }= 0.058Re^{0.2}. For a more thorough list of empirical skinfriction formulas, check out CFD Online. The output of the
assuming the desired y
^{+} of one for this scenario results in:
Though the y^{+ }is derived for turbulent boundary layers, I still commonly use this calculator for laminar flows as well. I simply swap the skinfriction coefficient to get a good first estimate of the required cell height. It is important to note this calculator is intended as a guide and the mesh should be run a few iterations and interrogated to ensure proper resolution by plotting y+, eddyviscosity contours and velocity vectors in your boundary layer to ensure adequate gradient resolution. This is a future topic of discussion, so for now, happy gridding!