facEdgeIntegrate.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 "facEdgeIntegrate.H"
29#include "faMesh.H"
31
32// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33
34namespace Foam
35{
36
37// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38
39namespace fac
40{
41
42// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43
44template<class Type>
45tmp<GeometricField<Type, faPatchField, areaMesh>>
47(
49)
50{
51 const faMesh& mesh = ssf.mesh();
52
54 (
56 (
58 (
59 "edgeIntegrate("+ssf.name()+')',
60 ssf.instance(),
61 ssf.db()
62 ),
63 mesh,
66 )
67 );
69
70
71 const labelUList& owner = mesh.owner();
72 const labelUList& neighbour = mesh.neighbour();
73
74 forAll(owner, faceI)
75 {
76 vf[owner[faceI]] += ssf[faceI];
77 vf[neighbour[faceI]] -= ssf[faceI];
78 }
79
80 forAll(mesh.boundary(), patchI)
81 {
82 const labelUList& pEdgeFaces =
83 mesh.boundary()[patchI].edgeFaces();
84
85 const faePatchField<Type>& pssf = ssf.boundaryField()[patchI];
86
87 forAll(mesh.boundary()[patchI], faceI)
88 {
89 vf[pEdgeFaces[faceI]] += pssf[faceI];
90 }
91 }
92
93 vf.primitiveFieldRef() /= mesh.S();
95
96 return tvf;
97}
98
99
100template<class Type>
103(
105)
106{
108 (
109 fac::edgeIntegrate(tssf())
110 );
111 tssf.clear();
112 return tvf;
113}
114
115
116template<class Type>
119(
121)
122{
123 const faMesh& mesh = ssf.mesh();
124
126 (
128 (
130 (
131 "edgeSum("+ssf.name()+')',
132 ssf.instance(),
133 ssf.db()
134 ),
135 mesh,
138 )
139 );
141
142 const labelUList& owner = mesh.owner();
143 const labelUList& neighbour = mesh.neighbour();
144
145 forAll(owner, faceI)
146 {
147 vf[owner[faceI]] += ssf[faceI];
148 vf[neighbour[faceI]] += ssf[faceI];
149 }
150
151 forAll(mesh.boundary(), patchI)
152 {
153 const labelUList& pEdgeFaces =
154 mesh.boundary()[patchI].edgeFaces();
155
156 const faePatchField<Type>& pssf = ssf.boundaryField()[patchI];
157
158 forAll(mesh.boundary()[patchI], faceI)
159 {
160 vf[pEdgeFaces[faceI]] += pssf[faceI];
161 }
162 }
163
165
166 return tvf;
167}
168
169
170template<class Type>
172(
174)
175{
177 edgeSum(tssf());
178 tssf.clear();
179 return tvf;
180}
181
182
183// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
184
185} // End namespace fac
186
187// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
188
189} // End namespace Foam
190
191// ************************************************************************* //
const dimensionSet & dimensions() const
Return dimensions.
const Mesh & mesh() const
Return mesh.
Generic GeometricField class.
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
void correctBoundaryConditions()
Correct boundary field.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
Definition: IOobject.H:170
const word & name() const noexcept
Return the object name.
Definition: IOobjectI.H:65
const objectRegistry & db() const noexcept
Return the local objectRegistry.
Definition: IOobject.C:500
const fileName & instance() const noexcept
Read access to instance path component.
Definition: IOobjectI.H:196
Generic dimensioned Type class.
Finite area mesh (used for 2-D non-Euclidian finite area method) defined using a patch of faces on a ...
Definition: faMesh.H:100
faePatchField<Type> abstract base class. This class gives a fat-interface to all derived classes cove...
Definition: faePatchField.H:83
A class for managing temporary objects.
Definition: tmp.H:65
T & ref() const
Definition: tmpI.H:227
Author Zeljko Tukovic, FMENA Hrvoje Jasak, Wikki Ltd.
dynamicFvMesh & mesh
Edge integrate edgeField creating a areaField. Edge sum a edgeField creating a areaField.
tmp< GeometricField< Type, faPatchField, areaMesh > > edgeSum(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
tmp< GeometricField< Type, faPatchField, areaMesh > > edgeIntegrate(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Namespace for OpenFOAM.
const dimensionSet dimArea(sqr(dimLength))
Definition: dimensionSets.H:59
static constexpr const zero Zero
Global zero (0)
Definition: zero.H:131
Calculate the second temporal derivative.
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333