Abstract: A machine-learned model enhances the accuracy of turbulence transport equations of RANS solver and applied for periodic hill test case. The accuracy is investigated in extrapolation modes. A parametric study is also performed to understand the effect of network hyperparameters on training and model accuracy and to quantify the uncertainty in model accuracy due to the non-deterministic nature of the neural network training. For any network, less than optimal mini-batch size results in overfitting, and larger than optimal reduces accuracy. Data clustering is an efficient approach to prevent the machine-learned model from over-training on more prevalent flow regimes, and results in a model with similar accuracy. Turbulence production is correlated with shear strain in the free-shear region, with shear strain and wall-distance and local velocity-based Reynolds number in the boundary layer regime, and with streamwise velocity gradient in the accelerating flow regime. The flow direction is key in identifying flow separation and reattachment regime. Machine-learned models perform poorly in extrapolation mode. A priori tests reveal model predictability improves as the hill dataset is partially added during training in a partial extrapolation model. These also provide better turbulent kinetic energy and shear stress predictions than RANS in a posteriori tests. Before a machine-learned model is applied for a posteriori tests, a priori tests should be performed.