We investigate the convergence of a relaxed version of the best reply numerical schemes (also known as best response or fictitious play) used to find Nash-mean field games equilibriums. This leads us to consider evolution equations in metric spaces similar to gradient flows except that the functional to be differentiated depends on the current point; these are called equilibrium flows. We give two definitions of solutions and prove, through the introduction of a specific index ϒ depending on the trajectory, that, as the time step tends to zero, the interpolated (à la de Giorgi) numerical curves converge to equilibrium flows. As a by-product we obtain a sufficient condition for the uniqueness of a mean field games equilibrium. We close with applications to congestion and vaccination mean field games.