Buggy twoPhaseEulerFoam turbulence Model
I am working on a multiphase Simulation with twoPhaseEulerFoam/compressibleTwoPhaseEulerFoam And have a problem with a buggy turbulence model. My case works fine in buoyantPimpleFoam, but when I try to use wall functions (i.e. epsilonWallFunction) for Epsilon, I always get this message:
When I replace the boundary conditions with "Dummy" to get shown available BC's, I see a lot possibilities. But only the simple ones like zeroGradient work.
It gets even more curious. When I calculate the Epsilon-value with the buoyantPimpleFoam-Case and copy them into the twoPhaseEulerFoam-Case with a 'fixedValue' BC, nothing changes at all.
In the solver/twoPhaseEulerFoam dictionary exist even special files called 'kEpsilon.H' or 'wallFunctions.H'.
But my case still doesn't work properly.
I have to add, I am calculating with kineticEnergy off. I did it because I try to make the dispersed Phase act more like a simple gas. Could this be part of the problem?
I am using OF220.
Thanks you for your attention or help.
I had a quick look at the tutorials for twoPhaseEulerFoam and they do not seem to use the conventional turbulence models. In addition, I checked to which libraries the solver links to, by running this command:
This explains why you're seeing that error message.
Problem is: why don't even the tutorials for this solver use the specific turbulence models, but still use the k-epsilon fields?
Specially, if we look at the output from this command:
But the thing is... why is twoPhaseEulerFoam able to compute the k-epsilon fields without these libraries in the first place? So lets try looking at the source code: https://github.com/OpenFOAM/OpenFOAM...PhaseEulerFoam
:eek: This feels like really old OpenFOAM code... the turbulence model is coded into the solver directly (have a look at the file "kEpsilon.H").
OK, from what I can figure out: this solver is too specifically coded for certain simulation scenarios. You cannot use the generic case preparation that most of the other solvers can handle, because this solver seems to be very restricted to what it was developed to simulate.
I had a quick look at the solver multiphaseEulerFoam and is relies on LES for the turbulence modelling, instead of having it coded directly: https://github.com/OpenFOAM/OpenFOAM...phaseEulerFoam
Of course the issue now is: which one is more suitable to your simulation?
I give up and start something else
Thanks wyldckat for your advice. I inserted those libraries in the controlDict.
I tried a second time bye changing a version of the solver itself bye looking on the buoyantPimpleFoam Solver and inserted the missing library entries in the main file.
I inserted in the EEqn.H the missing turbulence part
Both were missing in the compressibleTwoPhaseEulerFoam
But after adding them... i gain still got the same error message.
So I can give up now peacefully in trying further implementing the turbulence in this solver. I will start in creating an own solver.
|All times are GMT -4. The time now is 16:21.|