The design of low-pressure turbines (LPTs) must account for the losses generated by the unsteady interaction with the upstream blade row. The estimation of such unsteady wake-induced losses requires the accurate prediction of the incoming wake dynamics and decay. Existing linear turbulence closures (stress-strain relationships), however, do not offer an accurate prediction of the wake mixing. Therefore, machine-learnt, nonlinear turbulence closures (models) have been developed for LPT flows with unsteady inflow conditions using a zonal-based model development approach, with an aim to enhance the wake mixing prediction for unsteady Reynolds-averaged Navier-Stokes calculations. High-fidelity time-averaged and phase-lock averaged data at a realistic isentropic Reynolds number and two reduced frequencies, i.e., with discrete incoming wakes and with wake “fogging,” have been used as reference data for a machine learning algorithm based on gene expression programing to develop models. Models developed via phase-lock averaged data were able to capture the effect of certain prominent physical phenomena in LPTs such as wake-wake interactions, whereas models based on the time-averaged data could not. Correlations with the flow physics lead to a set of models that can effectively enhance the wake mixing prediction across the entire LPT domain for both cases. Based on a newly developed error metric, the developed models have reduced the a priori error over the Boussinesq approximation on average by 45%. This study thus aids blade designers in selecting the appropriate nonlinear closures capable of mimicking the physical mechanisms responsible for loss generation. Copyright © 2019 by ASME