Skip to main content Accessibility help

Extended Hydrodynamic Models and Multigrid Solver of a Silicon Diode Simulation


Extended hydrodynamic models for carrier transport are derived from the semiconductor Boltzmann equation with relaxation time approximation of the scattering term, by using the globally hyperbolic moment method and the moment-dependent relaxation time. Incorporating the microscopic relaxation time and the applied voltage bias, a formula is proposed to determine the relaxation time for each moment equation, which sets different relaxation rates for different moments such that higher moments damp faster. The resulting models would give more satisfactory results of macroscopic quantities of interest with a high-order convergence to those of the underlying Boltzmann equation as the involved moments increase, in comparison to the corresponding moment models using a single relaxation time. In order to simulate the steady states efficiently, a multigrid solver is developed for the derived moment models. Numerical simulations of an n +-n-n + silicon diode are carried out to demonstrate the validation of the presented moment models, and the robustness and efficiency of the designed multigrid solver.


Corresponding author

*Corresponding author. Email addresses: (Z. Hu), (R. Li), (Z. Qiao)


Hide All
[1] Anile, A.M. and Romano, V.. Hydrodynamical modeling of charge carrier transport in semiconductors. Meccanica, 35(3):249296, 2000.
[2] Anile, A.M., Romano, V., and Russo, G.. Extended hydrodynamical model of carrier transport in semiconductors. SIAM J. Appl. Math., 61(1):74101, 2000.
[3] Biegel, B. A. and Plummer, J. D.. Comparison of self-consistency iteration options for the Wigner function method of quantum device simulation. Phys. Rev. B, 54(11):80708082, 1996.
[4] Bløtekjær, K.. Transport equations for electrons in two-valley semiconductors. IEEE Trans. Electron Devices, 17(1):3847, Jan 1970.
[5] Brandt, A.. Multi-level adaptive solutions to boundary-value problems. Math. Comp., 31(138):333390, 1977.
[6] Brandt, A. and Livne, O. E.. Multigrid Techniques: 1984 Guide with Applications to Fluid Dynamics. Classics in Applied Mathematics. SIAM, revised edition, 2011.
[7] Cai, Z., Fan, Y., and Li, R.. Globally hyperbolic regularization of Grad's moment system in one dimensional space. Comm. Math. Sci., 11(2):547571, 2013.
[8] Cai, Z., Fan, Y., and Li, R.. Globally hyperbolic regularization of Grad's moment system. Comm. Pure Appl. Math., 67(3):464518, 2014.
[9] Cai, Z., Fan, Y., Li, R., Lu, T., and Wang, Y.. Quantum hydrodynamics models by moment closure of Wigner equation. J. Math. Phys., 53:103503, 2012.
[10] Cai, Z. and Li, R.. Numerical regularized moment method of arbitrary order for Boltzmann-BGK equation. SIAM J. Sci. Comput., 32(5):28752907, 2010.
[11] Cai, Z., Li, R., and Qiao, Z.. NRxx simulation of microflows with Shakhov model. SIAM J. Sci. Comput., 34(1):A339–A369, 2012.
[12] Cai, Z., Li, R., and Qiao, Z.. Globally hyperbolic regularized moment method with applications to microflow simulation. Computers and Fluids, 81:95109, 2013.
[13] Cai, Z., Li, R., and Wang, Y.. An efficient NRxx method for Boltzmann-BGK equation. J. Sci. Comput., 50(1):103119, 2012.
[14] Cai, Z., Li, R., and Wang, Y.. Numerical regularized moment method for high Mach number flow. Commun. Comput. Phys., 11(5):14151438, 2012.
[15] Carrillo, J.A., Gamba, I.M., and Shu, C.-W.. Computational macroscopic approximations to the one-dimensional relaxation-time kinetic system for semiconductors. Physica D, 146:289306, 2000.
[16] Carrillo, J.A., Gamba, I.M., Majorana, A., and Shu, C.-W.. 2D semiconductor device simulations by WENO-Boltzmann schemes: Efficiency, boundary conditions and comparison to Monte Carlo methods. J. Comput. Phys., 214(1):5580, 2006.
[17] Carrillo, J.A., Gamba, I.M., Muscato, O., and Shu, C.-W.. Comparison of Monte Carlo and deterministic simulations of a silicon diode. In IMA Volumes in Mathematics and Its Applications, pages 7584. Springer-Verlag, 2004.
[18] Cercignani, C., Gamba, I. M., Jerome, J. W., and Shu, C.-W.. Device benchmark comparisons via kinetic, hydrodynamic, and high-hield models. Comput. Methods Appl. Mech. Engrg., 181:381392, 2000.
[19] Cercignani, C., Gamba, I. M., and Levermore, C. D.. High field approximations to a Boltzmann-Poisson system and boundary conditions in a semiconductor. Appl. Math. Lett., 10(4):111117, 1997.
[20] Chen, S., E, W., Liu, Y., and Shu, C.-W.. A discontinuous Galerkin implementation of a domain decomposition method for kinetic-hydrodynamic coupling multiscale problems in gas dynamics and device simulations. J. Comput. Phys., 225(2):13141330, 2007.
[21] Grad, H.. On the kinetic theory of rarefied gases. Comm. Pure Appl. Math., 2(4):331407, 1949.
[22] Grasser, T., Tang, T.-W., Kosina, H., and Selberherr, S.. A review of hydrodynamic and energy-transport models for semiconductor device simulation. Proceedings of The IEEE, 91(2):251274, 2003.
[23] Gummel, H. K.. A self-consistent iterative scheme for one-dimensional steady state transistor calculations. IEEE Trans. Electron Devices, 11(10):455465, 1964.
[24] Hackbusch, W.. Multi-Grid Methods and Applications. Springer-Verlag, Berlin, 1985. second printing 2003.
[25] Hu, Z. and Li, R.. A nonlinear multigrid steady-state solver for 1D microflow. Computers and Fluids, 103:193203, 2014.
[26] Hu, Z., Li, R., Lu, T., Wang, Y., and Yao, W.. Simulation of an n +-n-n + diode by using globally-hyperbolically-closed high-order moment models. J. Sci. Comput., 59(3):761774, 2014.
[27] Imamura, T., Suzuki, K., Nakamura, T., and Yoshida, M.. Acceleration of steady-state lattice Boltzmann simulations on non-uniform mesh using local time step method. J. Comput. Phys., 202(2):645663, Jan 2005.
[28] Jüngel, A.. Transport Equations for Semiconductors. Number 773 in Lecture Notes in Physics. Springer-Verlag, Berlin, 2009.
[29] Kannan, R.. An implicit LU-SGS spectral volume method for the moment models in device simulations: Formulation in 1D and application to a p-multigrid algorithm. Int. J. Numer. Meth. Biomed. Engng., 27:13621375, 2011.
[30] Liu, Y. X. and Shu, C.-W.. Error analysis of the semi-discrete local discontinuous Galerkin method for semiconductor device simulation models. Sci. China Math., 53(12):32553278, 2010.
[31] Lu, T., Du, G., Liu, X., and Zhang, P.. A finite volume method for the multi subband Boltzmann equation with realistic 2D scattering in DG MOSFETs. Commun. Comput. Phys., 10:305338, 2011.
[32] Dal Maso, G., LeFloch, P. G., and Murat, F.. Definition and weak stability of nonconservative products. J. Math. Pures Appl., 74(6):483548, 1995.
[33] Mavriplis, D. J.. Multigrid solution of the steady-state lattice Boltzmann equation. Computers and Fluids, 35:793804, 2006.
[34] McClarren, R. G. and Hauck, C. D.. Robust and accurate filtered spherical harmonics expansions for radiative transfer. J. Comput. Phys., 229(16):55975614, 2010.
[35] Struchtrup, H.. Extendedmoment method for electrons in semiconductors. Physica A: Statistical Mechanics and its Applications, 275(12):229255, 2000.
[36] Tadmor, E.. Convergence of spectral methods for nonlinear conservation laws. SIAM J. Numer. Anal., 26(1):3044, 1989.
[37] Tölke, J., Krafczyk, M., and Rank, E.. A multigrid-solver for the discrete Boltzmann equation. J. Stat. Phys., 107(1/2):573591, 2002.
[38] Van Roosbroeck, W.. Theory of the flow of electrons and holes in germanium and other semiconductors. Bell System Technical Journal, 29(4):560607, 1950.
[39] Vasileska, D., Goodnick, S. M., and Klimeck, G.. Computational Electronics: Semiclassical and Quantum Device Modeling and Simulation. CRC Press, 2010.


MSC classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed