diff --git a/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H b/applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H index 926abe2a09..d5a5d20631 100644 --- 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); }