applications/solvers/basic/PODSolver/PODSolver.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 Application
  25     PODsolver
  26
  27 Description
  28     Proper orthogonal decomposition solver with run-time selectable
  29     physical model
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #include "fvCFD.H"
  34 #include "PODODE.H"
  35 #include "ODESolver.H"
  36
  37 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  38
  39 int main(int argc, char *argv[])
  40 {
  41
  42 #   include "setRootCase.H"
  43
  44 #   include "createTime.H"
  45 #   include "createMesh.H"
  46
  47     Info<< "Reading PODsolverDict\n" << endl;
  48
  49     IOdictionary PODsolverDict
  50     (
  51         IOobject
  52         (
  53             "PODsolverDict",
  54             runTime.system(),
  55             mesh,
  56             IOobject::MUST_READ,
  57             IOobject::NO_WRITE
  58         )
  59     );
  60
  61     // Create equation
  62     autoPtr<PODODE> equation
  63     (
  64         PODODE::New(mesh, PODsolverDict)
  65     );
  66
  67     // Create solver
  68     autoPtr<ODESolver> solver =
  69         ODESolver::New
  70         (
  71             PODsolverDict.lookup("solver"),
  72             equation()
  73         );
  74
  75     // Read required accuracy
  76     scalar eps = readScalar(PODsolverDict.lookup("eps"));
  77
  78     Info<< "\nStarting time loop\n" << endl;
  79
  80     for (runTime++; !runTime.end(); runTime++)
  81     {
  82         Info<< "Time = " << runTime.timeName() << nl << endl;
  83
  84         solver->solve
  85         (
  86             runTime.value(),
  87             runTime.value() + runTime.deltaT().value(),
  88             eps,
  89             runTime.deltaT().value()
  90         );
  91
  92         // Create fields and recalculate properties
  93         if (runTime.outputTime())
  94         {
  95             Info << "Writing fields" << endl;
  96             equation->write();
  97         }
  98     }
  99
 100     Info<< "End\n" << endl;
 101
 102     return(0);
 103 }
 104
 105
 106 // ************************************************************************* //