tutorials/basic/potentialFoam/cylinder/constant/polyMesh/blockMeshDict

   1 /*--------------------------------*- C++ -*----------------------------------*\
   2 | =========                 |                                                 |
   3 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
   4 |  \\    /   O peration     | Version:  2.0.0                                 |
   5 |   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
   6 |    \\/     M anipulation  |                                                 |
   7 \*---------------------------------------------------------------------------*/
   8 FoamFile
   9 {
  10     version     2.0;
  11     format      ascii;
  12     class       dictionary;
  13     object      blockMeshDict;
  14 }
  15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  16
  17 convertToMeters 1;
  18
  19 vertices #codeStream
  20 {
  21     codeInclude
  22     #{
  23         #include "pointField.H"
  24     #};
  25
  26     code
  27     #{
  28         pointField points(19);
  29         points[0]  = point(0.5, 0, -0.5);
  30         points[1]  = point(1, 0, -0.5);
  31         points[2]  = point(2, 0, -0.5);
  32         points[3]  = point(2, 0.707107, -0.5);
  33         points[4]  = point(0.707107, 0.707107, -0.5);
  34         points[5]  = point(0.353553, 0.353553, -0.5);
  35         points[6]  = point(2, 2, -0.5);
  36         points[7]  = point(0.707107, 2, -0.5);
  37         points[8]  = point(0, 2, -0.5);
  38         points[9]  = point(0, 1, -0.5);
  39         points[10] = point(0, 0.5, -0.5);
  40         points[11] = point(-0.5, 0, -0.5);
  41         points[12] = point(-1, 0, -0.5);
  42         points[13] = point(-2, 0, -0.5);
  43         points[14] = point(-2, 0.707107, -0.5);
  44         points[15] = point(-0.707107, 0.707107, -0.5);
  45         points[16] = point(-0.353553, 0.353553, -0.5);
  46         points[17] = point(-2, 2, -0.5);
  47         points[18] = point(-0.707107, 2, -0.5);
  48
  49         // Duplicate z points
  50         label sz = points.size();
  51         points.setSize(2*sz);
  52         for (label i = 0; i < sz; i++)
  53         {
  54             const point& pt = points[i];
  55             points[i+sz] = point(pt.x(), pt.y(), -pt.z());
  56         }
  57
  58         os  << points;
  59     #};
  60 };
  61
  62
  63 blocks
  64 (
  65     hex (5 4 9 10 24 23 28 29) (10 10 1) simpleGrading (1 1 1)
  66     hex (0 1 4 5 19 20 23 24) (10 10 1) simpleGrading (1 1 1)
  67     hex (1 2 3 4 20 21 22 23) (20 10 1) simpleGrading (1 1 1)
  68     hex (4 3 6 7 23 22 25 26) (20 20 1) simpleGrading (1 1 1)
  69     hex (9 4 7 8 28 23 26 27) (10 20 1) simpleGrading (1 1 1)
  70     hex (15 16 10 9 34 35 29 28) (10 10 1) simpleGrading (1 1 1)
  71     hex (12 11 16 15 31 30 35 34) (10 10 1) simpleGrading (1 1 1)
  72     hex (13 12 15 14 32 31 34 33) (20 10 1) simpleGrading (1 1 1)
  73     hex (14 15 18 17 33 34 37 36) (20 20 1) simpleGrading (1 1 1)
  74     hex (15 9 8 18 34 28 27 37) (10 20 1) simpleGrading (1 1 1)
  75 );
  76
  77 edges
  78 (
  79     arc 0 5 (0.469846 0.17101 -0.5)
  80     arc 5 10 (0.17101 0.469846 -0.5)
  81     arc 1 4 (0.939693 0.34202 -0.5)
  82     arc 4 9 (0.34202 0.939693 -0.5)
  83     arc 19 24 (0.469846 0.17101 0.5)
  84     arc 24 29 (0.17101 0.469846 0.5)
  85     arc 20 23 (0.939693 0.34202 0.5)
  86     arc 23 28 (0.34202 0.939693 0.5)
  87     arc 11 16 (-0.469846 0.17101 -0.5)
  88     arc 16 10 (-0.17101 0.469846 -0.5)
  89     arc 12 15 (-0.939693 0.34202 -0.5)
  90     arc 15 9 (-0.34202 0.939693 -0.5)
  91     arc 30 35 (-0.469846 0.17101 0.5)
  92     arc 35 29 (-0.17101 0.469846 0.5)
  93     arc 31 34 (-0.939693 0.34202 0.5)
  94     arc 34 28 (-0.34202 0.939693 0.5)
  95 );
  96
  97 boundary
  98 (
  99     down
 100     {
 101         type symmetryPlane;
 102         faces
 103         (
 104             (0 1 20 19)
 105             (1 2 21 20)
 106             (12 11 30 31)
 107             (13 12 31 32)
 108         );
 109     }
 110     right
 111     {
 112         type patch;
 113         faces
 114         (
 115             (2 3 22 21)
 116             (3 6 25 22)
 117         );
 118     }
 119     up
 120     {
 121         type symmetryPlane;
 122         faces
 123         (
 124             (7 8 27 26)
 125             (6 7 26 25)
 126             (8 18 37 27)
 127             (18 17 36 37)
 128         );
 129     }
 130     left
 131     {
 132         type patch;
 133         faces
 134         (
 135             (14 13 32 33)
 136             (17 14 33 36)
 137         );
 138     }
 139     cylinder
 140     {
 141         type symmetryPlane;
 142         faces
 143         (
 144             (10 5 24 29)
 145             (5 0 19 24)
 146             (16 10 29 35)
 147             (11 16 35 30)
 148         );
 149     }
 150 );
 151
 152 mergePatchPairs
 153 (
 154 );
 155
 156 // ************************************************************************* //