kineticTheoryModel: The frictional contribution nuFric is separated from the turbulent viscosity of the granular phase

Corrects granular phase heat transfer in cases where eddy diffusivity based
thermophysical transport models are used.

Patch contributed by VTT Technical Research Centre of Finland.

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());

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)