44 scalar backwardDdtScheme<Type>::deltaT_()
const
46 return mesh().time().deltaT().value();
51 scalar backwardDdtScheme<Type>::deltaT0_()
const
53 return mesh().time().deltaT0().value();
58 template<
class GeoField>
59 scalar backwardDdtScheme<Type>::deltaT0_(
const GeoField& vf)
const
61 if (vf.nOldTimes() < 2)
75 tmp<GeometricField<Type, fvPatchField, volMesh> >
86 mesh().time().timeName(),
90 scalar deltaT = deltaT_();
91 scalar deltaT0 = deltaT0_();
93 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
94 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
95 scalar coefft0 = coefft + coefft00;
114 tdtdt().internalField() = rDeltaT.
value()*dt.
value()*
116 coefft - (coefft0*
mesh().V0() - coefft00*
mesh().V00())/
mesh().V()
153 "ddt("+vf.
name()+
')',
154 mesh().time().timeName(),
158 scalar deltaT = deltaT_();
159 scalar deltaT0 = deltaT0_(vf);
161 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
162 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
163 scalar coefft0 = coefft + coefft00;
225 "ddt("+rho.
name()+
','+vf.
name()+
')',
226 mesh().time().timeName(),
230 scalar deltaT = deltaT_();
231 scalar deltaT0 = deltaT0_(vf);
233 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
234 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
235 scalar coefft0 = coefft + coefft00;
296 "ddt("+rho.
name()+
','+vf.
name()+
')',
297 mesh().time().timeName(),
301 scalar deltaT = deltaT_();
302 scalar deltaT0 = deltaT0_(vf);
304 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
305 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
306 scalar coefft0 = coefft + coefft00;
377 scalar rDeltaT = 1.0/deltaT_();
379 scalar deltaT = deltaT_();
380 scalar deltaT0 = deltaT0_(vf);
382 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
383 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
384 scalar coefft0 = coefft + coefft00;
386 fvm.
diag() = (coefft*rDeltaT)*
mesh().V();
428 scalar rDeltaT = 1.0/deltaT_();
430 scalar deltaT = deltaT_();
431 scalar deltaT0 = deltaT0_(vf);
433 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
434 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
435 scalar coefft0 = coefft + coefft00;
479 scalar rDeltaT = 1.0/deltaT_();
481 scalar deltaT = deltaT_();
482 scalar deltaT0 = deltaT0_(vf);
484 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
485 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
486 scalar coefft0 = coefft + coefft00;
528 "ddtPhiCorr(" + rA.
name() +
',' + U.
name() +
',' + phi.
name() +
')',
529 mesh().time().timeName(),
533 scalar deltaT = deltaT_();
534 scalar deltaT0 = deltaT0_(U);
536 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
537 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
538 scalar coefft0 = coefft + coefft00;
550 - coefft00*phi.
oldTime().oldTime()
586 + phiAbs.
name() +
')',
587 mesh().time().timeName(),
591 scalar deltaT = deltaT_();
592 scalar deltaT0 = deltaT0_(U);
594 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
595 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
596 scalar coefft0 = coefft + coefft00;
652 - coefft00*phiAbs.
oldTime().oldTime()
686 - coefft00*phiAbs.
oldTime().oldTime()
706 "backwardDdtScheme<Type>::fvcDdtPhiCorr"
707 ) <<
"dimensions of phiAbs are not correct"
710 return fluxFieldType::null();
721 scalar deltaT = deltaT_();
722 scalar deltaT0 = deltaT0_(vf);
725 scalar coefft0_00 = deltaT/(deltaT + deltaT0);
728 scalar coefftn_0 = 1 + coefft0_00;
730 return coefftn_0*
mesh().phi() - coefft0_00*
mesh().phi().oldTime();