From dbb8b04f5d2ac74f799571876cd3c822c28724d7 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 16 Nov 2016 16:27:11 -0800
Subject: [PATCH] ENH: Correcting Nu correlation for thermoCoupleProbes
---
.../utilities/thermoCoupleProbes/thermoCoupleProbes.C | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
index 91a6b26e2bd..740856a5a8d 100644
--- a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
+++ b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C
@@ -125,9 +125,9 @@ void Foam::functionObjects::thermoCoupleProbes::derivatives
muc = this->sample(thermo_.mu()());
Cpc = this->sample(thermo_.Cp()());
- scalarField Re(rhoc*Uc*d_/(muc + ROOTVSMALL));
- scalarField Pr(Cpc*muc/(kappac + ROOTVSMALL));
- //scalarField Nu(2.0 + 0.6*sqrt(Re)*cbrt(Pr));
+ scalarField Re(rhoc*Uc*d_/muc);
+ scalarField Pr(Cpc*muc/kappac);
+ Pr = max(ROOTVSMALL, Pr);
scalarField Nu(2.0 + (0.4*sqrt(Re) + 0.06*pow(Re, 2/3))*pow(Pr, 0.4));
scalarField htc(Nu*kappac/d_);
--
GitLab