Atmosphericmodels only works properly for default C1 and C2 values

It appears that the atm wall functions and both (r)kEpsilon and kOmegaSST are not in balance when non-default values for C1 and C2 are used.

For example, for higher levels of turbulence (e.g. C1 = 0, C2 = 10), k immediately decays near the wall. This can be reproduced in the atmDownstreamDevelopment tutorial case. Is this intended behaviour?

Can the wall functions for either Epsilon or Nut be modifed to include the appropriate corrections to allow for generation of k at the rate balanced with the specified inlet turbulence?

changed title from Atmospheric models only works for default C1 and C2 values to Atmosphericmodels only works for default C1 and C2 values

changed title from Atmosphericmodels only works for default C1 and C2 values to Atmosphericmodels only works properly for default C1 and C2 values

changed the description

added community label

Hello @crohr - we don't have any comments on the balance issue as we have followed the relevant resources as close as possible during the implementations. If you happen to have any suggestions or code patches we can try, please don't hesitate to share them.

Thank you for the reply! I believe the solution to this would require the implementation of a variable Cmu constant (usually just == 0.09) to allow for the equations to be in balance for any C1,C2. There's quite a bit of validation needed before that could be incorporated here though. I will loop back if I ever make progress on this.

closed