Optimization of flank wear width (VB) progression during face milling of Inconel 718 is challenging due to the synergistic effect of cutting parameters on the complex wear mechanisms and failure modes. The lack of quantitative understanding between VB and the cutting conditions limits the development of the tool life extension. In this study, a Gaussian kernel ridge regression was employed to develop the VB progression model for face milling of Inconel 718 using multi-layer physical vapor deposition-TiAlN/NbN-coated carbide inserts with the input feature of cutting speed, feed rate, axial depth of cut, and cutting length. The model showed a root mean square error of 30.9 (49.7) μm and R2 of 0.93 (0.81) in full fit (5-fold cross-validation test). The statistics along with the cross-plot analyses suggested that the model had a high predictive ability. A new promising condition at the cutting speed of 40 m/min, feed rate of 0.08 mm/tooth, and axial depth of cut of 0.9 mm was designed and experimentally validated. The measured and predicted VB agreed well with each other. This model is thus applicable for VB prediction and optimization in the real face milling operation of Inconel 718.