Biplane configuration of pure plunging airfoils is investigated in terms of vortex dynamics both experimentally and numerically by utilizing particle image velocimetry and unstructured finite volume solver of incompressible unsteady Navier–Stokes equations. Experiments are carried out to disclose the vortex shedding and interaction mechanisms for various values of frequency and amplitude of the plunging motion. For the non-dimensional plunge amplitude with respect to the chord of airfoil h = 0.2, the effect of the reduced circular frequency based on chord length and the free stream velocity k = 1 and 10 are considered, whereas for h = 0.3, k = 2, 4, 8, and 10 cases are examined. Influence of the plunge amplitude is studied for h = 0.25 at k = 2.5 and for h = 0.0875, 0.15, and 0.3 at k = 4. Numerical simulations are performed to investigate the effect of phase difference on vortex structures and propulsive characteristics, such as thrust and Froude efficiency. Two cases having the highest thrust and efficiency values k = 2.5, h = 0.25, and k = 4, h = 0.15 value couples are selected for the phase angle of [Formula: see text], and [Formula: see text]. Opposed plunge, [Formula: see text], was found as the most efficient amongst all phase angles that were investigated, where [Formula: see text] is beneficial in lift production. Additionally, three-dimensional simulations indicate no significant three dimensionalities for the parameters used herein.