SIMPLE gave this error message: "Illegal master cellID 0. Should be 0..0"
 

Hi!

I do a cylinder with gambit and then convert to OpenFOAM, but then when I run my case (tutorial buoyantBoussinesSimpleFoam) the program gives me an error:

--> FOAM FATAL IO ERROR:
Illegal master cellID 0. Should be 0..0

file: /home/xr12/Anna/buoyantBoussinesqSimpleFoam/gambit_rectangle/system/fvSolution::SIMPLE from line 39 to line 50.

From function void Foam::setRefCell
(
const volScalarField&,
const volScalarField&,
const dictionary&,
label& scalar&,
bool
)
in file cfdTools/general/findRefCell/findRefCell.C at line 66.

FOAM exiting


I understand that in file fvSolution, the part of SIMPLE i have to change the 0, but I don't know what i have to put... someone can help me??

my fvSolution:

solvers
{
p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 0.01;
}

"(U|T|k|epsilon|R)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;

residualControl
{
p_rgh 1e-2;
U 1e-4;
T 1e-2;

// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}

relaxationFactors
{
fields
{
p_rgh 0.7;
}
equations
{
U 0.3;
T 0.5;
"(k|epsilon|R)" 0.7;
}
}

Is not really a description of the flow case. It is a closed case. It is an open case? what kind of boundary conditions are you using. In order to get help in this forum you should go by the rules http://openfoamwiki.net/index.php/Ho..._Message_Board

Therefore, questions:

1. what kind of case is it?
2. what kind of boundary conditions are there?

and please consider using the code tags in order to insert longer parts of a property or case file, it reads much easier that way.

Dear libia
look at simple carefully.
you omit one } at the end of simple class:)

Nope, there is no "}" missing. It comes after residual control. But good point, this is what I meant with setting it into code blocks :)

only what i'm doing is the tutorial "buoyantBoussinesqSimpleFoam" but with different geometry.

First I change this geometry with GMSH and the program runs, but then I make the geometry with gambit and after convert to openfoam. But when I run the case appears this error, and I don't understand why, because changing the geometry with gmsh the program runs...

Does your geometry have an opening or is it also closed like the two tutorials? Try changing the number 0 for refCell to something different, since it will set a reference cell.

What you should also do is to expect your mesh with paraview (running paraFoam) and run checkMesh. See if there is something wrong. In paraview you can also select single cells and get their ID, which you can use to set for refCell instead of 0.

my geometry is closed like the tutorial. I try to change the number to 0.1 for example and gives me this error:

--> FOAM FATAL IO ERROR:
wrong token type - expected int, found on line 40 the doubleScalar 0.1

file: /home/xr12/Anna/buoyantBoussinesqSimpleFoam/gambit_rectangle/system/fvSolution::SIMPLE::pRefCell at line 40.

From function operator>>(Istream&, int&)
in file primitives/ints/int/intIO.C at line 68.

FOAM exiting

Then I try to put 1, and gives me the same error that when I put 0.

what are the outcomes from checkMesh? does the mesh looks good in paraview?

select a cell from paraview which you think is a good idea to set the reference pressure in?

I suspect that something is wrong with your mesh though, since you said that it works with gmsh.

when I star parafoam gives me another error:

--> FOAM Warning :
From function polyMesh(const IOobject&)
in file meshes/polyMesh/polyMesh.C at line 307
no points in mesh
--> FOAM Warning :
From function polyMesh(const IOobject&)
in file meshes/polyMesh/polyMesh.C at line 312
no cells in mesh


I'm new with gambit and what I do is to make the geometry with gambit, export as .msh and then copy the file into my case and convert to openfoam (gambitToFoam). What is what I'm doing wrong?

I would say you don't actually have a mesh. Since it tells you that there are no cells nor any points. what does checkMesh give you?

sorry but... how do you do the checkmesh??

sorry if I ask a lot, but I'm new with the programs and it's difficult to find some things.

gives me a longer error...

--> FOAM Serious Error :
From function bool zone::checkDefinition(const label maxSize, const bool report) const
in file meshes/polyMesh/zones/zone/zone.C at line 211
Zone internal20 contains invalid index label 19
Valid index labels are 0..-1
--> FOAM Serious Error :
From function bool zone::checkDefinition(const label maxSize, const bool report) const
in file meshes/polyMesh/zones/zone/zone.C at line 211
Zone internal20 contains invalid index label 20
Valid index labels are 0..-1
--> FOAM Serious Error :
From function bool zone::checkDefinition(const label maxSize, const bool report) const
in file meshes/polyMesh/zones/zone/zone.C at line 211
Zone internal20 contains invalid index label 21
Valid index labels are 0..-1
--> FOAM Serious Error :
From function bool zone::checkDefinition(const label maxSize, const bool report) const
in file meshes/polyMesh/zones/zone/zone.C at line 211
Zone internal20 contains invalid index label 22
Valid index labels are 0..-1
--> FOAM Serious Error :
From function bool zone::checkDefinition(const label maxSize, const bool report) const
in file meshes/polyMesh/zones/zone/zone.C at line 211
Zone internal20 contains invalid index label 23
Valid index labels are 0..-1

.
.
.
.
--> FOAM FATAL ERROR:
Too many errors

From function messageStream::operator OSstream&()
in file lnInclude/messageStream.C at line 200.

This just means that your mesh is not properly converted, or has not been properly set up from the beginning. I have never used gambit or gmsh, therefore I won't be able to help you with setting this up. I guess the best thing is to start a thread with a proper title and a detailed description of what you want to do, otherwise ppl might not look at the thread.

One tip, when posting error messages and terminal output, try to use the code block feature in the post (this is the little hash sign at the top of the editor when posting), this will get rid of smilies and will create a frame around long error messages.

Good luck with the meshing.

ok thanks for your help! :)

Hi

I am getting the same error for my mesh. please let me know how if you could resolve this error.

Thank you.

Saurabh

For anyone who may meet this problem, please check you 'constant/polymesh/owner' file or 'porcessorX/constant/polymesh/owner' (X should be every porcessors). This ERROR is caused by the number of cell is 0, due to the method of decomposePar.