|
[Sponsors] |
February 9, 2007, 14:43 |
There's a bug in the implement
|
#1 |
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 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |