Corrected calculation of iTp

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1456

Corrected calculation of iTp
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1456
bd2f16e2 · Henry · 11e3740f · bd2f16e2
Commit bd2f16e2 authored 10 years ago by Henry
--- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
+++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
@@ -339,7 +339,7 @@ Foam::scalar Foam::SprayParcel<ParcelType>::chi
            scalar TBoil = liq.pvInvert(p0);

            scalar hl = liq.hl(pAmb, TBoil);
-            scalar iTp = liq.h(pAmb, T0) - liq.rho(pAmb, T0);
+            scalar iTp = liq.h(pAmb, T0) - pAmb/liq.rho(pAmb, T0);
            scalar iTb = liq.h(pAmb, TBoil) - pAmb/liq.rho(pAmb, TBoil);

            chi += X[i]*(iTp - iTb)/hl;