Finite volume method was employed for analyzing the three-dimensional turbulent flow structures, species distributions, and mixing behaviors of combustion flows in a rotary kiln under various operation conditions. The modified £e-£`turbulence model together with wall functions was adopted. Devolatilization of solid wastes were simulated by gaseous methane (CH4) non-uniformly distributed along the kiln bed. Combustion process was considered as a two-step reaction when primary air entered and mixed with methane gas in the first combustion chamber. Mixing-controlled eddy-dissipation model was employed for predicting the reaction rates of CH4, O2, CO2, CO and H2O. Effects of inleakage air, kiln rotation speed and methane distribution along the kiln bed were also examined. Results show that 128% excess air will get the best combustion efficiency, above which the combustion efficiency will decrease. The temperature and species are not uniformly distributed and are vertically stratified on cross-sectional plane. The combustion efficiency will also be lowered if there is inleakage airflow. Results also show rotation speed and methane distributions have little effect on combustion efficiency.