famNDiv.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 -------------------------------------------------------------------------------
10 License
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 "famNDiv.H"
29 #include "faMesh.H"
30 #include "faMatrix.H"
31 #include "faConvectionScheme.H"
32 
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34 
35 namespace Foam
36 {
37 
38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
39 
40 namespace fam
41 {
42 
43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
44 
45 template<class Type>
46 tmp<faMatrix<Type>>
47 ndiv
48 (
49  const edgeScalarField& flux,
51  const word& name
52 )
53 {
55  (
56  vf.mesh(),
57  flux,
58  vf.mesh().divScheme(name)
59  ).ref().famDiv(flux, vf);//TODO calculate normal
60 }
61 
62 
63 template<class Type>
65 ndiv
66 (
67  const tmp<edgeScalarField>& tflux,
69  const word& name
70 )
71 {
72  tmp<faMatrix<Type>> Div(fam::ndiv(tflux(), vf, name));
73  tflux.clear();
74 
75  return Div;
76 }
77 
78 
79 template<class Type>
81 ndiv
82 (
83  const edgeScalarField& flux,
85 )
86 {
87  return fam::ndiv(flux, vf, "div("+flux.name()+','+vf.name()+')');
88 }
89 
90 
91 template<class Type>
93 ndiv
94 (
95  const tmp<edgeScalarField>& tflux,
97 )
98 {
99  tmp<faMatrix<Type>> Div(fam::ndiv(tflux(), vf));
100  tflux.clear();
101 
102  return Div;
103 }
104 
105 
106 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
107 
108 } // End namespace fam
109 
110 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
111 
112 } // End namespace Foam
113 
114 // ************************************************************************* //
Foam::fvc::flux
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Foam::word
A class for handling words, derived from Foam::string.
Definition: word.H:65
Foam::tmp::clear
void clear() const noexcept
Definition: tmpI.H:287
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:61
faConvectionScheme.H
faMesh.H
ref
rDeltaT ref()
faMatrix.H
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
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:37
Foam::fam::ndiv
tmp< faMatrix< Type > > ndiv(const edgeScalarField &flux, GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
Definition: famNDiv.C:48
Foam::name
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
Definition: exprTraits.C:59
Foam::flux
Definition: vector.H:56
fam
Calculate the matrix for the second temporal derivative.
Foam::GeometricField
Generic GeometricField class.
Definition: areaFieldsFwd.H:53
famNDiv.H
Calculate the matrix for the divergence of the given field and flux.