[bastil]

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

[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

[mattijs]

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]

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]

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]

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]

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]

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

Dragos

[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

[bastil]

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]

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]

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]

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]

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]

try simpleSRFFoam.you will find more

[antonio_ing]

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]

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

I hope this is helpful,
Dragos

[mixer]

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]

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]

Originally Posted by mixer

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