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
angularOscillatingVelocity
angularOscillatingVelocityPointPatchVectorField.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 "
angularOscillatingVelocityPointPatchVectorField.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
angularOscillatingVelocityPointPatchVectorField::
40
angularOscillatingVelocityPointPatchVectorField
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
angularOscillatingVelocityPointPatchVectorField::
57
angularOscillatingVelocityPointPatchVectorField
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
angularOscillatingVelocityPointPatchVectorField::
88
angularOscillatingVelocityPointPatchVectorField
89
(
90
const
angularOscillatingVelocityPointPatchVectorField
& 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_)
103
{}
104
105
106
angularOscillatingVelocityPointPatchVectorField::
107
angularOscillatingVelocityPointPatchVectorField
108
(
109
const
angularOscillatingVelocityPointPatchVectorField
& 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
angularOscillatingVelocityPointPatchVectorField::autoMap
126
(
127
const
pointPatchFieldMapper
& m
128
)
129
{
130
fixedValuePointPatchField<vector>::autoMap
(m);
131
132
p0_.autoMap(m);
133
}
134
135
136
void
angularOscillatingVelocityPointPatchVectorField::rmap
137
(
138
const
pointPatchField<vector>
& ptf,
139
const
labelList
& addr
140
)
141
{
142
const
angularOscillatingVelocityPointPatchVectorField
& aOVptf =
143
refCast<const angularOscillatingVelocityPointPatchVectorField>(ptf);
144
145
fixedValuePointPatchField<vector>::rmap
(aOVptf, addr);
146
147
p0_.rmap(aOVptf.p0_, addr);
148
}
149
150
151
void
angularOscillatingVelocityPointPatchVectorField::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
const
pointPatch
& p = this->
patch
();
161
162
scalar angle = angle0_ + amplitude_*
sin
(omega_*t.
value
());
163
vector
axisHat = axis_/
mag
(axis_);
164
vectorField
p0Rel = p0_ - origin_;
165
166
vectorField::operator
=
167
(
168
(
169
p0_
170
+ p0Rel*(
cos
(angle) - 1)
171
+ (axisHat ^ p0Rel*
sin
(angle))
172
+ (axisHat & p0Rel)*(1 -
cos
(angle))*axisHat
173
- p.
localPoints
()
174
)/t.
deltaTValue
()
175
);
176
177
fixedValuePointPatchField<vector>::updateCoeffs
();
178
}
179
180
181
void
angularOscillatingVelocityPointPatchVectorField::write
182
(
183
Ostream
& os
184
)
const
185
{
186
pointPatchField<vector>::write
(os);
187
os.
writeKeyword
(
"axis"
)
188
<< axis_ <<
token::END_STATEMENT
<<
nl
;
189
os.
writeKeyword
(
"origin"
)
190
<< origin_ <<
token::END_STATEMENT
<<
nl
;
191
os.
writeKeyword
(
"angle0"
)
192
<< angle0_ <<
token::END_STATEMENT
<<
nl
;
193
os.
writeKeyword
(
"amplitude"
)
194
<< amplitude_ <<
token::END_STATEMENT
<<
nl
;
195
os.
writeKeyword
(
"omega"
)
196
<< omega_ <<
token::END_STATEMENT
<<
nl
;
197
p0_.writeEntry(
"p0"
, os);
198
writeEntry(
"value"
, os);
199
}
200
201
202
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
203
204
makePointPatchTypeField
205
(
206
pointPatchVectorField
,
207
angularOscillatingVelocityPointPatchVectorField
208
);
209
210
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
211
212
}
// End namespace Foam
213
214
// ************************ vim: set sw=4 sts=4 et: ************************ //