oscillatingVelocityPointPatchVectorField.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) 2011 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 
27 #include "pointPatchFields.H"
29 #include "Time.H"
30 #include "polyMesh.H"
31 
32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 
37 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
38 
41 (
42  const pointPatch& p,
44 )
45 :
47  amplitude_(vector::zero),
48  omega_(0.0),
49  p0_(p.localPoints())
50 {}
51 
52 
55 (
56  const pointPatch& p,
58  const dictionary& dict
59 )
60 :
62  amplitude_(dict.lookup("amplitude")),
63  omega_(readScalar(dict.lookup("omega")))
64 {
65  if (!dict.found("value"))
66  {
67  updateCoeffs();
68  }
69 
70  if (dict.found("p0"))
71  {
72  p0_ = vectorField("p0", dict , p.size());
73  }
74  else
75  {
76  p0_ = p.localPoints();
77  }
78 }
79 
80 
83 (
85  const pointPatch& p,
87  const pointPatchFieldMapper& mapper
88 )
89 :
90  fixedValuePointPatchField<vector>(ptf, p, iF, mapper),
91  amplitude_(ptf.amplitude_),
92  omega_(ptf.omega_),
93  p0_(ptf.p0_, mapper)
94 {}
95 
96 
99 (
102 )
103 :
105  amplitude_(ptf.amplitude_),
106  omega_(ptf.omega_),
107  p0_(ptf.p0_)
108 {}
109 
110 
111 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
112 
114 (
115  const pointPatchFieldMapper& m
116 )
117 {
119 
120  p0_.autoMap(m);
121 }
122 
123 
125 (
126  const pointPatchField<vector>& ptf,
127  const labelList& addr
128 )
129 {
131  refCast<const oscillatingVelocityPointPatchVectorField>(ptf);
132 
134 
135  p0_.rmap(oVptf.p0_, addr);
136 }
137 
138 
140 {
141  if (this->updated())
142  {
143  return;
144  }
145 
146  const polyMesh& mesh = this->dimensionedInternalField().mesh()();
147  const Time& t = mesh.time();
148  const pointPatch& p = this->patch();
149 
151  (
152  (p0_ + amplitude_*sin(omega_*t.value()) - p.localPoints())
153  /t.deltaTValue()
154  );
155 
157 }
158 
159 
161 {
163  os.writeKeyword("amplitude")
164  << amplitude_ << token::END_STATEMENT << nl;
165  os.writeKeyword("omega")
166  << omega_ << token::END_STATEMENT << nl;
167  p0_.writeEntry("p0", os);
168  writeEntry("value", os);
169 }
170 
171 
172 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
173 
175 (
178 );
179 
180 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
181 
182 } // End namespace Foam
183 
184 // ************************************************************************* //
bool found(const word &, bool recursive=false, bool patternMatch=true) const
Search dictionary for given keyword.
Definition: dictionary.C:306
Basic pointPatch represents a set of points from the mesh.
Definition: pointPatch.H:56
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
virtual const vectorField & localPoints() const =0
Return mesh points.
dynamicFvMesh & mesh
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
Namespace for OpenFOAM.
virtual void autoMap(const pointPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
oscillatingVelocityPointPatchVectorField(const pointPatch &, const DimensionedField< vector, pointMesh > &)
Construct from patch and internal field.
makePointPatchTypeField(pointPatchVectorField, solidBodyMotionDisplacementPointPatchVectorField)
dictionary dict
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
Definition: Time.H:68
static const char nl
Definition: Ostream.H:260
scalar deltaTValue() const
Return time step value.
Definition: TimeState.H:100
volScalarField & p
Definition: createFields.H:51
const DimensionedField< vector, pointMesh > & dimensionedInternalField() const
Return dimensioned internal field reference.
Field< vector > vectorField
Specialisation of Field<T> for vector.
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if succesful.
Definition: doubleScalar.H:63
void writeEntry(const word &keyword, Ostream &os) const
Write the field as a dictionary entry.
Definition: Field.C:634
Macros for easy insertion into run-time selection tables.
virtual void rmap(const pointPatchField< Type > &, const labelList &)
Reverse map the given PointPatchField onto.
Pre-declare SubField and related Field type.
Definition: Field.H:57
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
Ostream & writeKeyword(const keyType &)
Write the keyword followed by an appropriate indentation.
Definition: Ostream.C:59
Foam::pointPatchFieldMapper.
Mesh consisting of general polyhedral cells.
Definition: polyMesh.H:74
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:452
virtual void autoMap(const pointPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
static const Vector zero
Definition: Vector.H:80
An Ostream is an abstract base class for all output systems (streams, files, token lists...
Definition: Ostream.H:53
const Time & time() const
Return time.
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
const pointPatch & patch() const
Return patch.
virtual void rmap(const pointPatchField< vector > &, const labelList &)
Reverse map the given pointPatchField onto this pointPatchField.
bool updated() const
Return true if the boundary condition has already been updated.
label size() const
Return size.
const Type & value() const
Return const reference to value.
virtual void write(Ostream &) const
Write.
dimensionedScalar sin(const dimensionedScalar &ds)