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
finiteVolume
fields
fvPatchFields
basic
basicSymmetry
basicSymmetryFvPatchField.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 "
basicSymmetryFvPatchField.H
"
27
#include <
OpenFOAM/symmTransformField.H
>
28
29
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
30
31
namespace
Foam
32
{
33
34
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
35
36
template
<
class
Type>
37
basicSymmetryFvPatchField<Type>::basicSymmetryFvPatchField
38
(
39
const
fvPatch
&
p
,
40
const
DimensionedField<Type, volMesh>
& iF
41
)
42
:
43
transformFvPatchField<Type>
(
p
, iF)
44
{}
45
46
47
template
<
class
Type>
48
basicSymmetryFvPatchField<Type>::basicSymmetryFvPatchField
49
(
50
const
basicSymmetryFvPatchField<Type>
& ptf,
51
const
fvPatch
& p,
52
const
DimensionedField<Type, volMesh>
& iF,
53
const
fvPatchFieldMapper
& mapper
54
)
55
:
56
transformFvPatchField<Type>
(ptf,
p
, iF, mapper)
57
{}
58
59
60
template
<
class
Type>
61
basicSymmetryFvPatchField<Type>::basicSymmetryFvPatchField
62
(
63
const
fvPatch
& p,
64
const
DimensionedField<Type, volMesh>
& iF,
65
const
dictionary
& dict
66
)
67
:
68
transformFvPatchField<Type>
(
p
, iF, dict)
69
{
70
this->evaluate();
71
}
72
73
74
template
<
class
Type>
75
basicSymmetryFvPatchField<Type>::basicSymmetryFvPatchField
76
(
77
const
basicSymmetryFvPatchField<Type>
& ptf
78
)
79
:
80
transformFvPatchField<Type>
(ptf)
81
{}
82
83
84
template
<
class
Type>
85
basicSymmetryFvPatchField<Type>::basicSymmetryFvPatchField
86
(
87
const
basicSymmetryFvPatchField<Type>
& ptf,
88
const
DimensionedField<Type, volMesh>
& iF
89
)
90
:
91
transformFvPatchField<Type>
(ptf, iF)
92
{}
93
94
95
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
96
97
// return gradient at boundary
98
template
<
class
Type>
99
tmp<Field<Type>
>
basicSymmetryFvPatchField<Type>::snGrad
()
const
100
{
101
vectorField
nHat = this->patch().nf();
102
return
103
(
104
transform
(
I
- 2.0*
sqr
(nHat), this->patchInternalField())
105
- this->patchInternalField()
106
)*(this->patch().deltaCoeffs()/2.0);
107
}
108
109
110
// Evaluate the field on the patch
111
template
<
class
Type>
112
void
basicSymmetryFvPatchField<Type>::evaluate
(
const
Pstream::commsTypes
)
113
{
114
if
(!this->updated())
115
{
116
this->updateCoeffs();
117
}
118
119
vectorField
nHat = this->patch().nf();
120
Field<Type>::operator
=
121
(
122
(
123
this->patchInternalField()
124
+
transform
(
I
- 2.0*
sqr
(nHat), this->patchInternalField())
125
)/2.0
126
);
127
128
transformFvPatchField<Type>::evaluate
();
129
}
130
131
132
// Return defining fields
133
template
<
class
Type>
134
tmp<Field<Type>
>
basicSymmetryFvPatchField<Type>::snGradTransformDiag
()
const
135
{
136
vectorField
nHat = this->patch().nf();
137
vectorField
diag
(nHat.
size
());
138
139
diag
.replace(
vector::X
,
mag
(nHat.
component
(
vector::X
)));
140
diag
.replace(
vector::Y
,
mag
(nHat.
component
(
vector::Y
)));
141
diag
.replace(
vector::Z
,
mag
(nHat.
component
(
vector::Z
)));
142
143
return
transformFieldMask<Type>(pow<vector, pTraits<Type>::rank>(
diag
));
144
}
145
146
147
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
148
149
}
// End namespace Foam
150
151
// ************************ vim: set sw=4 sts=4 et: ************************ //