Asymptotically preserving numerical methods for large reaction diffusion systems with applications to solar cell design