Fixed missing return for bad flux

src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C

@@ -488,7 +488,12 @@ EulerDdtScheme<Type>::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())

src/finiteVolume/finiteVolume/ddtSchemes/steadyStateDdtScheme/steadyStateDdtScheme.C

@@ -318,6 +318,16 @@ steadyStateDdtScheme<Type>::fvcDdtPhiCorr
             )
         );
     }
+    else
+    {
+        FatalErrorIn
+        (
+            "steadyStateDdtScheme<Type>::fvcDdtPhiCorr"
+        )   << "dimensions of phi are not correct"
+            << abort(FatalError);
+
+        return fluxFieldType::null();
+    }
 }