src/engine/engineTopoChangerMesh/simpleTwoStroke/simpleTwoStrokeCalculate.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM 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 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27
  28 #include "simpleTwoStroke.H"
  29 #include "slidingInterface.H"
  30 #include "layerAdditionRemoval.H"
  31 #include "surfaceFields.H"
  32 #include "regionSplit.H"
  33 #include "componentMixedTetPolyPatchVectorField.H"
  34
  35 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  36
  37 void Foam::simpleTwoStroke::correctVerticalMotion()
  38 {
  39 }
  40
  41
  42 void Foam::simpleTwoStroke::calcMovingMasks() const
  43 {
  44     if (debug)
  45     {
  46         Info<< "void movingSquaresTM::calcMovingMasks() const : "
  47             << "Calculating point and cell masks"
  48             << endl;
  49     }
  50
  51     if (movingPointsMaskPtr_)
  52     {
  53         FatalErrorIn("void movingSquaresTM::calcMovingMasks() const")
  54             << "point mask already calculated"
  55             << abort(FatalError);
  56     }
  57
  58     // Set the point mask
  59     movingPointsMaskPtr_ = new scalarField(allPoints().size(), 0);
  60     scalarField& movingPointsMask = *movingPointsMaskPtr_;
  61
  62     const cellList& c = cells();
  63     const faceList& f = allFaces();
  64
  65     const labelList& cellAddr =
  66         cellZones()[cellZones().findZoneID("movingCells")];
  67
  68     forAll (cellAddr, cellI)
  69     {
  70         const cell& curCell = c[cellAddr[cellI]];
  71
  72         forAll (curCell, faceI)
  73         {
  74             // Mark all the points as moving
  75             const face& curFace = f[curCell[faceI]];
  76
  77             forAll (curFace, pointI)
  78             {
  79                 movingPointsMask[curFace[pointI]] = 1;
  80             }
  81         }
  82     }
  83
  84
  85     if(foundScavPorts())
  86     {
  87
  88         const word innerScavZoneName
  89         (
  90             scavInCylPatchName_  + "Zone"
  91         );
  92
  93         const labelList& innerScavAddr =
  94             faceZones()[faceZones().findZoneID(innerScavZoneName)];
  95
  96         forAll (innerScavAddr, faceI)
  97         {
  98             const face& curFace = f[innerScavAddr[faceI]];
  99
 100             forAll (curFace, pointI)
 101             {
 102                 movingPointsMask[curFace[pointI]] = 1;
 103             }
 104         }
 105
 106         const word outerScavZoneName
 107         (
 108             scavInPortPatchName_ + "Zone"
 109         );
 110
 111         const labelList& outerScavAddr =
 112             faceZones()[faceZones().findZoneID(outerScavZoneName)];
 113
 114         forAll (outerScavAddr, faceI)
 115         {
 116             const face& curFace = f[outerScavAddr[faceI]];
 117
 118             forAll (curFace, pointI)
 119             {
 120                 movingPointsMask[curFace[pointI]] = 0;
 121             }
 122         }
 123
 124     }
 125
 126 }
 127
 128 // Return moving points mask.  Moving points marked with 1
 129 const Foam::scalarField& Foam::simpleTwoStroke::movingPointsMask() const
 130 {
 131     if(movingPointsMaskPtr_)
 132     {
 133         movingPointsMaskPtr_ = NULL;
 134     }
 135
 136     if (!movingPointsMaskPtr_)
 137     {
 138         calcMovingMasks();
 139     }
 140
 141     return *movingPointsMaskPtr_;
 142 }
 143
 144 void Foam::simpleTwoStroke::applyMovingMask()
 145 {
 146 }
 147