The objective of the present paper is to develop a methodology to inject strong wind gusts into the computational domain in order to efficiently simulate their effect on the fluid flow. The design of the methodology based on a source-term formulation takes the feedback effect of the resulting turbulent flow (and, if present, the impacted structure) on the wind gust itself into account. Since the injection of the wind gusts can be carried out close to the region of main interest, CPU-time intensive methods to ensure a proper transport of the gust through the flow field can be avoided. The methodology is mainly intended for the application within eddy-resolving simulations (e.g., LES), but it is not restricted to this class of simulation approaches. For the description of the gusts classical shape functions such as the Extreme Coherent Gust (ECG) and the Extreme Operating Gust (EOG) as well as a newly derived C2-”1-cosine” shape are applied. Two scenarios are taken into account to assess the proposed gust injection technique. On the one hand a (laminar) undisturbed flow field is considered and the effect of different time and length scales of the gusts on their evolution and propagation through the flow field is studied in detail. On the other hand a turbulent background flow is assumed demonstrating that the methodology suggested is also applicable for practically relevant turbulent flows.