src/immersedBoundary/immersedBoundaryDynamicMesh/movingImmersedBoundary/movingImmersedBoundary.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     |
   5     \\  /    A nd           | For copyright notice see file Copyright
   6      \\/     M anipulation  |
   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 \*---------------------------------------------------------------------------*/
  25
  26 #include "movingImmersedBoundary.H"
  27 #include "immersedBoundaryPolyPatch.H"
  28 #include "transformField.H"
  29
  30 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  31
  32 Foam::movingImmersedBoundary::movingImmersedBoundary
  33 (
  34     const word& name,
  35     const fvMesh& mesh,
  36     const dictionary& dict
  37 )
  38 :
  39     name_(name),
  40     mesh_(mesh),
  41     sbmfPtr_(solidBodyMotionFunction::New(dict, mesh.time())),
  42     refIbSurface_
  43     (
  44         IOobject
  45         (
  46             name  + ".ftr",
  47             mesh.time().constant(),      // instance
  48             "triSurface",                // local
  49             mesh,                        // registry
  50             IOobject::MUST_READ,
  51             IOobject::NO_WRITE
  52         )
  53     )
  54 {}
  55
  56
  57 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  58
  59 Foam::movingImmersedBoundary::~movingImmersedBoundary()
  60 {}
  61
  62
  63 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  64
  65 void Foam::movingImmersedBoundary::movePoints() const
  66 {
  67     // Get ibMesh from patch
  68     const label patchID = mesh().boundaryMesh().findPatchID(name());
  69
  70     if (patchID < 0)
  71     {
  72         FatalErrorIn
  73         (
  74             "void movingImmersedBoundary::movePoints() const"
  75         )   << "Patch " << name() << " not found.  Available patch names: "
  76             << mesh().boundaryMesh().names()
  77             << abort(FatalError);
  78     }
  79
  80     // Get non-const reference to velocity field
  81     volVectorField& U = const_cast<volVectorField&>
  82     (
  83         mesh().lookupObject<volVectorField>("U")
  84     );
  85
  86     // Get non-const reference to patch field
  87     immersedBoundaryFvPatchVectorField& ibPatchField =
  88         refCast<immersedBoundaryFvPatchVectorField>
  89         (
  90             U.boundaryField()[patchID]
  91         );
  92
  93
  94     const immersedBoundaryPolyPatch& cibPatch =
  95         refCast<const immersedBoundaryPolyPatch>
  96         (
  97             mesh().boundaryMesh()[patchID]
  98         );
  99
 100     // Get non-const reference to patch
 101     immersedBoundaryPolyPatch& ibPatch =
 102         const_cast<immersedBoundaryPolyPatch&>(cibPatch);
 103
 104     const vectorField oldIbPoints = ibPatch.ibMesh().coordinates();
 105
 106     // Move points
 107     ibPatch.moveTriSurfacePoints
 108     (
 109         transform(sbmfPtr_->transformation(), refIbSurface_.points())
 110     );
 111
 112     // Set refValue_ to moving boundary velocity
 113     ibPatchField.refValue() =
 114         (ibPatch.ibMesh().coordinates() - oldIbPoints)/
 115         mesh().time().deltaT().value();
 116 }
 117
 118
 119 // ************************************************************************* //