
[Sponsors] 
June 12, 2008, 07:30 
Sorry, I gave you the wrong li

#21 
Member
David Hora
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63
Rep Power: 9 
Sorry, I gave you the wrong link. But you have found what I meant :)


June 12, 2008, 07:40 
Hello,
Actually I am trying t

#22 
New Member
Anand
Join Date: Mar 2009
Posts: 6
Rep Power: 9 
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 : linux02 PID : 18936 Root : /home/ang/OpenFOAM/ang1.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 

October 19, 2008, 14:29 
Hi everybody!
I will validat

#23 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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. 

October 19, 2008, 14:34 
Hi everybody!
I will validat

#24 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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. 

October 19, 2008, 14:37 
Hi everybody!
I will validat

#25 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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. 

October 19, 2008, 16:01 
You can download this case fro

#26 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 

October 20, 2008, 08:56 
Hi Danielle,
I've downloaded

#27 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
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 

October 20, 2008, 09:23 
Thanks Dragos for your intere

#28 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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.41188e16, No Iterations 2 bounding k, min: 2.89079e+55 max: 2.92153e+127 average: 3.10906e+122 ExecutionTime = 103.27 s ClockTime = 180 s 

October 20, 2008, 10:00 
Hello again Danielle,
It is v

#29 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
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 OpenFOAM1.4.1dev. 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 

October 20, 2008, 11:04 
...and the pressure after 5000

#30 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
...and the pressure after 5000 iterations:


October 20, 2008, 12:48 
Hello again Dragos,
In the fi

#31 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
Hello again Dragos,
In the first it's a good idea to put this case on the wiki as a example :) I'm using OpenFOAM1.4.1 and OpenFOAM1.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. 

October 20, 2008, 21:49 
Now I try to recompile the ve

#32 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
Now I try to recompile the version 1.4.1dev. 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.1dev as mentionned above by Dragos??? Thanks 

October 21, 2008, 02:39 
Hello Danielle,
I'm not aware

#33 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
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 OpenFOAM1.4.1dev 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 

October 21, 2008, 09:31 
Hi Dragos,
with the same mod

#34 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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.1dev  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.76801e18, cumulative = 3.76801e18 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.03471e16, cumulative = 9.97035e17 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.56831e17, cumulative = 6.40204e17 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.36455e17, cumulative = 9.62512e18 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.07658e06, 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.85347e08, cumulative = 1.85347e08 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 

October 21, 2008, 10:37 
Hmmmm, you're right,
Try with

#35 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 

October 21, 2008, 11:43 
Thanks for your help
It worke

#36 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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. 

October 21, 2008, 21:55 
Hi Dragos,
I would like to th

#37 
Member
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 9 
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 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Cyclic vs ggi vs directMapped Patches  jens_klostermann  OpenFOAM Running, Solving & CFD  33  May 3, 2013 02:45 
MRFSimpleFoam  xdanielx  OpenFOAM Running, Solving & CFD  0  December 17, 2008 02:28 
Application of cyclic patches  klaus  OpenFOAM Meshing Format & General Technical  9  January 18, 2008 05:28 
Cyclic boundary conditions for non conformed Patches  turnow  OpenFOAM Running, Solving & CFD  1  October 19, 2007 01:17 
View cyclic boundaries using MRFSimpleFoam  gabriel  OpenFOAM Paraview & paraFoam  0  July 27, 2007 08:42 