src/engine/engineTopoChangerMesh/engineValveSliding/engineValveSlidingInitialize.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 Class
  25     verticalValvesGambit
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "engineValveSliding.H"
  30 #include "slidingInterface.H"
  31 #include "layerAdditionRemoval.H"
  32 #include "surfaceFields.H"
  33 #include "regionSplit.H"
  34
  35 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  36
  37 void Foam::engineValveSliding::checkAndCalculate()
  38 {
  39
  40     label pistonIndex = -1;
  41     bool foundPiston = false;
  42
  43     label linerIndex = -1;
  44     bool foundLiner = false;
  45
  46     label cylinderHeadIndex = -1;
  47     bool foundCylinderHead = false;
  48
  49
  50     forAll(boundary(), i)
  51     {
  52         if (boundary()[i].name() == "piston")
  53         {
  54             pistonIndex = i;
  55             foundPiston = true;
  56         }
  57         else if (boundary()[i].name() == "liner")
  58         {
  59             linerIndex = i;
  60             foundLiner = true;
  61         }
  62         else if (boundary()[i].name() == "cylinderHead")
  63         {
  64             cylinderHeadIndex = i;
  65             foundCylinderHead = true;
  66         }
  67     }
  68
  69     reduce(foundPiston, orOp<bool>());
  70     reduce(foundLiner, orOp<bool>());
  71     reduce(foundCylinderHead, orOp<bool>());
  72
  73     if (!foundPiston)
  74     {
  75         FatalErrorIn("Foam::verticalValvesGambit::checkAndCalculate()")
  76             << " : cannot find piston patch"
  77             << abort(FatalError);
  78     }
  79
  80     if (!foundLiner)
  81     {
  82         FatalErrorIn("Foam::verticalValvesGambit::checkAndCalculate()")
  83             << " : cannot find liner patch"
  84             << abort(FatalError);
  85     }
  86
  87     if (!foundCylinderHead)
  88     {
  89         FatalErrorIn("Foam::verticalValvesGambit::checkAndCalculate()")
  90             << " : cannot find cylinderHead patch"
  91             << exit(FatalError);
  92     }
  93
  94     {
  95         if (linerIndex != -1)
  96         {
  97             pistonPosition() =
  98                 max(boundary()[pistonIndex].patch().localPoints()).z();
  99         }
 100         reduce(pistonPosition(), minOp<scalar>());
 101
 102         if (cylinderHeadIndex != -1)
 103         {
 104             deckHeight() = min
 105             (
 106                 boundary()[cylinderHeadIndex].patch().localPoints()
 107             ).z();
 108
 109         }
 110         reduce(deckHeight(), minOp<scalar>());
 111
 112         Info<< "deckHeight: " << deckHeight() << nl
 113             << "piston position: " << pistonPosition() << endl;
 114     }
 115
 116
 117 }
 118
 119 void Foam::engineValveSliding::setVirtualPositions()
 120 {}