But again, regarding y+, the Shaumann BC involve applying a shear stress at the wall and it is function of the roughness lenght (z0), which mean its basically using this no? $$ u_\tau = \sqrt{\frac{\tau_w}{\rho}} $$