[daniel.almeida]

Hello,

I am trying to create a mesh of a wedge with a small void part representing a nozzle. I have the blockMeshDict, and p and U at 0. I created the mesh, but when I try to run checkMesh I get the following message:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.2-9240f8b967db
Exec : checkMesh
Date : Jul 31 2015
Time : 14:16:43
Host : "daniel-N450-P-BE55P1"
PID : 6298
Case : /home/daniel/OpenFOAM/daniel-2.2.2/mit/nozzle2D
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 polyMesh for time = 0

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:olyMesh::calcDirections() const at ??:?
#4 Foam:olyMesh:olyMesh(Foam::IOobject const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception



I could open the mesh on paraview without loading the pressure and velocity fields, as it occurred in other error:

--> FOAM FATAL IO ERROR:

patch type 'patch' not constraint type 'wedge'
for patch front of field U in file "/home/daniel/OpenFOAM/daniel-2.2.2/mit/nozzle2D/0/U"


Here is the blockMeshDict:

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

convertToMeters 0.001;

vertices
(
(0 0 0) //vertex number 0
(19.3937 0 0) //vertex number 1
(300 0 0) //vertex number 2
(0 -0.138496 3.1720) //vertex number 3
(0 0.138496 3.1720) //vertex number 4
(19.3937 -0.138496 3.1720) //vertex number 5
(19.3937 0.138496 3.1720) //vertex number 6
(300 -0.138496 3.1720) //vertex number 7
(300 0.138496 3.1720) //vertex number 8
(0 -0.207748 4.7581) //vertex number 9
(0 0.207748 4.7581) //vertex number 10
(6.4646 -0.207748 4.7581) //vertex number 11
(6.4646 0.207748 4.7581) //vertex number 12
(19.3937 -0.207748 4.7581) //vertex number 13
(19.3937 0.207748 4.7581) //vertex number 14
(300 -0.207748 4.7581) //vertex number 15
(300 0.207748 4.7581) //vertex number 16
(0 -2.1810 49.952) //vertex number 17
(0 2.1810 49.952) //vertex number 18
(6.4646 -2.1810 49.952) //vertex number 19
(6.4646 2.1810 49.952) //vertex number 20
(19.3937 -2.1810 49.952) //vertex number 21
(19.3937 2.1810 49.952) //vertex number 22
(300 -2.1810 49.952) //vertex number 23
(300 2.1810 49.952) //vertex number 24
);

blocks
(
hex (0 1 1 0 3 5 6 4) (20 1 6) simpleGrading (1 1 1) //0
hex (1 2 2 1 5 7 8 6) (280 1 6) simpleGrading (1 1 1) //1
hex (5 7 8 6 13 15 16 14) (280 1 3) simpleGrading (1 1 1) //2
hex (9 11 12 10 17 19 20 18) (6 1 100) simpleGrading (1 1 1) //3
hex (11 13 14 12 19 21 22 20) (14 1 100) simpleGrading (1 1 1) //4
hex (13 15 16 14 21 23 24 22) (280 1 100) simpleGrading (1 1 1) //5
hex (11 5 6 12 11 13 14 12) (14 1 3) simpleGrading (1 1 1) //6
);

edges
(
arc 3 4 (0 0 3.17502)
arc 5 6 (19.3937 0 3.17502)
arc 9 10 (0 0 4.76263)
arc 11 12 (6.4646 0 4.76263)
arc 17 18 (0 0 49.9996)
arc 23 24 (300 0 49.9996)

);

patches
(
patch
inlet
(
(0 3 4 0)
)


patch
outlet
(
(2 7 8 2)
(7 15 16 8)
(15 23 24 16)
)


wall
wall
(
(3 5 6 4)
(17 19 20 18)
(19 21 22 20)
(21 23 24 22)
(10 9 17 18)
(10 12 11 9)
(12 6 5 11)
)

wedge
front
(
(0 1 5 3)
(1 2 7 5)
(5 7 15 13)
(9 11 19 17)
(11 13 21 19)
(13 15 23 21)
(11 5 13 11)


)

wedge
back
(
(0 1 6 4)
(1 2 8 6)
(6 8 16 14)
(10 12 20 18)
(12 14 22 20)
(14 16 24 22)
(12 6 14 12)
)


);

mergePatchPairs //merge faces for patches created
(
);


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


p:

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

dimensions [0 2 -2 0 0 0 0]; //dimension [Mass Length Time Temperature Quantity Current Luminous intensity] IMPORTANT: OpenFOAM use rho normalized pressure, divide pressure by density!!

internalField uniform 0;

boundaryField
{

inlet
{
type uniformFixedValue;
uniformValue tableFile;
tableFileCoeffs
{
fileName "$FOAM_CASE/inletpressure"
outOfBounds clamp;
}

}

outlet
{
type fixedValue;
value uniform 0;
}

wall
{
type zeroGradient;
}

front
{
type wedge;
}

back
{
type wedge;
}
}

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


U:


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

dimensions [0 1 -1 0 0 0 0]; //dimension of the variable [Mass Length Time Temperature Quantity Current Luminous intensity]

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type zeroGradient;
}


outlet
{
type zeroGradient;
}

wall
{
type fixedValue;
value uniform (0 0 0); //no slip condition, 0m/s velocity
}

front
{
type wedge;
}

back
{
type wedge;
}
}

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

I also uploaded a picture of the detail of the wedge. The nozzle part.

Has anyone seen something similar before?


Thanks in advance,

Daniel