This work investigates the receptivity mechanisms of a NACA0008 airfoil to a level of free-stream turbulence (FST) through a direct numerical simulation (DNS) and an associated linearised simulation on the same mesh. By comparing velocity perturbation fields between the two simulations, the study reveals that the streaky structures that degenerate into turbulent spots are predominantly influenced by nonlinear convective terms, rather than the linear amplification of inflow perturbations around the laminar base flow. A power spectral analysis shows differences in the energy distribution between the DNS and linearised simulation, with the DNS containing more energy at higher wavenumbers, for structures located near the airfoil's leading edge. Representative wavenumbers are identified through modal analysis, revealing a dynamics dominated by streak-like structures. The study employs the Nek5000 numerical solver to distinguish between linear and nonlinear receptivity mechanisms over the NACA0008 airfoil, highlighting their respective contributions to the amplification of perturbations inside the boundary layer. In the high FST case studied, it is observed that the energy of the incoming turbulence is continuously transferred into the boundary layer along the length of the wing. The nonlinear interactions generate streaks with higher spanwise wavenumbers compared with those observed in purely linearised simulations. These thinner streaks align with the spanwise scales identified as susceptible to secondary instabilities. Finally, the procedures presented here generalise the workflow of previous works, allowing for the assessment of receptivity for simulations with arbitrary mesh geometries.