Purpose - The purpose of this paper is to develop a three-dimensional (31)) numerical model capable of predicting the vaporization rate of a liquid fuel droplet exposed to a convective turbulent airflow at ambient room temperature and atmospheric pressure conditions. Design/methodology/approach - The 3D Reynolds-Averaged Navier-Stokes equations, together with the mass, species, and energy conservation equations were solved in Cartesian coordinates. Closure for the turbulence stress terms for turbulent flow was accomplished by testing two different turbulence closure models; the low-Reynolds number (LRN) k-epsilon and shear-stress transport (SST). Numerical solution of the resulted set of equations was achieved by using blocked-off technique with finite volume method. Findings - The present predictions showed good agreement with published turbulent experimental data when using the SST turbulence closure model. However, the LRN k-epsilon model produced poor predictions. In addition, the simple numerical approach employed in the present code demonstrated its worth. Research limitations/implications - The present study is limited to ambient room temperature and atmospheric pressure conditions. However, in most practical spray flow applications droplets evaporate under ambient high-pressure and a hot turbulent environment. Therefore, an extension of this study to evaluate the effects of pressure and temperature will make it more practical. Originality/value - It is believed that the numerical code developed is of great importance to scientists and engineers working in the field of spray combustion. This paper also demonstrated for the first time that the simple blocked-off technique can be successfully used for treating a droplet in the flow calculation domain.