src/mesh/extrudeModel/wedge/wedge.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "wedge.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "mathematicalConstants.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34 namespace extrudeModels
  35 {
  36
  37 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  38
  39 defineTypeNameAndDebug(wedge, 0);
  40
  41 addToRunTimeSelectionTable(extrudeModel, wedge, dictionary);
  42
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 wedge::wedge(const dictionary& dict)
  47 :
  48     extrudeModel(typeName, dict),
  49     axisPt_(coeffDict_.lookup("axisPt")),
  50     axis_(coeffDict_.lookup("axis")),
  51     angle_
  52     (
  53         readScalar(coeffDict_.lookup("angle"))
  54        *mathematicalConstant::pi/180.0
  55     )
  56 {}
  57
  58
  59 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  60
  61 wedge::~wedge()
  62 {}
  63
  64
  65 // * * * * * * * * * * * * * * * * Operators * * * * * * * * * * * * * * * * //
  66
  67 point wedge::operator()
  68 (
  69     const point& surfacePoint,
  70     const vector& surfaceNormal,
  71     const label layer
  72 ) const
  73 {
  74     scalar sliceAngle;
  75     // For the case of a single layer extrusion assume a
  76     // symmetric wedge about the reference plane is required
  77     if (nLayers_ == 1)
  78     {
  79         if (layer == 0)
  80         {
  81             sliceAngle = -angle_/2.0;
  82         }
  83         else
  84         {
  85             sliceAngle = angle_/2.0;
  86         }
  87     }
  88     else
  89     {
  90         //sliceAngle = angle_*(layer + 1)/nLayers_;
  91         sliceAngle = angle_*layer/nLayers_;
  92     }
  93
  94     // Find projection onto axis (or rather decompose surfacePoint
  95     // into vector along edge (proj), vector normal to edge in plane
  96     // of surface point and surface normal.
  97     point d = surfacePoint - axisPt_;
  98
  99     d -= (axis_ & d)*axis_;
 100
 101     scalar dMag = mag(d);
 102
 103     point edgePt = surfacePoint - d;
 104
 105     // Rotate point around sliceAngle.
 106     point rotatedPoint = edgePt;
 107
 108     if (dMag > VSMALL)
 109     {
 110         vector n = (d/dMag) ^ axis_;
 111
 112         rotatedPoint +=
 113           + cos(sliceAngle)*d
 114           - sin(sliceAngle)*mag(d)*n; // Use either n or surfaceNormal
 115     }
 116
 117     return rotatedPoint;
 118 }
 119
 120
 121 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 122
 123 } // End namespace extrudeModels
 124 } // End namespace Foam
 125
 126 // ************************************************************************* //