In this paper, we construct a combined upwinding and mixed finite element m
ethod for the numerical solution of a two-dimensional mean field model of s
uperconducting vortices. An advantage of our method is that it works for an
y unstructured regular triangulation. A simple convergence analysis is give
n without resorting to the discrete maximum principle. Numerical examples a
re also presented. Mathematics Subject Classification. 65N99, 82D55.