A Higher-Order Stabilized Hybridized Discontinuous Galerkin Method for Simulating Semiconductor Devices

• Published in: IEEE journal on multiscale and multiphysics computational techniques Vol. 10; pp. 283 - 294
• Main Authors: Zhang, Nian-En, Zhao, Dongyan, Feng, Haoqiang, Wang, Yin-Da, Chen, Yanning, Wang, Qi-Chao, Du, Zheng-Wei, Liang, Yingzong, Liu, Fang, Xie, Hao, Zhan, Qiwei, Yin, Wen-Yan
• Format: Journal Article
• Language: English
• Published: Piscataway IEEE 2025; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Accuracy; Approximation algorithms; Carrier density; Carrier transport; Charge carrier density; Differential equations; Drift diffusion model; Electric potential; Electronic devices; Finite element analysis; Galerkin method; higher order method; hybridized discontinuous galerkin method; Junctions; Mathematical analysis; Mathematical models; Method of moments; Nonlinear equations; Numerical models; Numerical stability; Operators (mathematics); Semiconductor devices
• ISSN: 2379-8815; 2379-8815
• DOI: 10.1109/JMMCT.2025.3575845

The simulation of carrier transport in power electronic devices imposes stringent requirements on numerical stability, confining the previous methods to low-order schemes. To address this issue, a stabilized higher-order hybridized discontinuous Galerkin method (S-HDG) is proposed, where we decouple the exponentially varying carrier density from the differential operator and project it onto a lower-dimensional equation. Based on the numerical jumps as indicator, an adaptive artificial diffusion term is introduced to dynamically control oscillatory errors and over diffusion during the iterations for solving nonlinear equations. We validate the proposed method to abrupt junction models, demonstrating its high-order accuracy and robustness against severe mesh skewness and curvature. Furthermore, we apply the method to lateral double-diffused MOSFET (LDMOS), a class of typical power electronic devices, achieving good agreement with the industrial-standard FVSG solver in simulating electrical parameters. Notably, our method can offer higher-order convergence and better compatibility with unstructured meshes.