CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

MRFSimpleFoam amp cyclic patches

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 12, 2008, 07:30
Default Sorry, I gave you the wrong li
  #21
Member
 
David Hora
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63
Rep Power: 8
david is on a distinguished road
Sorry, I gave you the wrong link. But you have found what I meant :-)
david is offline   Reply With Quote

Old   June 12, 2008, 07:40
Default Hello, Actually I am trying t
  #22
New Member
 
Anand
Join Date: Mar 2009
Posts: 6
Rep Power: 8
aganesan is on a distinguished road
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
aganesan is offline   Reply With Quote

Old   October 19, 2008, 14:29
Default Hi everybody! I will validat
  #23
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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 is offline   Reply With Quote

Old   October 19, 2008, 14:34
Default Hi everybody! I will validat
  #24
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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 is offline   Reply With Quote

Old   October 19, 2008, 14:37
Default Hi everybody! I will validat
  #25
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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 is offline   Reply With Quote

Old   October 19, 2008, 16:01
Default You can download this case fro
  #26
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
You can download this case from:
http://www.megaupload.com/?d=LO8OJ5EK
thanks
danielle is offline   Reply With Quote

Old   October 20, 2008, 08:56
Default Hi Danielle, I've downloaded
  #27
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
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
dmoroian is offline   Reply With Quote

Old   October 20, 2008, 09:23
Default Thanks Dragos for your intere
  #28
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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
danielle is offline   Reply With Quote

Old   October 20, 2008, 10:00
Default Hello again Danielle, It is v
  #29
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
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 is offline   Reply With Quote

Old   October 20, 2008, 11:04
Default ...and the pressure after 5000
  #30
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
...and the pressure after 5000 iterations:

dmoroian is offline   Reply With Quote

Old   October 20, 2008, 12:48
Default Hello again Dragos, In the fi
  #31
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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 is offline   Reply With Quote

Old   October 20, 2008, 21:49
Default Now I try to recompile the ve
  #32
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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
danielle is offline   Reply With Quote

Old   October 21, 2008, 02:39
Default Hello Danielle, I'm not aware
  #33
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
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
dmoroian is offline   Reply With Quote

Old   October 21, 2008, 09:31
Default Hi Dragos, with the same mod
  #34
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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
danielle is offline   Reply With Quote

Old   October 21, 2008, 10:37
Default Hmmmm, you're right, Try with
  #35
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hmmmm, you're right,
Try with this:
fvSolution
dmoroian is offline   Reply With Quote

Old   October 21, 2008, 11:43
Default Thanks for your help It worke
  #36
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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 is offline   Reply With Quote

Old   October 21, 2008, 21:55
Default Hi Dragos, I would like to th
  #37
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 8
danielle is on a distinguished road
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
danielle is offline   Reply With Quote

Reply

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
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


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