<p>Waveguides are used to transmit electromagnetic signals.Their geometry is typically long and slender their particularshape can be used in the design of computational methods. Onlyspecial modes are transmitted and eigenvalue and eigenvectoranalysis becomes important.</p><p>We develop a .nite-element code for solving theelectromagnetic .eld problem in closed waveguides .lled withvarious materials. By discretizing the cross-section of thewaveguide into a number of triangles, an eigenvalue problem isderived. A general program based on Arnoldi’s method andARPACK has been written using node and edge elements toapproximate the .eld. A serious problem in the FEM was theoccurrence of spurious solution, that was due to impropermodeling of the null space of the curl operator. Therefore edgeelements has been chosen to remove non physical spurioussolutions that arises.</p><p>Numerical examples are given for homogeneous andinhomogeneous waveguides, in the homogeneous case the resultsare compared to analytical solutions and the right order ofconvergence is achieved. For the more complicated inhomogeneouswaveguides with and without striplines, comparison has beendone with results found in literature together with gridconvergence studies.</p><p>The code has been implemented to be used in an industrialenvironment, together with full 3-D time and frequency domainsolvers. The2-D simulations has been used as input for full3-D time domain simulations, and the results have been comparedto what an analytical input would give.</p>

Waveguides are used to transmit electromagnetic signals.Their geometry is typically long and slender their particularshape can be used in the design of computational methods. Onlyspecial modes are transmitted and eigenvalue and eigenvectoranalysis becomes important. We develop a .nite-element code for solving theelectromagnetic .eld problem in closed waveguides .lled withvarious materials. By discretizing the cross-section of thewaveguide into a number of triangles, an eigenvalue problem isderived. A general program based on Arnoldi’s method andARPACK has been written using node and edge elements toapproximate the .eld. A serious problem in the FEM was theoccurrence of spurious solution, that was due to impropermodeling of the null space of the curl operator. Therefore edgeelements has been chosen to remove non physical spurioussolutions that arises. Numerical examples are given for homogeneous andinhomogeneous waveguides, in the homogeneous case the resultsare compared to analytical solutions and the right order ofconvergence is achieved. For the more complicated inhomogeneouswaveguides with and without striplines, comparison has beendone with results found in literature together with gridconvergence studies. The code has been implemented to be used in an industrialenvironment, together with full 3-D time and frequency domainsolvers. The2-D simulations has been used as input for full3-D time domain simulations, and the results have been comparedto what an analytical input would give. / NR 20140805

In this thesis, quadratic hexahedral edge elements have been applied to the three dimensional for open region electromagnetic scattering problems. For this purpose, a semi-automatic all-hexahedral mesh generation algorithm is developed and implemented. Material properties inside the elements and along the edges are also determined and prescribed during the mesh generation phase in order to be used in the solution phase. Based on the condition number quality metric, the generated mesh is optimized by means of the Particle Swarm Optimization (PSO) technique. A framework implementing hierarchical hexahedral edge elements is implemented to investigate the performance of linear and quadratic hexahedral edge elements. Perfectly Matched Layers (PMLs), which are implemented by using a complex coordinate transformation, have been used for mesh truncation in the software. Sparse storage and relevant efficient matrix ordering are used for the representation of the system of equations. Both direct and indirect sparse matrix solution methods are implemented and used. Performance of quadratic hexahedral edge elements is deeply investigated over the radar cross-sections of several curved or flat objects with or without patches. Instead of the de-facto standard of 0.1 wavelength linear element size, 0.3-0.4 wavelength quadratic element size was observed to be a new potential criterion for electromagnetic scattering and radiation problems.

Ces dernières années, la Méthode Intégrale de Volume (MIV) a reçu une attention particulière pour lamodélisation des problèmes électromagnétiques en basse fréquence. Son intérêt principal est l’absencedu maillage de la région air, ce qui rend la méthode légère et rapide. Associée aux méthodes decompression matricielle la MIV devient aujourd'hui une alternative compétitive à la méthode deséléments finis pour la modélisation de dispositifs électromagnétiques ayant un volume d'airprépondérant.Ce rapport porte sur le développement de deux formulations intégrales de volume pour la résolution deproblèmes magnétostatiques avec prise en compte des matériaux non linéaires, des aimants, desbobines, des circuits magnétiques avec ou sans entrefer et des régions minces magnétiques. Lapremière est une formulation en flux de mailles indépendantes basée sur l'interpolation par éléments defacette. La deuxième est une formulation en potentiel vecteur magnétique basée sur l'interpolation paréléments d'arête. L'application de ces formulations permet d’une part d'obtenir des résultats précismême en présence d’un faible maillage et d’autre part de résoudre aisément des problèmes nonlinéaires. Des méthodes de calcul de la force magnétique globale ainsi que du flux magnétique dansles bobines ont été également mises en oeuvre. Les développements informatiques ont été réalisés dansla plateforme MIPSE et ont été validés sur des problèmes académiques ainsi que sur quelquesdispositifs industriels. / In recent years, the Volume Integral Method (VIM) has been received particular attention formodeling of low frequency electromagnetic problems. The main advantage of this method is thatinactive regions do not to be discretized, which makes it light and rapid. Associated with matrixcompression methods, the VIM is a competitive alternative to the finite element method for modelingelectromagnetic devices containing a predominant air volume.This PhD thesis focuses on the development of two volume integral formulations for solvingmagnetostatic problems, in the presence of nonlinear materials, magnets, coils, multiply connectedmagnetic regions, and the presence of magnetic shielding. The first one is a mesh magnetic fluxformulation based on the interpolation of facet elements and the second one is a magnetic vectorpotential formulation based on the interpolation of edge elements. The application of theseformulations provides accurate results even with coarse meshes and allows solving straightforwardnonlinear magnetostatic problems. Methods for computing global magnetic force and magnetic fluxthrough a coil were also implemented as part of this work. Developments performed in the MIPSEplatform were validated on academic case-tests as well as some industrial devices.

