oscillatingRotatingMotion.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) 2011-2023 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 
28 #include "mathematicalConstants.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 
53 (
54  const dictionary& SBMFCoeffs,
55  const Time& runTime
56 )
57 :
58  solidBodyMotionFunction(SBMFCoeffs, runTime)
59 {
60  read(SBMFCoeffs);
61 }
62 
63 
64 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
65 
68 {}
69 
70 
71 // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
72 
75 transformation() const
76 {
77  scalar t = time_.value();
78 
79  vector eulerAngles = amplitude_*sin(omega_*t);
80 
81  // Convert the rotational motion from deg to rad
82  eulerAngles *= pi/180.0;
83 
84  quaternion R(quaternion::XYZ, eulerAngles);
85  septernion TR(septernion(-origin_)*R*septernion(origin_));
86 
87  DebugInFunction << "Time = " << t << " transformation: " << TR << endl;
88 
89  return TR;
90 }
91 
92 
94 (
95  const dictionary& SBMFCoeffs
96 )
97 {
99 
100  SBMFCoeffs_.lookup("origin") >> origin_;
101  SBMFCoeffs_.lookup("amplitude") >> amplitude_;
102  SBMFCoeffs_.lookup("omega") >> omega_;
103 
104  return true;
105 }
106 
107 
108 // ************************************************************************* //
Macros for easy insertion into run-time selection tables.
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
Definition: Time.H:76
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
Quaternion class used to perform rotations in 3D space.
Definition: quaternion.H:61
Septernion class used to perform translations and rotations in 3D space.
Definition: septernion.H:66
Base class for defining solid-body motions.
virtual bool read(const dictionary &SBMFCoeffs)=0
Update properties from given dictionary.
SolidBodyMotionFvMesh 6DoF motion function. Oscillating rotation.
virtual septernion transformation() const
Return the solid-body motion transformation septernion.
oscillatingRotatingMotion(const dictionary &SBMFCoeffs, const Time &runTime)
Construct from components.
virtual bool read(const dictionary &SBMFCoeffs)
Update properties from given dictionary.
#define DebugInFunction
Report an information message using Foam::Info.
mathematical constants.
Namespace for OpenFOAM.
addToRunTimeSelectionTable(polyPatch, mergedCyclicPolyPatch, word)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
dimensionedScalar sin(const dimensionedScalar &ds)
defineTypeNameAndDebug(combustionModel, 0)
static scalar R(const scalar a, const scalar x)
Definition: invIncGamma.C:102