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