[Nicacio]

Hi, I'm Antonio.
I am pretty new with OpenFoam and I need to simulate a case. After reading many tutorials and problems of the rest of the forum, I can not solve mine.
The mesh is made in Salome ​​with NetGen. I defined the parameters through Openfoam Wizard and run OpenFOAM in the terminal with to many errors that I have not yet found how to fix them. Can anyone help me out, please?
Here is the code:

nicacio@CAE-Linux:~/Desktop/External_flow/Sphere$ ./ParallelRun.sh/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : decomposePar
Date   : Jun 09 2014
Time   : 14:17:43
Host   : "CAE-Linux"
PID    : 5074
Case   : /home/nicacio/Desktop/External_flow/Sphere
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh

Calculating distribution of cells
Selecting decompositionMethod scotch

Finished decomposition in 3.52 s

Calculating original mesh data

Distributing cells to processors

Distributing faces to processors

Distributing points to processors

Constructing processor meshes

Processor 0
    Number of cells = 77738
    Number of faces shared with processor 2 = 425
    Number of processor patches = 1
    Number of processor faces = 425
    Number of boundary faces = 11207

Processor 1
    Number of cells = 85490
    Number of faces shared with processor 2 = 958
    Number of faces shared with processor 3 = 1038
    Number of processor patches = 2
    Number of processor faces = 1996
    Number of boundary faces = 452

Processor 2
    Number of cells = 77694
    Number of faces shared with processor 0 = 425
    Number of faces shared with processor 1 = 958
    Number of faces shared with processor 3 = 66
    Number of processor patches = 3
    Number of processor faces = 1449
    Number of boundary faces = 4737

Processor 3
    Number of cells = 85747
    Number of faces shared with processor 1 = 1038
    Number of faces shared with processor 2 = 66
    Number of processor patches = 2
    Number of processor faces = 1104
    Number of boundary faces = 0

Number of processor faces = 2487
Max number of cells = 85747 (4.99558% above average 81667.2)
Max number of processor patches = 3 (50% above average 2)
Max number of faces between processors = 1996 (60.5147% above average 1243.5)

Time = 0


--> FOAM FATAL IO ERROR: 
keyword defaultFaces is undefined in dictionary "/home/nicacio/Desktop/External_flow/Sphere/0/p::boundaryField"

file: /home/nicacio/Desktop/External_flow/Sphere/0/p::boundaryField from line 25 to line 34.

    From function dictionary::subDict(const word& keyword) const
    in file db/dictionary/dictionary.C at line 461.

FOAM exiting

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : simpleFoam -parallel
Date   : Jun 09 2014
Time   : 14:17:52
Host   : "CAE-Linux"
PID    : 5076
Case   : /home/nicacio/Desktop/External_flow/Sphere
nProcs : 4
Slaves : 
3
(
"CAE-Linux.5077"
"CAE-Linux.5078"
"CAE-Linux.5079"
)

Pstream initialized with:
    floatTransfer     : 0
    nProcsSimpleSum   : 0
    commsType         : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] cannot find file
[0] 
[0] file: /home/nicacio/Desktop/External_flow/Sphere/processor0/0/p at line 0.
[0] 
[0]     From function regIOobject::readStream()
[0]     in file db/regIOobject/regIOobjectRead.C at line 73.
[0] 
FOAM parallel run exiting
[0] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[2] 
[2] 
[2] --> FOAM FATAL IO ERROR: 
[2] cannot find file
[2] 
[2] file: /home/nicacio/Desktop/External_flow/Sphere/processor2/0/p at line 0.
[2] 
[2]     From function regIOobject::readStream()
[2]     in file db/regIOobject/regIOobjectRead.C at line 73[3] 
[3] 
[3] --> FOAM FATAL IO ERROR: 
[3] cannot find file
[3] 
[3] file: /home/nicacio/Desktop/External_flow/Sphere/processor3/0/p at line 0.
[3] 
[3]     From function regIOobject::readStream()
[3]     in file db/regIOobject/regIOobjectRead.C at line 73.
[3] 
FOAM parallel run exiting
[3] 
.
[2] 
FOAM parallel run exiting
[2] 
[1] 
[1] 
[1] --> FOAM FATAL IO ERROR: 
[1] cannot find file
[1] 
[1] file: /home/nicacio/Desktop/External_flow/Sphere/processor1/0/p at line 0.
[1] 
[1]     From function regIOobject::readStream()
[1]     in file db/regIOobject/regIOobjectRead.C at line 73.
[1] 
FOAM parallel run exiting
[1] 
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 5076 on
node CAE-Linux exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------
[CAE-Linux:05075] 3 more processes have sent help message help-mpi-api.txt / mpi-abort
[CAE-Linux:05075] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : reconstructPar
Date   : Jun 09 2014
Time   : 14:17:54
Host   : "CAE-Linux"
PID    : 5080
Case   : /home/nicacio/Desktop/External_flow/Sphere
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time



--> FOAM FATAL ERROR: 
No times selected

    From function reconstructPar
    in file reconstructPar.C at line 139.

FOAM exiting

Thank you so much

[itchy]

Hi Antonio,

