From 8c3d09d83f979c602961cc9b0478c33b6afe6e5f Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 19 Aug 2013 11:07:22 +0100
Subject: [PATCH] BUG: constTransport - corrected Prandtl number evaluation for
operators
---
.../specie/transport/const/constTransportI.H | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/thermophysicalModels/specie/transport/const/constTransportI.H b/src/thermophysicalModels/specie/transport/const/constTransportI.H
index 9949fddae4e..c0c15372ebf 100644
--- a/src/thermophysicalModels/specie/transport/const/constTransportI.H
+++ b/src/thermophysicalModels/specie/transport/const/constTransportI.H
@@ -157,7 +157,7 @@ inline void Foam::constTransport<Thermo>::operator+=
scalar molr2 = st.nMoles()/this->nMoles();
mu_ = molr1*mu_ + molr2*st.mu_;
- rPr_ = molr1*rPr_ + molr2*st.rPr_;
+ rPr_ = 1.0/(molr1/rPr_ + molr2/st.rPr_);
}
@@ -175,7 +175,7 @@ inline void Foam::constTransport<Thermo>::operator-=
scalar molr2 = st.nMoles()/this->nMoles();
mu_ = molr1*mu_ - molr2*st.mu_;
- rPr_ = molr1*rPr_ - molr2*st.rPr_;
+ rPr_ = 1.0/(molr1/rPr_ - molr2/st.rPr_);
}
@@ -210,7 +210,7 @@ inline Foam::constTransport<Thermo> Foam::operator+
(
t,
molr1*ct1.mu_ + molr2*ct2.mu_,
- molr1*ct1.rPr_ + molr2*ct2.rPr_
+ 1.0/(molr1/ct1.rPr_ + molr2/ct2.rPr_)
);
}
@@ -234,7 +234,7 @@ inline Foam::constTransport<Thermo> Foam::operator-
(
t,
molr1*ct1.mu_ - molr2*ct2.mu_,
- molr1*ct1.rPr_ - molr2*ct2.rPr_
+ 1.0/(molr1/ct1.rPr_ - molr2/ct2.rPr_)
);
}
@@ -250,7 +250,7 @@ inline Foam::constTransport<Thermo> Foam::operator*
(
s*static_cast<const Thermo&>(ct),
ct.mu_,
- ct.rPr_
+ 1.0/ct.rPr_
);
}
--
GitLab