[Xabi]

Hi all!

I am reading through the wall functions code for OF and comparing it with fluent and I noticed some strange things:

1- For the value E instead of using 9.781 as fluent 9 is used.

2- When OF calculates the epsilon for the wall function uses the following formula:
epsilon_[faceCelli] +=Cmu75*pow(k_[faceCelli], 1.5)/(kappa_.value()*RASModel::y_[patchi][facei]);
In Fluent the same formula is used but what is bothering me is the faceCelli term. In fluent epsilon_ is calculated at the point (epsilonp), however, in OF the term epsilon_[faceCelli] seems to refer to the value at the wall.

Thank you in advance,
Xabi