A finite difference scheme for two-dimensional semiconductor device of heat conduction on composite triangular grids

The momentary state of a semiconductor device of heat conduction is described by a system of four nonlinear partial differential equations. One elliptic equation is for the electrostatic, two parabolic equations are for the electron concentration and the hole concentration, and one heat exchange equation is for the temperature. According to the necessary of practical numerical simulations and based on the balance equation, finite difference schemes for two-dimensional transient behavior of a semiconductor device of heat conduction on composite triangular grids are constructed. Studying their stability and convergence properties, the error estimate in the energy norm is obtained. Finally, a numerical example is given.