43 void SLTSDdtScheme<Type>::relaxedDiag
57 diag[owner[faceI]] += phi[faceI];
58 rD[neighbour[faceI]] += phi[faceI];
62 diag[neighbour[faceI]] -= phi[faceI];
63 rD[owner[faceI]] -= phi[faceI];
74 if (pphi[patchFacei] > 0.0)
76 diag[faceCells[patchFacei]] += pphi[patchFacei];
80 rD[faceCells[patchFacei]] -= pphi[patchFacei];
85 rD += (1.0/alpha_ - 2.0)*
diag;
90 tmp<volScalarField> SLTSDdtScheme<Type>::SLrDeltaT()
const
97 tmp<volScalarField> trDeltaT
109 zeroGradientFvPatchScalarField::typeName
115 relaxedDiag(rDeltaT, phi);
117 if (phi.dimensions() == dimensionSet(0, 3, -1, 0, 0))
119 rDeltaT.internalField() =
max
121 rDeltaT.internalField()/
mesh().V(),
122 scalar(1)/deltaT.value()
125 else if (phi.dimensions() == dimensionSet(1, 0, -1, 0, 0))
131 rDeltaT.internalField() =
max
133 rDeltaT.internalField()/(rho.internalField()*
mesh().V()),
134 scalar(1)/deltaT.value()
140 <<
"Incorrect dimensions of phi: " << phi.dimensions()
144 rDeltaT.correctBoundaryConditions();
151 tmp<GeometricField<Type, fvPatchField, volMesh> >
161 "ddt("+dt.
name()+
')',
162 mesh().time().timeName(),
183 tdtdt().internalField() =
220 "ddt("+vf.
name()+
')',
221 mesh().time().timeName(),
272 "ddt("+rho.
name()+
','+vf.
name()+
')',
273 mesh().time().timeName(),
305 rDeltaT*rho*(vf - vf.
oldTime())
324 "ddt("+rho.
name()+
','+vf.
name()+
')',
325 mesh().time().timeName(),
385 scalarField rDeltaT = SLrDeltaT()().internalField();
387 Info<<
"max/min rDeltaT " <<
max(rDeltaT) <<
" " <<
min(rDeltaT) <<
endl;
422 scalarField rDeltaT = SLrDeltaT()().internalField();
459 scalarField rDeltaT = SLrDeltaT()().internalField();
491 "ddtPhiCorr(" + rA.
name() +
',' + U.
name() +
',' + phi.
name() +
')',
492 mesh().time().timeName(),
547 mesh().time().timeName(),
646 "SLTSDdtScheme<Type>::fvcDdtPhiCorr"
647 ) <<
"dimensions of phi are not correct"
650 return fluxFieldType::null();