diff --git a/applications/solvers/combustion/XiFoam/bEqn.H b/applications/solvers/combustion/XiFoam/bEqn.H
index df37bf34c5..ff19598ba1 100644
--- a/applications/solvers/combustion/XiFoam/bEqn.H
+++ b/applications/solvers/combustion/XiFoam/bEqn.H
@@ -202,8 +202,11 @@ if (ign.ignited())
         volScalarField XiEq
         (
             scalar(1.001)
-          + (scalar(1) + (2*XiShapeCoef)*(scalar(0.5) - min(max(b, 0.0), 1.0)))
-            *(XiEqStar - scalar(1.001))
+          + (
+                scalar(1)
+              + (2*XiShapeCoef)
+               *(scalar(0.5) - min(max(b, scalar(0)), scalar(1)))
+            )*(XiEqStar - scalar(1.001))
         );
 
         volScalarField Gstar(0.28/tauEta);
diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
index 31c1023b8d..d733640149 100644
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/alphaEqn.H
@@ -86,7 +86,7 @@ surfaceScalarField alphaPhi2("alphaPhi2", phi2);
         alpha1Eqn.solve();
 
         //***HGW temporary boundedness-fix pending the introduction of MULES
-        alpha1 = max(min(alpha1, 1.0), 0.0);
+        alpha1 = max(min(alpha1, scalar(1)), scalar(0));
 
         #include "packingLimiter.H"
 
diff --git a/src/combustionModels/FSD/FSD.C b/src/combustionModels/FSD/FSD.C
index f092147a1d..f5b3d4e8d2 100644
--- a/src/combustionModels/FSD/FSD.C
+++ b/src/combustionModels/FSD/FSD.C
@@ -199,7 +199,7 @@ void FSD<CombThermoType, ThermoType>::calculateSourceNorm()
     );
 
     // Linear correlation between delta and flame thickness
-    volScalarField omegaF(max(deltaF*(4.0/3.0) + (2.0/3.0), 1.0));
+    volScalarField omegaF(max(deltaF*(4.0/3.0) + (2.0/3.0), scalar(1)));
 
     scalar deltaFt = 1.0/ftDim_;
 
@@ -319,9 +319,9 @@ void FSD<CombThermoType, ThermoType>::calculateSourceNorm()
         products += Yp;
     }
 
-    volScalarField c(max(scalar(1.0) - products/max(pc, 1e-5), 0.0));
+    volScalarField c(max(scalar(1) - products/max(pc, scalar(1e-5)), scalar(0)));
 
-    pc = min(C_*c, scalar(1.0));
+    pc = min(C_*c, scalar(1));
 
     const volScalarField fres(this->singleMixture_.fres(fuelI));
 
diff --git a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
index f1478a7db7..7dd7adb80e 100644
--- a/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
+++ b/src/turbulenceModels/incompressible/RAS/kkLOmega/kkLOmega.C
@@ -77,7 +77,7 @@ tmp<volScalarField> kkLOmega::Cmu(const volScalarField& S) const
 
 tmp<volScalarField> kkLOmega::BetaTS(const volScalarField& Rew) const
 {
-    return(scalar(1.0) - exp(-sqr(max(Rew - CtsCrit_, 0.0))/Ats_));
+    return(scalar(1) - exp(-sqr(max(Rew - CtsCrit_, scalar(0)))/Ats_));
 }
 
 
@@ -89,7 +89,7 @@ tmp<volScalarField> kkLOmega::fTaul
 {
     return
     (
-        scalar(1.0)
+        scalar(1)
       - exp
         (
             -CtauL_*ktL
@@ -130,7 +130,7 @@ tmp<volScalarField> kkLOmega::fOmega
 {
     return
     (
-        scalar(1.0)
+        scalar(1)
       - exp
         (
             -0.41
@@ -159,15 +159,12 @@ tmp<volScalarField> kkLOmega::gammaBP(const volScalarField& omega) const
         max
         (
             kt_/nu()
-          /
-            (
+          / (
                 omega
               + dimensionedScalar("ROTVSMALL", omega.dimensions(), ROOTVSMALL)
             )
-          -
-            CbpCrit_
-        ,
-            0.0
+          - CbpCrit_,
+            scalar(0)
         )
     );
 }
@@ -187,9 +184,8 @@ tmp<volScalarField> kkLOmega::gammaNAT
           - CnatCrit_
           / (
                 fNatCrit + dimensionedScalar("ROTVSMALL", dimless, ROOTVSMALL)
-            )
-          ,
-            0.0
+            ),
+            scalar(0)
         )
     );
 }