On the convergence of Gummel's numerical algorithm

Abstract The algorithm described by Gummel for computing carrier distributions in semiconductors is analyzed by a simple perturbation method. It is concluded that in one-dimensional problems, convergence failures may be expected under forward bias conditions when the ratio of the carrier diffusion lengths to the device dimension is sufficiently small.