×

Application of the shifted-Laplace preconditioner for iterative solution of a higher order finite element discretisation of the vector wave equation: first experiences. (English) Zbl 1203.78044

Summary: The analysis of electromagnetic scattering by electrically large one-side open-ended cavities remains a challenge. Finite element discretisation of the vector wave equation to solve for the electric field inside the cavity leads to ill-conditioned indefinite linear systems of large dimension, a result of the requirement of a fine nearly uniform discrete sampling in the computational domain. Direct methods based on frontal solution techniques to solve the resulting linear system have been used because efficient iterative methods were not available. With the arrival of the shifted-Laplace preconditioner of Erlangga, iterative solution of the indefinite system becomes tractable. This paper discusses the modifications required for application of the shifted-Laplace preconditioner to cavity scattering and some preliminary results of this approach. The shifted-Laplace preconditioner is shown to be very effective for improving the convergence rate of the iterative solution algorithm. However, to be able to handle problems with larger number of degrees of freedom, it is necessary to include a multigrid algorithm to solve the preconditioner system, as this will allow the use of short recurrence Krylov subspace methods, as opposed to the currently employed long-recurrence method, for which the storage requirements of the Krylov basis become unpractically large.

MSC:

78M10 Finite element, Galerkin and related methods applied to problems in optics and electromagnetic theory
78A45 Diffraction, scattering
65F10 Iterative numerical methods for linear systems
65F08 Preconditioners for iterative methods
35Q61 Maxwell equations
Full Text: DOI

References:

[1] Anastassiu, H., A review of electromagnetic scattering analysis for inlets, cavities and open ducts, IEEE Trans. Antennas and Propagation, 45, 6, 27-40 (2003)
[2] Coifman, R.; Rohklin, V.; Wandzura, S., The fast multipole method for the wave equation: A pedestrian description, IEEE Trans. Antennas and Propagation, 35, 3, 7-12 (1993)
[3] Duff, I.; Scott, J., The design of a new frontal code for solving sparse, unsymmetric systems, ACM Trans. Math. Software, 22, 1, 30-45 (1996) · Zbl 0884.65018
[4] Eisenstat, S.; Elman, H.; Schultz, M., Variational iterative methods for nonsymmetric systems of linear equations, SIAM J. Numer. Anal., 20, 345-357 (1983) · Zbl 0524.65019
[5] Erlangga, Y.; Oosterlee, C.; Vuik, C., A novel multigrid based preconditioner for heterogeneous Helmholtz problems, SIAM J. Sci. Comput., 27, 1471-1492 (2006) · Zbl 1095.65109
[6] Erlangga, Y.; Vuik, C.; Oosterlee, C., On a class of preconditioners for solving the Helmholtz equation, Appl. Numer. Math., 50, 409-425 (2004) · Zbl 1051.65101
[7] Graglia, R.; Wilton, D.; Peterson, A., Higher order interpolatory vector bases for computational electromagnetics, IEEE Trans. Antennas and Propagation, 45, 3, 329-342 (1997)
[8] Irons, R., A frontal method solution program for finite element analysis, Internat. J. Numer. Methods Engrg., 2, 5-32 (1970) · Zbl 0252.73050
[9] Jin, J., Electromagnetic scattering from large, deep and arbitrarily shaped open cavities, Electromagnetics, 18, 3-34 (1998)
[10] Jin, J.; Liu, J.; Lou, Z.; Liang, C., A fully high-order finite-element simulation of scattering by deep cavities, IEEE Trans. Antennas and Propagation, 51, 9, 2420-2429 (2003) · Zbl 1368.78139
[11] Liu, J.; Jin, J.-M., A special higher order finite-element method for scattering by deep cavities, IEEE Trans. Antennas and Propagation, 48, 5, 694-703 (2000) · Zbl 1113.78316
[12] Liu, J.; Jin, J.-M., A novel hybridization of higher order finite element and boundary integral methods for electromagnetic scattering and radiation problems, IEEE Trans. Antennas and Propagation, 49, 12, 1794-1806 (2001) · Zbl 1001.78021
[13] Liu, J.; Jin, J.-M., A highly effective preconditioner for solving the finite element-boundary integral matrix equation of 3-D scattering, IEEE Trans. Antennas and Propagation, 50, 9, 1212-1221 (2002)
[14] Liu, J.; Jin, J.-M., Scattering analysis of a large body with deep cavities, IEEE Trans. Antennas and Propagation, 51, 6, 1157-1167 (2003)
[15] Saad, Y., A flexible inner-outer preconditioned GMRES, SIAM J. Sci. Stat. Comput., 14, 461-469 (1993) · Zbl 0780.65022
[18] van der Vorst, H.; Vuik, C., GMRESR: a family of nested GMRES methods, Numer. Linear Algebra Appl., 1, 369-386 (1994) · Zbl 0839.65040
[19] van Gijzen, M. B.; Erlangga, Y. A.; Vuik, C., Spectral analysis of the discrete Helmholtz operator preconditioned with a shifted Laplacian, SIAM J. Sci. Comput., 29, 1942-1958 (2007) · Zbl 1155.65088
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.