From e2a5ffbd8f41ac8a27c874c88bb96f206d444b27 Mon Sep 17 00:00:00 2001
From: s126103 <t.j.m.nijssen@tue.nl>
Date: Mon, 11 May 2020 14:24:42 +0200
Subject: [PATCH] bugfix

---
 .../energyModel/wallHeatTransferYagi/wallHeatTransferYagi.C     | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/wallHeatTransferYagi/wallHeatTransferYagi.C b/src/lagrangian/cfdemParticle/subModels/energyModel/wallHeatTransferYagi/wallHeatTransferYagi.C
index e7a210dd..f6d2a158 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/wallHeatTransferYagi/wallHeatTransferYagi.C
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/wallHeatTransferYagi/wallHeatTransferYagi.C
@@ -265,7 +265,7 @@ void wallHeatTransferYagi::calcEnergyContribution()
                     // calculate H
                     scalar H;
                     if (voidfraction_[faceCelli]<=voidfractionMax_)
-                        H = 0.2087 * (pow(ReField_[faceCelli]+SMALL,-0.20)) * CpField_[faceCelli] * magG / (pow(PrField_[faceCelli],2/3) + SMALL);
+                        H = 0.2087 * pow(ReField_[faceCelli]+SMALL,-0.20) * CpField_[faceCelli] * magG * pow(PrField_[faceCelli]+SMALL,(-2./3.));
                     else
                         H = 0;