twoPhaseEulerFoam turbulence model
Hello,
I am trying to find out which equations are under the hood of twoPhaseEulerFoam. Most of them can be found in H. Rusches thesis, which is great because to me C++ source code is not self explanatory in the case of fluid dynamics. I have a question about the turbulence model. In H. Rusches thesis (page 105, eq. 3.14 and 3.15) the kequation reads as follows (forgive the LaTeX syntax) \frac{\partial k_b}{\partial t} + (U_b \cdot \nabla) k_b  and so on ... but the source code in kEpsilon.H is: Code:
fvm::ddt(beta, k) + fvm::div(phib, k) + and so on ... \frac{\partial \beta k_b}{\partial t} + \nabla \cdot (U_b k_b)  and so on ... The same is also the case for the epsilon equation. Also the right hand side of both equations seem to have been multiplied by beta. Code:
... == beta*G  fvm::Sp(beta*epsilon/k, k) 
Quote:
I also found this difference in twoPhaseEulerFoam. But when I run the code, it seems that no problem. May be you can modify it to standard ke. Good luck!:cool: 
Quote:

Hello,
you can find H. Rusches thesis amongst others at this page: http://www.foamcfd.org/resources/theses.html 
Quote:
by the way, can u recommend me some papers about twophaseeulerfoam? some of the parameters I dont know the exact mean. 
I haven't really found any research papers dealing with the basics of twoPhaseEulerFoam.
H. Rusches thesis is a good point to start digging into the solver. Also the article about bubbleFoam in the OpenFOAM Wiki helps to understand the solver http://openfoamwiki.net/index.php/BubbleFoam 
Quote:

All times are GMT 4. The time now is 21:24. 