Les équations de Maxwell en régime harmonique comportent plusieurs difficultés lorsque la fréquence est élevée. On peut notamment citer le fait que leur formulation variationnelle n’est pas définie positive et l’effet de pollution qui oblige à utiliser des maillages très fins, ce qui rend problématique la construction de solveurs itératifs. Nous proposons une stratégie de solution précise et rapide, qui associe une discrétisation par des éléments finis d’ordre élevé à des préconditionneurs de type décomposition de domaine. La conception, l’implémentation et l’analyse des deux méthodes sont assez difficiles pour les équations de Maxwell. Les éléments finis adaptés à l’approximation du champ électrique sont les éléments finis H(rot)-conformes ou d’arête. Ici nous revisitons les degrés de liberté classiques définis par Nédélec, afin d’obtenir une expression plus pratique par rapport aux fonctions de base d’ordre élevé choisies. De plus, nous proposons une technique pour restaurer la dualité entre les fonctions de base et les degrés de liberté. Nous décrivons explicitement une stratégie d’implémentation qui a été appliquée dans le langage open source FreeFem++. Ensuite, nous nous concentrons sur les techniques de préconditionnement du système linéaire résultant de la discrétisation par éléments finis. Nous commençons par la validation numérique d’un préconditionneur à un niveau, de type Schwarz avec recouvrement, avec des conditions de transmission d’impédance entre les sous-domaines. Enfin, nous étudions comment des préconditionneurs à deux niveaux, analysés récemment pour l’équation de Helmholtz, se comportent pour les équations de Maxwell, des points de vue théorique et numérique. Nous appliquons ces méthodes à un problème à grande échelle qui découle de la modélisation d’un système d’imagerie micro-onde, pour la détection et le suivi des accidents vasculaires cérébraux. La précision et la vitesse de calcul sont essentielles dans cette application. / The time-harmonic Maxwell’s equations present several difficulties when the frequency is large, such as the sign-indefiniteness of the variational formulation, the pollution effect and the problematic construction of iterative solvers. We propose a precise and efficient solution strategy that couples high order finite element (FE) discretizations with domain decomposition (DD) preconditioners. High order FE methods make it possible for a given precision to reduce significantly the number of unknowns of the linear system to be solved. DD methods are then used as preconditioners for the iterative solver: the problem defined on the global domain is decomposed into smaller problems on subdomains, which can be solved concurrently and using robust direct solvers. The design, implementation and analysis of both these methods are particularly challenging for Maxwell’s equations. FEs suited for the approximation of the electric field are the curl-conforming or edge finite elements. Here, we revisit the classical degrees of freedom (dofs) defined by Nédélec to obtain a new more friendly expression in terms of the chosen high order basis functions. Moreover, we propose a general technique to restore duality between dofs and basis functions. We explicitly describe an implementation strategy, which we embedded in the open source language FreeFem++. Then we focus on the preconditioning of the linear system, starting with a numerical validation of a one-level overlapping Schwarz preconditioner, with impedance transmission conditions between subdomains. Finally, we investigate how two-level preconditioners recently analyzed for the Helmholtz equation work in the Maxwell case, both from the theoretical and numerical points of view. We apply these methods to the large scale problem arising from the modeling of a microwave imaging system, for the detection and monitoring of brain strokes. In this application accuracy and computing speed are indeed of paramount importance.

Humanitarian clearance of minefields is most often carried out by hand, conventionally using a a metal detector and a probe. Detection is a very slow process, as every piece of detected metal must treated as if it were a landmine and carefully probed and excavated, while many of them are not. The process can be safely sped up by use of Ground-Penetrating Radar (GPR) to image the subsurface, to verify metal detection results and safely ignore any objects which could not possibly be a landmine. In this thesis, we explore the possibility of using Full Wave Inversion (FWI) to improve GPR imaging for landmine detection. Posing the imaging task as FWI means solving the large-scale, non-linear and ill-posed optimisation problem of determining the physical parameters of the subsurface (such as electrical permittivity) which would best reproduce the data. This thesis begins by giving an overview of all the mathematical and implementational aspects of FWI, so as to provide an informative text for both mathematicians (perhaps already familiar with other inverse problems) wanting to contribute to the mine detection problem, as well as a wider engineering audience (perhaps already working on GPR or mine detection) interested in the mathematical study of inverse problems and FWI.We present the first numerical 3D FWI results for GPR, and consider only surface measurements from small-scale arrays as these are suitable for our application. The FWI problem requires an accurate forward model to simulate GPR data, for which we use a hybrid finite-element boundary-integral solver utilising first order curl-conforming N\'d\'{e}lec (edge) elements. We present a novel `line search' type algorithm which prioritises inversion of some target parameters in a region of interest (ROI), with the update outside of the area defined implicitly as a function of the target parameters. This is particularly applicable to the mine detection problem, in which we wish to know more about some detected metallic objects, but are not interested in the surrounding medium. We may need to resolve the surrounding area though, in order to account for the target being obscured and multiple scattering in a highly cluttered subsurface. We focus particularly on spatial sensitivity of the inverse problem, using both a singular value decomposition to analyse the Jacobian matrix, as well as an asymptotic expansion involving polarization tensors describing the perturbation of electric field due to small objects. The latter allows us to extend the current theory of sensitivity in for acoustic FWI, based on the Born approximation, to better understand how polarization plays a role in the 3D electromagnetic inverse problem. Based on this asymptotic approximation, we derive a novel approximation to the diagonals of the Hessian matrix which can be used to pre-condition the GPR FWI problem.