BUG: Corrected mechanical energy in pEqn of compressibleTwoPhaseEulerFoam - mantis #624

2d419cc7 · andy · 0dae57fb · 2d419cc7
Commit 2d419cc7 authored 12 years ago by andy
--- a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H
@@ -178,7 +178,7 @@
    rho2 = rho20 + psi2*p;

    K1 = 0.5*magSqr(U1);
-    K2 = 0.5*magSqr(U1);
+    K2 = 0.5*magSqr(U2);

    dpdt = fvc::ddt(p);
 }