Reliable and Efficient Numerical Methods for Time Harmonic Electromagnetic Design Problems

Loading...
Thumbnail Image
Journal Title
Journal ISSN
Volume Title
School of Electrical Engineering | Doctoral thesis (article-based) | Defence date: 2014-10-24
Checking the digitized thesis and permission for publishing
Instructions for the author
Date
2014
Major/Subject
Mcode
Degree programme
Language
en
Pages
63 + app. 54
Series
Aalto University publication series DOCTORAL DISSERTATIONS, 155/2014
Abstract
This dissertation studies numerical methods for solving time-harmonic Maxwell's equations. In this work, boundary element and a least squares type finite element method for time-harmonic Maxwell's equations is analysed. The boundary element method is studied from the point of view of shape optimization and shape sensitivity analysis. The main result of the dissertation is the description how to calculate numerically the derivative of electromagnetic properties of a metallic scatterer with respect to variation of its shape. The motivation of the result is the adjoint variable method aimed at computing the derivative of a cost function of a generic optimization problem. The method has the property that, essentially, it suffices to differentiate the system matrix with respect to shape if the load functional or the cost functional do not depend on the shape explicitly. The boundary element method is also applied in as wide as possible characterisation of the shapes of wire dipoles and that way in shape optimisation. The least-squares finite-element method for time-harmonic Maxwell's equations results in a method where the system matrix is positive definite despite the underlying physics those of time-harmonic fields usually lead to strictly indefinite systems. The positive definiteness of the system matrix is inherited from the ellipticity of the sesquilinear form of the variational method. The implementation of a conformal discretization is not, however, always reasonable for the reason that in certain geometries the finite element solution vanishes on the boundary when only tangential, or normal, component is required to vanish. Therefore, the discrete solution space is widened in such a manner that on mesh dense enough the sesquilinear form is still elliptic and that the tangential, or normal, component vanishes in suitable weak sense.

Tämä työ käsittelee numeerisia menetelmiä aikaharmonisten sähkömagneettisten kenttien simuloinnissa. Työssä kehitetään reunaelementtimenetelmää sekä analysoidaan miten voidaan ratkaista Maxwellin aikaharmoniset yhtälöt duaalinormaaliyhtälöinä pienimmän neliösumman mielessä elementtimenetelmällä. Reunaelementtimenetelmää työssä tarkastellaan muodon optimoinnin ja herkkyysanalyysin näkökulmasta. Työn päätulos kertoo miten lasketaan numeerisesti metallisten sirottajien sähköisten ominaisuuksien derivaatta sirottajan geometrian muutoksen suhteen. Tuloksen kannustimena on yleisen optimointitehtävän kustannusfunktion derivaatan laskentaan kehitetty abstrakti dualimuuttujamenetelmä. Menetellä on ominaisuus, joka takaa, että olennaisesti riittää laskea systeemimatriisin alkioiden muotoderivaatta kun ollaan kiinnostuneita optimointitehtävän kustannusfunktion derivaatasta, kun lastifunktionaali ja kustannusfunktio eivät riipu eksplisiittisesti sirottajan muodosta. Reunaelementtimenetelmää sovelletaan työssä myös lankadipolien muodon mahdollisimman laajan parametriavaruuden karakterisointiin, ja sitä kautta muodon optimointiin. Pienimmän neliösumman elementtimenetelmä aikaharmonisille Maxwellin yhtälöille johtaa menetelmään, jossa systeemimatriisi on positiividefiniitti vaikka tavallisesti aikaharmoniset yhtälöt johtavat aidosti indefiniitteihin järjestelmiin. Systeemimatriisin positiividefiniittiys periytyy variaatiotehtävän seskvilineaarimuodon elliptisyydestä. Konformin diskretoinnin toteuttaminen ei kuitenkaan ole aina järkevää, sillä tietynlaisissa alueissa se johtaa elementtiratkaisuihin, jotka häviävät reunalla kokonaan kun vain tangentiaalisen tai normaalin osan tulisi hävitä. Työssä lavennetaan ratkaisuavaruuttaa sopivasti siten, että tarpeeksi tiheällä verkolla seskvilineaarimuoto on edelleen elliptinen ja elementtiratkaisu toteuttaa reunaehdon heikommassa mielessä.
Description
Supervising professor
Nikoskinen, Keijo, Prof., Aalto University, Department of Radio Science and Engineering, Finland
Thesis advisor
Ylä-Oijala, Pasi, PhD, Aalto University, Department of Radio Science and Engineering, Finland
Keywords
optimization, EFIE, boundary element method, finite element method, optimointi, reunaelementtimenetelmä, elementtimenetelmä
Other note
Parts
  • [Publication 1]: J. Kataja, and K. Nikoskinen. The parametric optimization of wire dipole antennas. IEEE Transactions on Antennas and Propagation, vol. 59, no. 2, pp. 350–356, Feb. 2011.
    DOI: 10.1109/TAP.2010.2096181 View at publisher
  • [Publication 2]: J. Kataja, A.G. Polimeridis, J.R. Mosig and P. Yla-Oijala. Analytical shape derivatives of the MFIE system matrix discretized with RWG functions. IEEE Transactions on Antennas and Propagation, vol. 61, no. 2, pp. 985–988, Feb. 2013.
    DOI: 10.1109/TAP.2012.2223447 View at publisher
  • [Publication 3]: J. Kataja and J.I. Toivanen. On shape differentiation of discretized electric field integral equation. Engineering Analysis with Boundary Elements, vol. 37, no. 9, pp. 1197–1203, Sept. 2013.
    DOI: 10.1016/j.enganabound.2013.05.006 View at publisher
  • [Publication 4]: J. Kataja, S. Jarvenpaa, Jukka I. Toivanen, Raino A.E. Makinen, Pasi Yla-Oijala. Shape sensitivity analysis and gradient based optimization of large structures using MLFMA. IEEE Transactions on Antennas and Propagation, accepted, 2014.
    DOI: 10.1109/TAP.2014.2347394 View at publisher
  • [Publication 5]: J. Kataja. A least squares finite element method for the extended Maxwell system. Progress in Electromagnetic Research M, vol. 33, pp. 137–151, 2013.
    DOI: 10.2528/PIERM13080702 View at publisher
Citation