[gns4566]

hello, I'd like to make hemisphere domain.


*************************/
R 38.1E-03;
deg1 4.0;
Lx #calc "3*$R";
x0 #calc "-$Lx";
x1 #calc "-$R";
x2 #calc "-$Lx*cos($deg1)";
x3 #calc "-$R*cos($deg1)";
x4 0;
y0 0;
y1 $R;
y2 $Lx;
y3 #calc "$R*sin($deg1)";
y4 #calc "$Lx*sin($deg1)";
z0 0;
z1 1E-03;

//point on curve-R
deg0 20.0;
deg2 2.0;
rx0 #calc "$R*cos($deg0)";
ry0 #calc "$R*sin($deg0)";
rx1 #calc "$Lx*cos($deg0)";
ry1 #calc "$Lx*sin($deg0)";
rx2 #calc "$R*cos($deg2)";
ry2 #calc "$R*sin($deg2)";
rx3 #calc "$Lx*cos($deg2)";
ry3 #calc "$Lx*sin($deg2)";

convertToMeters 1;

vertices
(
($x0 $y0 $z0)//0
($x1 $y0 $z0)//1
($x3 $y2 $z0)//2
($x4 $y1 $z0)//3
($x4 $y4 $z0)//4
($x2 $y3 $z0)//5


($x0 $y0 $z1)//6
($x1 $y0 $z1)//7
($x3 $y2 $z1)//8
($x4 $y1 $z1)//9
($x4 $y4 $z0)//10
($x2 $y3 $z0)//11


);

blocks
(
hex (2 3 4 5 8 9 10 11) (200 400 1) simpleGrading (1 1 1)
hex (1 2 5 0 7 8 11 6) (200 400 1) simpleGrading (1 1 1)


);


edges
(
arc 2 3 ($rx0 $ry0 $z0)
arc 8 9 ($rx0 $ry0 $z1)
arc 5 4 ($rx1 $ry1 $z0)
arc 11 10 ($rx1 $ry1 $z1)
arc 1 2 ($rx2 $ry2 $z0)
arc 7 8 ($rx2 $ry2 $z1)
arc 0 5 ($rx3 $ry3 $z0)
arc 6 11 ($rx3 $ry3 $z1)
);


boundary
(
axi
{
type symmetryPlane;
faces
(
(0 6 7 1)
);
}
wall
{
type wall;
faces
(
(1 7 8 2)
(2 8 9 3)
);
}
out
{
type patch;
faces
(
(10 4 3 9)
);
}
in
{
type patch;
faces
(
(0 5 11 6)
(5 4 10 11)
);
}
front
{
type empty;
faces
(
(6 7 8 11)
(11 8 9 10)
);
}
back
{
type empty;
faces
(
(0 1 2 5)
(5 2 3 4)
);
}
);

mergePatchPairs
(
);



but I don't know what is the real problem of this...???


I just dived one body into two part.



Pleas help this...

--> FOAM Warning :
From function Foam::vector Foam::cellModel::centre(const labelList&, const pointField&) const
in file meshes/meshShapes/cellModel/cellModel.C at line 70
zero or negative pyramid volume: -1.680368821e-06 for face 1
--> FOAM Warning :
From function Foam::vector Foam::cellModel::centre(const labelList&, const pointField&) const
in file meshes/meshShapes/cellModel/cellModel.C at line 70
zero or negative pyramid volume: -5.00112815e-07 for face 2
--> FOAM Warning :
From function Foam::vector Foam::cellModel::centre(const labelList&, const pointField&) const
in file meshes/meshShapes/cellModel/cellModel.C at line 70
zero or negative pyramid volume: -1.158958419e-06 for face 3
--> FOAM Warning :
From function Foam::vector Foam::cellModel::centre(const labelList&, const pointField&) const
in file meshes/meshShapes/cellModel/cellModel.C at line 70
zero or negative pyramid volume: -1.05676112e-06 for face 4
--> FOAM Warning :
From function Foam::vector Foam::cellModel::centre(const labelList&, const pointField&) const
in file meshes/meshShapes/cellModel/cellModel.C at line 70
zero or negative pyramid volume: -1.05676112e-06 for face 5






--> FOAM FATAL IO ERROR:
Block hex (1 2 5 0 7 8 11 6) (200 400 1) simpleGrading (1(1) 1(1) 1(1)) has inward-pointing faces
4(2 5 11 8) 4(1 2 8 7) 4(0 6 11 5) 4(1 0 5 2) 4(7 8 11 6)



FOAM exiting

[Krao]

Hey GNS,

you have forgotten to add the picture! if you add someone would help you

[Flowkersma]

Without more information, my best guess is that the vertices of the block are in wrong order.

[gns4566]

Originally Posted by Krao

Hey GNS,


you have forgotten to add the picture! if you add someone would help you

Thank you for your advice, I add my blockMesh