A numerical investigation on the stalled flow characteristics of a NACA0021 aerofoil with a sinusoidal wavy leading edge (WLE) at Re ∞ = 1.2×105 and α = 20• is presented in this paper. It is observed that laminar separation bubbles (LSBs) form at the trough areas of the WLE in a collocated fashion rather than uniformly/periodically distributed over the span. It is found that the distribution of LSBs and its influence on the aerodynamic forces is strongly dependent on the spanwise domain size of the simulation, i.e. the wavenumber of the WLE used. The creation of a pair of counter-rotating streamwise vortices from the WLE and their evolution as an interface/buffer between the LSBs and the adjacent fully separated shear layers are discussed in detail. The current simulation results confirm that an increased lift and a decreased drag are achieved by using the WLEs compared to the straight leading edge (SLE) case, as observed in previous experiments. Additionally, the WLE cases exhibit a significantly reduced level of unsteady fluctuations in aerodynamic forces at the frequency of periodic vortex shedding. The beneficial aerodynamic characteristics of the WLE cases are attributed to the following three major events observed in the current simulations: (1) the appearance of a large low-pressure zone near the leading edge created by the LSBs; (2) the reattachment of flow behind the LSBs resulting in a decreased volume of the rear wake; and, (3) the deterioration of von-Kármán (periodic) vortex shedding due to the breakdown of spanwise coherent structures.