Hello everyone, I'm trying to do some turbulence modelling and I found out that the standard or recommended way to calculate initial values for the internal mesh of the turbulent quantities doesn't converge, I'm now guessing the values but seems a bit random. I'm using Bounding Gauss Linear Upwind Grad(U) scheme for every quantity ( U, omega, k ) and sometimes I get good convergence, but the results show strange phenomena for high Reynolds Numbers, for instance, I'm looking for a trend on wake zone dimensions for a surface mounted bluff body to different Reynolds Numbers but for a semi-cylinder, the trend seems to slow or decrease at some point or oscillate for bigger Reynolds Numbers ( which doesn't make any sense to me), any help with this would be very appreciated, thanks in regards.