A Mixed-finite Volume Element Coupled with the Method of Characteristic Fractional Step Difference for Simulating Transient Behavior of Semiconductor Device of Heat Conductor And Its Numerical Analysis

1 Institute of Mathematics, Shandong University, Jinan 250100, China; 2 School of Mathematical Sciences, Shandong Normal University, Jinan 250014, China; 3 School of Economics, Shandong University, Jinan 250100, China

Abstract The mathematical system is formulated by four partial differential equations combined with initialboundary value conditions to describe transient behavior of three-dimensional semiconductor device with heat conduction. The first equation of an elliptic type is defined with respect to the electric potential, the successive two equations of convection dominated diffusion type are given to define the electron concentration and the hole concentration, and the fourth equation of heat conductor is for the temperature. The electric potential appears in the equations of electron concentration, hole concentration and the temperature in the formation of the intensity. A mass conservative numerical approximation of the electric potential is presented by using the mixed finite volume element, and the accuracy of computation of the electric intensity is improved one order. The method of characteristic fractional step difference is applied to discretize the other three equations, where the hyperbolic terms are approximated by a difference quotient in the characteristics and the diffusion terms are discretized by the method of fractional step difference. The computation of three-dimensional problem works efficiently by dividing it into three one-dimensional subproblems and every subproblem is solved by the method of speedup in parallel. Using a pair of different grids (coarse partition and refined partition), piecewise threefold quadratic interpolation, variation theory, multiplicative commutation rule of differential operators, mathematical induction and priori estimates theory and special technique of differential equations, we derive an optimal second order estimate in L^{2}-norm. This numerical method is valuable in the simulation of semiconductor device theoretically and actually, and gives a powerful tool to solve the international problem presented by J. Douglas, Jr.

Fund:Project supported by the National Natural Science Foundation of China (Grant Nos. 11101124 and 11271231), the National Tackling Key Problems Program for Science and Technology (Grant No. 20050200069), and the Doctorate Foundation of the Ministry of Education of China (Grant No. 20030422047).

Yi-rang YUAN,Qing YANG,Chang-feng LI等. A Mixed-finite Volume Element Coupled with the Method of Characteristic Fractional Step Difference for Simulating Transient Behavior of Semiconductor Device of Heat Conductor And Its Numerical Analysis[J]. Acta Mathematicae Applicatae Sinica, English Serie, 2017, 33(4): 1053-1072.

Bank, R.E., Fichtner, W.M., Rose, D.J., et. al. Transient simulation of sillcon devices and circuits. IEEE Computer-Aided Design, 6:436-451(1985)

[2]

Ciariet, P.G. The finite element methods for elliptic problem. North-Holland, Amsterdam, 1978

[3]

Douglas, Jr. J. Simulation of miscible displacement in porous media by a modified method of characteristics procedure. In:Numerical Analysis, Dundee, 1981, Lecture Note in Mathematics 912, Springer-Verlag, Berlin, 1982

[4]

Douglas, Jr. J., Ewing, R.E., Wheeler, M.F. Approximation of the pressure by a mixed method in the simulation of miscible displacement. RAIRO Anal. Numer., 17(1):17-33(1983)

[5]

Douglas, Jr.J., Yuan, Y.R. Numerical simulation of immiscible flow in porous media based on combining the method of characteristics with mixed finite element procedure. Numerical Simulation in Oil Recovery, 119-132, New York:Springer-Verlag, 1986

[6]

Douglas, Jr.J., Yuan, Y.R. Finite difference methods for transient behavior of a semiconductor device. Mat. Apli. Comp., 6(1):25-38(1987)

[7]

Gummel, H.K. A self-consistent iterative scheme for one-dimensional steady-state transistor calculation. IEEE Trans:Electron Device, 11:455-465(1964)

[8]

He Y., Wei T.F. Computer simulation method for semiconductor device. Scicence Press, Beijing, 1989

[9]

Jerome, J.W. Mathematical Theory and Approximation of Semiconductor Models. Philadelphia, SIAM, 1994

[10]

Jiang L.S., Pang Z.Y. Finite element method and its theory. People's Education Press, Beijing, 1979

[11]

