[Bernhard]

Originally Posted by aerogt3

I am a bit lost here with what looks like two conflicting suggestions. If I have a y+1 mesh, and I want to fully resolve the BL (no wall fcns), should I set k = 1e-XX and nut = 1e-XX, or k = kwRWallFunction, with nut set to nutUSpaldingWallFunction?

If you have a proper y+<1 mesh, you should not have problems setting nut=0 and k to zeroGradient. Physically, k goes to zero at the wall with y to some power, so dk/dn=0 or k=0 are equally valid (as long as your mesh is OK, one follows from the other or vice versa). You can check this a posteriori.

For omega I don't know, but I cannot imagine that you want to apply wall functions there. You should look up the limiting behavior of omega(y) for y to 0. I don't know whether OpenFOAM is shipped with these kind of boundary conditions.