diff --git a/applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/MichaelsBolger/MichaelsBolger.C b/applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/MichaelsBolger/MichaelsBolger.C
index 6cc4a62b62..2277bce35c 100644
--- a/applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/MichaelsBolger/MichaelsBolger.C
+++ b/applications/solvers/multiphase/driftFluxFoam/relativeVelocityModels/MichaelsBolger/MichaelsBolger.C
@@ -71,8 +71,10 @@ void Foam::relativeVelocityModels::MichaelsBolger::correct()
 {
     Udm_ =
         (mixture_.rhoc()/mixture_.rho())*Vc_
-       *(a0_ + pow(max(1 - mixture_.alphad()/mixture().alphaMax(), 0.0), a1_))
-       *acceleration();
+       *(
+           a0_
+         + pow(max(1 - mixture_.alphad()/mixture().alphaMax(), scalar(0)), a1_)
+        )*acceleration();
 }
 
 
diff --git a/src/twoPhaseModels/twoPhaseMixture/MPLIC/MPLICcell.C b/src/twoPhaseModels/twoPhaseMixture/MPLIC/MPLICcell.C
index 964439a2d3..c7c5387450 100644
--- a/src/twoPhaseModels/twoPhaseMixture/MPLIC/MPLICcell.C
+++ b/src/twoPhaseModels/twoPhaseMixture/MPLIC/MPLICcell.C
@@ -260,10 +260,10 @@ Foam::FixedList<Foam::scalar, 4> Foam::MPLICcell::solveVanderMatrix() const
 
     return FixedList<scalar, 4>
     {
-      - 4.5*b[0] + 13.5*b[1] - 13.5*b[2] + 4.5*b[3],
-        9.0*b[0] - 22.5*b[1] + 18.0*b[2] - 4.5*b[3],
-      - 5.5*b[0] +  9.0*b[1] -  4.5*b[2] + 1.0*b[3],
-        1.0*b[0]
+        scalar(- 4.5*b[0] + 13.5*b[1] - 13.5*b[2] + 4.5*b[3]),
+        scalar(9.0*b[0] - 22.5*b[1] + 18.0*b[2] - 4.5*b[3]),
+        scalar(- 5.5*b[0] +  9.0*b[1] -  4.5*b[2] + 1.0*b[3]),
+        scalar(1.0*b[0])
     };
 }