Numerical simulations for type II superconductors