src/engine/engineMesh/layeredEngineMesh/layeredEngineMesh.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "layeredEngineMesh.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "fvcMeshPhi.H"
  29 #include "surfaceInterpolate.H"
  30
  31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35     defineTypeNameAndDebug(layeredEngineMesh, 0);
  36     addToRunTimeSelectionTable(engineMesh, layeredEngineMesh, IOobject);
  37 }
  38
  39
  40 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  41
  42 Foam::layeredEngineMesh::layeredEngineMesh(const IOobject& io)
  43 :
  44     engineMesh(io),
  45     pistonLayers_("pistonLayers", dimLength, 0.0)
  46 {
  47     engineDB_.engineDict().readIfPresent("pistonLayers", pistonLayers_);
  48 }
  49
  50
  51 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  52
  53 Foam::layeredEngineMesh::~layeredEngineMesh()
  54 {}
  55
  56
  57 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  58
  59 void Foam::layeredEngineMesh::move()
  60 {
  61     scalar deltaZ = engineDB_.pistonDisplacement().value();
  62     Info<< "deltaZ = " << deltaZ << endl;
  63
  64     // Position of the top of the static mesh layers above the piston
  65     scalar pistonPlusLayers = pistonPosition_.value() + pistonLayers_.value();
  66
  67     pointField newPoints(points());
  68
  69     forAll(newPoints, pointi)
  70     {
  71         point& p = newPoints[pointi];
  72
  73         if (p.z() < pistonPlusLayers)           // In piston bowl
  74         {
  75             p.z() += deltaZ;
  76         }
  77         else if (p.z() < deckHeight_.value())   // In liner region
  78         {
  79             p.z() +=
  80                 deltaZ
  81                *(deckHeight_.value() - p.z())
  82                /(deckHeight_.value() - pistonPlusLayers);
  83         }
  84     }
  85
  86     if (engineDB_.foundObject<surfaceScalarField>("phi"))
  87     {
  88         surfaceScalarField& phi =
  89             const_cast<surfaceScalarField&>
  90             (engineDB_.lookupObject<surfaceScalarField>("phi"));
  91
  92         const volScalarField& rho =
  93             engineDB_.lookupObject<volScalarField>("rho");
  94
  95         const volVectorField& U =
  96             engineDB_.lookupObject<volVectorField>("U");
  97
  98         bool absolutePhi = false;
  99         if (moving())
 100         {
 101             phi += fvc::interpolate(rho)*fvc::meshPhi(rho, U);
 102             absolutePhi = true;
 103         }
 104
 105         movePoints(newPoints);
 106
 107         if (absolutePhi)
 108         {
 109             phi -= fvc::interpolate(rho)*fvc::meshPhi(rho, U);
 110         }
 111     }
 112     else
 113     {
 114         movePoints(newPoints);
 115     }
 116
 117     pistonPosition_.value() += deltaZ;
 118     scalar pistonSpeed = deltaZ/engineDB_.deltaTValue();
 119
 120     Info<< "clearance: " << deckHeight_.value() - pistonPosition_.value() << nl
 121         << "Piston speed = " << pistonSpeed << " m/s" << endl;
 122 }
 123
 124
 125 // ************************************************************************* //