steadyStateDdtScheme.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) 2011-2015 OpenFOAM Foundation
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 "steadyStateDdtScheme.H"
27 #include "fvcDiv.H"
28 #include "fvMatrices.H"
29 
30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
31 
32 namespace Foam
33 {
34 
35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
36 
37 namespace fv
38 {
39 
40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
41 
42 template<class Type>
43 tmp<GeometricField<Type, fvPatchField, volMesh> >
44 steadyStateDdtScheme<Type>::fvcDdt
45 (
46  const dimensioned<Type>& dt
47 )
48 {
49  return tmp<GeometricField<Type, fvPatchField, volMesh> >
50  (
51  new GeometricField<Type, fvPatchField, volMesh>
52  (
53  IOobject
54  (
55  "ddt("+dt.name()+')',
56  mesh().time().timeName(),
57  mesh()
58  ),
59  mesh(),
60  dimensioned<Type>
61  (
62  "0",
63  dt.dimensions()/dimTime,
64  pTraits<Type>::zero
65  )
66  )
67  );
68 }
69 
70 
71 template<class Type>
72 tmp<GeometricField<Type, fvPatchField, volMesh> >
73 steadyStateDdtScheme<Type>::fvcDdt
74 (
75  const GeometricField<Type, fvPatchField, volMesh>& vf
76 )
77 {
78  return tmp<GeometricField<Type, fvPatchField, volMesh> >
79  (
80  new GeometricField<Type, fvPatchField, volMesh>
81  (
82  IOobject
83  (
84  "ddt("+vf.name()+')',
85  mesh().time().timeName(),
86  mesh()
87  ),
88  mesh(),
89  dimensioned<Type>
90  (
91  "0",
92  vf.dimensions()/dimTime,
93  pTraits<Type>::zero
94  )
95  )
96  );
97 }
98 
99 
100 template<class Type>
101 tmp<GeometricField<Type, fvPatchField, volMesh> >
102 steadyStateDdtScheme<Type>::fvcDdt
103 (
104  const dimensionedScalar& rho,
105  const GeometricField<Type, fvPatchField, volMesh>& vf
106 )
107 {
108  return tmp<GeometricField<Type, fvPatchField, volMesh> >
109  (
110  new GeometricField<Type, fvPatchField, volMesh>
111  (
112  IOobject
113  (
114  "ddt("+rho.name()+','+vf.name()+')',
115  mesh().time().timeName(),
116  mesh()
117  ),
118  mesh(),
119  dimensioned<Type>
120  (
121  "0",
122  rho.dimensions()*vf.dimensions()/dimTime,
123  pTraits<Type>::zero
124  )
125  )
126  );
127 }
128 
129 
130 template<class Type>
131 tmp<GeometricField<Type, fvPatchField, volMesh> >
132 steadyStateDdtScheme<Type>::fvcDdt
133 (
134  const volScalarField& rho,
135  const GeometricField<Type, fvPatchField, volMesh>& vf
136 )
137 {
138  return tmp<GeometricField<Type, fvPatchField, volMesh> >
139  (
140  new GeometricField<Type, fvPatchField, volMesh>
141  (
142  IOobject
143  (
144  "ddt("+rho.name()+','+vf.name()+')',
145  mesh().time().timeName(),
146  mesh()
147  ),
148  mesh(),
149  dimensioned<Type>
150  (
151  "0",
152  rho.dimensions()*vf.dimensions()/dimTime,
153  pTraits<Type>::zero
154  )
155  )
156  );
157 }
158 
159 
160 template<class Type>
161 tmp<GeometricField<Type, fvPatchField, volMesh> >
162 steadyStateDdtScheme<Type>::fvcDdt
163 (
164  const volScalarField& alpha,
165  const volScalarField& rho,
166  const GeometricField<Type, fvPatchField, volMesh>& vf
167 )
168 {
169  return tmp<GeometricField<Type, fvPatchField, volMesh> >
170  (
171  new GeometricField<Type, fvPatchField, volMesh>
172  (
173  IOobject
174  (
175  "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
176  mesh().time().timeName(),
177  mesh()
178  ),
179  mesh(),
180  dimensioned<Type>
181  (
182  "0",
183  rho.dimensions()*vf.dimensions()/dimTime,
184  pTraits<Type>::zero
185  )
186  )
187  );
188 }
189 
190 
191 template<class Type>
192 tmp<fvMatrix<Type> >
193 steadyStateDdtScheme<Type>::fvmDdt
194 (
195  const GeometricField<Type, fvPatchField, volMesh>& vf
196 )
197 {
198  tmp<fvMatrix<Type> > tfvm
199  (
200  new fvMatrix<Type>
201  (
202  vf,
203  vf.dimensions()*dimVol/dimTime
204  )
205  );
206 
207  return tfvm;
208 }
209 
210 
211 template<class Type>
212 tmp<fvMatrix<Type> >
213 steadyStateDdtScheme<Type>::fvmDdt
214 (
215  const dimensionedScalar& rho,
216  const GeometricField<Type, fvPatchField, volMesh>& vf
217 )
218 {
219  tmp<fvMatrix<Type> > tfvm
220  (
221  new fvMatrix<Type>
222  (
223  vf,
224  rho.dimensions()*vf.dimensions()*dimVol/dimTime
225  )
226  );
227 
228  return tfvm;
229 }
230 
231 
232 template<class Type>
233 tmp<fvMatrix<Type> >
234 steadyStateDdtScheme<Type>::fvmDdt
235 (
236  const volScalarField& rho,
237  const GeometricField<Type, fvPatchField, volMesh>& vf
238 )
239 {
240  tmp<fvMatrix<Type> > tfvm
241  (
242  new fvMatrix<Type>
243  (
244  vf,
245  rho.dimensions()*vf.dimensions()*dimVol/dimTime
246  )
247  );
248 
249  return tfvm;
250 }
251 
252 
253 template<class Type>
254 tmp<fvMatrix<Type> >
255 steadyStateDdtScheme<Type>::fvmDdt
256 (
257  const volScalarField& alpha,
258  const volScalarField& rho,
259  const GeometricField<Type, fvPatchField, volMesh>& vf
260 )
261 {
262  tmp<fvMatrix<Type> > tfvm
263  (
264  new fvMatrix<Type>
265  (
266  vf,
267  alpha.dimensions()*rho.dimensions()*vf.dimensions()*dimVol/dimTime
268  )
269  );
270 
271  return tfvm;
272 }
273 
274 
275 template<class Type>
276 tmp<typename steadyStateDdtScheme<Type>::fluxFieldType>
277 steadyStateDdtScheme<Type>::fvcDdtUfCorr
278 (
279  const GeometricField<Type, fvPatchField, volMesh>& U,
280  const GeometricField<Type, fvsPatchField, surfaceMesh>& Uf
281 )
282 {
283  return tmp<fluxFieldType>
284  (
285  new fluxFieldType
286  (
287  IOobject
288  (
289  "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
290  mesh().time().timeName(),
291  mesh()
292  ),
293  mesh(),
294  dimensioned<typename flux<Type>::type>
295  (
296  "0",
297  Uf.dimensions()*dimArea/dimTime,
298  pTraits<typename flux<Type>::type>::zero
299  )
300  )
301  );
302 }
303 
304 
305 template<class Type>
306 tmp<typename steadyStateDdtScheme<Type>::fluxFieldType>
307 steadyStateDdtScheme<Type>::fvcDdtPhiCorr
308 (
309  const GeometricField<Type, fvPatchField, volMesh>& U,
310  const fluxFieldType& phi
311 )
312 {
313  return tmp<fluxFieldType>
314  (
315  new fluxFieldType
316  (
317  IOobject
318  (
319  "ddtCorr(" + U.name() + ',' + phi.name() + ')',
320  mesh().time().timeName(),
321  mesh()
322  ),
323  mesh(),
324  dimensioned<typename flux<Type>::type>
325  (
326  "0",
327  phi.dimensions()/dimTime,
328  pTraits<typename flux<Type>::type>::zero
329  )
330  )
331  );
332 }
333 
334 
335 template<class Type>
336 tmp<typename steadyStateDdtScheme<Type>::fluxFieldType>
337 steadyStateDdtScheme<Type>::fvcDdtUfCorr
338 (
339  const volScalarField& rho,
340  const GeometricField<Type, fvPatchField, volMesh>& U,
341  const GeometricField<Type, fvsPatchField, surfaceMesh>& Uf
342 )
343 {
344  return tmp<fluxFieldType>
345  (
346  new fluxFieldType
347  (
348  IOobject
349  (
350  "ddtCorr("
351  + rho.name()
352  + ',' + U.name() + ',' + Uf.name() + ')',
353  mesh().time().timeName(),
354  mesh()
355  ),
356  mesh(),
357  dimensioned<typename flux<Type>::type>
358  (
359  "0",
360  Uf.dimensions()*dimArea/dimTime,
361  pTraits<typename flux<Type>::type>::zero
362  )
363  )
364  );
365 }
366 
367 
368 template<class Type>
369 tmp<typename steadyStateDdtScheme<Type>::fluxFieldType>
370 steadyStateDdtScheme<Type>::fvcDdtPhiCorr
371 (
372  const volScalarField& rho,
373  const GeometricField<Type, fvPatchField, volMesh>& U,
374  const fluxFieldType& phi
375 )
376 {
377  return tmp<fluxFieldType>
378  (
379  new fluxFieldType
380  (
381  IOobject
382  (
383  "ddtCorr("
384  + rho.name()
385  + ',' + U.name() + ',' + phi.name() + ')',
386  mesh().time().timeName(),
387  mesh()
388  ),
389  mesh(),
390  dimensioned<typename flux<Type>::type>
391  (
392  "0",
393  phi.dimensions()/dimTime,
394  pTraits<typename flux<Type>::type>::zero
395  )
396  )
397  );
398 }
399 
400 
401 template<class Type>
402 tmp<surfaceScalarField> steadyStateDdtScheme<Type>::meshPhi
403 (
404  const GeometricField<Type, fvPatchField, volMesh>& vf
405 )
406 {
407  return tmp<surfaceScalarField>
408  (
409  new surfaceScalarField
410  (
411  IOobject
412  (
413  "meshPhi",
414  mesh().time().timeName(),
415  mesh(),
416  IOobject::NO_READ,
417  IOobject::NO_WRITE,
418  false
419  ),
420  mesh(),
421  dimensionedScalar("0", dimVolume/dimTime, 0.0)
422  )
423  );
424 }
425 
426 
427 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
428 
429 } // End namespace fv
430 
431 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
432 
433 } // End namespace Foam
434 
435 // ************************************************************************* //
Foam::fv::steadyStateDdtScheme::fvmDdt
tmp< fvMatrix< Type > > fvmDdt(const GeometricField< Type, fvPatchField, volMesh > &)
Definition: steadyStateDdtScheme.C:194
Foam::dimTime
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:51
Foam::dimVol
const dimensionSet dimVol(dimVolume)
Definition: dimensionSets.H:59
Foam::dimArea
const dimensionSet dimArea(sqr(dimLength))
Definition: dimensionSets.H:57
Foam::fvMatrix
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
Definition: fvPatchField.H:68
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
Foam
Namespace for OpenFOAM.
Definition: combustionModel.C:30
Foam::dimensioned::dimensions
const dimensionSet & dimensions() const
Return const reference to dimensions.
Definition: dimensionedType.C:250
Foam::dimensioned
Generic dimensioned Type class.
Definition: dimensionedScalarFwd.H:41
Foam::IOobject
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:91
Foam::dimensioned::name
const word & name() const
Return const reference to name.
Definition: dimensionedType.C:237
Foam::fv::steadyStateDdtScheme::fvcDdtPhiCorr
tmp< fluxFieldType > fvcDdtPhiCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const fluxFieldType &phi)
Definition: steadyStateDdtScheme.C:308
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:52
Foam::zero
A class representing the concept of 0 used to avoid unnecessary manipulations for objects that are kn...
Definition: zero.H:47
Foam::fv::steadyStateDdtScheme::fvcDdtUfCorr
tmp< fluxFieldType > fvcDdtUfCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Definition: steadyStateDdtScheme.C:278
Foam::DimensionedField::dimensions
const dimensionSet & dimensions() const
Return dimensions.
Definition: DimensionedFieldI.H:46
Foam::IOobject::name
const word & name() const
Return name.
Definition: IOobject.H:260
Foam::fv::steadyStateDdtScheme::fvcDdt
tmp< GeometricField< Type, fvPatchField, volMesh > > fvcDdt(const dimensioned< Type > &)
Definition: steadyStateDdtScheme.C:45
Foam::GeometricField
Generic GeometricField class.
Definition: surfaceFieldsFwd.H:52
fvcDiv.H
Calculate the divergence of the given field.
Foam::pTraits
Traits class for primitives.
Definition: pTraits.H:50
Foam::dimVolume
const dimensionSet dimVolume(pow3(dimLength))
Definition: dimensionSets.H:58
fv
labelList fv(nPoints)
Foam::fv::steadyStateDdtScheme::meshPhi
tmp< surfaceScalarField > meshPhi(const GeometricField< Type, fvPatchField, volMesh > &)
Definition: steadyStateDdtScheme.C:403
Foam::dimensionedScalar
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Definition: dimensionedScalarFwd.H:41
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:118
timeName
word timeName
Definition: getTimeIndex.H:3
fvMatrices.H
A special matrix type and solver, designed for finite volume solutions of scalar equations.
Generated by   doxygen 1.8.11