The flow-induced snap-through dynamics of side-by-side buckled flexible filaments in a uniform flow were investigated using the penalty immersed boundary method. We systematically examined the effects of gap distance and bending rigidity on the filament's motion and energy harvesting performance. Four distinct modes were observed as we varied the aforementioned parameters: the dual equilibrium mode, the dual streamwise oscillation mode, the single snap-through oscillation (STO) mode, and the dual STO mode. We analyzed the corresponding wake patterns associated with each of these modes. A phase transition from out-of-phase to inphase motion was observed at low gap distances and bending rigidities. At lower gap distances, the interaction between two filaments was pronounced, leading to a co-motion of the filaments and enhanced elastic energy. This interaction decreased as the gap distance increased. Filaments in the 2STO mode exhibited significantly higher strain energy compared to the other modes. The elastic energy was concentrated at the rear part of the filaments and increased with higher gap distances and bending rigidities. Considering both energy harvesting efficiency and space utilization rate, a gap distance equal to 1 exhibited the best performance.