applications/solvers/basic/sixDOFSolver/sixDOFSolver.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2006 H. Jasak All rights reserved
   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 Application
  26     sixDOFsolver
  27
  28 Author
  29     Dubravko Matijasevic, FSB Zagreb.  All rights reserved.
  30
  31 Description
  32     Six degrees of freedom solver for multiple bodies
  33
  34 \*---------------------------------------------------------------------------*/
  35
  36 #include "argList.H"
  37 #include "ODESolver.H"
  38 #include "sixDOFbodies.H"
  39 #include "OFstream.H"
  40
  41 using namespace Foam;
  42
  43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  44
  45 int main(int argc, char *argv[])
  46 {
  47 #   include "setRootCase.H"
  48 #   include "createTime.H"
  49
  50     sixDOFbodies structure(runTime);
  51     OFstream of(runTime.path()/"motion.dat");
  52
  53     Info<< "\nStarting time loop\n" << endl;
  54
  55     for (runTime++; !runTime.end(); runTime++)
  56     {
  57         Info << "Time = " << runTime.value() << endl;
  58
  59         structure.solve();
  60
  61         of << runTime.value() << tab;
  62
  63         forAll (structure.names(), bodyI)
  64         {
  65             Info<< nl << "Average velocity of " << structure.names()[bodyI]
  66                 << " in time step = "
  67                 << structure()[bodyI].Uaverage().value() << nl
  68                 << "Current velocity in time instant = "
  69                 << structure()[bodyI].U().value() << nl
  70                 << "Average omega of " << structure.names()[bodyI]
  71                 << " in time step = "
  72                 << structure()[bodyI].omegaAverage().value() << nl
  73                 << "Current omega in time instant = "
  74                 << structure()[bodyI].omega().value()  << nl
  75                 << "Average omegaAbs in time step = "
  76                 << structure()[bodyI].omegaAverageAbsolute().value() << nl
  77                 << endl;
  78
  79             of  << structure()[bodyI].X().value().x() << tab;
  80         }
  81
  82         of << endl;
  83
  84         runTime.write();
  85
  86         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
  87             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
  88             << nl << endl;
  89     }
  90
  91     Info<< "End\n" << endl;
  92
  93     return(0);
  94 }
  95
  96
  97 // ************************************************************************* //