Abstract: The Lattice Boltzmann Method (LBM) has been seen as an alternative tool for the computational simulation of fluid dynamics. In this study, we use the LBM with Single-Relaxation-Time (SRT) collision model to simulate two-dimensional (2D) laminar flow past a square cylinder. The main aim of the study is to systematically investigate the influences of the locations of the inflow, outflow and side walls boundaries, where Reynolds number is kept at 100 for all calculations. The side wall boundary locations will be analyzed using the periodic and symmetric boundary conditions. Analyses of the relaxation time parameter also have been investigated. Some physical quantities, such as the drag coefficient, mean drag coefficient, root mean square values of lift coefficient and the Strouhal number are examined for the purpose. We also examined the vortex shedding formation which provides an excellent means of visualizing the von Karman vortex street. We found that there is a certain range for inflow, outflow and side wall boundaries where physical quantities such as drag and mean drag coefficients, root mean square value of lift coefficient and the Strouhal number show some changes. Results also show that, there is a change for physical quantities when the relaxation time parameter is changed from certain range and also effect the computational time. The physical quantities are obtained and compared with other existing experimental and numerical results.