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

MRFSimpleFoam Tutorial

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 29, 2008, 01:25
Default Hi all, I try to get the MR
  #1
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Hi all,

I try to get the MRFSimpleFoam tutorial running. Compiling runs with both 1.4.1 and 1.5 without problems on my Core2Duo (using 64 Bit). However, solving fails:

1.4.1:
MRFSimpleFoam . mixerVessel2D
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / 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 . mixerVessel2D
Date : Aug 29 2008
Time : 07:05:21
Host : Basti-Notebook.lan
PID : 9432
Root : /home/basti/OpenFOAM/basti-1.4.1/tutorials
Case : mixerVessel2D
Nprocs : 1
Create time

Create mesh for time = 0

Reading field p

new cannot satisfy memory request.
This does not necessarily mean you have run out of virtual memory.
It could be due to a stack violation caused by e.g. bad use of pointers or an out of date shared library

1.5
MRFSimpleFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : MRFSimpleFoam
Date : Aug 29 2008
Time : 07:16:41
Host : Basti-Notebook.lan
PID : 11930
Case : /home/basti/OpenFOAM/basti-1.5/tutorials/mixerVessel2D
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 RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
}


Starting time loop

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0404075, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0403676, No Iterations 2
#0 Foam::error::printStack(Foam:stream&) in "/home/basti/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/home/basti/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::MRFZone::relativeFlux(Foam::GeometricField<d ouble,>&) const in "/home/basti/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::MRFZones::relativeFlux(Foam::GeometricField< double,>&) const in "/home/basti/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 main in "/home/basti/OpenFOAM/basti-1.5/applications/bin/linux64GccDPOpt/MRFSimpleFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 Foam::regIOobject::readIfModified() in "/home/basti/OpenFOAM/basti-1.5/applications/bin/linux64GccDPOpt/MRFSimpleFoam"

Any ideas?

BTW can MRFSimpleFOam be run in parallel? I have read somewhere it cannot, but it has a -parallel flag?

Regards BastiL
bastil is offline   Reply With Quote

Old   August 29, 2008, 03:00
Default HI BastiL try to check your
  #2
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
HI BastiL

try to check your boundary setting,itseems like the boundary value cause the iteration go divagenced

it can run in parallel way,look at the userguide of OF in appendix A

wayne
waynezw0618 is offline   Reply With Quote

Old   August 29, 2008, 04:37
Default Works fine for me in 1.5, also
  #3
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Works fine for me in 1.5, also under valgrind. Your 'new cannot satisfy ..' error is usually caused by inconsistent libraries. Make sure you compile everything with the same compiler.
mattijs is offline   Reply With Quote

Old   August 29, 2008, 11:39
Default I see: It seems like 1.5 start
  #4
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
I see: It seems like 1.5 starts solving and diverges in Iteration 1. I will check BC; even though I did not change them...
1.4.1 does not seem to start at all. I uses the compiler supplied with 1.5 and 1.4.1, respectively.
bastil is offline   Reply With Quote

Old   August 29, 2008, 15:42
Default I reviewed BC and verything lo
  #5
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
I reviewed BC and verything looks fine for me. Maybe I could upload the case here and somebody could try to run it.
How do I make sure to use only one compiler. I usually source bashrc file from appropriate version and run wmake afterwards... That is what I did for both 1.4.1 and 1.5
bastil is offline   Reply With Quote

Old   August 29, 2008, 15:43
Default I reviewed BC and everything l
  #6
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
I reviewed BC and everything looks fine for me. Maybe I could upload the case here and somebody could try to run it.
How do I make sure to use only one compiler. I usually source bashrc file from appropriate version and run wmake afterwards... That is what I did for both 1.4.1 and 1.5
bastil is offline   Reply With Quote

Old   August 30, 2008, 09:32
Default Thanks guys got it runnig. Obv
  #7
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Thanks guys got it runnig. Obviously something with compiling was wrong, don't know what exactly.

