src/dynamicMesh/meshMotion/RBFMotionSolver/RBFMotionSolver.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 Class
  26     RBFMotionSolver
  27
  28 Description
  29     Radial basis function motion solver
  30
  31 Author
  32     Frank Bos, TU Delft.  All rights reserved.
  33
  34 SourceFiles
  35     RBFMotionSolver.C
  36
  37 \*---------------------------------------------------------------------------*/
  38
  39 #ifndef RBFMotionSolver_H
  40 #define RBFMotionSolver_H
  41
  42 #include "motionSolver.H"
  43 #include "polyMesh.H"
  44 #include "RBFInterpolation.H"
  45
  46 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  47
  48 namespace Foam
  49 {
  50
  51 /*---------------------------------------------------------------------------*\
  52                       Class RBFMotionSolver Declaration
  53 \*---------------------------------------------------------------------------*/
  54
  55 class RBFMotionSolver
  56 :
  57     public motionSolver
  58 {
  59     // Private data
  60
  61         //- Moving patch names
  62         wordList movingPatches_;
  63
  64         //- Static patch names
  65         wordList staticPatches_;
  66
  67         //- Coarsening ratio
  68         label coarseningRatio_;
  69
  70         //- Include zero motion of static patches in RBF interpolation
  71         Switch includeStaticPatches_;
  72
  73         //- Frozen interpolation
  74         Switch frozenInterpolation_;
  75
  76         //- Moving point IDs
  77         labelList movingIDs_;
  78
  79         //- Moving points on the boundary
  80         mutable vectorField movingPoints_;
  81
  82         //- Static point IDs
  83         labelList staticIDs_;
  84
  85         //- Control point IDs
  86         labelList controlIDs_;
  87
  88         //- Control points on the boundary
  89         mutable vectorField controlPoints_;
  90
  91         //- Internal point IDs
  92         labelList internalIDs_;
  93
  94         //- Internal points
  95         vectorField internalPoints_;
  96
  97         //- Prescribed motion on control points
  98         vectorField motion_;
  99
 100         //- RBF interpolation
 101         RBFInterpolation interpolation_;
 102
 103
 104     // Private Member Functions
 105
 106         //- Disallow default bitwise copy construct
 107         RBFMotionSolver(const RBFMotionSolver&);
 108
 109         //- Disallow default bitwise assignment
 110         void operator=(const RBFMotionSolver&);
 111
 112
 113         //- Make control point IDs.  Constructor helper
 114         void makeControlIDs();
 115
 116         //- Set location of points
 117         void setMovingPoints() const;
 118
 119
 120 public:
 121
 122     //- Runtime type information
 123     TypeName("RBFMotionSolver");
 124
 125
 126     // Constructors
 127
 128         //- Construct from polyMesh
 129         RBFMotionSolver
 130         (
 131             const polyMesh& mesh,
 132             Istream& msData
 133         );
 134
 135
 136     // Destructor
 137
 138         virtual ~RBFMotionSolver();
 139
 140
 141     // Member Functions
 142
 143         //- Return moving point indices in mesh points
 144         const labelList& movingIDs() const
 145         {
 146             return movingIDs_;
 147         }
 148
 149         //- Return moving points
 150         const vectorField& movingPoints() const;
 151
 152         //- Set motion of moving points
 153         void setMotion(const vectorField&);
 154
 155
 156         //- Return point location obtained from the current motion field
 157         virtual tmp<pointField> curPoints() const;
 158
 159         //- Solve for motion
 160         virtual void solve();
 161
 162         //- Update the mesh corresponding to given map
 163         virtual void updateMesh(const mapPolyMesh&);
 164 };
 165
 166
 167 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 168
 169 } // End namespace Foam
 170
 171 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 172
 173 #endif
 174
 175 // ************************************************************************* //