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

This paper presents some preliminary work on the application of Schwarz alternating method, combined with a multigrid solver, for solving the steady 2D viscous incompressible Navier-Stokes equations. The difficulty caused by the violation of the compatibility condition imposed on the normal velocity at the boundary, which arises when subdomain problems are not solved exactly at each iteration, is surrounded by adding artificial compressiblity constants to the continuity equation each time when a subdomain problem is to be solved. These constants permit to maintain the compatibility condition of the subdomain problems and one can easily show that they tend to zero when the iteration number tends to infinity (in fact, they can be expressed explicitly as functions of the iteration number and some geometrical parameters of the computational domain). Numerical experiences have been made for flows in a channel with a facing step.