I can not find something about parallel in UG Appendix A, its about FoamX. Is MRFSimpleFoam fully parallel or does the MRF Fluid have to be on one partition? I read about developments for turbomachinery? is MRFSimpleFoam depreciated?

By the way I compiled PV3FoamReader but I do not understand what it is good for?
bastil is offline   Reply With Quote

Old   August 31, 2008, 05:21
Default Hello BastiL, Yes MRFSimpleFo
  #8
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hello BastiL,
Yes MRFSimpleFoam is fully parallel, and not at all obsolete!

Dragos
dmoroian is offline   Reply With Quote

Old   August 31, 2008, 08:39
Default Hi Dragos do you have tried t
  #9
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
Hi Dragos
do you have tried the simpleSRFFoam in OF1.5?i think it is similar to the MRFsimpleFoam in consumption of calculation.but when i tried a similar case of centrifugal pump impeller of same mesh and boundary condition.simpleSRFFoam in OF1.5 is too slow--only 10 steps for 10 hours.but hunderds steps for MRFSimpleFoam in OF1.4.that is why? how can i do?

wayne
waynezw0618 is offline   Reply With Quote

Old   August 31, 2008, 17:33
Default Thanks once more. I could run
  #10
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Thanks once more. I could run in parallel, too.However, everything runs in DP but not in SP. MAybe I did something wrong in compiling it in SP.

Regards.
bastil is offline   Reply With Quote

Old   September 1, 2008, 04:17
Default Hello Wayne, Unfortunately I
  #11
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hello Wayne,
Unfortunately I did not switch to 1.5 version, so I can only guess that some of the thresholds for the linear solvers are changed and too tight.

Dragos
dmoroian is offline   Reply With Quote

Old   September 22, 2008, 21:14
Default Hello every one, As of Differ
  #12
New Member
 
MR Amiralaei
Join Date: Mar 2009
Posts: 20
Rep Power: 8
reza is on a distinguished road
Hello every one,
As of Different ways to treat rotating geometries by Olivier Petit, the MRFSimpleFoam is on OpenCFD versions and as I get error in runnnig the solver in dev version I like to know if there is any body who is runnig this solver on dev version as well.
My second question is that if there is any adjustment except .bashrc when you want to run two versions of OpenFOAM on the same machine.
Thanks
reza is offline   Reply With Quote

Old   December 29, 2008, 06:41
Default hi guys i have tried to com
  #13
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
hi guys

i have tried to compile the MRFSimpleFoam.C
and run the mixerVessel2D case but i got this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5 |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : MRFSimpleFoam
Date : Dec 29 2008
Time : 11:39:15
Host : antonio-laptop
PID : 19346
Case : /home/cfduser/OpenFOAM/cfduser-1.5/tutorials/MRFSimpleFoam/mixerVessel2D
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 RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
}



cannot find MRF faceZone rotor

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

FOAM exiting


any suggestions?
antonio_ing is offline   Reply With Quote

Old   December 29, 2008, 06:59
Default strictly speaking if I have a
  #14
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
strictly speaking if I have a case that works with simpleFoam, What file should I add/modify to have a rotation?

thanks in advance
antonio_ing is offline   Reply With Quote

Old   January 7, 2009, 07:23
Default try simpleSRFFoam.you will fin
  #15
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
try simpleSRFFoam.you will find more
waynezw0618 is offline   Reply With Quote

Old   January 7, 2009, 12:28
Default I have been able to run the mi
  #16
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
I have been able to run the mixerVessel case but still i do not understand where to set the region which is rotating. Does anyone have a very simple case?
antonio_ing is offline   Reply With Quote

Old   January 8, 2009, 14:15
Default Hi Antonio, Have a look at th
  #17
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hi Antonio,
Have a look at this tutorial: MRFSimpleFoam Tutorial

I hope this is helpful,
Dragos
dmoroian is offline   Reply With Quote

Old   April 9, 2009, 04:06
Default
  #18
Member
 
