BUG: Using autoPtr for turbulence in interFoam and interIsoFoam

applications/solvers/multiphase/interFoam/UEqn.H

@@ -4,7 +4,7 @@
     (
         fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
       + MRF.DDt(rho, U)
-      + turbulence.divDevRhoReff(rho, U)
+      + turbulence->divDevRhoReff(rho, U)
      ==
         fvOptions(rho, U)
     );

applications/solvers/multiphase/interFoam/createFields.H

@@ -70,8 +70,18 @@ surfaceScalarField rhoPhi
     fvc::interpolate(rho)*phi
 );
 
-incompressibleInterPhaseTransportModel<immiscibleIncompressibleTwoPhaseMixture>
-    turbulence(rho, U, phi, rhoPhi, mixture);
+typedef incompressibleInterPhaseTransportModel
+    <
+        immiscibleIncompressibleTwoPhaseMixture
+    > transportModelType;
+
+autoPtr<transportModelType> turbulence
+(
+    new transportModelType
+    (
+        rho, U, phi, rhoPhi, mixture
+    )
+);
 
 #include "readGravitationalAcceleration.H"
 #include "readhRef.H"

applications/solvers/multiphase/interFoam/interFoam.C

@@ -168,7 +168,7 @@ int main(int argc, char *argv[])
 
             if (pimple.turbCorr())
             {
-                turbulence.correct();
+                turbulence->correct();
             }
         }
 

applications/solvers/multiphase/interIsoFoam/createFields.H

@@ -71,8 +71,18 @@ surfaceScalarField rhoPhi
 );
 
 
-incompressibleInterPhaseTransportModel<immiscibleIncompressibleTwoPhaseMixture>
-    turbulence(rho, U, phi, rhoPhi, mixture);
+typedef incompressibleInterPhaseTransportModel
+    <
+        immiscibleIncompressibleTwoPhaseMixture
+    > transportModelType;
+
+autoPtr<transportModelType> turbulence
+(
+    new transportModelType
+    (
+        rho, U, phi, rhoPhi, mixture
+    )
+);
 
 
 #include "readGravitationalAcceleration.H"

applications/solvers/multiphase/interIsoFoam/interIsoFoam.C

@@ -180,7 +180,7 @@ int main(int argc, char *argv[])
 
             if (pimple.turbCorr())
             {
-                turbulence.correct();
+                turbulence->correct();
             }
         }