CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

[BlockMesh] Problem with blockMesh internal faces

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 20, 2021, 06:10
Default [BlockMesh] Problem with blockMesh internal faces
  #1
New Member
 
parth
Join Date: Feb 2020
Posts: 23
Rep Power: 6
parthigcar is on a distinguished road
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?
parthigcar is offline   Reply With Quote

Reply

Tags
blockmesh, blockmeshdict, internal faces


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] Can't Shake Erros: patch type 'patch' not constraint type 'empty' BrendaEM OpenFOAM Meshing & Mesh Conversion 12 April 3, 2022 18:32
[snappyHexMesh] Help with Snappy: no layers growing GianF OpenFOAM Meshing & Mesh Conversion 2 September 23, 2020 08:26
SU2-7.0.1 on ubuntu 18.04 hyunko SU2 Installation 7 March 16, 2020 04:37
[snappyHexMesh] layer not added Rasmusiwersen OpenFOAM Meshing & Mesh Conversion 1 January 2, 2020 09:43
Cluster ID's not contiguous in compute-nodes domain. ??? Shogan FLUENT 1 May 28, 2014 15:03


All times are GMT -4. The time now is 18:16.