CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   MRFSimpleFoam Tutorial (http://www.cfd-online.com/Forums/openfoam-solving/58141-mrfsimplefoam-tutorial.html)

bastil August 29, 2008 01:25

Hi all, I try to get the MR
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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

waynezw0618 August 29, 2008 03:00

HI BastiL try to check your
 
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

mattijs August 29, 2008 04:37

Works fine for me in 1.5, also
 
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.

bastil August 29, 2008 11:39

I see: It seems like 1.5 start
 
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 August 29, 2008 15:42

I reviewed BC and verything lo
 
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 August 29, 2008 15:43

I reviewed BC and everything l
 
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 August 30, 2008 09:32

Thanks guys got it runnig. Obv
 
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?

dmoroian August 31, 2008 05:21

Hello BastiL, Yes MRFSimpleFo
 
Hello BastiL,
Yes MRFSimpleFoam is fully parallel, and not at all obsolete!

Dragos

waynezw0618 August 31, 2008 08:39

Hi Dragos do you have tried t
 
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

bastil August 31, 2008 17:33

Thanks once more. I could run
 
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.

dmoroian September 1, 2008 04:17

Hello Wayne, Unfortunately I
 
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

reza September 22, 2008 21:14

Hello every one, As of Differ
 
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

antonio_ing December 29, 2008 06:41

hi guys i have tried to com
 
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 December 29, 2008 06:59

strictly speaking if I have a
 
strictly speaking if I have a case that works with simpleFoam, What file should I add/modify to have a rotation?

thanks in advance

waynezw0618 January 7, 2009 07:23

try simpleSRFFoam.you will fin
 
try simpleSRFFoam.you will find more

antonio_ing January 7, 2009 12:28

I have been able to run the mi
 
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?

dmoroian January 8, 2009 14:15

Hi Antonio, Have a look at th
 
Hi Antonio,
Have a look at this tutorial: MRFSimpleFoam Tutorial

I hope this is helpful,
Dragos

mixer April 9, 2009 04:06

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.

maddalena April 9, 2009 09:33

Newer MRFsimpleFoam?
 
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 April 9, 2009 09:40

changed MRFZones?
 
Quote:

Originally Posted by mixer (Post 212431)
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


All times are GMT -4. The time now is 13:43.