From 3c0ce7b152ae00e21488d777c109ea33763e8f72 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Thu, 13 Dec 2018 13:28:17 +0000
Subject: [PATCH] reactingMultiphaseEulerFoam: Updated EEqns for stationary
 phase

Resolves bug report https://bugs.openfoam.org/view.php?id=3132
---
 .../reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H     | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
index 22d875dce7..26d4dc2283 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
@@ -12,8 +12,8 @@ for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
         phaseModel& phase = fluid.anisothermalPhases()[anisothermalPhasei];
 
         const volScalarField& alpha = phase;
-        const volScalarField& rho = phase.rho();
-        const volVectorField& U = phase.U();
+        tmp<volScalarField> rho = phase.rho();
+        tmp<volVectorField> U = phase.U();
 
         fvScalarMatrix EEqn
         (