src/engine/engineTopoChangerMesh/twoStrokeEngine/twoStrokeEngineInitialize.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
  28 #include "twoStrokeEngine.H"
  29 #include "slidingInterface.H"
  30 #include "layerAdditionRemoval.H"
  31 #include "surfaceFields.H"
  32 #include "regionSplit.H"
  33
  34 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  35
  36 void Foam::twoStrokeEngine::checkAndCalculate()
  37 {
  38     label pistonIndex = -1;
  39     bool foundPiston = false;
  40
  41     label linerIndex = -1;
  42     bool foundLiner = false;
  43
  44     label cylinderHeadIndex = -1;
  45     bool foundCylinderHead = false;
  46
  47     forAll(boundary(), i)
  48     {
  49         Info << boundary()[i].name() << endl;
  50         if (boundary()[i].name() == "piston")
  51         {
  52             pistonIndex = i;
  53             foundPiston = true;
  54         }
  55         else if (boundary()[i].name() == "liner")
  56         {
  57             linerIndex = i;
  58             foundLiner = true;
  59         }
  60         else if (boundary()[i].name() == "cylinderHead")
  61         {
  62             cylinderHeadIndex = i;
  63             foundCylinderHead = true;
  64         }
  65     }
  66
  67     reduce(foundPiston, orOp<bool>());
  68     reduce(foundLiner, orOp<bool>());
  69     reduce(foundCylinderHead, orOp<bool>());
  70
  71     if (!foundPiston)
  72     {
  73         FatalErrorIn("Foam::twoStrokeEngine::checkAndCalculate()")
  74             << " : cannot find piston patch"
  75             << abort(FatalError);
  76     }
  77
  78     if (!foundLiner)
  79     {
  80         FatalErrorIn("Foam::twoStrokeEngine::checkAndCalculate()")
  81             << " : cannot find liner patch"
  82             << abort(FatalError);
  83     }
  84
  85     if (!foundCylinderHead)
  86     {
  87         FatalErrorIn("Foam::twoStrokeEngine::checkAndCalculate()")
  88             << " : cannot find cylinderHead patch"
  89             << exit(FatalError);
  90     }
  91
  92     {
  93         if (linerIndex != -1)
  94         {
  95             pistonPosition() =
  96                 max(boundary()[pistonIndex].patch().localPoints()).z();
  97         }
  98         reduce(pistonPosition(), minOp<scalar>());
  99
 100         if (cylinderHeadIndex != -1)
 101         {
 102             deckHeight() = min
 103             (
 104                 boundary()[cylinderHeadIndex].patch().localPoints()
 105             ).z();
 106         }
 107         reduce(deckHeight(), minOp<scalar>());
 108
 109         Info<< "deckHeight: " << deckHeight() << nl
 110             << "piston position: " << pistonPosition() << endl;
 111     }
 112 }
 113
 114 void Foam::twoStrokeEngine::setVirtualPistonPosition()
 115 {
 116
 117     label pistonFaceIndex = faceZones().findZoneID("pistonLayerFaces");
 118
 119     bool foundPistonFace = (pistonFaceIndex != -1);
 120
 121     Info << "piston face index = " << pistonFaceIndex << endl;
 122
 123     if(!foundPistonFace)
 124     {
 125         FatalErrorIn("Foam::twoStrokeEngine::setVirtualPistonPosition()")
 126             << " : cannot find the pistonLayerFaces"
 127             << exit(FatalError);
 128     }
 129
 130     const labelList& pistonFaces = faceZones()[pistonFaceIndex];
 131     forAll(pistonFaces, i)
 132     {
 133         const face& f = faces()[pistonFaces[i]];
 134
 135         // should loop over facepoints...
 136         forAll(f, j)
 137         {
 138             virtualPistonPosition() =
 139                 Foam::max(virtualPistonPosition(), points()[f[j]].z());
 140         }
 141     }
 142
 143     reduce(virtualPistonPosition(), maxOp<scalar>());
 144 }