Abstract A two degree-of-freedom (2-DOF) galloping piezoelectric energy harvesting system where both oscillators are excited by the incoming flow is studied in this paper. A coupled lumped-parameter model is developed to simulate the electro-fluid-structural coupling behaviour assuming a quasi-steady aerodynamic flow field. The differential equations governing the dynamics of the lumped system are converted into nonlinear algebraic equations employing the harmonic balance method (HBM) and then solved using Newton’s method. The approximate analytical solutions are compared to the solutions obtained by numerical integration, and the results agree, both qualitatively and quantitatively. The solutions were subsequently used to investigate the effects of the bluff body dimension, mechanical, electrical, aerodynamic, and electromechanical parameters on the performance of the harvester and to illustrate how to optimise some of these parameters to reduce the cut-in wind speed and maximise the power output of the energy harvester. It will be shown that the energy harvesting performance could be significantly improved by introducing piezoelectric transducers on both oscillators and including both masses in the incoming flow. A comparative study between the proposed 2-DOF flow energy harvesting system, a single-degree-of-freedom (SDOF) galloping piezoelectric energy harvester (GPEH), and other configurations of 2-DOF GPEH shows that the present 2-DOF GPEH generates the largest power peak.