Numerical simulations of airless rotary spray coating process are presented using CFD-EDEM with JKR contact model. The JKR model parameters are suggested according to the virtual test of accumulation angles in the funnel device. The reasonable droplet-wall surface energy is recommended from a relation of accumulation angles. The gas turbulence model is evaluated based on the comparisons of velocities and droplet number fractions using RNG k-ε, realizable k-ε and transition SST models. The simulated film thicknesses agree with measurements in literature. The correlation of exit droplet velocity of sprayer holes is proposed with a discharge coefficient of 2.27. The velocities and film thickness are predicted with changing rotation speeds. The normal force is an order of magnitude larger than the tangential force of droplet-wall collisions, and the ratio of the normal to tangential energy losses of droplet-wall collisions varies from 2.0 to 5.0 in airless rotary spray coating processes.