A numerical study of the spray transfer processes in an electrostatic rotary bell applicator (ESRB) has been conducted utilizing code for a newly developed simulation code. This code consists of three modularized solvers: a fluid flow solver, a spray dynamics solver, and an electrostatic solver. The development of the code consisted of the following steps. First, the flow solver designed for an unsteady three-dimensional Navier-Stokes equation was developed to simulate the shaping airflow with the initial condition and the boundary condition supported by experimental data. Second, the particle trajectory solver, which interacts with the airflow by momentum coupling, was developed to apply the spray transport processes. Finally, the electrostatic solver was developed to calculate the electrostatic field within the two phase flow field. The integrated code created by combining those three solvers was then applied to simulate the paint spray transport processes according to the operating conditions of interest. The numerical results show that the spray shape is very sensitive to changes in the charge to mass ratio. The voltage setting is a dominant operating parameter affecting the numerical transfer efficiency. The voltage range studied was 0 kV to 90 kV. In addition, the transfer efficiency decreases as the shaping airflow rate increases. However, a high shaping airflow rate produces a more uniform distribution of spray mass on the target plane.