44 tmp<volScalarField> CoEulerDdtScheme<Type>::CorDeltaT()
const
48 tmp<volScalarField> tcorDeltaT
55 cofrDeltaT.instance(),
60 zeroGradientFvPatchScalarField::typeName
71 corDeltaT[owner[faceI]] =
72 max(corDeltaT[owner[faceI]], cofrDeltaT[faceI]);
74 corDeltaT[neighbour[faceI]] =
75 max(corDeltaT[neighbour[faceI]], cofrDeltaT[faceI]);
81 cofrDeltaT.boundaryField()[
patchi];
83 const fvPatch&
p = pcofrDeltaT.patch();
86 forAll(pcofrDeltaT, patchFacei)
88 corDeltaT[faceCells[patchFacei]] =
max
90 corDeltaT[faceCells[patchFacei]],
91 pcofrDeltaT[patchFacei]
96 corDeltaT.correctBoundaryConditions();
105 tmp<surfaceScalarField> CoEulerDdtScheme<Type>::CofrDeltaT()
const
110 static_cast<const objectRegistry&
>(
mesh())
111 .lookupObject<surfaceScalarField>(phiName_);
113 if (phi.dimensions() == dimensionSet(0, 3, -1, 0, 0))
122 return max(Co/maxCo_, scalar(1))/deltaT;
124 else if (phi.dimensions() == dimensionSet(1, 0, -1, 0, 0))
127 static_cast<const objectRegistry&
>(
mesh())
128 .lookupObject<volScalarField>(rhoName_).oldTime();
137 return max(Co/maxCo_, scalar(1))/deltaT;
141 FatalErrorIn(
"CoEulerDdtScheme<Type>::CofrDeltaT() const")
142 <<
"Incorrect dimensions of phi: " << phi.dimensions()
145 return tmp<surfaceScalarField>(NULL);
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 = CorDeltaT()().internalField();
420 scalarField rDeltaT = CorDeltaT()().internalField();
457 scalarField rDeltaT = CorDeltaT()().internalField();
489 "ddtPhiCorr(" + rA.
name() +
',' + U.
name() +
',' + phi.
name() +
')',
490 mesh().time().timeName(),
545 mesh().time().timeName(),
644 "CoEulerDdtScheme<Type>::fvcDdtPhiCorr"
645 ) <<
"dimensions of phi are not correct"
648 return fluxFieldType::null();