From ab7c7a7d1c2270d384cea8eda80d11fcc403e316 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 23 Feb 2015 18:56:39 +0000
Subject: [PATCH] heThermo: Correct kappaEff

---
 src/thermophysicalModels/basic/heThermo/heThermo.C | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/src/thermophysicalModels/basic/heThermo/heThermo.C b/src/thermophysicalModels/basic/heThermo/heThermo.C
index d1b26c5e..9bd7bccc 100644
--- a/src/thermophysicalModels/basic/heThermo/heThermo.C
+++ b/src/thermophysicalModels/basic/heThermo/heThermo.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -729,7 +729,7 @@ Foam::heThermo<BasicThermo, MixtureType>::kappaEff
     const volScalarField& alphat
 ) const
 {
-    tmp<Foam::volScalarField> kappaEff(Cp()*alphaEff(alphat));
+    tmp<Foam::volScalarField> kappaEff(Cp()*(this->alpha_ + alphat));
     kappaEff().rename("kappaEff");
     return kappaEff;
 }
@@ -749,7 +749,11 @@ Foam::heThermo<BasicThermo, MixtureType>::kappaEff
             this->p_.boundaryField()[patchi],
             this->T_.boundaryField()[patchi],
             patchi
-        )*alphaEff(alphat, patchi);
+        )
+       *(
+           this->alpha_.boundaryField()[patchi]
+         + alphat
+        );
 }