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