A computational analysis has been performed to study the flow instability of two-parallel wall motions in a Cuboidal enclosure incorporated by a cylinder under different radii sizes. A numerical methodology based on the Finite Volume Method (FVM) and a full Multigrid acceleration is utilized in this paper. Left and right parallel walls of the cavity are maintained driven and all the other walls completing the domain are motionless. Different radii sizes (R=0.075, 0.1, 0.125, 0.15 and 0.175) are employed encompassing descriptive Reynolds numbers that range three orders of magnitude 100, 400 and 800 for the steady state. The obtained results show positions R=0.15 and R=0.175 of the inner cylinder promote cell distortion. Also, when the radius equates to R=0.15, it may lead to the birth of tertiary cells at Re=400 which are more developed for Re=800. Thereafter, analysis of the flow evolution shows that with increasing Re beyond a certain critical value, the flow becomes unstable and undergoes a Hopf bifurcation. A nonuniform variation with the radius size of the inner cylinder is observed. Otherwise said, elongating the radius of the cylinder leads to decrease in the critical Reynolds number. Hence, the acceleration of the unsteadiness is realized. On the other hand, by further increasing Reynolds number more than the critical value from 1200 to 2100, we note that the kinetic energy is monotonously increasing with Reynolds number and a stronger motion in the velocity near the rear wall of the enclosure is observed. Furthermore, the symmetry of flow patterns observed in the steady state has been lost. Therefore, a systematic description of various effects illuminating the optimum geometrical parameters to achieve effective flow behavior in those systems has been successfully established through this paper.