Jones, J.E. A mixed finite volume method for accurate computation of fluid velocities in porous media. Ph. D. Thesis, University of Colorado, Donrer. Co., 1995

[12]

Lou, Y. On basic semiconductor equation with heat conduction. J. Partial. Diff. Eqs., 1:43-54(1995)

[13]

Nitsche, J. Linear splint-funktionen and die von Rits for elliptishce randwort probleme. Arch. for Rational Mech. and Anal., 36:348-355(1968)

[14]

Russell, T.F. Time stepping along characteristics with incomplete iteration for a Galerkin approximation of miscible displacement in porous media. SIAM J. Numer. Anal., 22(5):970-1013(1985)

[15]

Russell, T.F. Rigorous black-centered discritization on irregular grids:Improved simulation of complex reservoir systems. Project Report, Research Corporation, Tulsa, 1995

[16]

Shi, M. Physics of modern semiconductor device. Science Press, Beijing, 2002

[17]

Sun, C.W., Lu, Q.S., Fan, Z.X. Laser irradiation effect. National Defence Industry Press, Beijing, 2002

[18]

Sun, T.J., Yuan, Y.R. An approximation of incompressible miscible displacement in porous media by mixed finite element method and characteristics-mixed finite element method. J. Comp. Appl. Math., 228:391-411(2009)

[19]

Weiser, A., Wheeler, M.F. On convergence of block-centered finite difference for elliptic problems. SIAM J. Numer. Anal., 25(2):351-375(1988)

[20]

Yang, Q, Yuan, Y.R. An approximation of semiconductor device by mixed finite element method and characteristics-mixed finite element method. Applied Mathematics and Computation, 225:407-424(2013)

[21]

Yang, Q., Yuan, Y.R. An approximation of semiconductor device of heat conduction by mixed finite element method characteristics-mixed finite element method. Applied Numerical Mathematics, 70:42-57(2013)

[22]

Yang, Q., Yuan, Y.R. An approximation of three-dimensional semiconductor device by mixed finite element method and characteristics-mixed finite element method. Numerical Mathematics:Theory, Methods and Application, 8(3):356-382(2015)

[23]

Yuan, Y.R., Ding, L.Y., Yang, H. A new method and theoretical analysis of numerical analog of semiconductor. Chinese Science Bulletin, 27(7):790-795(1982)

[24]

Yuan, Y.R. Characteristics method with mixed finite element for transient behavior of semiconductor device. Chin. Sci. Bull., 36(17):1356-1357(1991)

[25]

Yuan Y.R. The approximation of the electronic potential by a mixed method in the simulation of semiconductor. J. Systems Sci. Math. Sci., 11(2):117-120(1991)

[26]

Yuan, Y.R. Time stepping along characteristics for the finite element approximation of compressible miscible displacement in porous media. Math. Numer. Sinica, 14(4):385-400(1992)

[27]

Yuan, Y.R. Finite difference methods for a compressible miscible displacement problem in porous media. Math. Numer. Sinica, 15(1):16-28(1993)

[28]

Yuan, Y.R. Finite element method and analysis of numerical simulation of semiconductor device. Acta Math. Sci., 13(3):241-251(1993)

[29]

Yuan, Y.R. Finite difference method and analysis for three-dimensional semiconductor device of heat conduction. Sci. China Math., 39(11):1140-1151(1996)

[30]

Yuan, Y.R. Characteristic finite difference fractional step methods for three-dimensional semiconductor device of heat conduction. Chin. Sci. Bull., 45(2):123-131(2000)

[31]

Yuan, Y.R. Finite difference fractional step method for transient behavior of a semiconductor device. Acta Mathematica Scientia (Series B), 3:427-438(2005)

[32]

Yuan, Y.R. Modification of upwind finite difference fractional step methods by the transient state of the semiconductor device. Numer. Methods Partial Differential, 24:400-417(2008)

[33]

Yuan, Y.R. Recent progress in numerical methods for semiconductor devices. Chinese J. Computational Physics, 26(3):317-324(2009)

[34]

Yuan, Y.R. Theory and application of reservoir numerical simulation, Chapter 7, numercial method of transient behavior of semiconductor device. Science Press, Beijing, 2013