In this work, an atomizer with a de Laval-type nozzle is designed and studied by commercial computational fluid dynamics (CFD) software, and the secondary breakup process during atomization is simulated by two-way coupling and the discrete particle model (DPM) using the Euler-Lagrange method. The simulation result demonstrates that the gas flow patterns greatly change with the introduction of liquid droplets, which clearly indicates that the mass loading effect is quite significant as a result of the gas-droplet interactions. An hourglass shape of the cloud of disintegrating molten metal particles is observed by using a stochastic tracking model. Finally, this simulation approach is used for the quantitative evaluation of the effects of altering the atomizing process conditions (gas-to-melt ratio, operating pressure P, and operating gas temperature T) and nozzle geometry (protrusion length h, half-taper angle α, and gas slit nozzle diameter D) on the particle size distribution of the powders produced.