This paper presents the stationary Bingham fluid flow simulations past a circular cylinder placed in a channel. The governing equations of motion are discretized using the mixed finite element method. The biquadratic element Q 2 is used to approximate the velocity space and P 1 disc for the pressure space. This LBB-stable finite element pair Q 2 =P disc 1 leads to an accuracy of third and second order in L 2-norm for velocity and pressure approximations, respectively. The discrete version of the nonlinear system is linearized by the Newton's method, and the linear subproblems as an inner loop are coped with UMFPACK solver. Code validation and grid independence study is also presented to confirm the reliability of the results. Simulations are carried out for various values of dimensionless Bingham number Bn, and its impact is investigated with great detail by demonstrating the velocity, viscosity and pressure contours. The benchmark quantities like drag coefficient (C D) and lift coefficient (C L) are also discussed graphically and quantitatively. Moreover, two different locations of the circular cylinder are considered for the analysis of hydrodynamic forces.