A comprehensive three-dimensional model of droplet-gas flow was presented to study the evolution of spray in the effervescent atomization spray with an impinging plate. For gas phase, the N-S equation with the k-epsilon turbulence model was solved, considering two-way coupling interaction between droplets and gas phase. Dispersed droplet phase is modeled as Lagrangian entities, accounting for the physics of droplet generation from primary and secondary breakup, droplet collision and coalescence, droplet momentum and heat transfer. The mean size and statistical distribution of atomized droplets at various nozzle-to-plate distances were calculated. Some simulation results were compared well with experimental data. The results show that the existence of the impinging plate has a pronounced influence on the droplet mean size, size distribution and the droplet spatial distribution. The air-to-liquid ratio has obvious effects on the droplet size and distribution.