Flow past a single rotating cylinder has been studied, both numerically and experimentally, by many authors. In contrast is the flow past rotating cylinders in a side-by-side arrangement, which has only very limited numerical and experimental data available, despite the pioneering theoretical study by Jeffery 1 as early as 1922. In this paper, we discuss a numerical investigation of the steady laminar viscous flow past two infinite rotating cylinders in a side-by-side configuration. The solution of the compressible two-dimensional Navier-Stokes equations is determined numerically using the high-order spectral difference scheme over an unstructured qudralateral grid. Third order accurate results in both time and space were obtained and compared with existing data. In addition to obtaining a highorder accurate result for flow past two rotating cylinders, we extend the current numerical effort to investigate the effects of Reynolds number, compressibility, and high rotation speed, which have not been comprehensively studied in the past.