A new volume of fluid (VOF) based two-phase flow solver has been developed that has been parallelized for general purpose graphical processor units (GPGPU). Incompressible Navier-Stokes equation is explicitly solved using finite volume method on a 2-D Cartesian grid. The pressure Poisson equation is solved on the GPU in conjunction with the geometric multigrid method. This solver performance is established by validating it against standard benchmark test cases. Furthermore, the solver acceleration is determined by comparing it with the serial version of the code. The solver was then used to study the primary jet breakup of a 2-D planar liquid jet due to a co-flowing gas stream. The effect of liquid to gas phase velocity ratio on jet breakup was investigated. © 2015 Elsevier Ltd.