Loading...
Searching...
No Matches
fvcDdt.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) 2011-2017 OpenFOAM Foundation
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 "fvcDdt.H"
29#include "fvMesh.H"
30#include "ddtScheme.H"
31
32// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33
34namespace Foam
35{
36
37// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38
39namespace fvc
40{
41
42// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
43
44template<class Type>
45tmp<GeometricField<Type, fvPatchField, volMesh>>
46ddt
47(
48 const dimensioned<Type> dt,
49 const fvMesh& mesh
50)
51{
52 return fv::ddtScheme<Type>::New
53 (
54 mesh,
55 mesh.ddtScheme("ddt(" + dt.name() + ')')
56 ).ref().fvcDdt(dt);
57}
58
59
60template<class Type>
61tmp<GeometricField<Type, fvPatchField, volMesh>>
62ddt
63(
64 const GeometricField<Type, fvPatchField, volMesh>& vf
65)
66{
67 return fv::ddtScheme<Type>::New
68 (
69 vf.mesh(),
70 vf.mesh().ddtScheme("ddt(" + vf.name() + ')')
71 ).ref().fvcDdt(vf);
72}
73
74
75template<class Type>
76tmp<GeometricField<Type, fvPatchField, volMesh>>
77ddt
78(
79 const dimensionedScalar& rho,
80 const GeometricField<Type, fvPatchField, volMesh>& vf
81)
82{
83 return fv::ddtScheme<Type>::New
84 (
85 vf.mesh(),
86 vf.mesh().ddtScheme("ddt(" + rho.name() + ',' + vf.name() + ')')
87 ).ref().fvcDdt(rho, vf);
88}
89
90
91template<class Type>
92tmp<GeometricField<Type, fvPatchField, volMesh>>
93ddt
94(
95 const volScalarField& rho,
96 const GeometricField<Type, fvPatchField, volMesh>& vf
97)
98{
99 return fv::ddtScheme<Type>::New
100 (
101 vf.mesh(),
102 vf.mesh().ddtScheme("ddt(" + rho.name() + ',' + vf.name() + ')')
103 ).ref().fvcDdt(rho, vf);
104}
105
106
107template<class Type>
108tmp<GeometricField<Type, fvPatchField, volMesh>>
109ddt
110(
111 const volScalarField& alpha,
112 const volScalarField& rho,
113 const GeometricField<Type, fvPatchField, volMesh>& vf
114)
115{
116 return fv::ddtScheme<Type>::New
117 (
118 vf.mesh(),
119 vf.mesh().ddtScheme
120 (
121 "ddt("
122 + alpha.name() + ','
123 + rho.name() + ','
124 + vf.name() + ')'
125 )
126 ).ref().fvcDdt(alpha, rho, vf);
127}
128
129
130template<class Type>
131tmp<GeometricField<Type, fvsPatchField, surfaceMesh>>
132ddt
133(
134 const GeometricField<Type, fvsPatchField, surfaceMesh>& sf
135)
136{
137 return fv::ddtScheme<Type>::New
138 (
139 sf.mesh(),
140 sf.mesh().ddtScheme("ddt(" + sf.name() + ')')
141 ).ref().fvcDdt(sf);
142}
143
144
145template<class Type>
146tmp<GeometricField<Type, fvPatchField, volMesh>>
147ddt
148(
149 const Foam::one,
150 const GeometricField<Type, fvPatchField, volMesh>& vf
151)
152{
153 return ddt(vf);
154}
155
156
157template<class Type>
158tmp<GeometricField<Type, fvPatchField, volMesh>>
159ddt
160(
161 const GeometricField<Type, fvPatchField, volMesh>& vf,
162 const Foam::one
163)
176{
177 return fv::ddtScheme<Type>::New
178 (
179 U.mesh(),
180 U.mesh().ddtScheme("ddt(" + U.name() + ')')
181 ).ref().fvcDdtUfCorr(U, Uf);
182}
183
184
185template<class Type>
186tmp<GeometricField<typename flux<Type>::type, fvsPatchField, surfaceMesh>>
187ddtCorr
188(
189 const GeometricField<Type, fvPatchField, volMesh>& U,
190 const GeometricField
191 <
192 typename flux<Type>::type,
193 fvsPatchField,
194 surfaceMesh
195 >& phi
196)
197{
198 return fv::ddtScheme<Type>::New
199 (
200 U.mesh(),
201 U.mesh().ddtScheme("ddt(" + U.name() + ')')
202 ).ref().fvcDdtPhiCorr(U, phi);
203}
204
205
206template<class Type>
207tmp<GeometricField<typename flux<Type>::type, fvsPatchField, surfaceMesh>>
208ddtCorr
209(
210 const GeometricField<Type, fvPatchField, volMesh>& U,
211 const GeometricField
212 <
213 typename flux<Type>::type,
214 fvsPatchField,
215 surfaceMesh
216 >& phi,
217 const autoPtr<GeometricField<Type, fvsPatchField, surfaceMesh>>& Uf
218)
219{
220 if (U.mesh().dynamic())
221 {
222 return ddtCorr(U, Uf());
223 }
224 else
225 {
226 return ddtCorr(U, phi);
227 }
228}
229
230
231template<class Type>
232tmp<GeometricField<typename flux<Type>::type, fvsPatchField, surfaceMesh>>
233ddtCorr
234(
235 const volScalarField& rho,
236 const GeometricField<Type, fvPatchField, volMesh>& U,
237 const GeometricField<Type, fvsPatchField, surfaceMesh>& Uf
238)
239{
240 return fv::ddtScheme<Type>::New
241 (
242 U.mesh(),
243 U.mesh().ddtScheme("ddt(" + U.name() + ')')
244 ).ref().fvcDdtUfCorr(rho, U, Uf);
245}
246
247
248template<class Type>
249tmp<GeometricField<typename flux<Type>::type, fvsPatchField, surfaceMesh>>
250ddtCorr
251(
252 const volScalarField& rho,
253 const GeometricField<Type, fvPatchField, volMesh>& U,
254 const GeometricField
255 <
256 typename flux<Type>::type,
257 fvsPatchField,
258 surfaceMesh
259 >& phi
260)
261{
262 return fv::ddtScheme<Type>::New
263 (
264 U.mesh(),
265 U.mesh().ddtScheme("ddt(" + rho.name() + ',' + U.name() + ')')
266 ).ref().fvcDdtPhiCorr(rho, U, phi);
267}
268
269
270template<class Type>
271tmp<GeometricField<typename flux<Type>::type, fvsPatchField, surfaceMesh>>
272ddtCorr
273(
274 const volScalarField& rho,
275 const GeometricField<Type, fvPatchField, volMesh>& U,
276 const GeometricField
277 <
278 typename flux<Type>::type,
279 fvsPatchField,
280 surfaceMesh
281 >& phi,
282 const autoPtr<GeometricField<Type, fvsPatchField, surfaceMesh>>& Uf
283)
284{
285 if (U.mesh().dynamic())
286 {
287 return ddtCorr(rho, U, Uf());
288 }
289 else
290 {
291 return ddtCorr(rho, U, phi);
292 }
293}
294
295
296// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
297
298} // End namespace fvc
299
300// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
301
302} // End namespace Foam
303
304// ************************************************************************* //
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::IOobject::name
const word & name() const noexcept
Return the object name.
Definition IOobjectI.H:205
Foam::autoPtr
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
Definition autoPtr.H:65
Foam::dimensioned
Generic dimensioned Type class.
Definition dimensionedType.H:70
Foam::dimensioned::name
const word & name() const noexcept
Return const reference to name.
Definition dimensionedType.H:352
Foam::fvMesh
Mesh data needed to do the Finite Volume discretisation.
Definition fvMesh.H:85
Foam::fv::ddtScheme::New
static tmp< ddtScheme< Type > > New(const fvMesh &mesh, Istream &schemeData)
Return a pointer to a new ddtScheme created on freestore.
Definition ddtScheme.C:43
Foam::fvsPatchField
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
Definition fvsPatchField.H:298
Foam::innerProduct< vector, Type >::type
typeOfRank< typenamepTraits< vector >::cmptType, direction(pTraits< vector >::rank)+direction(pTraits< Type >::rank) -2 >::type type
Definition products.H:155
Foam::one
A class representing the concept of 1 (one) that can be used to avoid manipulating objects known to b...
Definition one.H:57
Foam::surfaceMesh
Mesh data needed to do the Finite Volume discretisation.
Definition surfaceMesh.H:47
Foam::tmp
A class for managing temporary objects.
Definition tmp.H:75
Foam::tmp::ref
T & ref() const
Return non-const reference to the contents of a non-null managed pointer.
Definition tmpI.H:235
U
U
Definition pEqn.H:72
mesh
dynamicFvMesh & mesh
Definition createDynamicFvMesh.H:6
Uf
autoPtr< surfaceVectorField > Uf
Definition createUfIfPresent.H:27
fvcDdt.H
Calculate the first temporal derivative.
Foam::fvc
Namespace of functions to calculate explicit derivatives.
Foam::fvc::ddt
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition fvcDdt.C:40
Foam::fvc::ddtCorr
tmp< GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > > ddtCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Definition fvcDdt.C:165
Foam
Namespace for OpenFOAM.
Definition atmBoundaryLayer.C:27
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition volFieldsFwd.H:72
Foam::dimensionedScalar
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Definition dimensionedScalarFwd.H:35
alpha
volScalarField & alpha
Definition readThermalProperties.H:212