axisRotationMotion.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2012-2018 OpenFOAM Foundation
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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your 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, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "axisRotationMotion.H"
28 #include "unitConversion.H"
29 
30 using namespace Foam::constant::mathematical;
31 
32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 namespace solidBodyMotionFunctions
37 {
40  (
44  );
45 }
46 }
47 
48 
49 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
50 
52 (
53  const dictionary& SBMFCoeffs,
54  const Time& runTime
55 )
56 :
57  solidBodyMotionFunction(SBMFCoeffs, runTime)
58 {
59  read(SBMFCoeffs);
60 }
61 
62 
63 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
64 
66 {}
67 
68 
69 // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
70 
73 {
74  scalar t = time_.value();
75 
76  // Rotation origin (in radians)
77  vector omega
78  (
79  t*degToRad(radialVelocity_.x()),
80  t*degToRad(radialVelocity_.y()),
81  t*degToRad(radialVelocity_.z())
82  );
83 
84  scalar magOmega = mag(omega);
85  quaternion R(omega/magOmega, magOmega);
86  septernion TR(septernion(-origin_)*R*septernion(origin_));
87 
88  DebugInFunction << "Time = " << t << " transformation: " << TR << endl;
89 
90  return TR;
91 }
92 
93 
95 (
96  const dictionary& SBMFCoeffs
97 )
98 {
100 
101  SBMFCoeffs_.lookup("origin") >> origin_;
102  SBMFCoeffs_.lookup("radialVelocity") >> radialVelocity_;
103 
104  return true;
105 }
106 
107 
108 // ************************************************************************* //
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
Unit conversion functions.
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
axisRotationMotion(const dictionary &SBMFCoeffs, const Time &runTime)
Construct from components.
Septernion class used to perform translations and rotations in 3D space.
Definition: septernion.H:65
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
Definition: Time.H:68
Macros for easy insertion into run-time selection tables.
scalar degToRad(const scalar deg)
Conversion from degrees to radians.
bool read(const char *, int32_t &)
Definition: int32IO.C:85
Base class for defining solid-body motions.
mathematical constants.
#define DebugInFunction
Report an information message using Foam::Info.
virtual bool read(const dictionary &SBMFCoeffs)=0
Update properties from given dictionary.
Quaternion class used to perform rotations in 3D space.
Definition: quaternion.H:60
virtual bool read(const dictionary &SBMFCoeffs)
Update properties from given dictionary.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
Constant velocity rotation around CoG. Similar to rotatingMotion but motion specified as rotation vec...
defineTypeNameAndDebug(combustionModel, 0)
#define R(A, B, C, D, E, F, K, M)
virtual septernion transformation() const
Return the solid-body motion transformation septernion.
dimensioned< scalar > mag(const dimensioned< Type > &)
Namespace for OpenFOAM.