src/engine/engineTopoChangerMesh/simpleTwoStroke/simpleTwoStroke.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 "simpleTwoStroke.H"
  28 #include "layerAdditionRemoval.H"
  29 #include "componentMixedTetPolyPatchVectorField.H"
  30 #include "mapPolyMesh.H"
  31 #include "polyTopoChange.H"
  32 #include "addToRunTimeSelectionTable.H"
  33 #include "GeometricField.H"
  34 #include "volMesh.H"
  35 #include "engineTime.H"
  36 #include "pointField.H"
  37 #include "fvPatchField.H"
  38
  39 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  40
  41 namespace Foam
  42 {
  43     defineTypeNameAndDebug(simpleTwoStroke, 0);
  44     addToRunTimeSelectionTable(engineTopoChangerMesh, simpleTwoStroke, IOobject);
  45 }
  46
  47
  48 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  49
  50
  51
  52
  53 bool Foam::simpleTwoStroke::realDeformation() const
  54 {
  55
  56
  57     if (virtualPistonPosition()+ engTime().pistonDisplacement().value() > deckHeight()-engTime().clearance().value()-SMALL)
  58     {
  59         return true;
  60     }
  61     else
  62     {
  63         return deformation();
  64     }
  65
  66 }
  67
  68 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  69
  70 // Construct from components
  71 Foam::simpleTwoStroke::simpleTwoStroke
  72 (
  73     const IOobject& io
  74 )
  75 :
  76     engineTopoChangerMesh(io),
  77     piston_(*this, engTime().engineDict().subDict("piston")),
  78     scavInCylPatchName_(engTime().engineDict().lookup("scavInCylPatch")),
  79     scavInPortPatchName_(engTime().engineDict().lookup("scavInPortPatch")),
  80     movingPointsMaskPtr_(NULL),
  81     deformSwitch_(readScalar(engTime().engineDict().lookup("deformAngle"))),
  82     delta_(readScalar(engTime().engineDict().lookup("delta"))),
  83     offSet_(readScalar(engTime().engineDict().lookup("offSet"))),
  84     pistonPosition_(-GREAT),
  85     virtualPistonPosition_(-GREAT),
  86     deckHeight_(GREAT),
  87     csPtr_
  88     (
  89         coordinateSystem::New
  90         (
  91             "coordinateSystem",
  92             engTime().engineDict().subDict("coordinateSystem")
  93         )
  94     ),
  95     foundScavPorts_
  96     (
  97         boundaryMesh().findPatchID(scavInCylPatchName_) != -1
  98         &&
  99         boundaryMesh().findPatchID(scavInPortPatchName_) != -1
 100     ),
 101     scavPortsTol_(readScalar(engTime().engineDict().lookup("scavPortsTol")))
 102 {
 103     // Add zones and modifiers if not already there.
 104     addZonesAndModifiers();
 105 }
 106
 107
 108 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 109
 110
 111 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 112
 113
 114 void Foam::simpleTwoStroke::setBoundaryVelocity(volVectorField& U)
 115 {
 116     vector pistonVel = piston().cs().axis()*engTime().pistonSpeed().value();
 117
 118     // On the piston movingWallVelocity is used. There is no need to update
 119     // the piston velocity
 120
 121 //    U.boundaryField()[piston().patchID().index()] = pistonVel;
 122
 123     Info << "setting the boundary velocity" << endl;
 124     U.boundaryField()[boundaryMesh().findPatchID(scavInPortPatchName_)] ==
 125         pistonVel;
 126
 127 }
 128
 129 bool Foam::simpleTwoStroke::portsOpened() const
 130 {
 131
 132     scalar maxScavPortsZ = max
 133        (
 134            boundary()
 135            [
 136                boundaryMesh().findPatchID(scavInPortPatchName_)
 137            ].patch().localPoints()
 138        ).z();
 139
 140
 141     scalar minLinerZ = min
 142        (
 143            boundary()
 144            [
 145                boundaryMesh().findPatchID(scavInCylPatchName_)
 146            ].patch().localPoints()
 147        ).z();
 148
 149
 150     if (maxScavPortsZ - scavPortsTol_ < minLinerZ)
 151     {
 152         return false;
 153     }
 154     else
 155     {
 156         return true;
 157     }
 158 }
 159
 160
 161 // ************************************************************************* //