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. Alberto |
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. Dave |
All times are GMT -4. The time now is 06:03. |