src/engine/engineTopoChangerMesh/layerSmooth/addLayerSmoothMeshModifiers.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 #include "layerSmooth.H"
  28 #include "slidingInterface.H"
  29 #include "layerAdditionRemoval.H"
  30 #include "surfaceFields.H"
  31 #include "regionSplit.H"
  32 #include "attachDetach.H"
  33 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  34
  35 void Foam::layerSmooth::addZonesAndModifiers()
  36 {
  37     // Add the zones and mesh modifiers to operate piston motion
  38
  39     if
  40     (
  41 /*
  42         pointZones().size() > 0
  43      || faceZones().size() > 0
  44      || cellZones().size() > 0
  45      || topoChanger_.size() > 0
  46 */
  47         pointZones().size() > 0
  48      && faceZones().size() > 0
  49      && topoChanger_.size() > 0
  50     )
  51     {
  52         Info<< "Time = " << engTime().theta() << endl;
  53         Info<< "void Foam::layerSmooth::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     else
  67     {
  68         pointZones().setSize(0);
  69         faceZones().setSize(0);
  70         cellZones().setSize(0);
  71         topoChanger_.setSize(0);
  72     }
  73
  74     if
  75     (
  76         engTime().engineDict().found("zOffsetGambit")
  77      && engTime().engineDict().found("zDisplGambit")
  78     )
  79     {
  80         Info << "Assembling the cylinder mesh" << endl;
  81
  82         scalar zOffset
  83         (
  84             readScalar(engTime().engineDict().lookup("zOffsetGambit"))
  85         );
  86
  87         scalar zDispl
  88         (
  89             readScalar(engTime().engineDict().lookup("zDisplGambit"))
  90         );
  91
  92         const pointField& ap = allPoints();
  93
  94         pointField pDispl = allPoints();
  95
  96         forAll(ap, pointI)
  97         {
  98             const point p = ap[pointI];
  99
 100             if (p.z() >= zOffset)
 101             {
 102                 pDispl[pointI].z() -= zDispl;
 103             }
 104         }
 105
 106
 107         movePoints(pDispl);
 108         write();
 109         resetMotion();
 110
 111         Info << "Cylinder mesh assembled" << endl;
 112     }
 113
 114
 115
 116     Info << "checkAndCalculate()" << endl;
 117     checkAndCalculate();
 118
 119     Info<< "Time = " << engTime().theta() << endl
 120         << "Adding zones to the engine mesh" << endl;
 121
 122 /*
 123     Point zones
 124     1) Piston points
 125 */
 126
 127     DynamicList<pointZone*> pz;
 128
 129 /*
 130     Face zones
 131     1) Piston layer faces
 132
 133 */
 134     DynamicList<faceZone*> fz;
 135
 136 /*
 137 */
 138
 139     DynamicList<cellZone*> cz;
 140
 141     label nPointZones = 0;
 142     label nFaceZones = 0;
 143     label nCellZones = 0;
 144
 145 /*
 146
 147     Adding the following faces zones:
 148     1:  pistonLayerFaces
 149
 150     Adding the following point zones:
 151     1:  pistonPoints
 152
 153 */
 154
 155 #   include "addPistonFacesPointZonesLayerSmooth.H"
 156
 157
 158 #   include "addAttachDetachFacesLayerSmooth.H"
 159
 160     Info<< "Adding " << nPointZones << " point zones, "
 161         << nFaceZones << " face zones and "
 162         << nCellZones << " cell zones" << endl;
 163
 164     pz.setSize(nPointZones);
 165     Info << "setSize pz" << endl;
 166     fz.setSize(nFaceZones);
 167     Info << "setSize fz" << endl;
 168     cz.setSize(nCellZones);
 169     Info << "setSize cz" << endl;
 170
 171     addZones(pz, fz, cz);
 172
 173 #   include "addMeshModifiersLayerSmooth.H"
 174
 175     // Calculating the virtual positions of piston and valves
 176
 177     setVirtualPositions();
 178
 179     // Write mesh
 180     write();
 181
 182     Info << "virtualPistonPosition = " << virtualPistonPosition() << endl;
 183     Info << "piston position = " << pistonPosition() << endl;
 184 }
 185
 186
 187 // ************************************************************************* //