[yousefaz]

Hello,

I am trying to simulate a homogeneous atmospheric boundary layer in FLUENT using steady state RANS. I am using a Standard k-e model with inlet UDF functions for velocity, k, and epsilon as specified by Richards and Hoxey (1993).

For the most part, I seem to have set up the model correctly. Here is a comparison of the inlet and outlet profiles (which would collapse for a perfectly homogeneous ABL):



However, there seems to be a divergence issue at the bottom corner of the domain near the outlet (right hand side), which particularly affects the kinetic energy throughout the entire domain. It seems to occur in the boundary layer. The solution is fully converged with residual values dropping to 10E-7.

Here is the velocity contour (note the blue (=0 velocity) on the bottom right):



There is a concentration of negative pressure at that corner. Interestingly, no warnings of reverse flow appeared:



The turbulent kinetic energy contour shows values blowing up near that corner:




Zoomed pics with mesh shown:

Pressure:



TKE:



Epsilon:




Any ideas on how to resolve this issue would be greatly appreciated.
Thanks in advance.