From fc2e2b54da38b3f38cb2dac246dec94a0dbb1c75 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Mon, 16 Sep 2019 11:33:48 +0100
Subject: [PATCH] TwoResistanceHeatTransferPhaseSystem: Restored asymmetric
 stabilisation term

See commit 7e799e9f for details.
---
 .../TwoResistanceHeatTransferPhaseSystem.C                    | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/TwoResistanceHeatTransferPhaseSystem/TwoResistanceHeatTransferPhaseSystem.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/TwoResistanceHeatTransferPhaseSystem/TwoResistanceHeatTransferPhaseSystem.C
index d1af6ad165..cda1aca9b6 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/TwoResistanceHeatTransferPhaseSystem/TwoResistanceHeatTransferPhaseSystem.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/TwoResistanceHeatTransferPhaseSystem/TwoResistanceHeatTransferPhaseSystem.C
@@ -479,10 +479,10 @@ heatTransfer() const
         const volScalarField HEff(H1*H2/(H1 + H2));
         *eqns[phase1.name()] +=
             HEff*(phase2.thermo().T() - phase1.thermo().T())
-          + HEff/Cpv1*he1 - fvm::Sp(HEff/Cpv1, he1);
+          + H1/Cpv1*he1 - fvm::Sp(H1/Cpv1, he1);
         *eqns[phase2.name()] +=
             HEff*(phase1.thermo().T() - phase2.thermo().T())
-          + HEff/Cpv2*he2 - fvm::Sp(HEff/Cpv2, he2);
+          + H2/Cpv2*he2 - fvm::Sp(H2/Cpv2, he2);
     }
 
     return eqnsPtr;