From 2176098a3ed7280c0df0dc4740aa5a4ecfb27d74 Mon Sep 17 00:00:00 2001 From: Henry Date: Fri, 1 Mar 2013 17:45:42 +0000 Subject: [PATCH] compressibleTwoPhaseEulerFoam: select the phase energy source terms correctly --- .../solvers/multiphase/compressibleTwoPhaseEulerFoam/EEqns.H | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/EEqns.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/EEqns.H index b9a47a917d..1f12fc9a12 100644 --- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/EEqns.H +++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/EEqns.H @@ -27,7 +27,7 @@ - (fvc::ddt(alpha1) + fvc::div(alphaPhi1))*K1 + ( - he1.name() == "e" + he1.name() == thermo1.phasePropertyName("e") ? fvc::div(alphaPhi1, p) : -dalpha1pdt )/rho1 @@ -49,7 +49,7 @@ - (fvc::ddt(alpha2) + fvc::div(alphaPhi2))*K2 + ( - he2.name() == "e" + he2.name() == thermo2.phasePropertyName("e") ? fvc::div(alphaPhi2, p) : -dalpha2pdt )/rho2