From 3ac2b9388afd088a38483560c3bad4a0e52d496a Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Wed, 21 Oct 2009 17:50:55 +0100
Subject: [PATCH] updates to enthalpy retained by particle due to phase change
---
.../Templates/ReactingParcel/ReactingParcel.C | 14 ++++++++++----
1 file changed, 10 insertions(+), 4 deletions(-)
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index f5975720e46..cadaed0d368 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -462,15 +462,21 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
{
const label idc =
td.cloud().composition().localToGlobalCarrierId(idPhase, i);
- const scalar hv = td.cloud().mcCarrierThermo().speciesData()[idc].H(T);
-
const label idl = td.cloud().composition().globalIds(idPhase)[i];
+
+ const scalar hv = td.cloud().mcCarrierThermo().speciesData()[idc].H(Ts);
const scalar hl =
- td.cloud().composition().liquids().properties()[idl].h(pc_, T);
+ td.cloud().composition().liquids().properties()[idl].h(pc_, Ts);
- // Enthalphy transfer to carrier phase
+ // Enthalphy transfer to carrier phase - method 1 using enthalpy diff
Sh += dMassPC[i]*(hl - hv)/dt;
+ // Enthalphy transfer to carrier phase - method 2 using latent heat
+// const scalar hl =
+// td.cloud().composition().liquids().properties()[idl].hl(pc_, Ts);
+// Sh -= dMassPC[i]*hl/dt;
+
+ // Update particle surface thermo properties
const scalar Dab =
td.cloud().composition().liquids().properties()[idl].D(pc_, Ts, Wc);
--
GitLab