diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
index 6ae60e42dc..2224030613 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
@@ -629,8 +629,9 @@ Foam::multiphaseSystem::dragCoeffs() const
             (
                 max
                 (
-                    fvc::average(dm.phase1())*fvc::average(dm.phase2()),
-                    //dm.phase1()*dm.phase2(),
+                    //fvc::average(dm.phase1()*dm.phase2()),
+                    //fvc::average(dm.phase1())*fvc::average(dm.phase2()),
+                    dm.phase1()*dm.phase2(),
                     dm.residualPhaseFraction()
                 )
                *dm.K