The aerodynamic performance of wind turbine airfoils is crucial for the efficiency and reliability of wind energy systems, with leading-edge roughness significantly impacting blade performance. This study conducts numerical simulations on the DU 00-W-401 airfoil to investigate the effects of leading-edge roughness. Results reveal that the rough airfoil exhibits a distinctive “N”-shaped lift coefficient curve. The formation mechanism of this nonlinear lift curve is primarily attributed to the development of the trailing-edge separation vortex and variations in the adverse pressure gradient from the maximum thickness position to the trailing-edge confluence. The impact of different roughness heights is further investigated. It is discovered that when the roughness height is higher than 0.3 mm, the boundary layer can be considered fully turbulent, and the lift curve shows the “N” shape stably. When the roughness height is between 0.07 mm and 0.1 mm, a transitional state can be observed, with several saltation points in the lift curve. The main characteristics of different flow regimes based on different lift curve segments are summarized. This research enhances the understanding of the effects of leading-edge roughness on the aerodynamic performance of a thick wind turbine airfoil, and the simulation method for considering the effect of leading-edge roughness is practical to be applied on large-scale wind turbine blade to estimate the aerodynamic performance under rough leading-edge conditions, thereby supporting advancements in wind turbine technology and promoting the broader adoption of renewable energy.