期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 163, 期 1, 页码 150-171出版社
ACADEMIC PRESS INC
DOI: 10.1006/jcph.2000.6545
关键词
vector potential; Helmholtz decomposition; Coulomb gauge; Maxwell's equations; solution discontinuities; finite volume; Krylov space methods; preconditioning
We consider solving three-dimensional electromagnetic problems in parameter regimes where the quasi-static approximation applies, the permeability is constant, the conductivity may vary significantly, and the range of frequencies is moderate. The difficulties encountered include handling solution discontinuities across interfaces and accelerating convergence of traditional iterative methods for the solution of the linear systems of algebraic equations that arise when discretizing Maxwell's equations in the frequency domain. We use a potential-current formulation (A, phi, (J) over circle) with a Coulomb gauge. The potentials A and phi decompose the electric field E into components in the active and null spaces of the del x operator. We develop a finite volume discretization on a staggered grid that naturally employs harmonic averages for the conductivity at cell faces. After discretization, we eliminate the current and the resulting large, sparse, linear system of equations has a block structure that is diagonally dominant, allowing an efficient solution with preconditioned Krylov space methods. A particularly efficient algorithm results from the combination of BICGSTAB and an incomplete LU-decomposition. We demonstrate the efficacy of our method in several numerical experiments. (C) 2000 Academic Press.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据