diff --git a/applications/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C b/applications/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
index 0403e693d5..3436d613bc 100644
--- a/applications/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
+++ b/applications/modules/multiphaseEuler/multiphaseEuler/thermophysicalPredictor.C
@@ -85,7 +85,8 @@ void Foam::solvers::multiphaseEuler::energyPredictor()
 
         const volScalarField& alpha = phase;
         const volScalarField& rho = phase.rho();
-        const volVectorField& U = phase.URef();
+        const tmp<volVectorField> tU(phase.U());
+        const volVectorField& U(tU());
 
         fvScalarMatrix EEqn
         (