diff --git a/applications/solvers/cfdemSolverIBRO/UEqn.H b/applications/solvers/cfdemSolverIBRO/UEqn.H
index af4edf22..a20a5986 100644
--- a/applications/solvers/cfdemSolverIBRO/UEqn.H
+++ b/applications/solvers/cfdemSolverIBRO/UEqn.H
@@ -5,7 +5,7 @@ fvVectorMatrix UEqn
   + turbulence->divDevReff(U)
  ==
     fvOptions(U)
- +  (lambda*(1-voidfractionNext)/U.mesh().time().deltaT())*(Us - fvm::Sp(1,U))
+ +  (lambda*(1-voidfractionNext)/U.mesh().time().deltaT())*(fvc::Sp(1,Us)-fvm::Sp(1,U))
 );
 
 UEqn.relax();
diff --git a/src/lagrangian/cfdemParticle/derived/cfdemCloudIBmodified/cfdemCloudIBmodified.C b/src/lagrangian/cfdemParticle/derived/cfdemCloudIBmodified/cfdemCloudIBmodified.C
index 6bbec0f3..c120228f 100644
--- a/src/lagrangian/cfdemParticle/derived/cfdemCloudIBmodified/cfdemCloudIBmodified.C
+++ b/src/lagrangian/cfdemParticle/derived/cfdemCloudIBmodified/cfdemCloudIBmodified.C
@@ -177,7 +177,7 @@ void cfdemCloudIBmodified::calcForcingTerm(volVectorField& Us)
                     // calc particle velocity
                     //for(int i=0;i<3;i++) rVec[i]=Us.mesh().C()[cell][i]-position(index)[i];
                     for(int i=0;i<3;i++) uP[i] = velocities()[index][i];
-                    Us[cell] = uP;//(1-voidfractions_[index][subCell])*uP;
+                    Us[cell] = (1-voidfractions_[index][subCell])*uP;
                 }
             }
         }