14.01.2013 Views

Abstracts - Dipartimento di Elettronica Applicata

Abstracts - Dipartimento di Elettronica Applicata

Abstracts - Dipartimento di Elettronica Applicata

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!