Research of Parallel Computing Electromagnetic Field Based on the Domain Decomposition Method

LI Tian-wei,RUAN Jiang-jun,ZHANG Yu,HUANG Dao-chun,YU Shi-feng
DOI: https://doi.org/10.3969/j.issn.1003-6520.2007.05.015
2007-01-01
Abstract:The calculation velocity and memory storage of single PC are often limited in large-scale electromagnetic analysis by finite element method(FEM),parallel processing is an important means to overcome such problem.The domain decomposition method(DDM) which decomposes the whole domain into some sub-domains suiting well for parallel processing was described first.By DDM,local system of linear equations were formed directly on each sub-domain instead of forming the global system of linear equations first and then distributing it.So the parallel efficiency was improved greatly.After the nodes in a sub-domain were assigned to one processor,they are classified into 3 types,internal,border and external to deal with the update in the sub-domain and communications among neighbor sub-domains.A 2-dimensional electrostatic model was built and decomposed by the DDM,the system of linear equations formed by FEM were solved with parallel conjugate gradient method(CG) on the distributed parallel system composed of 6 PC.CG method for system of linear equations is superior in storage and suitable for the sparseness of the obtained high-order matrix which would reduce the arithmetic operations per iteration.But when the condition number of coefficient matrix A is larger than 102,iteration speed will slow down.So it's necessary to decrease condition number of coefficient matrix for increasing the iteration speed by transforming Ax=b into [C-1A(CT)-1]CTx=C-1b,where A was incompletely factorized to A≈CCT and C is sparse the same as A and need no extra memory space to allocate for it.Through analyzing the 2-dimensional simplified capacity model,the effective speed up reaching 97.5% was satisfying.For large-scale simulation which consists of more than millions of freedoms,the parallel processing reduces computing time and increases the computing velocity greatly.It's the base on which large-scale 3D electromagnetic parallel computing.
What problem does this paper attempt to address?