diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H
index 2c2f1281f8..a2802f8564 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/UEqns.H
@@ -3,7 +3,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
 
 {
     {
-        volTensorField gradU1T = fvc::grad(U1)().T();
+        volTensorField gradU1T(fvc::grad(U1)().T());
 
         if (kineticTheory.on())
         {
@@ -45,7 +45,7 @@ fvVectorMatrix U2Eqn(U2, U2.dimensions()*dimVol/dimTime);
     }
 
     {
-        volTensorField gradU2T = fvc::grad(U2)().T();
+        volTensorField gradU2T(fvc::grad(U2)().T());
         volTensorField Rc2
         (
             "Rc2",
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
index f0ebde6332..31c1023b8d 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
@@ -10,8 +10,8 @@ surfaceScalarField alphaPhi2("alphaPhi2", phi2);
 
     if (g0.value() > 0.0)
     {
-        surfaceScalarField alpha1f = fvc::interpolate(alpha1);
-        surfaceScalarField phipp = ppMagf*fvc::snGrad(alpha1)*mesh.magSf();
+        surfaceScalarField alpha1f(fvc::interpolate(alpha1));
+        surfaceScalarField phipp(ppMagf*fvc::snGrad(alpha1)*mesh.magSf());
         phir += phipp;
         phic += fvc::interpolate(alpha1)*phipp;
     }
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/createFields.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/createFields.H
index b3eca666c2..4642e01330 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/createFields.H
@@ -228,15 +228,19 @@
 
     Info<< "Calculating field DDtU1 and DDtU2\n" << endl;
 
-    volVectorField DDtU1 =
+    volVectorField DDtU1
+    (
         fvc::ddt(U1)
       + fvc::div(phi1, U1)
-      - fvc::div(phi1)*U1;
+      - fvc::div(phi1)*U1
+    );
 
-    volVectorField DDtU2 =
+    volVectorField DDtU2
+    (
         fvc::ddt(U2)
       + fvc::div(phi2, U2)
-      - fvc::div(phi2)*U2;
+      - fvc::div(phi2)*U2
+    );
 
 
     Info<< "Calculating field g.h\n" << endl;
@@ -354,8 +358,10 @@
     setRefCell(p, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);
 
 
-    volScalarField dgdt =
-        pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001));
+    volScalarField dgdt
+    (
+        pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001))
+    );
 
 
     Info<< "Creating field dpdt\n" << endl;
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialCoeffs.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialCoeffs.H
index 436b68a2f2..3d32463050 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialCoeffs.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialCoeffs.H
@@ -35,8 +35,8 @@ volScalarField heatTransferCoeff
 );
 
 {
-    volVectorField Ur = U1 - U2;
-    volScalarField magUr = mag(Ur);
+    volVectorField Ur(U1 - U2);
+    volScalarField magUr(mag(Ur));
 
     if (dispersedPhase == "1")
     {
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
index cb9399921a..bac588c738 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
@@ -72,7 +72,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::GidaspowErgunWenYu::K
 ) const
 {
     volScalarField beta(max(scalar(1) - alpha_, scalar(1.0e-6)));
-    volScalarField d = phase1_.d();
+    volScalarField d(phase1_.d());
     volScalarField bp(pow(beta, -2.65));
     volScalarField Re(max(Ur*d/phase2_.nu(), scalar(1.0e-3)));
 
@@ -87,7 +87,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::GidaspowErgunWenYu::K
     }
 
     // Wen and Yu (1966)
-    tmp<volScalarField> tKWenYu = 0.75*Cds*phase2_.rho()*Ur*bp/d;
+    tmp<volScalarField> tKWenYu(0.75*Cds*phase2_.rho()*Ur*bp/d);
     volScalarField& KWenYu = tKWenYu();
 
     // Ergun
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
index 34a945690a..cc53c67d76 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
@@ -71,10 +71,10 @@ Foam::tmp<Foam::volScalarField> Foam::heatTransferModels::RanzMarshall::K
     const volScalarField& Ur
 ) const
 {
-    volScalarField Re = max(Ur*phase1_.d()/phase2_.nu(), scalar(1.0e-3));
+    volScalarField Re(max(Ur*phase1_.d()/phase2_.nu(), scalar(1.0e-3)));
     dimensionedScalar Prb =
         phase2_.rho()*phase2_.nu()*phase2_.Cp()/phase2_.kappa();
-    volScalarField Nu = scalar(2) + 0.6*sqrt(Re)*cbrt(Prb);
+    volScalarField Nu(scalar(2) + 0.6*sqrt(Re)*cbrt(Prb));
 
     return 6.0*phase2_.kappa()*Nu/sqr(phase1_.d());
 }
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
index 6e2fc44b08..34a9a2600a 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
@@ -24,14 +24,14 @@
           - fvc::Sp(fvc::div(phid2), p);
     }
 
-    surfaceScalarField alpha1f = fvc::interpolate(alpha1);
-    surfaceScalarField alpha2f = scalar(1) - alpha1f;
+    surfaceScalarField alpha1f(fvc::interpolate(alpha1));
+    surfaceScalarField alpha2f(scalar(1) - alpha1f);
 
     volVectorField U10 = U1;
     volVectorField U20 = U2;
 
