Home
Downloads
Documentation
Installation
User Guide
man-pages
API Documentation
README
Release Notes
Changes
License
Support
SourceForge Project
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
src
fvMotionSolver
pointPatchFields
derived
angularOscillatingDisplacement
angularOscillatingDisplacementPointPatchVectorField.C
Go to the documentation of this file.
1
/*---------------------------------------------------------------------------*\
2
========= |
3
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4
\\ / O peration |
5
\\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
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 "
angularOscillatingDisplacementPointPatchVectorField.H
"
27
#include <
OpenFOAM/pointPatchFields.H
>
28
#include <
OpenFOAM/addToRunTimeSelectionTable.H
>
29
#include <
OpenFOAM/Time.H
>
30
#include <
OpenFOAM/polyMesh.H
>
31
32
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33
34
namespace
Foam
35
{
36
37
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
38
39
angularOscillatingDisplacementPointPatchVectorField::
40
angularOscillatingDisplacementPointPatchVectorField
41
(
42
const
pointPatch
&
p
,
43
const
DimensionedField<vector, pointMesh>
& iF
44
)
45
:
46
fixedValuePointPatchField<vector>
(
p
, iF),
47
axis_(
vector::zero
),
48
origin_(
vector::zero
),
49
angle0_(0.0),
50
amplitude_(0.0),
51
omega_(0.0),
52
p0_(p.
localPoints
())
53
{}
54
55
56
angularOscillatingDisplacementPointPatchVectorField::
57
angularOscillatingDisplacementPointPatchVectorField
58
(
59
const
pointPatch
& p,
60
const
DimensionedField<vector, pointMesh>
& iF,
61
const
dictionary
& dict
62
)
63
:
64
fixedValuePointPatchField<vector>
(
p
, iF, dict),
65
axis_(dict.
lookup
(
"axis"
)),
66
origin_(dict.
lookup
(
"origin"
)),
67
angle0_(
readScalar
(dict.
lookup
(
"angle0"
))),
68
amplitude_(
readScalar
(dict.
lookup
(
"amplitude"
))),
69
omega_(
readScalar
(dict.
lookup
(
"omega"
)))
70
{
71
if
(!dict.
found
(
"value"
))
72
{
73
updateCoeffs();
74
}
75
76
if
(dict.
found
(
"p0"
))
77
{
78
p0_ =
vectorField
(
"p0"
, dict , p.
size
());
79
}
80
else
81
{
82
p0_ = p.
localPoints
();
83
}
84
}
85
86
87
angularOscillatingDisplacementPointPatchVectorField::
88
angularOscillatingDisplacementPointPatchVectorField
89
(
90
const
angularOscillatingDisplacementPointPatchVectorField
& ptf,
91
const
pointPatch
& p,
92
const
DimensionedField<vector, pointMesh>
& iF,
93
const
pointPatchFieldMapper
& mapper
94
)
95
:
96
fixedValuePointPatchField<vector>
(ptf,
p
, iF, mapper),
97
axis_(ptf.axis_),
98
origin_(ptf.origin_),
99
angle0_(ptf.angle0_),
100
amplitude_(ptf.amplitude_),
101
omega_(ptf.omega_),
102
p0_(ptf.p0_, mapper)
103
{}
104
105
106
angularOscillatingDisplacementPointPatchVectorField::
107
angularOscillatingDisplacementPointPatchVectorField
108
(
109
const
angularOscillatingDisplacementPointPatchVectorField
& ptf,
110
const
DimensionedField<vector, pointMesh>
& iF
111
)
112
:
113
fixedValuePointPatchField<vector>
(ptf, iF),
114
axis_(ptf.axis_),
115
origin_(ptf.origin_),
116
angle0_(ptf.angle0_),
117
amplitude_(ptf.amplitude_),
118
omega_(ptf.omega_),
119
p0_(ptf.p0_)
120
{}
121
122
123
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
124
125
void
angularOscillatingDisplacementPointPatchVectorField::autoMap
126
(
127
const
pointPatchFieldMapper
& m
128
)
129
{
130
fixedValuePointPatchField<vector>::autoMap
(m);
131
132
p0_.autoMap(m);
133
}
134
135
136
void
angularOscillatingDisplacementPointPatchVectorField::rmap
137
(
138
const
pointPatchField<vector>
& ptf,
139
const
labelList
& addr
140
)
141
{
142
const
angularOscillatingDisplacementPointPatchVectorField
& aODptf =
143
refCast<const angularOscillatingDisplacementPointPatchVectorField>(ptf);
144
145
fixedValuePointPatchField<vector>::rmap
(aODptf, addr);
146
147
p0_.rmap(aODptf.p0_, addr);
148
}
149
150
151
void
angularOscillatingDisplacementPointPatchVectorField::updateCoeffs
()
152
{
153
if
(this->
updated
())
154
{
155
return
;
156
}
157
158
const
polyMesh
&
mesh
= this->
dimensionedInternalField
().
mesh
()();
159
const
Time
& t = mesh.
time
();
160
161
scalar angle = angle0_ + amplitude_*
sin
(omega_*t.
value
());
162
vector
axisHat = axis_/
mag
(axis_);
163
vectorField
p0Rel = p0_ - origin_;
164
165
vectorField::operator
=
166
(
167
p0Rel*(
cos
(angle) - 1)
168
+ (axisHat ^ p0Rel*
sin
(angle))
169
+ (axisHat & p0Rel)*(1 -
cos
(angle))*axisHat
170
);
171
172
fixedValuePointPatchField<vector>::updateCoeffs
();
173
}
174
175
176
void
angularOscillatingDisplacementPointPatchVectorField::write
177
(
178
Ostream
& os
179
)
const
180
{
181
pointPatchField<vector>::write
(os);
182
os.
writeKeyword
(
"axis"
)
183
<< axis_ <<
token::END_STATEMENT
<<
nl
;
184
os.
writeKeyword
(
"origin"
)
185
<< origin_ <<
token::END_STATEMENT
<<
nl
;
186
os.
writeKeyword
(
"angle0"
)
187
<< angle0_ <<
token::END_STATEMENT
<<
nl
;
188
os.
writeKeyword
(
"amplitude"
)
189
<< amplitude_ <<
token::END_STATEMENT
<<
nl
;
190
os.
writeKeyword
(
"omega"
)
191
<< omega_ <<
token::END_STATEMENT
<<
nl
;
192
p0_.writeEntry(
"p0"
, os);
193
writeEntry(
"value"
, os);
194
}
195
196
197
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
198
199
makePointPatchTypeField
200
(
201
pointPatchVectorField
,
202
angularOscillatingDisplacementPointPatchVectorField
203
);
204
205
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
206
207
}
// End namespace Foam
208
209
// ************************ vim: set sw=4 sts=4 et: ************************ //