diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
index 0b899dc22f..f2ba60a720 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
@@ -156,14 +156,6 @@ void Foam::JohnsonJacksonParticleSlipFvPatchVectorField::updateCoeffs()
         )
     );
 
-    const scalarField nuFric
-    (
-        patch().lookupPatchField<volScalarField, scalar>
-        (
-            IOobject::groupName("nuFric", phase.name())
-        )
-    );
-
     word ThetaName(IOobject::groupName("Theta", phase.name()));
 
     const fvPatchScalarField& Theta
@@ -181,7 +173,7 @@ void Foam::JohnsonJacksonParticleSlipFvPatchVectorField::updateCoeffs()
        *gs0
        *specularityCoefficient_.value()
        *sqrt(3*Theta)
-       /max(6*(nu - nuFric)*phase.alphaMax(), small)
+       /max(6*nu*phase.alphaMax(), small)
     );
 
     this->valueFraction() = c/(c + patch().deltaCoeffs());
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
index 354bdfaa66..c301657a0c 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
@@ -293,7 +293,7 @@ Foam::RASModels::kineticTheoryModel::sigma() const
         volSymmTensorField::New
         (
             IOobject::groupName("R", U_.group()),
-          - (nut_)*dev(twoSymm(fvc::grad(U_)))
+          - (nut_ + nuFric_)*dev(twoSymm(fvc::grad(U_)))
           - (lambda_*fvc::div(phi_))*symmTensor::I
         )
     );
@@ -371,7 +371,7 @@ Foam::RASModels::kineticTheoryModel::devTau() const
         volSymmTensorField::New
         (
             IOobject::groupName("devTau", U_.group()),
-          - (rho_*nut_)
+          - (rho_*(nut_ + nuFric_))
            *dev(twoSymm(fvc::grad(U_)))
           - ((rho_*lambda_)*fvc::div(phi_))*symmTensor::I
         )
@@ -387,10 +387,10 @@ Foam::RASModels::kineticTheoryModel::divDevTau
 {
     return
     (
-      - fvm::laplacian(rho_*nut_, U)
+      - fvm::laplacian(rho_*(nut_ + nuFric_), U)
       - fvc::div
         (
-            (rho_*nut_)*dev2(T(fvc::grad(U)))
+            (rho_*(nut_ + nuFric_))*dev2(T(fvc::grad(U)))
           + ((rho_*lambda_)*fvc::div(phi_))
            *dimensioned<symmTensor>("I", dimless, symmTensor::I),
             "divDevTau(" + U_.name() + ')'
@@ -615,9 +615,6 @@ void Foam::RASModels::kineticTheoryModel::correct()
             ),
             maxNut_ - nut_
         );
-
-        // Add frictional viscosity
-        nut_ += nuFric_;
     }
 
     if (debug)