Turek-hron FSI model continuity convergence
 

Hello, I've been trying to verify this Turek-Hron benchmark for FSI simulation, The velocity convergence in x,y and z directions is going well but it is giving me huge values in terms of continuity residuals. How can I reduce continuity residuals for this model? I have tried to refine the mesh and change the solution methods (SIMPLE,SIMPLEC,PISO...) but still having same results. I'm attaching some pictures of residuals and mesh here.

I would appreciate any comments.

thanks.