From 8f8612fd0c0ab95fa598f82bfa01993f58b5b442 Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Fri, 22 Apr 2022 09:51:29 -0700
Subject: [PATCH] BUG: multiphaseInterSystem: fix kappaEff functions

---
 .../multiphaseInterSystem.C                      | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/multiphaseInterSystem/multiphaseInterSystem.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/multiphaseInterSystem/multiphaseInterSystem.C
index dbe7501a11..952032240c 100644
--- a/src/phaseSystemModels/multiphaseInter/phasesSystem/multiphaseInterSystem/multiphaseInterSystem.C
+++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/multiphaseInterSystem/multiphaseInterSystem.C
@@ -926,15 +926,25 @@ Foam::tmp<Foam::volScalarField> Foam::multiphaseInterSystem::nuEff() const
 
 Foam::tmp<Foam::volScalarField> Foam::multiphaseInterSystem::kappaEff() const
 {
-    return(this->kappa() + this->Cp()*turb_->mut()/Prt_);
+    const volScalarField kappaEff
+    (
+        this->kappa() + this->Cp()*turb_->mut()/Prt_
+    );
+
+    return tmp<volScalarField>::New(kappaEff);
 }
 
 
 Foam::tmp<Foam::scalarField>
 Foam::multiphaseInterSystem::kappaEff(const label patchi) const
 {
-    tmp<scalarField> tCp(this->Cp()().boundaryField()[patchi]);
-    return this->kappa(patchi) + tCp()*turb_->mut(patchi)/Prt_.value();
+    const scalarField Cp(this->Cp()().boundaryField()[patchi]);
+    const scalarField kappaEffp
+    (
+         this->kappa(patchi) + Cp*turb_->mut(patchi)/Prt_.value()
+    );
+
+    return tmp<scalarField>::New(kappaEffp);
 }