From f412d39174d6e80407f3414277425742f95e0fd0 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Sun, 25 Sep 2022 21:49:37 +0100 Subject: [PATCH] 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. --- .../JohnsonJacksonParticleSlipFvPatchVectorField.C | 10 +--------- .../kineticTheoryModel/kineticTheoryModel.C | 11 ++++------- 2 files changed, 5 insertions(+), 16 deletions(-) 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 - ( - 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("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)