The present work aims to investigate the recirculation and incipient mixing zones in a channel flow supplied with obstacles. The main objective is to develop a new technique to control these recirculation zones by setting a variable roughness. For the purpose of varying that roughness, 4 small bars of heights 0.25H, 0.5H, 0.75H and H were placed downstream of the obstacle; H is the height of the obstacle. For this, a three-dimensional numerical approach was carried out using the ANSYS CFX computer code. In addition, the governing equations were solved using the finite volume method. The K-ω shear-stress transport (SST) turbulence model was utilized to model the turbulent stresses. In the end, we presented the time-averaged simulation results of the contours of the current lines (3D time-averaged streamlines, trace-lines), three components of the velocities: <u> (velocity u contour), <v> (velocity v contour) and <w> (velocity w contour), trace-lines, stream ribbons and mean Q-criterion iso-surface.