Investigations of vortex dynamics about two circular cylinders in a side-by-side arrangement help the understanding of flows around more complex structures, which are found to have many engineering applications. These applications involve offshore structures, power generation, micro-turbine engines, cooling towers, and paper machine forming fabrics, etc. Therefore, two-dimensional compressible laminar flows over two cylinders in side-by-side arrangement are numerically investigated at low Reynolds number. The high-order discontinuous Galerkin method is employed to simulate the flow, which combines the advantages associated with finite element and finite volume methods. As in classical finite element method, the spatial accuracy can be obtained by the high-order polynomial approximation within an element rather than by stencils as in finite volume method. The curved triangle is used to represent the wall boundary of cylinder to maintain the high-order accurate simulation. Then the characteristics of the wake flow are identified by capturing the vortex structure. After verifying the rationality of the method, the influences of gap spacing on vortex shedding and mechanical characteristics are analyzed. The results reveal that the flow depends to a large extent on the gap spacing between the two cylinders, which can change the vortex shedding pattern. At the gap spacing S*=1.1, wake flow pattern resembles the vortex street of a single bluff body. The flow in the gap is too weak to affect the wake pattern, leading to the complete suppression of vortices shed on the gap sides of both cylinders. At the gap spacing S*=1.4, the results reveal that the gap flow is deflected from one cylinder to another. Meanwhile, the wakes represent randomly flip-flopping between two states of the gap flow direction, which is called the flip-flopping wake pattern. The flow is no longer periodic but becomes drastically unsteady. Anti-symmetric flow pattern is predicted for gap spacing S*=2.5, indicating that two parallel vortex streets are anti-symmetric with respect to the centerline. With further increasing the gap spacing to S*=4, the symmetric flow pattern is observed. Furthermore, the flow preserves its structure very far downstream without any distortion. With the increase of the cylinder spacing, the average drag coefficients are declined significantly, and the absolute value of average lift coefficient decreases simultaneously. The Reynolds number has a little influence on the average drag coefficient. As the Reynolds number increases, the average lift coefficient decreases, while the vortex shedding frequency increases.