43 tmp<GeometricField<Type, fvPatchField, volMesh> >
50 4.0/
sqr(
mesh().time().deltaT() +
mesh().time().deltaT0());
54 "d2dt2("+vf.
name()+
')',
55 mesh().time().timeName(),
61 scalar deltaT =
mesh().time().deltaT().value();
62 scalar deltaT0 =
mesh().time().deltaT0().value();
64 scalar coefft = (deltaT + deltaT0)/(2*deltaT);
65 scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
66 scalar coefft0 = coefft + coefft00;
70 scalar halfRdeltaT2 = rDeltaT2.
value()/2.0;
86 - (coefft*VV0 + coefft00*V0V00)
128 4.0/
sqr(
mesh().time().deltaT() +
mesh().time().deltaT0());
132 "d2dt2("+rho.
name()+
','+vf.
name()+
')',
133 mesh().time().timeName(),
139 scalar deltaT =
mesh().time().deltaT().value();
140 scalar deltaT0 =
mesh().time().deltaT0().value();
142 scalar coefft = (deltaT + deltaT0)/(2*deltaT);
143 scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
147 scalar halfRdeltaT2 = 0.5*rDeltaT2.
value();
148 scalar quarterRdeltaT2 = 0.25*rDeltaT2.
value();
172 - (coefft*VV0rhoRho0 + coefft00*V0V00rho0Rho00)
175 + (coefft00*V0V00rho0Rho00)
221 - (coefft*rhoRho0 + coefft00*rho0Rho00)*vf.
oldTime()
248 scalar deltaT =
mesh().time().deltaT().value();
249 scalar deltaT0 =
mesh().time().deltaT0().value();
251 scalar coefft = (deltaT + deltaT0)/(2*deltaT);
252 scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
253 scalar coefft0 = coefft + coefft00;
255 scalar rDeltaT2 = 4.0/
sqr(deltaT + deltaT0);
259 scalar halfRdeltaT2 = rDeltaT2/2.0;
264 fvm.
diag() = (coefft*halfRdeltaT2)*VV0;
266 fvm.
source() = halfRdeltaT2*
268 (coefft*VV0 + coefft00*V0V00)
276 fvm.
diag() = (coefft*rDeltaT2)*
mesh().V();
309 scalar deltaT =
mesh().time().deltaT().value();
310 scalar deltaT0 =
mesh().time().deltaT0().value();
312 scalar coefft = (deltaT + deltaT0)/(2*deltaT);
313 scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
315 scalar rDeltaT2 = 4.0/
sqr(deltaT + deltaT0);
319 scalar halfRdeltaT2 = 0.5*rDeltaT2;
325 fvm.
diag() = rho.
value()*(coefft*halfRdeltaT2)*VV0;
329 (coefft*VV0 + coefft00*V0V00)
370 scalar deltaT =
mesh().time().deltaT().value();
371 scalar deltaT0 =
mesh().time().deltaT0().value();
373 scalar coefft = (deltaT + deltaT0)/(2*deltaT);
374 scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
376 scalar rDeltaT2 = 4.0/
sqr(deltaT + deltaT0);
380 scalar quarterRdeltaT2 = 0.25*rDeltaT2;
393 fvm.
diag() = (coefft*quarterRdeltaT2)*VV0rhoRho0;
395 fvm.
source() = quarterRdeltaT2*
397 (coefft*VV0rhoRho0 + coefft00*V0V00rho0Rho00)
400 - (coefft00*V0V00rho0Rho00)
406 scalar halfRdeltaT2 = 0.5*rDeltaT2;
417 fvm.
diag() = (coefft*halfRdeltaT2)*
mesh().V()*rhoRho0;
421 (coefft*rhoRho0 + coefft00*rho0Rho00)
424 - (coefft00*rho0Rho00)