In this article, we study numerically a diagnostic model, based on mass conservation, to recover solenoidal vector fields from experimental data. Based on a reformulation of the mathematical model as a saddle-point problem, we introduce an iterative preconditioned conjugate gradient algorithm, applied to an associated operator equation of elliptic type, to solve the problem. To obtain a stable algorithm, we use a secondorder mixed finite element approximation for discretization. We show, using synthetic vector fields, that this new approach, yields very accurate solutions at a low computational cost compared to traditional methods with the same order of approximation.