The problem of flapping filaments in an uniform incoming flow is tackled using a Lattice Boltzmann—Immersed Boundary method. The fluid momentum equations are solved on a Cartesian uniform lattice while the beating filaments are tracked through a series of markers, whose dynamics are functions of the forces exerted by the fluid, the filament flexural rigidity and the tension. The instantaneous wall conditions on the filament are imposed via a system of singular body forces, consistently discretised on the lattice of the Boltzmann equation. We first consider the case of a single beating filament, and then the case of multiple beating filaments in a side-by-side configuration, focussing on the modal behaviour of the whole dynamical systems.