pherhaps there is going smth wrong when you convert the mesh into OpenFoam. Can you please post your /0/p file and /constant/polyMesh/boundary file. In my opinion you have not declare all the faces in the mesh as boundary faces. So if you convert the mesh into OF you get additional boundaries called "default faces". This should be the problem.

I have never used this wizard or Salome. In my opinion the best way is to work without a Wizard (only terminal) and gridgeneration with GMSH (GMSH is free). Or any commercial gridgenerator instead of GMSH. That can save you so much time.

kind regards
Florian

[Nicacio]

Yes, you are right. I get always the "default faces error". I'm trying to simulate a sphere inside a domain, and I have declared as boundary inlet, outlet, wall. I guess openfoam adds the sphere volume as a boundary as well and this is the error. I have tried in different ways but always the same... should I declare the sphere in salome in a mesh volume group separately or like faces??

Or do you have any tutorial for external flow (for dummies) from the geometry?? or let me know where I can find a good one because the tutorials I have found are more "advanced"

[itchy]

Hi Antonio,

At the moment I have no tutorial for that. There are good tutorials at Chalmers university:
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2013/

In my opinion it is a problem in Salome with the boundarys or by converting it into OpenFOAM.

First do this:

type in your terminal in the case-directory (where the /system, /constant, /0 folders are):
checkMesh

if checkMesh says everthing is right, then go into paraView (type paraFoam in your case-directory):
there you check where the defaultFaces are. Regard that no fields are selected in paraView. If you select field p then paraView will crash.

Please post the checkMesh output if smth. is wrong.
And please post /0/p and /constant/polyMesh/boundary
without these two files nobody can help you with your problem.

kind regards
Florian

How many cellzones do you have?? You can find the answer in your /constant/polyMesh/cellZones
If you have two then pherhaps you have a mesh for the fluid and for the sphere. That could be wrong.

[Nicacio]

Thank you very very much Itchy.
The mesh is Ok.

/*--------------------------------*- C++ -*----------------------------------*| =========                 |                                                 |
| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \    /   O peration     | Version:  2.0.1                                 |
|   \  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
Inlet
{
       type            zeroGradient;
}
Outlet
{
        type            fixedValue; 
        value           uniform 0;
}
Wall
{
        type            zeroGradient;
}
}

And here are my /0/p

/*--------------------------------*- C++ -*----------------------------------*| =========                 |                                                 |
| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \    /   O peration     | Version:  2.0.1                                 |
|   \  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
Inlet
{
       type            zeroGradient;
}
Outlet
{
        type            fixedValue; 
        value           uniform 0;
}
Wall
{
        type            zeroGradient;
}
}

and the /constant/polyMesh/boundary

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
    Inlet
    {
        type            patch;
        nFaces          18;
        startFace       645140;
    }
    Outlet
    {
        type            patch;
        nFaces          18;
        startFace       645158;
    }
    Wall
    {
        type            wall;
        nFaces          416;
        startFace       645176;
    }
    defaultFaces
    {
        type            patch;
        nFaces          15944;
        startFace       645592;
    }
)

// ************************************************************************* //

[adambarfi]

you should specify some boundary conditions for 'defaultFaces'.

[itchy]

ok we have the problem.

It is in /constant/polyMesh/boundary

here you see your boundaryfaces of the mesh. You have Inlet, Outlet and Wall. These are ok. But what is about defaultFaces??
This comes from the converting of your mesh in OpenFOAM.
In my opinion, there are a few faces at the boundary you have not set in Salome. But for a first shot this isn't a problem. In the most cases, these faces are walls. So what you can do is:

add the following lines into p file after Wall{...}
...
defaultFaces
{
type fixedValue;
value uniform 0;
}
...

and the same for U :

...
defaultFaces
{
type fixedValue;
value uniform (0 0 0);
}
...

do you have any other files in /0 ?? If yes change it in the same way!!

Or these defaultFaces could be internal faces, than use zeroGradient for them!!

Or it is a 2d case, then use empty for the faces.

Than go to paraView and look where these faces are. If these faces are walls, everything is ok. Then you change the type of defaultFaces from patch into wall (/constant/polyMesh/boundary). If not we have to think about a proper type.

If this will not work you can delete defaultFaces from the /constant/polyMesh/boundary file. And then export the meshfile into some other format (for example foamToFluent) and then convert this file back into openFoam (fluentToFoam). That's a bad way but it should work

So in /constant/polyMesh/boundary you can find all the boundarys of your mesh, what openfoam has converted. These boundarys have to be equal to the boundary in your e.g. /0/U or /0/p
If you find boundarys there you have not set in Salome, then you have done something wrong in Salome.

I hope you can fix your problem.

kind regards
Florian

[Nicacio]

I got it!!! Thank you very much Itchy for your help!!
I figured it out where was the problem. I defined the sphere inside the domain like a boundary as well, so then, I had inlet, outlet, wall and default faces or sphere... I changed it, having just inlet, outlet and wall (including the sphere and the domain walls) and it worked... This is not going to happen again (newbie problems)

Thank you very much and sorry for all the inconveniences that may I have caused
Kind Regards
Antonio