[david]

Sorry, I gave you the wrong link. But you have found what I meant :-)

[aganesan]

Hello,
Actually I am trying to run a Fluent Generated Mesh case using MRFSimpleFoam (my application is impeller in a tank with baffles very similar to mixerVessel2D but the outer cylinder's dimension is bigger than the inner one). But when I try to compile it, I am getting the following errors (very similar to the one that has been reported in this discussion forum):
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : MRFSimpleFoam ./ impeller
Date : Jun 12 2008
Time : 15:43:01
Host : linux-02
PID : 18936
Root : /home/ang/OpenFOAM/ang-1.4.1/run/tutorials/MRFSimpleFoam/
Case : impeller
Nprocs : 1
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model laminar


--> FOAM FATAL ERROR : cannot find MRF patch rotor

From function Foam::MRFZone::MRFZone(const fvMesh& , const dictionary&)
in file cfdTools/general/MRF/MRFZone.C at line 86.

FOAM exiting

Actually I tried doing cellSet (successful) and setstoZones (unsuccessful, it basically complains of the some missing file called wallcells...).
It would be great if I can get help to run this case successfully.
Thanks.

Anand

[danielle]

Hi everybody!
I will validate my case with MRFsimpleFoam,
but I get a strange. I have one nonOrthogonal element...
I added the gravity force in the solver but the result still very bad. The same mesh work good with Fluent6.3.
you find my case here and thanks for your help.

[danielle]

Hi everybody!
I will validate my case with MRFsimpleFoam,
but I get a strange. I have one nonOrthogonal element...
I added the gravity force in the solver but the result still very bad. The same mesh work good with Fluent6.3.
you find my case here and thanks for your help.

[danielle]

Hi everybody!
I will validate my case with MRFsimpleFoam,
but I get a strange. I have one nonOrthogonal element...
I added the gravity force in the solver but the result still very bad. The same mesh work good with Fluent6.3.
you find my case here and thanks for your help.

[danielle]

You can download this case from:
http://www.megaupload.com/?d=LO8OJ5EK
thanks

[dmoroian]

Hi Danielle,
I've downloaded your case and the only thing that I can object is that although you ask for kEpsilon turbulence model, you set turbulence off. In my opinion, this is why you get an asymmetric velocity distribution in your solution. If you set turbulence on, then this is how the solution looks after 400 iterations:


I hope this is helpful,
Dragos

[danielle]

Thanks Dragos for your interest,
I already tried with turbulence on, but it diverge after 9 iterations !

have you change other things?
----------------------
Time = 9

smoothSolver: Solving for Ux, Initial residual = 0.414047, Final residual = 0.0125807, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.774418, Final residual = 0.00930591, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.349594, Final residual = 0.00896896, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 31.5247, No Iterations 1000
GAMG: Solving for p, Initial residual = 0.60274, Final residual = 0.0233148, No Iterations 2
GAMG: Solving for p, Initial residual = 0.999192, Final residual = 0.000139628, No Iterations 1
GAMG: Solving for p, Initial residual = 0.918478, Final residual = 0.0121946, No Iterations 1
GAMG: Solving for p, Initial residual = 0.792618, Final residual = 0.00782516, No Iterations 1
time step continuity errors : sum local = 3.82222e+89, global = 1.26894e+74, cumulative = 1.26894e+74
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 4.84359e+30, No Iterations 1000
smoothSolver: Solving for k, Initial residual = 1, Final residual = 1.41188e-16, No Iterations 2
bounding k, min: -2.89079e+55 max: 2.92153e+127 average: 3.10906e+122
ExecutionTime = 103.27 s ClockTime = 180 s

[dmoroian]

Hello again Danielle,
It is very strange that it diverges so fast. It still runs on my computer even after 2800 iterations:

I'm using OpenFOAM-1.4.1-dev. What are you using?
Maybe after I finish all the 5000 iterations I will put it on the wiki as a nice example... if you agree.

Dragos

[dmoroian]

...and the pressure after 5000 iterations:

[danielle]

Hello again Dragos,
In the first it's a good idea to put this case on the wiki as a example :-)

I'm using OpenFOAM-1.4.1 and OpenFOAM-1.5 download from the http://www.opencfd.co.uk/openfoam/download.html
what do you mean by -dev ??
is it another version ?
how can get it ?
NB: I like the stream line !! I can't do it.

[danielle]

Now I try to recompile the version 1.4.1-dev. It cost a lot of time ...
Can you tel me what's the difference for the MRF implementation between this version and 1.5. It's very strange that my case diverge with 1.4.1 and 1.5 but not with 1.4.1-dev as mentionned above by Dragos???
Thanks

[dmoroian]

Hello Danielle,
I'm not aware of any difference in the MRFSimpleFoam solver, between those versions. However, I did some small modifications to your setup files. One for instance is switching on the turbulence in constant/turbublenceProperties. Some other are in system/fvSolution. All the differences in the fvSolution file come from the fact that in OpenFOAM-1.4.1-dev there are no implicit defined values (according to Hrv), so you have to specify all the necessary options. This is why I got a running case and you didn't.
I tried the 1.4.1 version with your case, but having this fvSolution instead:
fvSolution
It worked at least the first 18 iterations without any problems.

I hope this is helpful,
Dragos

[danielle]

Hi Dragos,
with the same modification it run only the 4 iterations and diverge. if I mahke turbulence off it work but the result is very bad !!
I hope to resolve this problem!
I have not succeded to compile 1.4.1-dev
---------
version 1.4.1
Starting time loop

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0219844, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0224696, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0211405, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0468348, No Iterations 5
GAMG: Solving for p, Initial residual = 0.673131, Final residual = 0.0102878, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0760988, Final residual = 0.00309761, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0186959, Final residual = 0.00069648, No Iterations 3
GAMG: Solving for p, Initial residual = 0.00616283, Final residual = 0.000306048, No Iterations 2
time step continuity errors : sum local = 0.00134537, global = 3.76801e-18, cumulative = 3.76801e-18
smoothSolver: Solving for epsilon, Initial residual = 0.050018, Final residual = 0.00229752, No Iterations 2
bounding epsilon, min: -31.614 max: 11145.1 average: 39.4265
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0926464, No Iterations 2
ExecutionTime = 4.83 s ClockTime = 5 s

Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.554542, Final residual = 0.0478937, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.569828, Final residual = 0.0495053, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.57678, Final residual = 0.0482794, No Iterations 2
GAMG: Solving for p, Initial residual = 0.851186, Final residual = 0.0365313, No Iterations 2
GAMG: Solving for p, Initial residual = 0.633769, Final residual = 0.0229999, No Iterations 2
GAMG: Solving for p, Initial residual = 0.359745, Final residual = 0.0177159, No Iterations 2
GAMG: Solving for p, Initial residual = 0.106486, Final residual = 0.00306774, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0431512, Final residual = 0.0021308, No Iterations 2
time step continuity errors : sum local = 0.0322882, global = -1.03471e-16, cumulative = -9.97035e-17
smoothSolver: Solving for epsilon, Initial residual = 0.0161207, Final residual = 0.000744167, No Iterations 2
bounding epsilon, min: -6.24393 max: 21495.4 average: 76.4378
smoothSolver: Solving for k, Initial residual = 0.91095, Final residual = 0.0454208, No Iterations 2
ExecutionTime = 7.96 s ClockTime = 9 s

Time = 3

smoothSolver: Solving for Ux, Initial residual = 0.251846, Final residual = 0.0204955, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.250968, Final residual = 0.0204949, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.123233, Final residual = 0.00961137, No Iterations 2
GAMG: Solving for p, Initial residual = 0.894017, Final residual = 0.0364758, No Iterations 3
GAMG: Solving for p, Initial residual = 0.385103, Final residual = 0.0130242, No Iterations 2
GAMG: Solving for p, Initial residual = 0.105777, Final residual = 0.00524001, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0356128, Final residual = 0.00101145, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0102751, Final residual = 0.000496442, No Iterations 3
time step continuity errors : sum local = 0.0179467, global = 3.56831e-17, cumulative = -6.40204e-17
smoothSolver: Solving for epsilon, Initial residual = 0.0245249, Final residual = 0.00132025, No Iterations 2
bounding epsilon, min: -7.51668 max: 54431.7 average: 136.557
smoothSolver: Solving for k, Initial residual = 0.0502388, Final residual = 0.00403868, No Iterations 2
ExecutionTime = 11.28 s ClockTime = 12 s

Time = 4

smoothSolver: Solving for Ux, Initial residual = 0.136687, Final residual = 0.00884054, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.139049, Final residual = 0.00936171, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0763141, Final residual = 0.0049499, No Iterations 2
GAMG: Solving for p, Initial residual = 0.659499, Final residual = 0.0324916, No Iterations 3
GAMG: Solving for p, Initial residual = 0.295368, Final residual = 0.0110679, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0693061, Final residual = 0.00207471, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0215861, Final residual = 0.000810994, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0058221, Final residual = 0.000244605, No Iterations 3
time step continuity errors : sum local = 0.00615604, global = 7.36455e-17, cumulative = 9.62512e-18
smoothSolver: Solving for epsilon, Initial residual = 0.00728371, Final residual = 0.000434512, No Iterations 2
bounding epsilon, min: -15.4742 max: 71860.5 average: 211.005
smoothSolver: Solving for k, Initial residual = 0.0141166, Final residual = 0.000618886, No Iterations 2
ExecutionTime = 14.61 s ClockTime = 16 s

Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.0613413, Final residual = 4.07658e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.368884, Final residual = 0.0109739, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.40229, Final residual = 0.0119589, No Iterations 2
GAMG: Solving for p, Initial residual = 0.502816, Final residual = 0.0193722, No Iterations 3
GAMG: Solving for p, Initial residual = 0.994737, Final residual = 0.0316806, No Iterations 3
GAMG: Solving for p, Initial residual = 0.231567, Final residual = 0.00704651, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0409088, Final residual = 0.00166323, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0155019, Final residual = 0.000753901, No Iterations 2
time step continuity errors : sum local = 5.40986e+07, global = -1.85347e-08, cumulative = -1.85347e-08
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0118739, No Iterations 4
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00940916, No Iterations 4
ExecutionTime = 17.83 s ClockTime = 19 s

[dmoroian]

Hmmmm, you're right,
Try with this:
fvSolution

[danielle]

Thanks for your help
It worke :-)
now 50 iter
51 iter
52 ...
I will see after that the result and tell you if it's the same thing with Fluent...
many thanks for your precious time.

[danielle]

Hi Dragos,
I would like to thank you for your help.
It gives a good results.
Now I can try with my solver twoPhaseEulerFoamMRF.
good day