problems couplying pbe-turbulence model
 

I am trying to implement a coalescence model for modeling of multiphase systems where turbulence model is coupled with coalescence model, the solver is based on twoPhaseEulerFoam..my problem is the following: when I run a simulation using kepsilon, the solver doesn't calculate k and neither epsilon :S, it just calculate k-e when I delete ke corrector :
if (pimple.turbCorr())

{
#include "kEpsilon.H"
}
to putting only

#include "kEpsilon.H"

inside the pimple routine.. the problem is that i can't obtain convergence when run my simulations without "if" ..
anybody can help me please??