src/engine/engineTopoChangerMesh/accordionEngineMesh/addAccordionEngineMeshModifiers.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 "accordionEngineMesh.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::accordionEngineMesh::addZonesAndModifiers()
  36 {
  37     // Add the zones and mesh modifiers to operate piston motion
  38
  39     if
  40     (
  41         pointZones().size() > 0
  42 /*
  43      && faceZones().size() > 0
  44      && cellZones().size() > 0
  45      && topoChanger_.size() > 0
  46 */
  47     )
  48     {
  49         Info<< "Time = " << engTime().theta() << endl;
  50         Info<< "void Foam::accordionEngineMesh::addZonesAndModifiers() : "
  51             << "Zones and modifiers already present.  Skipping."
  52             << endl;
  53
  54      //   setVirtualPositions();
  55         checkAndCalculate();
  56
  57         Info << "Point zones found = " << pointZones().size() << endl;
  58         Info << "Face zones found = " << faceZones().size() << endl;
  59         Info << "Cell zones found = " << cellZones().size() << endl;
  60
  61         return;
  62
  63     }
  64     else
  65     {
  66         pointZones().setSize(0);
  67         faceZones().setSize(0);
  68         cellZones().setSize(0);
  69         topoChanger_.setSize(0);
  70     }
  71
  72     Info << "checkAndCalculate()" << endl;
  73     checkAndCalculate();
  74
  75     Info<< "Time = " << engTime().theta() << endl
  76         << "Adding zones to the engine mesh" << endl;
  77
  78 /*
  79     Point zones
  80     1) Piston points
  81 */
  82
  83     DynamicList<pointZone*> pz;
  84
  85 /*
  86     Face zones
  87     1) Piston layer faces
  88
  89 */
  90     DynamicList<faceZone*> fz;
  91
  92 /*
  93 */
  94
  95     DynamicList<cellZone*> cz;
  96
  97     label nPointZones = 0;
  98     label nFaceZones = 0;
  99     label nCellZones = 0;
 100
 101 /*
 102
 103     Adding the following faces zones:
 104     1:  pistonLayerFaces
 105     nv: movingFaces
 106     nv: staticFaces
 107     nv: detachFaces
 108
 109     Adding the following point zones:
 110     1:  pistonPoints
 111     nv:  movingValvePoints
 112     nv:  staticValvePoints
 113
 114 */
 115
 116 #   include "addPistonFacesPointZonesAccordionEngineMesh.H"
 117
 118 #   include "addAttachDetachFacesAccordionEngineMesh.H"
 119 #   include "addValveFaceZonesAccordionEngineMesh.H"
 120
 121 #   include "addOutputCellsAccordionEngineMesh.H"
 122     Info<< "Adding " << nPointZones << " point, "
 123         << nFaceZones << " face zones and "
 124         << nCellZones << " cell zones" << endl;
 125
 126
 127     pz.setSize(nPointZones);
 128     Info << "setSize pz" << endl;
 129     fz.setSize(nFaceZones);
 130     Info << "setSize fz" << endl;
 131     cz.setSize(nCellZones);
 132     Info << "setSize cz" << endl;
 133
 134     addZones(pz, fz, cz);
 135
 136 #   include "addMeshModifiersAccordionEngineMesh.H"
 137
 138     // Calculating the virtual positions of piston and valves
 139
 140  //   setVirtualPositions();
 141
 142     // Write mesh and modifiers
 143     topoChanger_.writeOpt() = IOobject::AUTO_WRITE;
 144     write();
 145
 146     Info << "virtualPistonPosition = " << virtualPistonPosition() << endl;
 147     Info << "piston position = " << pistonPosition() << endl;
 148 }
 149
 150
 151 // ************************************************************************* //