We discuss the basics of modeling semiconductors and semiconductor devices. Drift-diffusion equations and the implementation issues associated with their nonlinear solvers will be overviewed. Difficulties in modeling interfaces with offsets in primary variables will be mentioned. The talk relates to our research within theĀ larger effort to improve solar efficiency through a hypothesized process known as Heterojunction Assisted Impact Ionization.