Hydrodynamics of the traditional wedges, cylinders, and cones entering water have been heavily studied. However, the hydrodynamic characteristics of airfoils are practically unknown; the problem is inspired by the recently widely developed fixed-wing unmanned aerial-underwater vehicle, which can fly in the air, cruise in water and frequently enters and exits free water surface. In this study, the hydrodynamic characteristics and fully nonlinear free surface deformation for the water entry of airfoils have been investigated by the time domain boundary element method. First, convergence study based on the NACA0012 airfoil is considered. The time step and element size are adjusted, and the results show that the free surface elevation and pressure distribution agree well with those obtained through Reynolds-averaged Navier–Stokes equations. The results show that the higher the entry velocity, the greater the peak hydrodynamic force. Furthermore, the higher thickness ratio causes greater added mass. Gravity and velocity are found to significantly influence the hydrodynamic behavior, including the time-varying force, while thickness ratio has a great impact on added mass and pressure.