In order to reduce the gust response level of the aircraft wing aeroelastic systems in supersonic flows, the nonlinear energy sink (NES) is introduced into the aeroelastic system, and the effect of the nonlinear energy sink on the aeroelastic gust response is studied. The gust response dynamic model without and with NES is established by Lagrange equation, considering the cubic nonlinear characteristics of both plunge and pitch stiffness. The piston theory model as well as the '1-cos' model are employed to describe the unsteady aerodynamic load and the gust velocity, respectively. Besides, the comparison study of the gust response amplitude and energy characteristics of the wing aeroelastic systems with and without NES are analyzed. For systems with NES, the amplitude of gust dynamic response decreases rapidly, while vibration can be transmitted to the NES mechanism and causing its vibration. From an energy perspective, about 72.15 % of the original vibration energy has been transferred into NES, and further consumed through its damping. The influence of NES parameters, such as frequency ratio, mass ratio, damping ratio, and installation position on gust response, is discussed. Numerical simulations demonstrate that: as the frequency ratio, damping ratio and mass ratio increase, the dynamic response amplitude of the aeroelastic system decays faster, and the energy transferred from the wing system also increases. In addition, an optimal relative distance for the installation of NES can be found with best gust alleviation effect. The influence of nonlinear energy sink on the targeted transfer of vibration energy in aeroelastic systems are analyzed under different working conditions. The NES shows great gust alleviation effect and good robustness under different flight speeds, gust intensities, and gust lengths.