In this study, a thermal energy storage system is modeled as an axisymmetric ventilated cavity partially filled with a porous medium that is subject to turbulent flow. Local thermal non-equilibrium is considered to analyze the heat transport and the turbulent k−ε model is used to account for the high inlet flow velocities. Finite volume method is employed for discretization of the equations that are relaxed with the SIMPLE method. The investigation is concerned with the turbulence field, heat fluxes and pressure drop in the system relative to variations in Reynolds number (from 8.3 × 103 to 5 × 104), porosity (ϕ from 0.6 to 0.8) and Da number (from 4x10−5 to 4x10−7). Results indicate that porosity effects slight changes in the turbulence field while the temperatures increase significantly faster for higher porosity cases. Turbulence in higher porosity cases was lower, accompanying the increase in thermal efficiencies. Variations in Da number with fixed porosity showed that, for lower Da number, recirculation in the clear region increased while for higher Da number porous region turbulence increased. Also, higher Da decreased heat exchange between fluid and solid phases. Finally, an increase in thermal efficiency for lower Da number flows was followed by increased average turbulent kinetic energy and relative pressure drop.