一种可扩展的三维半导体器件并行数值模拟算法

计算物理 ›› 2012, Vol. 29 ›› Issue (3): 439-448.

一种可扩展的三维半导体器件并行数值模拟算法

成杰, 张林波

中国科学院数学与系统科学研究院科学与工程计算国家重点实验室, 北京, 100190

收稿日期:2011-07-28 修回日期:2011-12-02 出版日期:2012-05-25 发布日期:2012-05-25
作者简介:成杰(1984-),男,北京,博士生,从事并行有限元算法研究,E-mail:chengjie@lsec.cc.ac.cn
基金资助:
国家973计划(2011CB309703);国家自然科学基金(60873177,11171334);创新群体基金(11021101)资助项目

A Scalable Parallel Algorithm for Three-dimensional Semiconductor Device Simulation

CHENG Jie, ZHANG Linbo

State Key Laboratory of Scientific and Engineering Computing, Academy of Mathematics and Systems Science, CAS, Beijing 100190, China

Received:2011-07-28 Revised:2011-12-02 Online:2012-05-25 Published:2012-05-25

摘要/Abstract

摘要： 在基于漂移-扩散模型的三维半导体器件数值模拟中,通过有限体积法进行数值离散,采用完全耦合的牛顿迭代求解非线性代数方程组,并使用基于代数多重网格预条件子的GMRES方法求解牛顿迭代中的线性方程组,构造一种稳健且高度可扩展的非结构四面体网格上求解半导体方程的并行算法.基于PHG平台实现该算法的并行计算程序,并对PN结和MOS场效应晶体管等问题进行了最大网格规模达到5亿单元、最大并行规模达到1 024进程的大规模数值模拟实验,结果表明,该算法计算效率高,可扩展性好.

关键词: 三维半导体器件, 漂移-扩散模型, 并行数值模拟, 有限体积方法, 代数多重网格方法

Abstract: A scalable parallel algorithm for three-dimensional semiconductor device simulation on unstructured tetrahedral meshes using drift-diffusion model is proposed,which is characterized by finite volume discretization,fully coupled Newton iterations for discretized nonlinear equations,and GMRES iterations using algebraic multigrid(AMG) preconditioner for linear equations in Newton iterations.The algorithm was implemented using a parallel adaptive finite element toolbox PHG.Large scale parallel numerical experiments with several problems,including PN diode and MOSFET,were carried out.In numerical simulations the largest mesh size are 500 million elements and the largest number of MPI processes used are 1 024.It shows that the proposed algorithm is efficient,robust and highly scalable.

Key words: semiconductor device, drift-diffusion, finite volume, 3D parallel simulation, algebraic multigrid

中图分类号:

O241.82

成杰, 张林波. 一种可扩展的三维半导体器件并行数值模拟算法[J]. 计算物理, 2012, 29(3): 439-448.

CHENG Jie, ZHANG Linbo. A Scalable Parallel Algorithm for Three-dimensional Semiconductor Device Simulation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 29(3): 439-448.

[1]	张淑娴, 杭旭登. 二维扩散方程的局部二阶线性保极值节点计算方法及应用[J]. 计算物理, 2022, 39(1): 17-32.
[2]	王宏亮, 田宙, 郭永辉. 高温高压条件下堆积颗粒的吸能效应[J]. 计算物理, 2014, 31(6): 706-712.
[3]	贾祖朋, 蔚喜军. 基于近似Riemann解的有限体积ALE方法[J]. 计算物理, 2007, 24(5): 543-549.
[4]	高智, 向华, 申义庆. 摄动有限体积法重构近似高精度的意义[J]. 计算物理, 2004, 21(2): 131-136.
[5]	侯慧清, 张莉, 金希卓. 无结构网格有限体积方法及在热对流中的应用研究[J]. 计算物理, 2000, 17(4): 381-387.