The leading-edge bulges along the extension direction are designed on the marine wingsail. The height and the spanwise wavelength of the protuberances are 0.1c and 0.25c, respectively. At Reynolds number Re=5×10⁵, the Reynolds Averaged Navier-Stokes equations are applied to the simulation of the wingsail with the bulges thanks to ANSYS Fluent finite-volume solver based on the SST K-ω models. The grid independence analysis is carried out with the lift and drag coefficients of the wingsail at AOA = 8° and AOA=20°. The results show that while the efficiency of the wingsail is reduced by devising the leading-edge bulges before stall, the bulges help to improve the lift coefficient of the wingsail when stalling. At AOA=22° under the action of the leading-edge tubercles, a convective vortex is formed on the suction surface of the modified wingsail, which reduces the flow loss. So the bulges of the wingsail can delay the stall.