Simulating microstructure evolution during forging of superalloys is of great interest for manufacturer of critical components, particularly in aerospace applications. In this research, a phenomenological approach using a mathematical model was employed into a commercial finite element code, i.e. Abaqus, to study recrystallization and grain growth of IN718 supperalloys during forging. Model validation was performed by compression testing. Results show that the measured recrystallized volume fraction and the grain size were in good agreement with the model predictions. Subsequently, the model was used to simulate the evolution of microstructure in a turbine disc. Finally, a technique based on Tagushi method was used to determine the influence of forging parameters such as forging temperature, ingot height to diameter ratio, and top die velocity on the resultant microstructure. Uniform and refined microstructure of final forging was considered into the objective function. Optimum as well as the poorest parameters combination was attained through analysis. As well the contribution of each parameter on microstructure development was determined through variance analysis.