An analysis is performed to study the two-phase natural convection flow of nanofluid along a vertical wavy surface. The model includes equations expressing conservation of total mass, momentum and thermal energy for two-phase nanofluid. Primitive variable formulations (PVF) are used to transform the dimensionless boundary layer equations into a convenient coordinate system and the resulting equations are integrated numerically via implicit finite difference iterative scheme. The effect of controlling parameters on the dimensionless quantities such as skin friction coefficient, rate of heat transfer and rate of mass transfer is explored. It is concluded from the present analysis, that the diffusivity ratio parameter, NA and particle-density increment number, NB have pronounced influence on the reduction of heat transfer rate. (C) 2017 Elsevier Ltd. All rights reserved.