Spray formation in a liquid jet in air crossflow is investigated numerically using a hybrid approach for two cases: We = 11 and We = 53, where We is the air Weber number. The VOF method was used to solve the interaction between the liquid jet column and the air flow up to the primary breakup. The Eulerian liquid portion is converted to discrete Lagrangian drops by the use of source terms of mass and momentum and according to empirical correlations. The effects of two secondary breakup models, as well as two-way coupling with droplet-to-droplet collisions were evaluated in this work, in order to establish an improved hybrid model suitable to solve liquid jet in crossflow at low computational cost. The results were very satisfactory relating droplet velocity and mass fraction distribution for the lower Weber number case. For the higher Weber number case, the deviation of the numerical mass fraction distribution in comparison with the experimental one was less than 14 x 10(-2). Improvement to the droplet size distribution was observed by using a new secondary breakup model, specially for the higher Weber number case. (C) 2019 Elsevier Ltd. All rights reserved.