src/dynamicMesh/motionSolver/motionSolver.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 \*---------------------------------------------------------------------------*/
  26
  27 #include "motionSolver.H"
  28 #include "Time.H"
  29 #include "polyMesh.H"
  30 #include "dlLibraryTable.H"
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36     defineTypeNameAndDebug(motionSolver, 0);
  37     defineRunTimeSelectionTable(motionSolver, dictionary);
  38 }
  39
  40 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  41
  42 Foam::motionSolver::motionSolver(const polyMesh& mesh)
  43 :
  44     IOdictionary
  45     (
  46         IOobject
  47         (
  48             "dynamicMeshDict",
  49             mesh.time().constant(),
  50             mesh,
  51             IOobject::MUST_READ,
  52             IOobject::NO_WRITE
  53         )
  54     ),
  55     mesh_(mesh),
  56     twoDPointCorrector_(mesh)
  57 {}
  58
  59
  60 // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * //
  61
  62 Foam::autoPtr<Foam::motionSolver> Foam::motionSolver::New(const polyMesh& mesh)
  63 {
  64     IOdictionary solverDict
  65     (
  66         IOobject
  67         (
  68             "dynamicMeshDict",
  69             mesh.time().constant(),
  70             mesh,
  71             IOobject::MUST_READ,
  72             IOobject::NO_WRITE,
  73             false
  74         )
  75     );
  76
  77     Istream& msData = solverDict.lookup("solver");
  78
  79     word solverTypeName(msData);
  80
  81     Info << "Selecting motion solver: " << solverTypeName << endl;
  82
  83     dlLibraryTable::open
  84     (
  85         solverDict,
  86         "motionSolverLibs",
  87         dictionaryConstructorTablePtr_
  88     );
  89
  90     if (!dictionaryConstructorTablePtr_)
  91     {
  92         FatalErrorIn
  93         (
  94             "motionSolver::New(const polyMesh& mesh)"
  95         )   << "solver table is empty"
  96             << exit(FatalError);
  97     }
  98
  99     dictionaryConstructorTable::iterator cstrIter =
 100         dictionaryConstructorTablePtr_->find(solverTypeName);
 101
 102     if (cstrIter == dictionaryConstructorTablePtr_->end())
 103     {
 104         FatalErrorIn
 105         (
 106             "motionSolver::New(const polyMesh& mesh)"
 107         )   << "Unknown solver type " << solverTypeName
 108             << endl << endl
 109             << "Valid solver types are: " << endl
 110             << dictionaryConstructorTablePtr_->toc()
 111             << exit(FatalError);
 112     }
 113
 114     return autoPtr<motionSolver>(cstrIter()(mesh, msData));
 115 }
 116
 117
 118 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 119
 120 Foam::motionSolver::~motionSolver()
 121 {}
 122
 123
 124 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 125
 126 Foam::tmp<Foam::pointField> Foam::motionSolver::newPoints()
 127 {
 128     solve();
 129     return curPoints();
 130 }
 131
 132
 133 void Foam::motionSolver::twoDCorrectPoints(pointField& p) const
 134 {
 135     twoDPointCorrector_.correctPoints(p);
 136 }
 137
 138
 139 void Foam::motionSolver::updateMesh(const mapPolyMesh& mpm)
 140 {
 141     twoDPointCorrector_.updateMesh();
 142 }
 143
 144
 145 // ************************************************************************* //