Abstracts - Dipartimento di Elettronica Applicata
Abstracts - Dipartimento di Elettronica Applicata
Abstracts - Dipartimento di Elettronica Applicata
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Meta 2010 & FEM 2010 – Rome, 13-15 December 2010<br />
Some comments on the solution of the linear<br />
algebraic systems defined by the finite element<br />
method when applied to electromagnetic<br />
problems involving bianisotropic me<strong>di</strong>a<br />
Paolo Fernandes (1) , Marina Ottonello (2) and Mirco Raffetto (2)<br />
(1) Istituto <strong>di</strong> Matematica <strong>Applicata</strong> e Tecnologie Informatiche del<br />
Consiglio Nazionale Delle Ricerche,<br />
Via De Marini 6, I 16149 Genoa, Italy – E-mail: fernandes@ge.imati.cnr.it<br />
(2) Department of Biophysical and Electronic Engineering,<br />
University of Genoa, Via Opera Pia 11a,<br />
I 16145, Genoa, Italy – E-mail: ottonello@<strong>di</strong>be.unige.it,<br />
raffetto@<strong>di</strong>be.unige.it<br />
It is well known that the finite element (FE) method, when applied to the<br />
solution of time-harmonic electromagnetic boundary value problems involving<br />
linear me<strong>di</strong>a, requires the solution of linear algebraic systems of equations<br />
characterized by very sparse matrices [1]. For this task <strong>di</strong>fferent techniques of<br />
solution can be considered: <strong>di</strong>rect solvers and iterative solvers [1]. Usually<br />
iterative solvers [2] work very well and are used worldwide in FE simulations<br />
of electromagnetic boundary value problems. Many types of iterative solvers<br />
have been developed [2] and, among these, some solvers like the conjugate<br />
gra<strong>di</strong>ent [1], [2], the biconjugate gra<strong>di</strong>ent [1], [2] and the Conjugate<br />
Orthogonal Conjugate Gra<strong>di</strong>ent [3], which have received a particular attention<br />
in the research community working on the FE method, assume some kind of<br />
symmetry of the FE matrices.<br />
When the time-harmonic electromagnetic boundary value problem of interest<br />
involves innovative and linear me<strong>di</strong>a, like the so-called bianisotropic materials<br />
[4], [5], the usual symmetry of the FE matrices can be lost and for the most<br />
widely used algebraic iterative algorithms convergence is not guaranteed<br />
anymore. This is the reason why in [6] the authors considered the me<strong>di</strong>um<br />
analyzed in [5] setting a parameter to zero, so reducing the bianisotropic<br />
me<strong>di</strong>um to a biisotropic one. In this contribution we analyze what can be done<br />
to overcome this <strong>di</strong>fficulty.<br />
References<br />
[1] J. Jin, The finite element method in electromagnetics, John Wiley & Sons, 1993.<br />
[2] R. Barrett, M. Berry, T. F. Chan, J. Demmel, J. Donato, J. Dongarra, V. Eijkhout, R.<br />
Pozo, C. Romine and H. Van der Vorst, Templates for the Solution of Linear Systems:<br />
Buil<strong>di</strong>ng Blocks for Iterative Methods, 2nd E<strong>di</strong>tion, SIAM, 1994.<br />
[3] H. A. van der Vorst and J. B. M. Melissen, “A Petrov-Galerkin type method for solving<br />
Ax=b, where A is symmetric complex,” IEEE Trans. Magnetics, vol. 26, pp. 706-708,<br />
1990.<br />
[4] J. A. Kong, Theory of Electromagnetic Waves, Wiley, 1975.<br />
[5] S. Maruyama and M. Koshiba, “A vector finite element formulation for general<br />
bianisotropic waveguides,” IEEE Transactions on Magnetics, vol. 33, pp. 1528- 1531,<br />
1997.<br />
[6] P. Fernandes and M. Raffetto, “Well posedness and finite element approximability of<br />
time-harmonic electromagnetic boundary value problems involving bianisotropic<br />
materials and metamaterials, Mathematical Models and Methods in Applied Sciences,<br />
vol. 19, pp. 2299-2335, December 2009.<br />
62