tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/patchifyObstacles

   1 #!/bin/bash
   2 #--------------------------------*- C++ -*----------------------------------#
   3 # =========                 |                                               #
   4 # \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox         #
   5 #  \\    /   O peration     | Version:  2.0.0                               #
   6 #   \\  /    A nd           | Web:      www.OpenFOAM.org                    #
   7 #    \\/     M anipulation  |                                               #
   8 #---------------------------------------------------------------------------#
   9 cd ${0%/*} || exit 1    # run from this directory
  10
  11 x0=0.4
  12 y0=0.1
  13 z0=0.1
  14 l=0.2
  15 offset=0.3
  16 tol=0.001
  17
  18 nx=1
  19 ny=2
  20 nz=2
  21
  22
  23 tmpSetSet='patchify.setSet'
  24 rm -f patchify.setSet >/dev/null 2>&1
  25
  26
  27 initialiseFaceSet()
  28 {
  29     echo "faceSet cubeFaces clear" >> $tmpSetSet;
  30 }
  31
  32 addToFaceSet()
  33 {
  34     name=$1
  35     x1=$2
  36     x2=$3
  37     y1=$4
  38     y2=$5
  39     z1=$6
  40     z2=$7
  41
  42     echo "faceSet cubeFaces add boxToFace ($x1 $y1 $z1) ($x2 $y2 $z2)" >> $tmpSetSet
  43     echo "   $name: (x1=$x1 y1=$y1 z1=$z1) (x2=$x2 y2=$y2 z2=$z2)"
  44 }
  45
  46 createSetsAndZone()
  47 {
  48     name=$1
  49     x1=$2
  50     x2=$3
  51     y1=$4
  52     y2=$5
  53     z1=$6
  54     z2=$7
  55
  56     echo "faceSet ${name}Faces new boxToFace ($x1 $y1 $z1) ($x2 $y2 $z2)" >> $tmpSetSet
  57     echo "   $name: (x1=$x1 y1=$y1 z1=$z1) (x2=$x2 y2=$y2 z2=$z2)"
  58 }
  59
  60 initialiseFaceSet
  61
  62
  63 # loop over cubes
  64 n=1
  65 x=$x0
  66 for xi in $(seq 1 1 $nx); do
  67    y=$y0
  68    for yi in $(seq 1 1 $ny); do
  69        z=$z0
  70        for zi in $(seq 1 1 $nz); do
  71            echo "processing cube $n"
  72
  73            pad=""
  74            if [ $n -lt 10 ]; then
  75                pad="0"
  76            fi
  77
  78            # side 1
  79            side=1
  80            x1=`echo $x $tol | awk '{print $1 - $2}'`
  81            x2=`echo $x $l $tol | awk '{print $1 + $2 + $3}'`
  82            y1=`echo $y $tol | awk '{print $1 - $2}'`
  83            y2=`echo $y $l $tol | awk '{print $1 + $2 + $3}'`
  84            z1=`echo $z $l $tol | awk '{print $1 + $2 - $3}'`
  85            z2=`echo $z $l $tol | awk '{print $1 + $2 + $3}'`
  86            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
  87
  88            #  side 2
  89            side=2
  90            x1=`echo $x $l $tol | awk '{print $1 + $2 - $3}'`
  91            x2=`echo $x $l $tol | awk '{print $1 + $2 + $3}'`
  92            y1=`echo $y $tol | awk '{print $1 - $2}'`
  93            y2=`echo $y $l $tol | awk '{print $1 + $2 + $3}'`
  94            z1=`echo $z $tol  | awk '{print $1 - $2}'`
  95            z2=`echo $z $l $tol | awk '{print $1 + $2 + $3}'`
  96            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
  97
  98            # side 3
  99            side=3
 100            x1=`echo $x $tol | awk '{print $1 - $2}'`
 101            x2=`echo $x $l $tol | awk '{print $1 + $2 + $3}'`
 102            y1=`echo $y $tol | awk '{print $1 - $2}'`
 103            y2=`echo $y $l $tol | awk '{print $1 + $2 + $3}'`
 104            z1=`echo $z $tol | awk '{print $1 - $2}'`
 105            z2=`echo $z $tol | awk '{print $1 + $2}'`
 106            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
 107
 108            # side 4
 109            side=4
 110            x1=`echo $x $tol | awk '{print $1 - $2}'`
 111            x2=`echo $x $tol | awk '{print $1 + $2}'`
 112            y1=`echo $y $tol | awk '{print $1 - $2}'`
 113            y2=`echo $y $l $tol | awk '{print $1 + $2 + $3}'`
 114            z1=`echo $z $tol | awk '{print $1 - $2}'`
 115            z2=`echo $z $l $tol | awk '{print $1 + $2 + $3}'`
 116            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
 117
 118            # side 5
 119            side=5
 120            x1=`echo $x $tol | awk '{print $1 - $2}'`
 121            x2=`echo $x $l $tol | awk '{print $1 + $2 + $3}'`
 122            y1=`echo $y $l $tol | awk '{print $1 + $2 - $3}'`
 123            y2=`echo $y $l $tol | awk '{print $1 + $2 + $3}'`
 124            z1=`echo $z $tol | awk '{print $1 - $2}'`
 125            z2=`echo $z $l $tol | awk '{print $1 + $2 + $3}'`
 126            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
 127
 128            # side 6
 129            side=6
 130            x1=`echo $x $tol | awk '{print $1 - $2}'`
 131            x2=`echo $x $l $tol | awk '{print $1 + $2 + $3}'`
 132            y1=`echo $y $tol | awk '{print $1 - $2}'`
 133            y2=`echo $y $tol | awk '{print $1 + $2}'`
 134            z1=`echo $z $tol | awk '{print $1 - $2}'`
 135            z2=`echo $z $l $tol | awk '{print $1 + $2 + $3}'`
 136            addToFaceSet cube${pad}${n}_side${side} $x1 $x2 $y1 $y2 $z1 $z2
 137
 138            let n+=1
 139
 140            z=`echo $z $offset | awk '{print $1 + $2}'`
 141        done
 142
 143        y=`echo $y $offset | awk '{print $1 + $2}'`
 144    done
 145
 146    x=`echo $x $offset | awk '{print $1 + $2}'`
 147 done
 148
 149 echo "cellSet cubeFacesCells new faceToCell cubeFaces owner" >> $tmpSetSet
 150 echo "faceZoneSet cubeFaces new setsToFaceZone cubeFaces cubeFacesCells" >> $tmpSetSet
 151
 152 echo "processing floor"
 153 floorMax=`echo 1 $tol | awk '{print $1 + $2}'`
 154 createSetsAndZone floor -$tol $floorMax -$tol $floorMax -$tol $tol
 155 echo "cellSet floorCells new faceToCell floorFaces owner" >> $tmpSetSet
 156 echo "faceZoneSet floorFaces new setsToFaceZone floorFaces floorCells" >> $tmpSetSet
 157
 158 setSet -batch $tmpSetSet > log.setSet.patchifyObstacles >/dev/null 2>&1
 159
 160
 161 # *************************************************************************
 162