Numerical simulations are carried out to look at the primary atomization of a 2-D planar liquid jet. A finite volume method based solver is developed and interface capturing is done by volume of fluid (VOF) method. The solver uses a projection algorithm to solve the governing equations. Preconditioned conjugate gradient method is used to solve the pressure poisson equation. This part of the solver is ported on to graphics processing unit (GPU) to meet the computational demand required. The solver is validated against standard benchmark test cases. Initially the parallelized version on GPU is compared with the serial version on single CPU core to estimate the speed up achieved. Effect of liquid inlet velocity on jet disintegration is studied. © 2015 The Authors. Published by Elsevier B.V.