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-2021 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  const scalar t = time_.value() - time_.beginTime().value();
75 
76  if (t > small)
77  {
78  // Rotational position (in radians)
79  const vector omega
80  (
81  t*degToRad(radialVelocity_.x()),
82  t*degToRad(radialVelocity_.y()),
83  t*degToRad(radialVelocity_.z())
84  );
85 
86  const scalar magOmega = mag(omega);
87  const quaternion R(omega/magOmega, magOmega);
88  const septernion TR(septernion(-origin_)*R*septernion(origin_));
89 
90  DebugInFunction << "Time = " << t << " transformation: " << TR << endl;
91 
92  return TR;
93  }
94  else
95  {
96  return septernion::I;
97  }
98 }
99 
100 
102 (
103  const dictionary& SBMFCoeffs
104 )
105 {
106  solidBodyMotionFunction::read(SBMFCoeffs);
107 
108  SBMFCoeffs_.lookup("origin") >> origin_;
109  SBMFCoeffs_.lookup("radialVelocity") >> radialVelocity_;
110 
111  return true;
112 }
113 
114 
115 // ************************************************************************* //
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:156
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:69
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
static const septernion I
Definition: septernion.H:83
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.