CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Pre-Processing

Buggy twoPhaseEulerFoam turbulence Model

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   December 1, 2013, 05:52
Default Buggy twoPhaseEulerFoam turbulence Model
New Member
Join Date: Jul 2013
Posts: 14
Rep Power: 5
Kenna is on a distinguished road
Hello everyone,

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:


    request for turbulenceModel turbulenceModel from objectRegistry region0 failed
    available objects of type turbulenceModel are


    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/andreas/OpenFOAM/OpenFOAM-2.2.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 136.

FOAM aborting

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.

Kenna is offline   Reply With Quote

Old   December 1, 2013, 10:03
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,748
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Greetings Andreas,

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:
ldd $(which twoPhaseEulerFoam)
and it indicates that it does not link to any of the turbulence related libraries.
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:
ldd $(which interFoam)
it does indicate that it links to the libraries needed for these models, namely:
Now, the hackish way to add these libraries to the application is to edit the file "system/controlDict" and add this piece of code:
libs (
This will tell the solver to load these libraries, before it starts to try to solve anything .

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:
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:
Of course the issue now is: which one is more suitable to your simulation?

Best regards,
wyldckat is offline   Reply With Quote

Old   December 2, 2013, 08:43
Default I give up and start something else
New Member
Join Date: Jul 2013
Posts: 14
Rep Power: 5
Kenna is on a distinguished road
Thanks wyldckat for your advice. I inserted those libraries in the controlDict.

libs ( );
But got the same error.

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
      - fvm::laplacian(turbulence->alphaEff(), he)
and the same for the velocity UEqn.H
      + turbulence->divDevRhoReff(U)
and adjusted a phi into it, to make it fit with the rest of the code.

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.
Kenna is offline   Reply With Quote


turbulence model, twophaseeulerfoam, wall function

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam case with SpalartAllmaras turbulence model implemented nedved OpenFOAM Running, Solving & CFD 2 November 30, 2014 23:43
turbulence model for Refrigerator cicatrix Main CFD Forum 0 October 3, 2012 01:47
Wrong calculation of nut in the kOmegaSST turbulence model FelixL OpenFOAM Bugs 27 March 27, 2012 09:02
Low Reynolds k-epsilon model YJZ ANSYS 1 August 20, 2010 13:57
KOmega Turbulence model from wwwopenFOAMWikinet philippose OpenFOAM Running, Solving & CFD 30 August 4, 2010 10:26

All times are GMT -4. The time now is 07:25.