src/engine/engineTopoChangerMesh/pistonSliding/addPistonSlidingMeshModifiers.C

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