The study of innovative energy systems often involves complex fluid flows problems and the Computational Fluid-Dynamics (CFD) is one of the main tools of analysis. It is important to put in evidence that in several energy systems the flow field experiences the laminar-to-turbulent transition. Direct Numerical Simulations (DNS) or Large Eddy Simulation (LES) are able to predict the flow transition but they are still inapplicable to the study of real problems due to the significant computational resources requirements. Differently standard Reynolds Averaged Navier Stokes (RANS) approaches are not always reliable since they assume a fully turbulent regime. In order to overcome this drawback in the recent years some locally formulated transition RANS models have been developed. In this work, we present a local correlation–based transition approach adding two equations that control the laminar-toturbulent transition process –γ and – to the well–known Spalart–Allmaras (SA) turbulence model. The new model was implemented within OpenFOAM code. The energy equation is also implemented in order to evaluate the model performance in thermal–fluid dynamics applications. In all the considered cases a very good agreement between numerical and experimental data was observed.