src/engine/engineTopoChangerMesh/engineValveSliding/addEngineValveSlidingMeshModifiers.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 "engineValveSliding.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::engineValveSliding::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     Info << "checkAndCalculate()" << endl;
  67
  68     Info<< "Time = " << engTime().theta() << endl
  69         << "Adding zones to the engine mesh" << endl;
  70
  71 /*
  72     Point zones
  73     1) Piston points
  74     1) Cut point zone for liner in head
  75
  76     nValves*
  77     1) cutPointsV
  78     2) valveTopPoints
  79     3) valveBottomPoints
  80 */
  81
  82     DynamicList<pointZone*> pz;
  83
  84 /*
  85     Face zones
  86     1) Piston layer faces
  87
  88     nValves*
  89     1) valveTopLayerFaces
  90     2) valveBottomLayerFaces
  91     3) valveCurtainPort
  92     4) valveCurtainCyl
  93     5) cutFaceV
  94 */
  95     DynamicList<faceZone*> fz;
  96
  97 /*
  98     cell zones
  99     1) moving cells inside piston
 100
 101     nValves*
 102     1) moving cells in the top of the valve
 103     2) moving cells in the bottom of the valve
 104 */
 105
 106     DynamicList<cellZone*> cz;
 107
 108     label nPointZones = 0;
 109     label nFaceZones = 0;
 110     label nCellZones = 0;
 111
 112 /*
 113     Adding the following faces zones:
 114     1:  pistonLayerFaces
 115     nV: pistonLayerFacesV
 116
 117     Adding the following cell zones:
 118     1:  movingCellsPiston
 119     nV:  movingCellsPistonV
 120
 121     Adding the following point zones:
 122     1: pistonPoints
 123     nV: valvePistonPointsV
 124
 125 */
 126
 127 #   include "addPistonFacesPointZonesEngineValveSliding.H"
 128
 129 /*
 130     Adding the following face zones:
 131
 132     nV: curtainCylZoneV
 133     nV: curtainPortZoneV
 134     nV: cutFaceZoneV
 135     nV: poppetZoneV
 136     nV: bottomZoneV
 137
 138     Adding the following point zones:
 139
 140     nV: cutPointsV
 141
 142 */
 143
 144 #   include "addValvesFacesPointZonesEngineValveSliding.H"
 145
 146 /*
 147
 148     Adding the following point zones:
 149
 150     nV: valveTopPointsV
 151     nV: valveBottomPointsV
 152
 153     Adding the following cell zones:
 154
 155     nV: movingCellsTopV
 156     nV: movingCellsBotV
 157
 158 */
 159
 160 #   include "addOutputCellsEngineValveSliding.H"
 161
 162     Info<< "Adding " << nPointZones << " point, "
 163         << nFaceZones << " face zones and "
 164         << nCellZones << " cell zones" << endl;
 165
 166     pz.setSize(nPointZones);
 167     Info << "setSize pz" << endl;
 168     fz.setSize(nFaceZones);
 169     Info << "setSize fz" << endl;
 170     cz.setSize(nCellZones);
 171     Info << "setSize cz" << endl;
 172
 173     addZones(pz, fz, cz);
 174
 175 #   include "addMeshModifiersEngineValveSliding.H"
 176
 177     // Calculating the virtual positions of piston and valves
 178
 179     setVirtualPositions();
 180
 181     Info << " Write mesh " << endl;
 182     // Write mesh and modifiers
 183     topoChanger_.writeOpt() = IOobject::AUTO_WRITE;
 184     write();
 185
 186     Info << " Mesh written " << endl;
 187
 188 //    Info << "virtualPistonPosition = " << virtualPistonPosition() << endl;
 189 //    Info << "piston position = " << pistonPosition() << endl;
 190
 191 }
 192
 193
 194 // ************************************************************************* //