localEulerDdtScheme.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 "localEulerDdtScheme.H"
27 #include "surfaceInterpolate.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 const volScalarField& localEulerDdtScheme<Type>::localRDeltaT() const
44 {
45  return localEulerDdt::localRDeltaT(mesh());
46 }
47 
48 
49 template<class Type>
50 tmp<GeometricField<Type, fvPatchField, volMesh> >
51 localEulerDdtScheme<Type>::fvcDdt
52 (
53  const dimensioned<Type>& dt
54 )
55 {
56  const volScalarField& rDeltaT = localRDeltaT();
57 
58  IOobject ddtIOobject
59  (
60  "ddt(" + dt.name() + ')',
61  mesh().time().timeName(),
62  mesh()
63  );
64 
65  if (mesh().moving())
66  {
67  tmp<GeometricField<Type, fvPatchField, volMesh> > tdtdt
68  (
69  new GeometricField<Type, fvPatchField, volMesh>
70  (
71  ddtIOobject,
72  mesh(),
73  dimensioned<Type>
74  (
75  "0",
76  dt.dimensions()/dimTime,
77  pTraits<Type>::zero
78  )
79  )
80  );
81 
82  tdtdt().internalField() =
83  rDeltaT.internalField()*dt.value()
84  *(1.0 - mesh().Vsc0()/mesh().Vsc());
85 
86  return tdtdt;
87  }
88  else
89  {
90  return tmp<GeometricField<Type, fvPatchField, volMesh> >
91  (
92  new GeometricField<Type, fvPatchField, volMesh>
93  (
94  ddtIOobject,
95  mesh(),
96  dimensioned<Type>
97  (
98  "0",
99  dt.dimensions()/dimTime,
100  pTraits<Type>::zero
101  ),
102  calculatedFvPatchField<Type>::typeName
103  )
104  );
105  }
106 }
107 
108 
109 template<class Type>
110 tmp<GeometricField<Type, fvPatchField, volMesh> >
111 localEulerDdtScheme<Type>::fvcDdt
112 (
113  const GeometricField<Type, fvPatchField, volMesh>& vf
114 )
115 {
116  const volScalarField& rDeltaT = localRDeltaT();
117 
118  IOobject ddtIOobject
119  (
120  "ddt(" + vf.name() + ')',
121  mesh().time().timeName(),
122  mesh()
123  );
124 
125  if (mesh().moving())
126  {
127  return tmp<GeometricField<Type, fvPatchField, volMesh> >
128  (
129  new GeometricField<Type, fvPatchField, volMesh>
130  (
131  ddtIOobject,
132  mesh(),
133  rDeltaT.dimensions()*vf.dimensions(),
134  rDeltaT.internalField()*
135  (
136  vf.internalField()
137  - vf.oldTime().internalField()*mesh().Vsc0()/mesh().Vsc()
138  ),
139  rDeltaT.boundaryField()*
140  (
141  vf.boundaryField() - vf.oldTime().boundaryField()
142  )
143  )
144  );
145  }
146  else
147  {
148  return tmp<GeometricField<Type, fvPatchField, volMesh> >
149  (
150  new GeometricField<Type, fvPatchField, volMesh>
151  (
152  ddtIOobject,
153  rDeltaT*(vf - vf.oldTime())
154  )
155  );
156  }
157 }
158 
159 
160 template<class Type>
161 tmp<GeometricField<Type, fvPatchField, volMesh> >
162 localEulerDdtScheme<Type>::fvcDdt
163 (
164  const dimensionedScalar& rho,
165  const GeometricField<Type, fvPatchField, volMesh>& vf
166 )
167 {
168  const volScalarField& rDeltaT = localRDeltaT();
169 
170  IOobject ddtIOobject
171  (
172  "ddt(" + rho.name() + ',' + vf.name() + ')',
173  mesh().time().timeName(),
174  mesh()
175  );
176 
177  if (mesh().moving())
178  {
179  return tmp<GeometricField<Type, fvPatchField, volMesh> >
180  (
181  new GeometricField<Type, fvPatchField, volMesh>
182  (
183  ddtIOobject,
184  mesh(),
185  rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
186  rDeltaT.internalField()*rho.value()*
187  (
188  vf.internalField()
189  - vf.oldTime().internalField()*mesh().Vsc0()/mesh().Vsc()
190  ),
191  rDeltaT.boundaryField()*rho.value()*
192  (
193  vf.boundaryField() - vf.oldTime().boundaryField()
194  )
195  )
196  );
197  }
198  else
199  {
200  return tmp<GeometricField<Type, fvPatchField, volMesh> >
201  (
202  new GeometricField<Type, fvPatchField, volMesh>
203  (
204  ddtIOobject,
205  rDeltaT*rho*(vf - vf.oldTime())
206  )
207  );
208  }
209 }
210 
211 
212 template<class Type>
213 tmp<GeometricField<Type, fvPatchField, volMesh> >
214 localEulerDdtScheme<Type>::fvcDdt
215 (
216  const volScalarField& rho,
217  const GeometricField<Type, fvPatchField, volMesh>& vf
218 )
219 {
220  const volScalarField& rDeltaT = localRDeltaT();
221 
222  IOobject ddtIOobject
223  (
224  "ddt(" + rho.name() + ',' + vf.name() + ')',
225  mesh().time().timeName(),
226  mesh()
227  );
228 
229  if (mesh().moving())
230  {
231  return tmp<GeometricField<Type, fvPatchField, volMesh> >
232  (
233  new GeometricField<Type, fvPatchField, volMesh>
234  (
235  ddtIOobject,
236  mesh(),
237  rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
238  rDeltaT.internalField()*
239  (
240  rho.internalField()*vf.internalField()
241  - rho.oldTime().internalField()
242  *vf.oldTime().internalField()*mesh().Vsc0()/mesh().Vsc()
243  ),
244  rDeltaT.boundaryField()*
245  (
246  rho.boundaryField()*vf.boundaryField()
247  - rho.oldTime().boundaryField()
248  *vf.oldTime().boundaryField()
249  )
250  )
251  );
252  }
253  else
254  {
255  return tmp<GeometricField<Type, fvPatchField, volMesh> >
256  (
257  new GeometricField<Type, fvPatchField, volMesh>
258  (
259  ddtIOobject,
260  rDeltaT*(rho*vf - rho.oldTime()*vf.oldTime())
261  )
262  );
263  }
264 }
265 
266 
267 template<class Type>
268 tmp<GeometricField<Type, fvPatchField, volMesh> >
269 localEulerDdtScheme<Type>::fvcDdt
270 (
271  const volScalarField& alpha,
272  const volScalarField& rho,
273  const GeometricField<Type, fvPatchField, volMesh>& vf
274 )
275 {
276  const volScalarField& rDeltaT = localRDeltaT();
277 
278  IOobject ddtIOobject
279  (
280  "ddt("+alpha.name()+','+rho.name()+','+vf.name()+')',
281  mesh().time().timeName(),
282  mesh()
283  );
284 
285  if (mesh().moving())
286  {
287  return tmp<GeometricField<Type, fvPatchField, volMesh> >
288  (
289  new GeometricField<Type, fvPatchField, volMesh>
290  (
291  ddtIOobject,
292  mesh(),
293  rDeltaT.dimensions()
294  *alpha.dimensions()*rho.dimensions()*vf.dimensions(),
295  rDeltaT.internalField()*
296  (
297  alpha.internalField()
298  *rho.internalField()
299  *vf.internalField()
300 
301  - alpha.oldTime().internalField()
302  *rho.oldTime().internalField()
303  *vf.oldTime().internalField()*mesh().Vsc0()/mesh().Vsc()
304  ),
305  rDeltaT.boundaryField()*
306  (
307  alpha.boundaryField()
308  *rho.boundaryField()
309  *vf.boundaryField()
310 
311  - alpha.oldTime().boundaryField()
312  *rho.oldTime().boundaryField()
313  *vf.oldTime().boundaryField()
314  )
315  )
316  );
317  }
318  else
319  {
320  return tmp<GeometricField<Type, fvPatchField, volMesh> >
321  (
322  new GeometricField<Type, fvPatchField, volMesh>
323  (
324  ddtIOobject,
325  rDeltaT
326  *(
327  alpha*rho*vf
328  - alpha.oldTime()*rho.oldTime()*vf.oldTime()
329  )
330  )
331  );
332  }
333 }
334 
335 
336 template<class Type>
337 tmp<fvMatrix<Type> >
338 localEulerDdtScheme<Type>::fvmDdt
339 (
340  const GeometricField<Type, fvPatchField, volMesh>& vf
341 )
342 {
343  tmp<fvMatrix<Type> > tfvm
344  (
345  new fvMatrix<Type>
346  (
347  vf,
348  vf.dimensions()*dimVol/dimTime
349  )
350  );
351 
352  fvMatrix<Type>& fvm = tfvm();
353 
354  const scalarField& rDeltaT = localRDeltaT().internalField();
355 
356  fvm.diag() = rDeltaT*mesh().Vsc();
357 
358  if (mesh().moving())
359  {
360  fvm.source() = rDeltaT*vf.oldTime().internalField()*mesh().Vsc0();
361  }
362  else
363  {
364  fvm.source() = rDeltaT*vf.oldTime().internalField()*mesh().Vsc();
365  }
366 
367  return tfvm;
368 }
369 
370 
371 template<class Type>
372 tmp<fvMatrix<Type> >
373 localEulerDdtScheme<Type>::fvmDdt
374 (
375  const dimensionedScalar& rho,
376  const GeometricField<Type, fvPatchField, volMesh>& vf
377 )
378 {
379  tmp<fvMatrix<Type> > tfvm
380  (
381  new fvMatrix<Type>
382  (
383  vf,
384  rho.dimensions()*vf.dimensions()*dimVol/dimTime
385  )
386  );
387  fvMatrix<Type>& fvm = tfvm();
388 
389  const scalarField& rDeltaT = localRDeltaT().internalField();
390 
391  fvm.diag() = rDeltaT*rho.value()*mesh().Vsc();
392 
393  if (mesh().moving())
394  {
395  fvm.source() = rDeltaT
396  *rho.value()*vf.oldTime().internalField()*mesh().Vsc0();
397  }
398  else
399  {
400  fvm.source() = rDeltaT
401  *rho.value()*vf.oldTime().internalField()*mesh().Vsc();
402  }
403 
404  return tfvm;
405 }
406 
407 
408 template<class Type>
409 tmp<fvMatrix<Type> >
410 localEulerDdtScheme<Type>::fvmDdt
411 (
412  const volScalarField& rho,
413  const GeometricField<Type, fvPatchField, volMesh>& vf
414 )
415 {
416  tmp<fvMatrix<Type> > tfvm
417  (
418  new fvMatrix<Type>
419  (
420  vf,
421  rho.dimensions()*vf.dimensions()*dimVol/dimTime
422  )
423  );
424  fvMatrix<Type>& fvm = tfvm();
425 
426  const scalarField& rDeltaT = localRDeltaT().internalField();
427 
428  fvm.diag() = rDeltaT*rho.internalField()*mesh().Vsc();
429 
430  if (mesh().moving())
431  {
432  fvm.source() = rDeltaT
433  *rho.oldTime().internalField()
434  *vf.oldTime().internalField()*mesh().Vsc0();
435  }
436  else
437  {
438  fvm.source() = rDeltaT
439  *rho.oldTime().internalField()
440  *vf.oldTime().internalField()*mesh().Vsc();
441  }
442 
443  return tfvm;
444 }
445 
446 
447 template<class Type>
448 tmp<fvMatrix<Type> >
449 localEulerDdtScheme<Type>::fvmDdt
450 (
451  const volScalarField& alpha,
452  const volScalarField& rho,
453  const GeometricField<Type, fvPatchField, volMesh>& vf
454 )
455 {
456  tmp<fvMatrix<Type> > tfvm
457  (
458  new fvMatrix<Type>
459  (
460  vf,
461  alpha.dimensions()*rho.dimensions()*vf.dimensions()*dimVol/dimTime
462  )
463  );
464  fvMatrix<Type>& fvm = tfvm();
465 
466  const scalarField& rDeltaT = localRDeltaT().internalField();
467 
468  fvm.diag() = rDeltaT*alpha.internalField()*rho.internalField()*mesh().Vsc();
469 
470  if (mesh().moving())
471  {
472  fvm.source() = rDeltaT
473  *alpha.oldTime().internalField()
474  *rho.oldTime().internalField()
475  *vf.oldTime().internalField()*mesh().Vsc0();
476  }
477  else
478  {
479  fvm.source() = rDeltaT
480  *alpha.oldTime().internalField()
481  *rho.oldTime().internalField()
482  *vf.oldTime().internalField()*mesh().Vsc();
483  }
484 
485  return tfvm;
486 }
487 
488 
489 template<class Type>
490 tmp<typename localEulerDdtScheme<Type>::fluxFieldType>
491 localEulerDdtScheme<Type>::fvcDdtUfCorr
492 (
493  const GeometricField<Type, fvPatchField, volMesh>& U,
494  const GeometricField<Type, fvsPatchField, surfaceMesh>& Uf
495 )
496 {
497  const surfaceScalarField rDeltaT(fvc::interpolate(localRDeltaT()));
498 
499  fluxFieldType phiCorr
500  (
501  mesh().Sf() & (Uf.oldTime() - fvc::interpolate(U.oldTime()))
502  );
503 
504  return tmp<fluxFieldType>
505  (
506  new fluxFieldType
507  (
508  IOobject
509  (
510  "ddtCorr(" + U.name() + ',' + Uf.name() + ')',
511  mesh().time().timeName(),
512  mesh()
513  ),
514  this->fvcDdtPhiCoeff
515  (
516  U.oldTime(),
517  (mesh().Sf() & Uf.oldTime()),
518  phiCorr
519  )
520  *rDeltaT*phiCorr
521  )
522  );
523 }
524 
525 
526 template<class Type>
527 tmp<typename localEulerDdtScheme<Type>::fluxFieldType>
528 localEulerDdtScheme<Type>::fvcDdtPhiCorr
529 (
530  const GeometricField<Type, fvPatchField, volMesh>& U,
531  const fluxFieldType& phi
532 )
533 {
534  const surfaceScalarField rDeltaT(fvc::interpolate(localRDeltaT()));
535 
536  fluxFieldType phiCorr
537  (
538  phi.oldTime() - (mesh().Sf() & fvc::interpolate(U.oldTime()))
539  );
540 
541  return tmp<fluxFieldType>
542  (
543  new fluxFieldType
544  (
545  IOobject
546  (
547  "ddtCorr(" + U.name() + ',' + phi.name() + ')',
548  mesh().time().timeName(),
549  mesh()
550  ),
551  this->fvcDdtPhiCoeff(U.oldTime(), phi.oldTime(), phiCorr)
552  *rDeltaT*phiCorr
553  )
554  );
555 }
556 
557 
558 template<class Type>
559 tmp<typename localEulerDdtScheme<Type>::fluxFieldType>
560 localEulerDdtScheme<Type>::fvcDdtUfCorr
561 (
562  const volScalarField& rho,
563  const GeometricField<Type, fvPatchField, volMesh>& U,
564  const GeometricField<Type, fvsPatchField, surfaceMesh>& Uf
565 )
566 {
567  if
568  (
569  U.dimensions() == dimVelocity
570  && Uf.dimensions() == dimDensity*dimVelocity
571  )
572  {
573  const surfaceScalarField rDeltaT(fvc::interpolate(localRDeltaT()));
574 
575  GeometricField<Type, fvPatchField, volMesh> rhoU0
576  (
577  rho.oldTime()*U.oldTime()
578  );
579 
580  fluxFieldType phiCorr
581  (
582  mesh().Sf() & (Uf.oldTime() - fvc::interpolate(rhoU0))
583  );
584 
585  return tmp<fluxFieldType>
586  (
587  new fluxFieldType
588  (
589  IOobject
590  (
591  "ddtCorr("
592  + rho.name() + ',' + U.name() + ',' + Uf.name() + ')',
593  mesh().time().timeName(),
594  mesh()
595  ),
596  this->fvcDdtPhiCoeff
597  (
598  rhoU0,
599  mesh().Sf() & Uf.oldTime(),
600  phiCorr
601  )
602  *rDeltaT*phiCorr
603  )
604  );
605  }
606  else if
607  (
608  U.dimensions() == dimDensity*dimVelocity
609  && Uf.dimensions() == dimDensity*dimVelocity
610  )
611  {
612  return fvcDdtUfCorr(U, Uf);
613  }
614  else
615  {
616  FatalErrorIn
617  (
618  "localEulerDdtScheme<Type>::fvcDdtPhiCorr"
619  ) << "dimensions of Uf are not correct"
620  << abort(FatalError);
621 
622  return fluxFieldType::null();
623  }
624 }
625 
626 
627 template<class Type>
628 tmp<typename localEulerDdtScheme<Type>::fluxFieldType>
629 localEulerDdtScheme<Type>::fvcDdtPhiCorr
630 (
631  const volScalarField& rho,
632  const GeometricField<Type, fvPatchField, volMesh>& U,
633  const fluxFieldType& phi
634 )
635 {
636  if
637  (
638  U.dimensions() == dimVelocity
639  && phi.dimensions() == rho.dimensions()*dimVelocity*dimArea
640  )
641  {
642  dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
643 
644  GeometricField<Type, fvPatchField, volMesh> rhoU0
645  (
646  rho.oldTime()*U.oldTime()
647  );
648 
649  fluxFieldType phiCorr
650  (
651  phi.oldTime() - (mesh().Sf() & fvc::interpolate(rhoU0))
652  );
653 
654  return tmp<fluxFieldType>
655  (
656  new fluxFieldType
657  (
658  IOobject
659  (
660  "ddtCorr("
661  + rho.name() + ',' + U.name() + ',' + phi.name() + ')',
662  mesh().time().timeName(),
663  mesh()
664  ),
665  this->fvcDdtPhiCoeff(rhoU0, phi.oldTime(), phiCorr)
666  *rDeltaT*phiCorr
667  )
668  );
669  }
670  else if
671  (
672  U.dimensions() == rho.dimensions()*dimVelocity
673  && phi.dimensions() == rho.dimensions()*dimVelocity*dimArea
674  )
675  {
676  return fvcDdtPhiCorr(U, phi);
677  }
678  else
679  {
680  FatalErrorIn
681  (
682  "localEulerDdtScheme<Type>::fvcDdtPhiCorr"
683  ) << "dimensions of phi are not correct"
684  << abort(FatalError);
685 
686  return fluxFieldType::null();
687  }
688 }
689 
690 
691 template<class Type>
692 tmp<surfaceScalarField> localEulerDdtScheme<Type>::meshPhi
693 (
694  const GeometricField<Type, fvPatchField, volMesh>&
695 )
696 {
697  return tmp<surfaceScalarField>
698  (
699  new surfaceScalarField
700  (
701  IOobject
702  (
703  "meshPhi",
704  mesh().time().timeName(),
705  mesh(),
706  IOobject::NO_READ,
707  IOobject::NO_WRITE,
708  false
709  ),
710  mesh(),
711  dimensionedScalar("0", dimVolume/dimTime, 0.0)
712  )
713  );
714 }
715 
716 
717 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
718 
719 } // End namespace fv
720 
721 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
722 
723 } // End namespace Foam
724 
725 // ************************************************************************* //
Foam::fvc::interpolate
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &vf, const surfaceScalarField &faceFlux, Istream &schemeData)
Definition: surfaceInterpolate.C:110
Foam::calculatedFvPatchField
This boundary condition is not designed to be evaluated; it is assmued that the value is assigned via...
Definition: calculatedFvPatchField.H:65
Foam::GeometricField::boundaryField
GeometricBoundaryField & boundaryField()
Return reference to GeometricBoundaryField.
Definition: GeometricField.C:699
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::GeometricField::internalField
InternalField & internalField()
Return internal field.
Definition: GeometricField.C:688
surfaceInterpolate.H
Surface Interpolation.
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
Foam::fv::localEulerDdtScheme::meshPhi
tmp< surfaceScalarField > meshPhi(const GeometricField< Type, fvPatchField, volMesh > &)
Definition: localEulerDdtScheme.C:693
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
Foam::dimDensity
const dimensionSet dimDensity
Foam::fv::localEulerDdtScheme::fvcDdt
tmp< GeometricField< Type, fvPatchField, volMesh > > fvcDdt(const dimensioned< Type > &)
Definition: localEulerDdtScheme.C:52
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::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
Foam::dimensioned::name
const word & name() const
Return const reference to name.
Definition: dimensionedType.C:237
Foam::fv::localEulerDdt::localRDeltaT
static const volScalarField & localRDeltaT(const fvMesh &mesh)
Return the reciprocal of the local time-step.
Definition: localEulerDdt.C:45
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:52
Foam::DimensionedField::dimensions
const dimensionSet & dimensions() const
Return dimensions.
Definition: DimensionedFieldI.H:46
Foam::abort
errorManip< error > abort(error &err)
Definition: errorManip.H:131
Foam::IOobject::name
const word & name() const
Return name.
Definition: IOobject.H:260
Foam::lduMatrix::diag
scalarField & diag()
Definition: lduMatrix.C:183
FatalErrorIn
#define FatalErrorIn(functionName)
Report an error message using Foam::FatalError.
Definition: error.H:314
Foam::pTraits
Traits class for primitives.
Definition: pTraits.H:50
Foam::GeometricField::oldTime
const GeometricField< Type, PatchField, GeoMesh > & oldTime() const
Return old time field.
Definition: GeometricField.C:768
Foam::FatalError
error FatalError
Foam::fv::localEulerDdtScheme::fvcDdtUfCorr
tmp< fluxFieldType > fvcDdtUfCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Definition: localEulerDdtScheme.C:492
Foam::dimVolume
const dimensionSet dimVolume(pow3(dimLength))
Definition: dimensionSets.H:58
fv
labelList fv(nPoints)
Foam::fvMatrix::source
Field< Type > & source()
Definition: fvMatrix.H:291
Foam::dimVelocity
const dimensionSet dimVelocity
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
Foam::fv::localEulerDdtScheme::fvcDdtPhiCorr
tmp< fluxFieldType > fvcDdtPhiCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const fluxFieldType &phi)
Definition: localEulerDdtScheme.C:529
Foam::fv::localEulerDdtScheme::fvmDdt
tmp< fvMatrix< Type > > fvmDdt(const GeometricField< Type, fvPatchField, volMesh > &)
Definition: localEulerDdtScheme.C:339
Foam::dimensioned::value
const Type & value() const
Return const reference to value.
Definition: dimensionedType.C:263
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