tutorials/basic/PODSolver/1DPODsin/setTcos/setTcos.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 Description
  25     Set up the gamma field
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "fvCFD.H"
  30 #include "mathematicalConstants.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 int main(int argc, char *argv[])
  35 {
  36
  37 #   include "setRootCase.H"
  38 #   include "createTime.H"
  39 #   include "createMesh.H"
  40
  41     Info<< "Time = " << runTime.value() << endl;
  42
  43     Info<< "    Reading field T" << endl;
  44     volScalarField T
  45     (
  46         IOobject
  47         (
  48             "T",
  49             runTime.timeName(),
  50             mesh,
  51             IOobject::MUST_READ,
  52             IOobject::NO_WRITE
  53         ),
  54         mesh
  55     );
  56
  57     // Do cells
  58     const volVectorField& centres = mesh.C();
  59
  60     Info << pos(0.02 - centres.internalField().component(vector::X))*
  61             centres.internalField().component(vector::X)/0.02 << endl;
  62     T.internalField() =
  63         pos(0.02 - centres.internalField().component(vector::X))*
  64         (
  65             Foam::cos
  66             (
  67                 centres.internalField().component(vector::X)/0.02*
  68                 mathematicalConstant::pi
  69             )/2.0 + 0.5
  70         );
  71
  72     T.write();
  73
  74     Info<< "End\n" << endl;
  75
  76     return(0);
  77 }
  78
  79
  80 // ************************************************************************* //