Fish schools and bird flock formation is an interesting phenomenon that could emerge from specific hydrodynamic interactions. Aerodynamics for self-propulsion of a single heaving NACA0017 airfoil has been numerically investigated to understand these interactions. Simulations using Ansys Fluent have been performed to examine the effect of heaving frequency variation on the maximum velocity downstream of the airfoil and the wavelength of fluctuating wake velocity. Thereafter, tandem heaving airfoil simulations are performed, and stable equilibrium positions between the airfoils are derived. The stable equilibrium position is found to be a function of the heaving frequency and amplitude.