-    volScalarField rAU1 = 1.0/U1Eqn.A();
-    volScalarField rAU2 = 1.0/U2Eqn.A();
+    volScalarField rAU1(1.0/U1Eqn.A());
+    volScalarField rAU2(1.0/U2Eqn.A());
 
     surfaceScalarField rAlphaAU1f = fvc::interpolate(alpha1*rAU1);
     surfaceScalarField rAlphaAU2f = fvc::interpolate(alpha2*rAU2);
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/turbulenceModel/kEpsilon.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/turbulenceModel/kEpsilon.H
index cb19fba0b7..1ed20444bc 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/turbulenceModel/kEpsilon.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/turbulenceModel/kEpsilon.H
@@ -1,12 +1,12 @@
-if(turbulence)
+if (turbulence)
 {
     if (mesh.changing())
     {
         y.correct();
     }
 
-    tmp<volTensorField> tgradU2 = fvc::grad(U2);
-    volScalarField G = 2*nut2*(tgradU2() && dev(symm(tgradU2())));
+    tmp<volTensorField> tgradU2(fvc::grad(U2));
+    volScalarField G(2*nut2*(tgradU2() && dev(symm(tgradU2()))));
     tgradU2.clear();
 
     #include "wallFunctions.H"
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C b/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
index 34a945690a..cc53c67d76 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
@@ -71,10 +71,10 @@ Foam::tmp<Foam::volScalarField> Foam::heatTransferModels::RanzMarshall::K
     const volScalarField& Ur
 ) const
 {
-    volScalarField Re = max(Ur*phase1_.d()/phase2_.nu(), scalar(1.0e-3));
+    volScalarField Re(max(Ur*phase1_.d()/phase2_.nu(), scalar(1.0e-3)));
     dimensionedScalar Prb =
         phase2_.rho()*phase2_.nu()*phase2_.Cp()/phase2_.kappa();
-    volScalarField Nu = scalar(2) + 0.6*sqrt(Re)*cbrt(Prb);
+    volScalarField Nu(scalar(2) + 0.6*sqrt(Re)*cbrt(Prb));
 
     return 6.0*phase2_.kappa()*Nu/sqr(phase1_.d());
 }
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H b/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H
index 98fa7b139e..744faab968 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H
@@ -172,7 +172,7 @@
 
         if (pimple.finalNonOrthogonalIter())
         {
-            surfaceScalarField mSfGradp = pEqnIncomp.flux()/Dp;
+            surfaceScalarField mSfGradp(pEqnIncomp.flux()/Dp);
 
             phasei = 0;
             phi = dimensionedScalar("phi", phi.dimensions(), 0);
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/readMechanicalProperties.H b/applications/solvers/stressAnalysis/solidDisplacementFoam/readMechanicalProperties.H
index 762245abf4..6cfd4dbcc9 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/readMechanicalProperties.H
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/readMechanicalProperties.H
@@ -185,7 +185,7 @@
     volScalarField& nu = nuPtr();
 
     Info<< "Normalising E : E/rho\n" << endl;
-    volScalarField E = rhoE/rho;
+    volScalarField E(rhoE/rho);
 
     Info<< "Calculating Lame's coefficients\n" << endl;
 
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermalProperties.H b/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermalProperties.H
index 54ac7d6db9..c8817c5975 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermalProperties.H
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermalProperties.H
@@ -222,7 +222,7 @@ if (thermalStress)
     volScalarField& alpha = alphaPtr();
 
     Info<< "Normalising k : k/rho\n" << endl;
-    volScalarField k = rhoK/rho;
+    volScalarField k(rhoK/rho);
 
     Info<< "Calculating thermal coefficients\n" << endl;
 
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/tractionDisplacement/tractionDisplacementFvPatchVectorField.C b/applications/solvers/stressAnalysis/solidDisplacementFoam/tractionDisplacement/tractionDisplacementFvPatchVectorField.C
index b24239c0f8..a0056d3c9a 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/tractionDisplacement/tractionDisplacementFvPatchVectorField.C
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/tractionDisplacement/tractionDisplacementFvPatchVectorField.C
@@ -159,10 +159,10 @@ void tractionDisplacementFvPatchVectorField::updateCoeffs()
     const fvPatchField<scalar>& nu =
         patch().lookupPatchField<volScalarField, scalar>("nu");
 
-    scalarField E = rhoE/rho;
-    scalarField mu = E/(2.0*(1.0 + nu));
-    scalarField lambda = nu*E/((1.0 + nu)*(1.0 - 2.0*nu));
-    scalarField threeK = E/(1.0 - 2.0*nu);
+    scalarField E(rhoE/rho);
+    scalarField mu(E/(2.0*(1.0 + nu)));
+    scalarField lambda(nu*E/((1.0 + nu)*(1.0 - 2.0*nu)));
+    scalarField threeK(E/(1.0 - 2.0*nu));
 
     Switch planeStress(mechanicalProperties.lookup("planeStress"));
 
@@ -172,7 +172,7 @@ void tractionDisplacementFvPatchVectorField::updateCoeffs()
         threeK = E/(1.0 - nu);
     }
 
-    scalarField twoMuLambda = (2*mu + lambda);
+    scalarField twoMuLambda(2*mu + lambda);
 
     vectorField n(patch().nf());