From 3b65f924ec4dfc77e7c0d449098fe0acf0aebcef Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Sat, 18 Apr 2015 21:40:54 +0100
Subject: [PATCH] twoPhaseEulerFoam: Added fvOptions correction to the phase
 velocities after matrix construction

---
 applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
index 6c765ec66b..abd3c8b646 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@@ -31,6 +31,7 @@ volScalarField dragCoeff(fluid.dragCoeff());
         U1Eqn += fvm::Sp(dragCoeff, U1);
         fvOptions.constrain(U1Eqn);
         U1.correctBoundaryConditions();
+        fvOptions.correct(U1);
     }
 
     {
@@ -54,5 +55,6 @@ volScalarField dragCoeff(fluid.dragCoeff());
         U2Eqn += fvm::Sp(dragCoeff, U2);
         fvOptions.constrain(U2Eqn);
         U2.correctBoundaryConditions();
+        fvOptions.correct(U2);
     }
 }