There's a bug in the implement
There's a bug in the implementation of wall-functions implemented in twoPhaseEulerFoam.
In wallFunctions.H, the dissipation epsilon is not set to zero as the production (line 20) and then it's used with a += operator (line 52).
The correct implementation is provided in bubbleFoam:
This bug was reported to me by David Huckaby @ NETL/DOE. Thank you David.
All the credit should go to my
All the credit should go to my colleague Sofiane Benyahia.
While your in looking at this routine you might also check the (yPlus > 11.6) statement. The value 11.6 might not be consistent with the values of kappa and E. The single phase turbulence models include a function to calculate this.
|All times are GMT -4. The time now is 23:33.|