Temperature distribution in a energy cell significantly affects the performance and effectiveness of the energy cell system. Particularly, in low temperature energy cells, enhancement of the system requires proper thermal operation, which indicates the need for developing accurate thermal models. In this study, a 3D numerical thermal model is presented to dissect the heat transfer and prognosticate the temperature distribution in air cooled proton exchange membrane energy cells (PEMFC). In the modeled energy cell mound, forced air flux inventories oxidant as well as cooling. Conservation equations of mass, instigation, and energy are answered in the oxidant channel, while energy equation is answered in the entire sphere, including the gas prolixity layers and division plates, which play a significant part in heat transfer. Parametric studies are performed to probe the goods of colorful parcels and operating conditions on the maximum cell temperature. The present results are farther validated with trial. This model provides a theoretical foundation for thermal analysis of air- cooled PEMFC heaps, where temperature non- uniformity is high and thermal operation and mound cooling is a significant challenge.