An experimentally validated computer simulation model has been developed for the analysis of gas-phase and droplet characteristics of isothermal sprays generated by pressure jet atomizers. Employing a coupled Euler-Lagrange approach for the gas-droplet flow, secondary droplet breakup (based on the ETAB model), was assumed to be dominant and the k-ε model was selected for simulating the gas flow. Specifically, transient spray formation in terms of turbulent gas flow as well as droplet velocities and size distributions are provided for different back pressures. Clearly, two-way coupling of the phases is important because of the impact of significant gas entrainment, droplet momentum transfer, and turbulent dispersion. Several spray phenomena are discussed in light of low back-pressure (1atm) and high back-pressure (30atm) environments. At low back-pressure, sprays have long thin geometric features and penetrate faster and deeper than at high back-pressures because of the measurable change in air density and hence drag force. Away from the nozzle exit under relatively high back pressures, there is no distinct droplet size difference between peripheral and core regions because of the high droplet Weber numbers, leading to very small droplets which move randomly. In contrast to transient spray developments, under steady-state conditions droplets are subject to smaller drag forces due to the fully-developed gas entrainment velocities which reduce gas-liquid slip. Turbulent dispersion influences droplet trajectories significantly because of the impact of random gas-phase fluctuations.