Mandar
Join Date: Mar 2009
Posts: 39
Rep Power: 8
mixer is on a distinguished road
Hello Dragos , Wayne and people who have solved this problem

I need help from you guys. I am facing similar problem: "cannot find MRF faceZone rotor". I tried but i am unsuccessful.

I have done the following steps:

1. fluentMeshToFoam mixer.msh -writeSets -writeZones

It seems to have been able to generate the cellZones and FaceZones, and shows three domains that are there:

TANK domain, IMPELLER domain and a BOTTOM domain (where inlet is).

IMPELLER domain is the rotating domain.

2. So, i felt that i do not need to do cellSet, setsToZones, faceSet, setsToZones routine. Am I right?

3. then, I corrected the cellSetDict file, see below, is this right?

name rotor;
action new;
topoSetSources
(
// Cells in cell zone
zoneToCell
{
name IMPELLER; // name of cellZone
}


as IMPELLER is the rotating domain.

4. Now, on running ./Allrun from the problem directory

i get the error, "cannot find MRF faceZone rotor".

so, i did

cellset and
setsToZones -noFlipMap

It worked upto this.


5. But error cropped up when i did faceSet,

Reading faceSetDict
unexpected class name cellSet expected faceSet
while reading object rotor
file: /usr6/tab01c/OpenFOAM/tab01c-1.5/run/tutorials/MRFSimpleFoam/mixermandar/constant/polyMesh/sets/rotor at line 15.
From function regIOobject::readStream(const word&)
in file db/regIOobject/regIOobjectRead.C at line 114.
FOAM exiting


Can anyone please guide me on this. I have three domains, impeller domian being the rotating domain. Is it problem with ./Allrun file?. How i should go about solving this. thanks for your help.

Please let me know your suggestions, i wil ltry implementing this. Enjoy Easter by the way. Thanks.
mixer is offline   Reply With Quote

Old   April 9, 2009, 09:33
Default Newer MRFsimpleFoam?
  #19
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Hello Dragos,
I have a short question about your MRFsimpleFoam tutorial. Since I want to implement a similar case (a ducted fan), I started to run your case to learn how to set up mine. However, after no more than 15 iterations, I had a Floating Point exception error, due to the fast increase in epsilon and k values.
Since to make it run I had to copy an old RASproperties file into /costant , I guess that there is a newer version of the solver. Indeed trasportProperties file has the same function of RASProperties but with a slightly different formulation. I tried to install the solver files that I got when downloaded your files, but it didn't work. So I think it is due to the new OF 1.5.dev version. Am I right? If so, where can I get it? Thanks a lot!

Maddalena
maddalena is offline   Reply With Quote

Old   April 9, 2009, 09:40
Default changed MRFZones?
  #20
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 12
maddalena is on a distinguished road
Quote:
Originally Posted by mixer View Post
2. So, i felt that i do not need to do cellSet, setsToZones, faceSet, setsToZones routine. Am I right?
Yes, you are right. You do not need to do any cellSet, setsToZones, faceSet, setsToZones routines, neither ./Allrun, that does the same. After doing fluentMeshToFoam mixer.msh -writeSets -writeZones, cells and zones are ok. Btw, did you change MRFZones file to match your rotating patch and domain?
Maddalena
maddalena 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
Equations in the MRFsimpleFOAM waynezw0618 OpenFOAM Running, Solving & CFD 5 May 7, 2015 04:43
Convergence with MRFSimpleFoam grugg OpenFOAM Running, Solving & CFD 7 March 28, 2014 05:56
CellZons and MRFSimpleFoam vinz OpenFOAM Running, Solving & CFD 18 July 30, 2012 08:11
MRFSimpleFoam PropellerMixer tino_boelke OpenFOAM Running, Solving & CFD 0 December 17, 2008 11:25
MRFSimpleFoam xdanielx OpenFOAM Running, Solving & CFD 0 December 17, 2008 02:28


All times are GMT -4. The time now is 14:22.