CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Bugs (
-   -   Bug in twoPhaseEulerFoam wallfunctions (

alberto February 9, 2007 14:43

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:

~/OpenFOAM/OpenFOAM-1.3/applications/solvers/multiphase/bubbleFoam/wallFunctions .H

This bug was reported to me by David Huckaby @ NETL/DOE. Thank you David.


david_h February 9, 2007 15:15

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 19:46.