44 scalar boundedBackwardFaDdtScheme::deltaT_()
const
50 scalar boundedBackwardFaDdtScheme::deltaT0_()
const
70 mesh()().time().timeName(),
76 scalar deltaT = deltaT_();
77 scalar deltaT0 = deltaT0_();
79 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
80 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
81 scalar coefft0 = coefft + coefft00;
97 coefft - (coefft0*
mesh().S0() - coefft00*
mesh().S00())/
mesh().S()
108 calculatedFaPatchScalarField::typeName
124 "ddt("+dt.
name()+
')',
125 mesh()().time().timeName(),
131 scalar deltaT = deltaT_();
132 scalar deltaT0 = deltaT0_();
134 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
135 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
136 scalar coefft0 = coefft + coefft00;
144 -rDeltaT*(coefft0 - coefft00)*dt
150 tdtdt0.
ref().primitiveFieldRef() = (-rDeltaT.
value()*dt.value())*
152 (coefft0*
mesh().S0() - coefft00*
mesh().S00())/
mesh().S()
169 "ddt("+vf.name()+
')',
170 mesh()().time().timeName(),
176 scalar deltaT = deltaT_();
177 scalar deltaT0 = deltaT0_(vf);
266 "ddt0("+vf.name()+
')',
267 mesh()().time().timeName(),
273 scalar deltaT = deltaT_();
274 scalar deltaT0 = deltaT0_(vf);
361 "ddt("+
rho.name()+
','+vf.name()+
')',
362 mesh()().time().timeName(),
368 scalar deltaT = deltaT_();
369 scalar deltaT0 = deltaT0_(vf);
458 "ddt0("+
rho.name()+
','+vf.name()+
')',
459 mesh()().time().timeName(),
465 scalar deltaT = deltaT_();
466 scalar deltaT0 = deltaT0_(vf);
553 "ddt("+
rho.name()+
','+vf.name()+
')',
554 mesh()().time().timeName(),
560 scalar deltaT = deltaT_();
561 scalar deltaT0 = deltaT0_(vf);
603 rho.oldTime().primitiveField()*
606 rho.oldTime().oldTime().primitiveField()
615 rho.oldTime().boundaryField()*
618 rho.oldTime().oldTime().boundaryField()*
654 "ddt0("+
rho.name()+
','+vf.name()+
')',
655 mesh()().time().timeName(),
661 scalar deltaT = deltaT_();
662 scalar deltaT0 = deltaT0_(vf);
703 rho.oldTime().primitiveField()*
706 rho.oldTime().oldTime().primitiveField()*
714 rho.oldTime().boundaryField()*
717 rho.oldTime().oldTime().boundaryField()*
758 scalar rDeltaT = 1.0/deltaT_();
760 scalar deltaT = deltaT_();
761 scalar deltaT0 = deltaT0_(vf);
789 fam.diag() = (coefft*rDeltaT)*
mesh().S();
793 fam.source() = rDeltaT*
802 fam.source() = rDeltaT*
mesh().S()*
829 scalar rDeltaT = 1.0/deltaT_();
831 scalar deltaT = deltaT_();
832 scalar deltaT0 = deltaT0_(vf);
860 fam.diag() = (coefft*rDeltaT*
rho.value())*
mesh().S();
864 fam.source() = rDeltaT*
rho.value()*
873 fam.source() = rDeltaT*
mesh().S()*
rho.value()*
900 scalar rDeltaT = 1.0/deltaT_();
902 scalar deltaT = deltaT_();
903 scalar deltaT0 = deltaT0_(vf);
912 rho.oldTime().oldTime().internalField()*
914 -
rho.oldTime().oldTime().oldTime().internalField()*
920 rho.oldTime().internalField()*
922 -
rho.oldTime().oldTime().internalField()*
935 fam.diag() = (coefft*rDeltaT)*
rho.primitiveField()*
mesh().S();
939 fam.source() = rDeltaT*
941 coefft0*
rho.oldTime().primitiveField()
943 - coefft00*
rho.oldTime().oldTime().primitiveField()
949 fam.source() = rDeltaT*
mesh().S()*
951 coefft0*
rho.oldTime().primitiveField()
953 - coefft00*
rho.oldTime().oldTime().primitiveField()