The low-Reynolds number flow around two square cylinders placed side-by-side is investigated using the lattice Boltzmann method (LBM). The effects of the gap ratio s/d (s is the separation between the cylinders and d is the characteristic dimension) on the flow are studied. These simulations reveal the existence of regimes with either synchronized or non-synchronized vortex-shedding, with transition occurring at s/d approximate to 2, which is larger than for circular cylinders. Detailed results are presented at Re = 73 for s/d = 2.5 and 0.7 corresponding to the synchronized and flip-flop regimes, respectively. Vortex-shedding from the cylinder occurs either in-phase or in-antiphase in the synchronized regime. However, linear stochastic estimate (LSE) calculations show that in-phase locking is the predominant mode. LSE is also employed to educe the underlying modes in the flip-flop regime, where evidence for both in-phase and anti-phase locked vortices is found, indicating that this regime is in a quasi-stable state between these two modes. The merging of the wakes, which is gradual for the synchronized regime, occurs rapidly in the flip-flop regime. The mean pressure on the upstream surface is symmetric and asymmetric for the synchronized and flip-flop regimes, respectively. Differences in results between the two regimes are interpreted in terms of the interaction of the jet formed between the cylinders with the adjoining wakes, the strength of this interaction depending on the spacing. (c) 2005 Elsevier Ltd. All rights reserved.