Application of Helmholtz/Hodge Decomposition to Finite Element Methods for Two-Dimensional Maxwell\u27s Equations

In this work we apply the two-dimensional Helmholtz/Hodge decomposition to develop new finite element schemes for two-dimensional Maxwell\u27s equations. We begin with the introduction of Maxwell\u27s equations and a brief survey of finite element methods for Maxwell\u27s equations. Then we review the related fundamentals in Chapter 2. In Chapter 3, we discuss the related vector function spaces and the Helmholtz/Hodge decomposition which are used in Chapter 4 and 5. The new results in this dissertation are presented in Chapter 4 and Chapter 5. In Chapter 4, we propose a new numerical approach for two-dimensional Maxwell\u27s equations that is based on the Helmholtz/Hodge decomposition for divergence-free vector fields. In this approach an approximate solution for Maxwell\u27s equations can be obtained by solving standard second order scalar elliptic boundary value problems. This new approach is illustrated by a P1 finite element method. In Chapter 5, we further extend the new approach described in Chapter 4 to the interface problem for Maxwell\u27s equations. We use the extraction formulas and multigrid method to overcome the low regularity of the solution for the Maxwell interface problem. The theoretical results obtained in this dissertation are confirmed by numerical experiments

Experimental and theoretical study of an integrated silicon Mach-Zehnder interferometer for chemical sensing applications

This thesis involves the design, fabrication and characterization of an integrated optical waveguide sensor. Prior to fabrication, design parameters of the waveguide need to be determined and optimized. The waveguide parameters such as waveguide dimension and the refractive index of the core and cladding are obtained from the single-mode cutoff frequency calculated using either analytical or numerical methods. In this thesis, details of analytical calculations to determine the cutoff frequency in terms of the waveguide parameters will be presented. The method discussed here is Marcatili\u27s approximation. The purpose is to solve the scalar wave equation derived from Maxwell\u27s equations because it describes the mode properties inside the waveguides. The Finite Element Method is used to simulate the electric and magnetic fields inside the waveguides and to determine the propagation characteristics in optical waveguides. This method is suited for problems involving complicated geometries and variable index of refraction. Fabrication of the Integrated Mach-Zehnder Interferometer sensor involves several important standard processes such as Chemical Vapor Deposition (CVD) for thin film fabrication, photolithography for mask transfer, and etching for ridge waveguide formation. The detailed fabrication procedures of the tested Mach-Zehnder Interferometer sensors are discussed. After completion of the sensor fabrication processes, the characterizations were carried out for the thin film of Si02 and PSG, the waveguides and the Y-junction separately. The waveguides were analyzed to make sure that the sensors are working as expected. The experimental testing on the separated waveguide portions of the first batch Integrated Mach-Zehnder Interferometer (MZI) sensors are described. These testing procedures were also performed for the subsequent fabricated batches of the integrated MZI sensors until optimum performance is achieved. A new concept has been proposed for chemical sensing applications. The novelty of the approach is mainly based on utilizing the multi -wavelength or broadband source instead of single wavelength input to the integrated MZI. The shifting of output spectra resulting from the interference has shown the ability of the MZI to analyze the different concentrations of a chemical analyte. The sensitivity of the sensor is also determined from the plot of intensity versus concentration, which is around 0.013 (%ml)-1 and 0.007 (%ml)-1 for the white light source and the 1.5 ~tm broadband source, respectively, while the lowest detectable concentration of ethanol for the sensor detection is around 8% using a intensity variation method and 0.6% using a peak wavelength variation method

Advanced techniques in scientific computing: application to electromagnetics

