A novel two-stage multiphase model is developed for close-coupled gas atomization by combining formulations available in the literature. Primary atomization is simulated using an Eulerian atomization model, and the outputsa reused as input so fa Lagrangian particle tracking approach to predict the particle size distribution resulting from second aryatomization. Theresults given by the primary atomization mode lare validated with published Direct Numerical Simulations (DNS) values and by comparisonwith experimental images of the spray, while the particle size distributions obtained are accurately fitted toa log-normal distribution, and with powder meandiameters showing good agreemen twith experimental data. The variation of powder characteristic diameters 𝑑10,𝑑50,and𝑑 90 as a function of thegas-to-meltmass flow rate ratio follows correct trends, and the value sofpowder mean diameters 𝑑50 are correctly predicted by the model.