Loading...
Searching...
No Matches
facDiv.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 | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2016-2017 Wikki Ltd
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26\*---------------------------------------------------------------------------*/
27
28#include "facDiv.H"
29#include "faMesh.H"
30#include "facEdgeIntegrate.H"
31#include "faDivScheme.H"
32#include "faConvectionScheme.H"
33#include "transformField.H"
34
35// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
36
37namespace Foam
38{
39
40// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
41
42namespace fac
43{
44
45// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
46
47template<class Type>
48tmp<GeometricField<Type, faPatchField, areaMesh>>
49div
50(
51 const GeometricField<Type, faePatchField, edgeMesh>& ssf
52)
53{
54 const areaVectorField& n = ssf.mesh().faceAreaNormals();
55
56 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv =
57 fac::edgeIntegrate(ssf);
58
59 GeometricField<Type, faPatchField, areaMesh>& Div = tDiv.ref();
60
61 Div.primitiveFieldRef() =
62 transform(tensor::I - sqr(n), Div.internalField());
63 Div.correctBoundaryConditions();
64
65 return tDiv;
66}
67
68
69template<class Type>
70tmp<GeometricField<Type, faPatchField, areaMesh>>
71div
72(
73 const tmp<GeometricField<Type, faePatchField, edgeMesh>>& tssf
74)
75{
76 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv(fac::div(tssf()));
77 tssf.clear();
78 return tDiv;
79}
80
81
82template<class Type>
83tmp
84<
85 GeometricField
86 <
87 typename innerProduct<vector, Type>::type, faPatchField, areaMesh
88 >
89>
90div
91(
92 const GeometricField<Type, faPatchField, areaMesh>& vf,
93 const word& name
94)
95{
96 const areaVectorField& n = vf.mesh().faceAreaNormals();
97
98 tmp
99 <
100 GeometricField
101 <
102 typename innerProduct<vector, Type>::type,
103 faPatchField,
104 areaMesh
105 >
106 > tDiv
107 (
108 fa::divScheme<Type>::New
109 (
110 vf.mesh(), vf.mesh().divScheme(name)
111 ).ref().facDiv(vf)
112 );
113 GeometricField
114 <
115 typename innerProduct<vector, Type>::type,
116 faPatchField,
117 areaMesh
118 >& Div = tDiv.ref();
119
120 Div.primitiveFieldRef() =
121 transform(tensor::I - sqr(n), Div.internalField());
122
123 Div.correctBoundaryConditions();
124
125 return tDiv;
126}
127
128
129template<class Type>
130tmp
131<
132 GeometricField
133 <
134 typename innerProduct<vector, Type>::type, faPatchField, areaMesh
135 >
136>
137div
138(
139 const tmp<GeometricField<Type, faPatchField, areaMesh>>& tvvf,
140 const word& name
141)
142{
143 typedef typename innerProduct<vector, Type>::type DivType;
144 tmp<GeometricField<DivType, faPatchField, areaMesh>> tDiv
145 (
146 fac::div(tvvf(), name)
147 );
148 tvvf.clear();
149 return tDiv;
150}
151
152template<class Type>
160div
161(
162 const GeometricField<Type, faPatchField, areaMesh>& vf
163)
164{
165 return fac::div(vf, "div("+vf.name()+')');
166}
167
168
169template<class Type>
184 (
185 fac::div(tvvf())
186 );
187
188 tvvf.clear();
189 return tDiv;
190}
191
192
193template<class Type>
194tmp<GeometricField<Type, faPatchField, areaMesh>>
195div
196(
197 const edgeScalarField& flux,
198 const GeometricField<Type, faPatchField, areaMesh>& vf,
199 const word& name
200)
201{
202 const areaVectorField& n = vf.mesh().faceAreaNormals();
203
204 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
205 (
206 fa::convectionScheme<Type>::New
207 (
208 vf.mesh(),
209 flux,
210 vf.mesh().divScheme(name)
211 ).ref().facDiv(flux, vf)
212 );
213 GeometricField<Type, faPatchField, areaMesh>& Div = tDiv.ref();
214
215 Div.primitiveFieldRef() = transform(tensor::I - sqr(n), Div.internalField());
216 Div.correctBoundaryConditions();
217
218 return tDiv;
219
220}
221
222
223template<class Type>
224tmp<GeometricField<Type, faPatchField, areaMesh>>
225div
226(
227 const tmp<edgeScalarField>& tflux,
228 const GeometricField<Type, faPatchField, areaMesh>& vf,
229 const word& name
230)
231{
232 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
233 (
234 fac::div(tflux(), vf, name)
235 );
236 tflux.clear();
237 return tDiv;
238}
239
240
241template<class Type>
242tmp<GeometricField<Type, faPatchField, areaMesh>>
243div
244(
245 const edgeScalarField& flux,
246 const tmp<GeometricField<Type, faPatchField, areaMesh>>& tvf,
247 const word& name
248)
249{
250 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
251 (
252 fac::div(flux, tvf(), name)
253 );
254 tvf.clear();
255 return tDiv;
256}
257
258
259template<class Type>
260tmp<GeometricField<Type, faPatchField, areaMesh>>
261div
262(
263 const tmp<edgeScalarField>& tflux,
264 const tmp<GeometricField<Type, faPatchField, areaMesh>>& tvf,
265 const word& name
266)
267{
268 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
269 (
270 fac::div(tflux(), tvf(), name)
271 );
272 tflux.clear();
273 tvf.clear();
274 return tDiv;
275}
276
277
278template<class Type>
279tmp<GeometricField<Type, faPatchField, areaMesh>>
280div
281(
282 const edgeScalarField& flux,
283 const GeometricField<Type, faPatchField, areaMesh>& vf
284)
285{
286 return fac::div
287 (
288 flux, vf, "div("+flux.name()+','+vf.name()+')'
289 );
290}
291
292
293template<class Type>
294tmp<GeometricField<Type, faPatchField, areaMesh>>
295div
296(
297 const tmp<edgeScalarField>& tflux,
298 const GeometricField<Type, faPatchField, areaMesh>& vf
299)
300{
301 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
302 (
303 fac::div(tflux(), vf)
304 );
305 tflux.clear();
306 return tDiv;
307}
308
309
310template<class Type>
311tmp<GeometricField<Type, faPatchField, areaMesh>>
312div
313(
314 const edgeScalarField& flux,
315 const tmp<GeometricField<Type, faPatchField, areaMesh>>& tvf
316)
317{
318 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
319 (
320 fac::div(flux, tvf())
321 );
322 tvf.clear();
323 return tDiv;
324}
325
326
327template<class Type>
328tmp<GeometricField<Type, faPatchField, areaMesh>>
329div
330(
331 const tmp<edgeScalarField>& tflux,
332 const tmp<GeometricField<Type, faPatchField, areaMesh>>& tvf
333)
334{
335 tmp<GeometricField<Type, faPatchField, areaMesh>> tDiv
336 (
337 fac::div(tflux(), tvf())
338 );
339 tflux.clear();
340 tvf.clear();
341 return tDiv;
342}
343
344
345// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
346
347} // End namespace fac
348
349// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
350
351} // End namespace Foam
352
353// ************************************************************************* //
n
label n
Definition TABSMDCalcMethod2.H:31
Foam::DimensionedField::mesh
const Mesh & mesh() const noexcept
Return const reference to mesh.
Definition DimensionedField.H:688
Foam::GeometricField
Generic GeometricField class.
Definition GeometricField.H:75
Foam::GeometricField::primitiveFieldRef
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field values.
Definition GeometricField.C:881
Foam::GeometricField::internalField
const Internal & internalField() const noexcept
Return a const-reference to the dimensioned internal field.
Definition GeometricFieldI.H:27
Foam::GeometricField::correctBoundaryConditions
void correctBoundaryConditions()
Correct boundary field.
Definition GeometricField.C:1059
Foam::IOobject::name
const word & name() const noexcept
Return the object name.
Definition IOobjectI.H:205
Foam::Tensor< scalar >::I
static const Tensor I
Definition Tensor.H:81
Foam::areaMesh
Mesh data needed to do the Finite Area discretisation.
Definition areaFaMesh.H:50
Foam::faPatchField
faPatchField<Type> abstract base class. This class gives a fat-interface to all derived classes cover...
Definition faPatchField.H:321
Foam::fa::convectionScheme::New
static tmp< convectionScheme< Type > > New(const faMesh &mesh, const edgeScalarField &faceFlux, Istream &schemeData)
Return a pointer to a new convectionScheme created on freestore.
Definition faConvectionScheme.C:30
Foam::fa::divScheme::New
static tmp< divScheme< Type > > New(const faMesh &mesh, Istream &schemeData)
Return a pointer to a new divScheme created on freestore.
Definition faDivScheme.C:35
Foam::innerProduct::type
typeOfRank< typenamepTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) -2 >::type type
Definition products.H:155
Foam::tmp
A class for managing temporary objects.
Definition tmp.H:75
Foam::tmp::clear
void clear() const noexcept
If object pointer points to valid object: delete object and set pointer to nullptr.
Definition tmpI.H:289
Foam::tmp::ref
T & ref() const
Return non-const reference to the contents of a non-null managed pointer.
Definition tmpI.H:235
Foam::word
A class for handling words, derived from Foam::string.
Definition word.H:66
facDiv.H
Calculate the divergence of the given field.
facEdgeIntegrate.H
Edge integrate edgeField creating a areaField. Edge sum a edgeField creating a areaField.
Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition facDiv.C:43
Foam::fac::edgeIntegrate
tmp< GeometricField< Type, faPatchField, areaMesh > > edgeIntegrate(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition facEdgeIntegrate.C:40
Foam
Namespace for OpenFOAM.
Definition atmBoundaryLayer.C:27
Foam::sqr
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Definition dimensionedSymmTensor.C:44
Foam::transform
refinementData transform(const tensor &, const refinementData val)
No-op rotational transform for base types.
Definition refinementData.C:26
Foam::areaVectorField
GeometricField< vector, faPatchField, areaMesh > areaVectorField
Definition areaFieldsFwd.H:76
Foam::edgeScalarField
GeometricField< scalar, faePatchField, edgeMesh > edgeScalarField
Definition edgeFieldsFwd.H:56
Foam::name
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
Definition exprTraits.C:127
fac
Calculate the second temporal derivative.
transformField.H
Spatial transformation functions for primitive fields.