Mención Internacional en el título de doctorDurante los últimos años, los componentes de radiofrecuencia que forman parte de un sistema de comunicaciones necesitan simulaciones cada vez más exigentes desde el punto de vista de recursos computacionales. Para ello, se han desarrollado diferentes técnicas con el método de los elementos finitos (FEM) como la conocida como adaptatividad hp, que consiste en estimar el error en el problema electromagnético para generar mallas de elementos adecuadas al problema que obtienen una aproximación de forma más efectiva que las mallas estándar; o métodos de descomposición de dominios (DDM), basado en la división del problema original en problemas más pequeños que se pueden resolver en paralelo. El principal problema de las técnicas de adaptatividad es que ofrecen buenas prestaciones para problemas bidimensionales, mientras que en tres dimensiones el tiempo de generación de las mallas adaptadas es prohibitivo. Por otra parte, DDM se ha utilizado satisfactoriamente para la simulación de problemas eléctricamente muy grandes y de gran complejidad, convirtiéndose en uno de los temas más actuales en la comunidad de electromagnetismo computacional. El principal objetivo de este trabajo es estudiar la viabilidad de algoritmos escalables (en términos de paralelización) combinando DDM no conformes y adaptatividad automática en tres dimensiones. Esto permitir ía la ejecución de algoritmos de adaptatividad independiente en cada subdominio de DDM. En este trabajo se presenta y discute un prototipo que combina técnicas de adaptatividad y DDM, que aún no se han tratado en detalle en la comunidad científica. Para ello, se implementan tres bloques fundamentales: i) funciones de base para los elementos finitos que permitan órdenes variables dentro de la misma malla; ii) DDM no conforme y sin solapamiento; y iii) algoritmos de adaptatividad en tres dimensiones. Estos tres bloques se han implementado satisfactoriamente en un código FEM mediante un método sistemático basado en el método de las soluciones manufacturadas (MMS). Además, se ha llevado a cabo una paralelización a tres niveles: a nivel de algoritmo, con DDM; a nivel de proceso, con MPI (Message Passing Interface); y a nivel de hebra, con OpenMP; todo en un código modular que facilita el mantenimiento y la introducción de nuevas características. Con respecto al primer bloque fundamental, se ha desarrollado una familia de funciones base con un enfoque sistemático que permite la expansión correcta del espacio de funciones. Por otra parte, se han introducido funciones de base jerárquicas de otros autores (con los que el grupo al que pertenece el autor de la tesis ha colaborado estrechamente en los últimos años) para facilitar la introducción de diferentes órdenes de aproximación en el mismo mallado. En lo relativo a DDM, se ha realizado un estudio cuantitativo del error generado por las disconformidades en la interfaz entre subdominios, incluidas las discontinuidades generadas por un algoritmo de adaptatividad. Este estudio es fundamental para el correcto funcionamiento de la adaptatividad, y no ha sido evaluado con detalle en la comunidad científica. Además, se ha desarrollado un algoritmo de adaptatividad con prismas triangulares, haciendo especial énfasis en las peculiaridades debidas a la elección de este elemento. Finalmente, estos tres bloques básicos se han utilizado para desarrollar, y discutir, un prototipo que une las técnicas de adaptatividad y DDM.In the last years, more and more accurate and demanding simulations of radiofrequency components in a system of communications are requested by the community. To address this need, some techniques have been introduced in finite element methods (FEM), such as hp adaptivity (which estimates the error in the problem and generates tailored meshes to achieve more accuracy with less unknowns than in the case of uniformly refined meshes) or domain decomposition methods (DDM, consisting of dividing the whole problem into more manageable subdomains which can be solved in parallel). The performance of the adaptivity techniques is good up to two dimensions, whereas for three dimensions the generation time of the adapted meshes may be prohibitive. On the other hand, large scale simulations have been reported with DDM becoming a hot topic in the computational electromagnetics community. The main objective of this dissertation is to study the viability of scalable (in terms of parallel performance) algorithms combining nonconformal DDM and automatic adaptivity in three dimensions. Specifically, the adaptivity algorithms might be run in each subdomain independently. This combination has not been detailed in the literature and a proof of concept is discussed in this work. Thus, three building blocks must be introduced: i) basis functions for the finite elements which support non-uniform approximation orders p; ii) non-conformal and non-overlapping DDM; and iii) adaptivity algorithms in 3D. In this work, these three building blocks have been successfully introduced in a FEM code with a systematic procedure based on the method of manufactured solutions (MMS). Moreover, a three-level parallelization (at the algorithm level, with DDM; at the process level, with message passing interface (MPI), and at the thread level, with OpenMP) has been developed using the paradigm of modular programming which eases the software maintenance and the introduction of new features. Regarding first building block, a family of basis functions which follows a sound mathematical approach to expand the correct space of functions is developed and particularized for triangular prisms. Also, to ease the introduction of different approximation orders in the same mesh, hierarchical basis functions from other authors are used as a black box. With respect to DDM, a thorough study of the error introduced by the non-conformal interfaces between subdomains is required for the adaptivity algorithm. Thus, a quantitative analysis is detailed including non-conformalities generated by independent refinements in neighbor subdomains. This error has not been assessed with detail in the literature and it is a key factor for the adaptivity algorithm to perform properly. An adaptivity algorithm with triangular prisms is also developed and special considerations for the implementation are explained. Finally, on top of these three building blocks, the proof of concept of adaptivity with DDM is discussed.Programa Oficial de Doctorado en Multimedia y ComunicacionesPresidente: Daniel Segovia Vargas.- Secretario: David Pardo Zubiaur.- Vocal: Romanus Dyczij-Edlinge

Modeling EMI Resulting from a Signal Via Transition Through Power/Ground Layers

Signal transitioning through layers on vias are very common in multi-layer printed circuit board (PCB) design. For a signal via transitioning through the internal power and ground planes, the return current must switch from one reference plane to another reference plane. The discontinuity of the return current at the via excites the power and ground planes, and results in noise on the power bus that can lead to signal integrity, as well as EMI problems. Numerical methods, such as the finite-difference time-domain (FDTD), Moment of Methods (MoM), and partial element equivalent circuit (PEEC) method, were employed herein to study this problem. The modeled results are supported by measurements. In addition, a common EMI mitigation approach of adding a decoupling capacitor was investigated with the FDTD method

Finite element analysis and design optimisation of shaded pole induction motors

SIGLEAvailable from British Library Document Supply Centre-DSC:DX212729 / BLDSC - British Library Document Supply CentreGBUnited Kingdo

Investigations of the analysis and modelling of magnetotelluric data

SIGLEAvailable from British Library Document Supply Centre- DSC:D75707/87 / BLDSC - British Library Document Supply CentreGBUnited Kingdo