tutorials/compressible/rhoPorousSimpleFoam/angledDuctImplicit/constant/polyMesh/blockMeshDict.m4

   1 /*--------------------------------*- C++ -*----------------------------------*\
   2 | =========                 |                                                 |
   3 | \\      /  F ield         | foam-extend: Open Source CFD                    |
   4 |  \\    /   O peration     | Version:     3.0                                |
   5 |   \\  /    A nd           | Web:         http://www.extend-project.de       |
   6 |    \\/     M anipulation  |                                                 |
   7 \*---------------------------------------------------------------------------*/
   8 FoamFile
   9 {
  10     version     2.0;
  11     format      ;
  12     class       dictionary;
  13     object      blockMeshDict;
  14 }
  15 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  16 // block definition for a porosity with an angled inlet/outlet
  17 // the porosity is not aligned with the main axes
  18 //
  19 dnl> -----------------------------------------------------------------
  20 dnl> <STANDARD DEFINTIONS>
  21 dnl>
  22 changecom(//)changequote([,]) dnl>
  23 define(calc, [esyscmd(perl -e 'print ($1)')]) dnl>
  24 define(VCOUNT, 0)  dnl>
  25 define(vlabel, [[// ]pt VCOUNT ($1) define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))])  dnl>
  26 dnl>
  27 define(hex2D, hex ($1b $2b $3b $4b $1f $2f $3f $4f)) dnl>
  28 define(quad2D, ($1f $1b $2b $2f))  dnl>
  29 define(frontQuad, ($1f $2f $3f $4f)) dnl>
  30 define(backQuad, ($4b $3b $2b $1b)) dnl>
  31 dnl>
  32 dnl> </STANDARD DEFINTIONS>
  33 dnl> -----------------------------------------------------------------
  34 dnl>
  35 define(ncells, 20) dnl>
  36 define(ninlet, 15) dnl>
  37 define(nporo, 20) dnl>
  38 define(noutlet, 20) dnl>
  39 dnl>
  40 define(x0,0) dnl>
  41 define(y0,0) dnl>
  42 define(y0,0) dnl>
  43 define(Cos,0.7071067812) dnl>   == cos(45)
  44 define(Sin,0.7071067812) dnl>   == sin(45)
  45 dnl>
  46 define(width,50) dnl>
  47 define(zBack,calc(-width/2)) dnl>
  48 define(zFront,calc(width/2)) dnl>
  49 define(leninlet,150)dnl>
  50 define(lenporo,100)dnl>
  51 define(lenoutlet,100)dnl>
  52 dnl>
  53 define(xhyp,calc(Sin*width)) dnl>
  54 define(yhyp,calc(Cos*width)) dnl>
  55 define(xinlet,leninlet)dnl>
  56 define(xporo,calc(Cos*lenporo)) dnl>
  57 define(yporo,calc(Sin*lenporo)) dnl>
  58 define(xoutlet,calc(xporo + Cos*lenoutlet)) dnl>
  59 define(youtlet,calc(yporo + Sin*lenoutlet)) dnl>
  60 dnl>
  61
  62 convertToMeters 0.001;
  63
  64 vertices
  65 (
  66     // inlet region
  67     ( -xinlet  y0  zBack )  vlabel(in1b)
  68     ( -xinlet yhyp  zBack ) vlabel(in2b)
  69     ( -xinlet  y0  zFront )  vlabel(in1f)
  70     ( -xinlet yhyp  zFront ) vlabel(in2f)
  71
  72     // join inlet->outlet
  73     (  x0 y0  zBack )    vlabel(join1b)
  74     ( -xhyp   yhyp  zBack ) vlabel(join2b)
  75     (  x0 y0  zFront )    vlabel(join1f)
  76     ( -xhyp   yhyp  zFront ) vlabel(join2f)
  77
  78     // porosity ends ->outlet
  79     ( xporo yporo  zBack )  vlabel(poro1b)
  80     ( calc(xporo - xhyp) calc(yporo + yhyp)  zBack )  vlabel(poro2b)
  81     ( xporo yporo  zFront )  vlabel(poro1f)
  82     ( calc(xporo - xhyp) calc(yporo + yhyp)  zFront )  vlabel(poro2f)
  83
  84     // outlet
  85     ( xoutlet youtlet zBack ) vlabel(out1b)
  86     ( calc(xoutlet - xhyp) calc(youtlet + yhyp) zBack ) vlabel(out2b)
  87     ( xoutlet youtlet zFront ) vlabel(out1f)
  88     ( calc(xoutlet - xhyp) calc(youtlet + yhyp) zFront ) vlabel(out2f)
  89 );
  90
  91 blocks
  92 (
  93     // inlet block
  94     hex2D(in1, join1, join2, in2)
  95     inlet ( ninlet ncells ncells ) simpleGrading (1 1 1)
  96
  97     // porosity block
  98     hex2D(join1, poro1, poro2, join2)
  99     porosity ( nporo ncells ncells ) simpleGrading (1 1 1)
 100
 101     // outlet block
 102     hex2D(poro1, out1, out2, poro2)
 103     outlet ( noutlet ncells ncells )  simpleGrading (1 1 1)
 104 );
 105
 106 edges
 107 (
 108 );
 109
 110 patches
 111 (
 112     // is there no way of defining all my 'defaultFaces' to be 'wall'?
 113     wall front
 114     (
 115     // inlet block
 116     frontQuad(in1, join1, join2, in2)
 117     // outlet block
 118     frontQuad(poro1, out1, out2, poro2)
 119     )
 120
 121     wall back
 122     (
 123     // inlet block
 124     backQuad(in1, join1, join2, in2)
 125     // outlet block
 126     backQuad(poro1, out1, out2, poro2)
 127     )
 128
 129     wall wall
 130     (
 131     // inlet block
 132     quad2D(in1, join1)
 133     quad2D(join2, in2)
 134     // outlet block
 135     quad2D(poro1, out1)
 136     quad2D(out2, poro2)
 137     )
 138
 139     wall porosityWall
 140     (
 141     // porosity block
 142     frontQuad(join1, poro1, poro2, join2)
 143     // porosity block
 144     backQuad(join1, poro1, poro2, join2)
 145     // porosity block
 146     quad2D(join1, poro1)
 147     quad2D(poro2, join2)
 148     )
 149
 150     patch inlet
 151     (
 152     quad2D(in2, in1)
 153     )
 154
 155     patch outlet
 156     (
 157     quad2D(out2, out1)
 158     )
 159 );
 160
 161 mergePatchPairs
 162 (
 163 );
 164
 165 // ************************************************************************* //