Bug in twoPhaseEulerFoam wallfunctions

alberto · February 9, 2007, 14:43

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

david_h · February 9, 2007, 15:15

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

February 9, 2007, 14:43	There's a bug in the implement	#1
alberto Senior Member Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36	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 __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

February 9, 2007, 15:15	All the credit should go to my	#2
david_h Member E. David Huckaby Join Date: Mar 2009 Posts: 57 Rep Power: 17	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

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Wallfunctions onoff	gjesing	OpenFOAM Running, Solving & CFD	5	February 11, 2010 00:55
TwoPhaseEulerFoam	sara	OpenFOAM Running, Solving & CFD	2	November 6, 2008 20:26
Bug in twoPhaseEulerFoam	alberto	OpenFOAM Bugs	2	May 20, 2008 22:25
TwoPhaseEulerFoam Bug	alondono	OpenFOAM Bugs	1	February 19, 2008 21:01
TwoPhaseEulerFoam	newbee	OpenFOAM	0	March 27, 2006 09:41