From 46df3ab4561efe663a1096905c5e90d78d65417c Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Fri, 19 Feb 2016 13:26:41 +0000 Subject: [PATCH] kineticTheoryModel: Corrected handling of the SuSp term in the Theta equation Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=2002 --- .../kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C | 2 +- .../kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C index 43ad071962..9f436cd5a0 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C @@ -434,7 +434,7 @@ void Foam::RASModels::kineticTheoryModel::correct() ) - fvm::laplacian(kappa_, Theta_, "laplacian(kappa,Theta)") == - fvm::SuSp(-((PsCoeff*I) && gradU), Theta_) + - fvm::SuSp((PsCoeff*I) && gradU, Theta_) + (tau && gradU) + fvm::Sp(-gammaCoeff, Theta_) + fvm::Sp(-J1, Theta_) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C index f07624c001..de35e7900b 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C @@ -437,7 +437,7 @@ void Foam::RASModels::kineticTheoryModel::correct() ) - fvm::laplacian(kappa_, Theta_, "laplacian(kappa,Theta)") == - fvm::SuSp(-((PsCoeff*I) && gradU), Theta_) + - fvm::SuSp((PsCoeff*I) && gradU, Theta_) + (tau && gradU) + fvm::Sp(-gammaCoeff, Theta_) + fvm::Sp(-J1, Theta_)