CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   MRFSimpleFoam amp cyclic patches (https://www.cfd-online.com/Forums/openfoam-solving/58445-mrfsimplefoam-amp-cyclic-patches.html)

david June 12, 2008 07:30

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

aganesan June 12, 2008 07:40

Hello, Actually I am trying t
 
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 October 19, 2008 14:29

Hi everybody! I will validat
 
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 October 19, 2008 14:34

Hi everybody! I will validat
 
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 October 19, 2008 14:37

Hi everybody! I will validat
 
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 October 19, 2008 16:01

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

dmoroian October 20, 2008 08:56

Hi Danielle, I've downloaded
 
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:
http://www.cfd-online.com/OpenFOAM_D...ges/1/9614.jpg

I hope this is helpful,
Dragos

danielle October 20, 2008 09:23

Thanks Dragos for your intere
 
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 October 20, 2008 10:00

Hello again Danielle, It is v
 
Hello again Danielle,
It is very strange that it diverges so fast. It still runs on my computer even after 2800 iterations:
http://www.cfd-online.com/OpenFOAM_D...ges/1/9619.jpg
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 October 20, 2008 11:04

...and the pressure after 5000
 
...and the pressure after 5000 iterations:
http://www.cfd-online.com/OpenFOAM_D...ges/1/9623.jpg

danielle October 20, 2008 12:48

Hello again Dragos, In the fi
 
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 October 20, 2008 21:49

Now I try to recompile the ve
 
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 October 21, 2008 02:39

Hello Danielle, I'm not aware
 
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:
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif fvSolution
It worked at least the first 18 iterations without any problems.

I hope this is helpful,
Dragos

danielle October 21, 2008 09:31

Hi Dragos, with the same mod
 
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 http://www.cfd-online.com/OpenFOAM_D...lipart/sad.gif
---------
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 October 21, 2008 10:37

Hmmmm, you're right, Try with
 
Hmmmm, you're right,
Try with this:
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif fvSolution

danielle October 21, 2008 11:43

Thanks for your help It worke
 
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 October 21, 2008 21:55

Hi Dragos, I would like to th
 
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


All times are GMT -4. The time now is 09:47.