Hostname: page-component-74656f9884-8r8nv Total loading time: 0 Render date: 2026-11-10T04:32:05.176Z Has data issue: false hasContentIssue false
  • English
  • Français

Extended Hydrodynamic Models and Multigrid Solver of a Silicon Diode Simulation

Part of: Partial differential equations, boundary value problems Equations of mathematical physics and other areas of application Applications to specific types of physical systems Hyperbolic equations and systems

Published online by Cambridge University Press:  31 August 2016

Zhicheng Hu ,
Ruo Li  and
Zhonghua Qiao
Zhicheng Hu*
Affiliation:
Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Hong Kong
Ruo Li*
Affiliation:
HEDPS & CAPT, LMAM & School of Mathematical Sciences, Peking University, Beijing, China
Zhonghua Qiao*
Affiliation:
Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Hong Kong
*
*Corresponding author. Email addresses:huzhicheng1986@gmail.com (Z. Hu), rli@math.pku.edu.cn (R. Li), zqiao@polyu.edu.hk (Z. Qiao)
*Corresponding author. Email addresses:huzhicheng1986@gmail.com (Z. Hu), rli@math.pku.edu.cn (R. Li), zqiao@polyu.edu.hk (Z. Qiao)
*Corresponding author. Email addresses:huzhicheng1986@gmail.com (Z. Hu), rli@math.pku.edu.cn (R. Li), zqiao@polyu.edu.hk (Z. Qiao)

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the 'Save PDF' action button.

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.

Keywords

Boltzmann transport equationextended hydrodynamic modelmoment-dependent relaxation timemultigridsemiconductor device simulation

MSC classification

Secondary: 35L75: Nonlinear higher-order hyperbolic equations 35Q20: Boltzmann equations 65N55: Multigrid methods; domain decomposition 82D37: Semiconductors

Information

Type
Research Article
Information
Communications in Computational Physics , Volume 20 , Issue 3 , September 2016 , pp. 551 - 582
Copyright
Copyright © Global-Science Press 2016 

References

[1] Anile, A.M. and Romano, V.. Hydrodynamical modeling of charge carrier transport in semiconductors. Meccanica, 35(3):249–296, 2000.CrossRefGoogle Scholar
[2] Anile, A.M., Romano, V., and Russo, G.. Extended hydrodynamical model of carrier transport in semiconductors. SIAM J. Appl. Math., 61(1):74–101, 2000.Google Scholar
[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):8070–8082, 1996.CrossRefGoogle ScholarPubMed
[4] Bløtekjær, K.. Transport equations for electrons in two-valley semiconductors. IEEE Trans. Electron Devices, 17(1):38–47, Jan 1970.CrossRefGoogle Scholar
[5] Brandt, A.. Multi-level adaptive solutions to boundary-value problems. Math. Comp., 31(138):333–390, 1977.Google Scholar
[6] Brandt, A. and Livne, O. E.. Multigrid Techniques: 1984 Guide with Applications to Fluid Dynamics. Classics in Applied Mathematics. SIAM, revised edition, 2011.Google Scholar
[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):547–571, 2013.Google Scholar
[8] Cai, Z., Fan, Y., and Li, R.. Globally hyperbolic regularization of Grad's moment system. Comm. Pure Appl. Math., 67(3):464–518, 2014.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[10] Cai, Z. and Li, R.. Numerical regularized moment method of arbitrary order for Boltzmann-BGK equation. SIAM J. Sci. Comput., 32(5):2875–2907, 2010.CrossRefGoogle Scholar
[11] Cai, Z., Li, R., and Qiao, Z.. NRxx simulation of microflows with Shakhov model. SIAM J. Sci. Comput., 34(1):A339–A369, 2012.CrossRefGoogle Scholar
[12] Cai, Z., Li, R., and Qiao, Z.. Globally hyperbolic regularized moment method with applications to microflow simulation. Computers and Fluids, 81:95–109, 2013.CrossRefGoogle Scholar
[13] Cai, Z., Li, R., and Wang, Y.. An efficient NRxx method for Boltzmann-BGK equation. J. Sci. Comput., 50(1):103–119, 2012.Google Scholar
[14] Cai, Z., Li, R., and Wang, Y.. Numerical regularized moment method for high Mach number flow. Commun. Comput. Phys., 11(5):1415–1438, 2012.CrossRefGoogle Scholar
[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:289–306, 2000.CrossRefGoogle Scholar
[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):55–80, 2006.CrossRefGoogle Scholar
[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 75–84. Springer-Verlag, 2004.Google Scholar
[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:381–392, 2000.Google Scholar
[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):111–117, 1997.Google Scholar
[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):1314–1330, 2007.CrossRefGoogle Scholar
[21] Grad, H.. On the kinetic theory of rarefied gases. Comm. Pure Appl. Math., 2(4):331–407, 1949.Google Scholar
[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):251–274, 2003.Google Scholar
[23] Gummel, H. K.. A self-consistent iterative scheme for one-dimensional steady state transistor calculations. IEEE Trans. Electron Devices, 11(10):455–465, 1964.CrossRefGoogle Scholar
[24] Hackbusch, W.. Multi-Grid Methods and Applications. Springer-Verlag, Berlin, 1985. second printing 2003.Google Scholar
[25] Hu, Z. and Li, R.. A nonlinear multigrid steady-state solver for 1D microflow. Computers and Fluids, 103:193–203, 2014.CrossRefGoogle Scholar
[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):761–774, 2014.CrossRefGoogle Scholar
[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.Google Scholar
[28] Jüngel, A.. Transport Equations for Semiconductors. Number 773 in Lecture Notes in Physics. Springer-Verlag, Berlin, 2009.Google Scholar
[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:1362–1375, 2011.CrossRefGoogle Scholar
[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.CrossRefGoogle Scholar
[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:305–338, 2011.CrossRefGoogle Scholar
[32] Dal Maso, G., LeFloch, P. G., and Murat, F.. Definition and weak stability of nonconservative products. J. Math. Pures Appl., 74(6):483–548, 1995.Google Scholar
[33] Mavriplis, D. J.. Multigrid solution of the steady-state lattice Boltzmann equation. Computers and Fluids, 35:793–804, 2006.CrossRefGoogle Scholar
[34] McClarren, R. G. and Hauck, C. D.. Robust and accurate filtered spherical harmonics expansions for radiative transfer. J. Comput. Phys., 229(16):5597–5614, 2010.CrossRefGoogle Scholar
[35] Struchtrup, H.. Extendedmoment method for electrons in semiconductors. Physica A: Statistical Mechanics and its Applications, 275(12):229–255, 2000.Google Scholar
[36] Tadmor, E.. Convergence of spectral methods for nonlinear conservation laws. SIAM J. Numer. Anal., 26(1):30–44, 1989.Google Scholar
[37] Tölke, J., Krafczyk, M., and Rank, E.. A multigrid-solver for the discrete Boltzmann equation. J. Stat. Phys., 107(1/2):573–591, 2002.Google Scholar
[38] Van Roosbroeck, W.. Theory of the flow of electrons and holes in germanium and other semiconductors. Bell System Technical Journal, 29(4):560–607, 1950.CrossRefGoogle Scholar
[39] Vasileska, D., Goodnick, S. M., and Klimeck, G.. Computational Electronics: Semiclassical and Quantum Device Modeling and Simulation. CRC Press, 2010.Google Scholar