PurposeGrowing application of micro aerial vehicle (MAV) sets in demand for accurate computations of low Reynolds number flows past their wings. The purpose of this study is to investigate the effect of unsteady freestream velocity or wind gust on a harmonically plunging symmetric NACA0012 airfoil at Re = 1,000. The influence of unsteady parameters, such as reduced frequency of plunging motion (0.25 < k < 1.5), non-dimensional plunging amplitude (ho = 0.2) and non-dimensional amplitude of wind gust (0.1 = λ = 0.4) has been studied.Design/methodology/approachComputations have been carried out using commercial software ANSYS Fluent 16.0. To incorporate the plunging motion, the entire reference frame is oscillating, and thereby, a source term is added in the Navier–Stokes equation.FindingsThe results have been presented in the form of streamlines, vorticity contours, lift and drag signals and their spectra. It is observed that the ratio of plunging frequency to gust frequency (f/fg) has strong influence on periodic characteristics of unsteady wake. It has also been observed that for a fixed plunging amplitude, an increase in value of k results into a change from positive drag to thrust.Practical implicationsThe research has implications in the development of MAV.Originality/valueThis study is intended to get a better understanding of unsteady parameters associated with gusty flow in flapping wing applications and possible ways to alleviate its adverse effect on it.