diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C index 0e3015729..f179cf95c 100644 --- a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C +++ b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C @@ -488,7 +488,12 @@ EulerDdtScheme::fvcDdtPhiCorr ( ddtIOobject, rDeltaT - *this->fvcDdtPhiCoeff(rho.oldTime(), U.oldTime(), phiAbs.oldTime()) + *this->fvcDdtPhiCoeff + ( + rho.oldTime(), + U.oldTime(), + phiAbs.oldTime() + ) *( fvc::interpolate(rA)*phiAbs.oldTime() - (fvc::interpolate(rA*U.oldTime()) & mesh().Sf()) diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C index 9f18dffa2..7f3b91a40 100644 --- a/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C +++ b/src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C @@ -318,6 +318,16 @@ steadyStateDdtScheme::fvcDdtPhiCorr ) ); } + else + { + FatalErrorIn + ( + "steadyStateDdtScheme::fvcDdtPhiCorr" + ) << "dimensions of phi are not correct" + << abort(FatalError); + + return fluxFieldType::null(); + } }