The Ensemble Monte Carlo method has been applied to the study of the relaxation dynamics of coupled free-carriers and excitons in bulk Wurtzite GaN, after non-resonant optical excitation. The semiclassical Boltzmann equations for electron and hole populations are completed by an additional equation for the exciton distribution, and are coupled by reaction terms, describing exciton formation and dissociation mediated by phonons. The temporal evolution in the range from 1 to 400 ps show that the exciton formation process and the temporal evolution of the resulting population are sensitive to the excitation energy. The rates for the acoustic piezoelectric interaction have been calculated, including screening effects. Results show that piezo-acoustic rates are bigger in the wurtzite phase, respect to the cubic one, and they are really sensitive to the background doping of the sample.