From 2602523b6c419b3181e794b16a6aae19e17e9bf0 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Wed, 13 Apr 2011 11:08:18 +0100
Subject: [PATCH] ENH: Updated momentum coupling for reacting parcels

---
 .../ReactingMultiphaseParcel/ReactingMultiphaseParcel.C        | 3 +--
 .../parcels/Templates/ReactingParcel/ReactingParcel.C          | 3 +--
 2 files changed, 2 insertions(+), 4 deletions(-)

diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
index bec3481a0bb..0d605da15d3 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
@@ -327,7 +327,6 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
 
 
 
-
     // Heat transfer
     // ~~~~~~~~~~~~~
 
@@ -369,7 +368,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
             d0,
             U0,
             rho0,
-            mass0,
+            0.5*(mass0 + mass1),
             Su,
             dUTrans,
             Spu
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 5fe67854f85..b5a27067326 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -379,13 +379,12 @@ void Foam::ReactingParcel<ParcelType>::calc
             d0,
             U0,
             rho0,
-            mass0,
+            0.5*(mass0 + mass1),
             Su,
             dUTrans,
             Spu
         );
 
-    dUTrans += 0.5*(mass0 - mass1)*(U0 + U1);
 
     // Accumulate carrier phase source terms
     // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-- 
GitLab