The present study provides a numerical method for the atomization and evaporation in liquid jet flows with Volume of Fluid Method (VOF) and Ghost Fluid Method (GFM). Based on mass, momentum and energy conservation equations, the heat and mass transfer on the gas-liquid interface are calculated. The motion of the sharp interface is simulated by the VOF method. The GFM method solves the discretizing problem of the interface with velocity jump conditions. In order to verify the accuracy of the numerical method, the numerical and theoretical results for the droplet evaporation are compared. The method is employed to simulate the atomization and evaporation in liquid jet flows. According to numerical results, the breakup process is analyzed for four different stages. The discussion for the effect of the phase change on the jet atomization is given. The results show that the increase of the ambient gas temperature will help to stabilize the motion of the interface and delay the breakup of the liquid jet. (C) 2019 Elsevier Ltd. All rights reserved.