January 20, 2021, 07:10
|
[BlockMesh] Problem with blockMesh internal faces
|
#1
|
|
New Member
parth
Join Date: Feb 2020
Posts: 24
Rep Power: 7
|
Dear all,
I am new to OpenFOAM. My blockMeshDict is looking as follows:
Code:
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
scale 1;
rin 2.09;
rout 2.1;
h0 0;
h1 1.8;
delh 1.0005;
sin30 0.5;
cos30 0.866;
cos0 1;
sin0 0;
dr1 0.01;
dr2 0.06;
// (#calc "$rin*$cos30" #calc "$rin*-$sin30" $h0) //0
// (#calc "$rout*$cos30" #calc "$rout*-$sin30" $h0) //1
// (#calc "$rout*$cos30" #calc "$rout*$sin30" $h0) //2
// (#calc "$rin*$cos30" #calc "$rin*$sin30" $h0) //3
vertices
(
(#calc "$rin*$cos30" #calc "$rin*-$sin30" $h0) //0
(#calc "$rout*$cos30" #calc "$rout*-$sin30" $h0) //1
(#calc "$rout*$cos30" #calc "$rout*$sin30" $h0) //2
(#calc "$rin*$cos30" #calc "$rin*$sin30" $h0) //3
(#calc "$rin*$cos30" #calc "$rin*-$sin30" $h1) //4
(#calc "$rout*$cos30" #calc "$rout*-$sin30" $h1) //5
(#calc "$rout*$cos30" #calc "$rout*$sin30" $h1) //6
(#calc "$rin*$cos30" #calc "$rin*$sin30" $h1) //7
(#calc "$rin*$cos30" #calc "$rin*-$sin30" $h1) //8
(#calc "$rout*$cos30" #calc "$rout*-$sin30" $h1) //9
(#calc "$rout*$cos30" #calc "$rout*$sin30" $h1) //10
(#calc "$rin*$cos30" #calc "$rin*$sin30" $h1) //11
(#calc "$rin*$cos30" #calc "$rin*-$sin30" #calc "$h1+$delh") //12
(#calc "$rout*$cos30" #calc "$rout*-$sin30" #calc "$h1+$delh") //13
(#calc "$rout*$cos30" #calc "$rout*$sin30" #calc "$h1+$delh") //14
(#calc "$rin*$cos30" #calc "$rin*$sin30" #calc "$h1+$delh") //15
(#calc "($rin+$dr1)*$cos30" #calc "($rin+$dr1)*-$sin30" $h1) //16
(#calc "($rout+$dr2)*$cos30" #calc "($rout+$dr2)*-$sin30" $h1) //17
(#calc "($rout+$dr2)*$cos30" #calc "($rout+$dr2)*$sin30" $h1) //18
(#calc "($rin+$dr1)*$cos30" #calc "($rin+$dr1)*$sin30" $h1) //19
(#calc "($rin+$dr1)*$cos30" #calc "($rin+$dr1)*-$sin30" #calc "$h1+$delh") //20
(#calc "($rout+$dr2)*$cos30" #calc "($rout+$dr2)*-$sin30" #calc "$h1+$delh") //21
(#calc "($rout+$dr2)*$cos30" #calc "($rout+$dr2)*$sin30" #calc "$h1+$delh") //22
(#calc "($rin+$dr1)*$cos30" #calc "($rin+$dr1)*$sin30" #calc "$h1+$delh") //23
);
blocks
(
//
hex (0 1 2 3 4 5 6 7) ( 30 300 300) simpleGrading (1 1 1) // bird
hex (8 9 10 11 12 13 14 15) ( 30 300 20) simpleGrading (1 1 1) //beak
hex (16 17 18 19 20 21 22 23) ( 70 300 10) simpleGrading (3 1 1) //beak
);
edges
(
arc 0 3 ($rin 0 $h0)
arc 1 2 ($rout 0 $h0)
arc 4 7 ($rin 0 $h1)
arc 5 6 ($rout 0 $h1)
arc 8 11 ($rin 0 #calc "$h1")
arc 9 10 ($rout 0 #calc "$h1")
arc 12 15 ($rin 0 #calc "$h1+$delh")
arc 13 14 ($rout 0 #calc "$h1+$delh")
arc 16 19 (#calc "$rin+$dr1" 0 #calc "$h1")
arc 17 18 (#calc "$rout+$dr2" 0 #calc "$h1")
arc 20 23 (#calc "$rin+$dr1" 0 #calc "$h1+$delh")
arc 21 22 (#calc "$rout+$dr2" 0 #calc "$h1+$delh")
);
boundary
(
inlet
{
type patch;
faces
(
(0 3 2 1)
);
}
walls
{
type wall;
faces
(
(7 3 0 4)
(12 15 11 8)
(2 6 5 1)
);
}
wedge1
{
type symmetry;
faces
(
(5 4 0 1)
(8 9 13 12)
(16 17 21 20)
);
}
wedge2
{
type symmetry;
faces
(
(3 7 6 2)
(15 14 10 11)
(23 22 18 19)
);
}
top_walls
{
type wall;
faces
(
(16 19 18 17)
(12 13 14 15)
(20 21 22 23)
);
}
outlet
{
type patch;
faces
(
(21 17 18 22)
);
}
defaultFaces
{
type empty;
faces
(
(4 5 6 7)
(8 9 10 11)
(23 20 16 19)
(14 13 9 10)
);
}
);
mergePatchPairs
(
);
It generates some internal faces which blocks the flow during simulation.
Is my blockMesh is correct?
|
|
|