src/engine/engineTopoChangerMesh/thoboisSliding/addThoboisSlidingMeshModifiers.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 Class
  25     verticalValvesGambit
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "thoboisSliding.H"
  30 #include "slidingInterface.H"
  31 #include "layerAdditionRemoval.H"
  32 #include "surfaceFields.H"
  33 #include "regionSplit.H"
  34 #include "attachDetach.H"
  35
  36 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  37
  38 void Foam::thoboisSliding::addZonesAndModifiers()
  39 {
  40     // Add the zones and mesh modifiers to operate piston motion
  41
  42     if
  43     (
  44         pointZones().size() > 0
  45      || faceZones().size() > 0
  46      || cellZones().size() > 0
  47      || topoChanger_.size() > 0
  48     )
  49     {
  50         Info<< "Time = " << engTime().theta() << endl;
  51         Info<< "void Foam::verticalValvesGambit::addZonesAndModifiers() : "
  52             << "Zones and modifiers already present.  Skipping."
  53             << endl;
  54
  55         setVirtualPositions();
  56 //        checkAndCalculate();
  57
  58         Info << "Point zones found = " << pointZones().size() << endl;
  59         Info << "Face zones found = " << faceZones().size() << endl;
  60         Info << "Cell zones found = " << cellZones().size() << endl;
  61
  62         return;
  63
  64     }
  65
  66     if
  67     (
  68         engTime().engineDict().found("zOffsetGambit")
  69      && engTime().engineDict().found("zDisplGambit")
  70     )
  71     {
  72         Info << "Assembling the cylinder mesh" << endl;
  73
  74         scalar zOffset
  75         (
  76             readScalar(engTime().engineDict().lookup("zOffsetGambit"))
  77         );
  78
  79         scalar zDispl
  80         (
  81             readScalar(engTime().engineDict().lookup("zDisplGambit"))
  82         );
  83
  84         pointField pDispl = points();
  85
  86         forAll(points(), pointI)
  87         {
  88             const point p = points()[pointI];
  89
  90             if(p.z() >= zOffset)
  91             {
  92                 pDispl[pointI].z() -= zDispl;
  93             }
  94         }
  95
  96         movePoints(pDispl);
  97         write();
  98         resetMotion();
  99
 100         Info << "Cylinder mesh assembled" << endl;
 101     }
 102
 103     Info << "checkAndCalculate()" << endl;
 104 //    checkAndCalculate();
 105
 106     Info<< "Time = " << engTime().theta() << endl
 107         << "Adding zones to the engine mesh" << endl;
 108
 109 /*
 110     Point zones
 111     1) Piston points
 112     1) Cut point zone for liner in head
 113
 114     nValves*
 115     1) cutPointsV
 116     2) valveTopPoints
 117     3) valveBottomPoints
 118 */
 119     DynamicList<pointZone*> pz;
 120
 121 /*
 122     Face zones
 123     1) Piston layer faces
 124
 125     nValves*
 126     1) valveTopLayerFaces
 127     2) valveBottomLayerFaces
 128     3) valveCurtainPort
 129     4) valveCurtainCyl
 130     5) cutFaceV
 131 */
 132     DynamicList<faceZone*> fz;
 133
 134 /*
 135     cell zones
 136     1) moving cells inside piston
 137
 138     nValves*
 139     1) moving cells in the top of the valve
 140     2) moving cells in the bottom of the valve
 141 */
 142
 143     DynamicList<cellZone*> cz;
 144
 145     label nPointZones = 0;
 146     label nFaceZones = 0;
 147     label nCellZones = 0;
 148
 149 /*
 150     Adding the following faces zones:
 151     1:  pistonLayerFaces
 152     nV: pistonLayerFacesV
 153
 154     Adding the following cell zones:
 155     1:  movingCellsPiston
 156     nV:  movingCellsPistonV
 157
 158     Adding the following point zones:
 159     1: pistonPoints
 160     nV: valvePistonPointsV
 161
 162 */
 163
 164 //#   include "addPistonFacesPointZonesThoboisSliding.H"
 165
 166 /*
 167     Adding the following face zones:
 168
 169     nV: curtainCylZoneV
 170     nV: curtainPortZoneV
 171     nV: cutFaceZoneV
 172     nV: poppetZoneV
 173     nV: bottomZoneV
 174
 175     Adding the following point zones:
 176
 177     nV: cutPointsV
 178
 179 */
 180
 181
 182 #   include "addValvesFacesPointZonesThoboisSliding.H"
 183
 184 /*
 185
 186     Adding the following point zones:
 187
 188     nV: valveTopPointsV
 189     nV: valveBottomPointsV
 190
 191     Adding the following cell zones:
 192
 193     nV: movingCellsTopV
 194     nV: movingCellsBotV
 195
 196 */
 197
 198 //#   include "addValvePistonCellZonesThoboisSliding.H"
 199
 200 //#   include "addAttachDetachFacesThoboisSliding.H"
 201
 202     Info<< "Adding " << nPointZones << " point, "
 203         << nFaceZones << " face zones and "
 204         << nCellZones << " cell zones" << endl;
 205
 206     pz.setSize(nPointZones);
 207     Info << "setSize pz" << endl;
 208     fz.setSize(nFaceZones);
 209     Info << "setSize fz" << endl;
 210     cz.setSize(nCellZones);
 211     Info << "setSize cz" << endl;
 212
 213     addZones(pz, fz, cz);
 214
 215 #   include "addMeshModifiersThoboisSliding.H"
 216
 217     // Calculating the virtual positions of piston and valves
 218
 219     setVirtualPositions();
 220
 221     Info << " Write mesh " << endl;
 222     // Write mesh and modifiers
 223     topoChanger_.writeOpt() = IOobject::AUTO_WRITE;
 224     write();
 225     Info << " Mesh written " << endl;
 226
 227 //    Info << "virtualPistonPosition = " << virtualPistonPosition() << endl;
 228 //    Info << "piston position = " << pistonPosition() << endl;
 229 }
 230
 231
 232 // ************************************************************************* //