Three-dimensional time-harmonic eddy current problems solved by the geometric multigrid preconditioned conjugate gradient method

The focus of this study is on the efficient solution for three-dimensional (3D) time-harmonic eddy current problems discretised by the finite-element method (FEM) with thin elements. The systems of equations arising from the finite-element formulations are solved by the geometric multigrid preconditioned conjugate gradient (MGCG) method. Numerical examples show that the MGCG method has stable convergence, which does not deteriorate for elements with high aspect ratios and its computation time is much less than that of the conjugate gradient method with incomplete Cholesky factorisation as the preconditioner (ICCG).