From f7996e45a062ffb62c179a261109fc8d9315b28a Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Mon, 7 Nov 2016 23:14:13 +0000 Subject: [PATCH] twoPhaseEulerFoam::EEqns: Updated pressure-work See commit b5206472b5b5f738352cd035a23d73642e888d1b --- applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H index 25b97c17a1..bd676b55be 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/EEqns.H @@ -16,8 +16,8 @@ - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") - ? fvc::ddt(alpha1)*p + fvc::div(alphaPhi1, p) - : -alpha1*dpdt + ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + : -fvc::ddt(alpha1, p) - alpha1*(dpdt - fvc::ddt(p)) ) - fvm::laplacian @@ -48,8 +48,8 @@ - contErr2*K2 + ( he2.name() == thermo2.phasePropertyName("e") - ? fvc::ddt(alpha2)*p + fvc::div(alphaPhi2, p) - : -alpha2*dpdt + ? fvc::div(fvc::absolute(alphaPhi2, alpha2, U2), p) + : -fvc::ddt(alpha2, p) - alpha2*(dpdt - fvc::ddt(p)) ) - fvm::laplacian