The principle of coordinate invariance states that all physical laws must be formulated in a mathematical form that is independent of the geometrical properties of any particular coordinate system. Embracing this principle is the key to understand how to systematically incorporate curved material interfaces into a numerical solution of Maxwell’s equations. This dissertation describes how to generate a coordinate invariant representation of Maxwell’s equations in differential form, and it demonstrates why employing such representation is crucial to the development of robust finite-difference discretisations with consistent global error properties. As part of this process, two original contributions are presented that address the issue of constructing finite-difference approximations at the locations of material interfaces. The first contribution is a domain-decomposition procedure to enforce the tangential field continuity conditions with a second-order local truncation error that can be applied in 2-D or 3-D. The second contribution is a similar domain-decomposition procedure that enforces the tangential field continuity conditions with a local truncation of order 2L—where L is an integer greater or equal to one—but that can only be applied in 1-D. To conclude, the dissertation also describes the interesting connection that exists between the use of a coordinate invariant representation of Maxwell’s equations to design artificial materials and the use of the same representation to model curved material interfaces in a finite-difference discretisation.

First, one of the purposes of this thesis is to estimate the complete crosstalk effects including the package and the pads on the surface acoustic wave (SAW) substrate. A new approach based on finite-difference time-domain (FDTD) with equivalent current source method is applied. Two kinds of patterns of one-port SAW resonators with the same package structure and inter-digital transducer (IDT) design are studied. Verification with the measurement results shows that our method is able to obtain good agreement and be used to observe the influence from the SAW pattern. Second, the equivalent current source method is extended to model the excitation of human-body¡¦s electrostatic discharge (ESD) situations. The efficiencies of sacrificial electrodes are also discussed. Finally, a novel sacrificial electrode with fractal to protect SAW devices from ESD break is proposed. Comparing with traditional electrode, the simulation results show that fractal can improve the protective efficiency greatly. Finally, a novel analysis model that can be used to analyze and optimize the impedance of an RFID transponder integrated circuit (IC) which uses backscatter encoding based on simultaneously maintaining the BER of the reader and maximizing the received power of the transponder IC is proposed. The analysis method utilizes mapping from signal constellation of the backscattered signal to the Smith chart to relate the two parameters. Given the system specification and characteristics of the reader and transponder antennas, the optimum impedances of transponder IC for binary communication system can be easily designed by using this model.

Ziel dieser Arbeit war die Berechnung elektromagnetischer Nahfelder, die wesentlich sind für die Charakterisierung von Strukturen im Submikrometerbereich. Diese Aufgabenstellung wurde im Rahmen der klassischen Elektrodynamik unter Vernachlässigung von quantenmechanischen und relativistischen Effekten durchgeführt. Die untersuchten Modellsysteme bestanden aus stückweise homogenen Medien. Eine Beschränkung auf eine harmonische Zeitabhängigkeit der Felder fand nicht statt. Zum Einsatz kamen analytische (Mie-Theorie und deren Erweiterungen), semi-analytische (MMP) und rein numerische Methoden (FDTD). Besonders umfassend wurden die Eigenschaften evaneszenter Wellen untersucht. Entgegen der oft üblichen Vorgehensweise wurden die Beschränkung auf 2-dimensionale Modelle vermieden und Metalle nicht idealisiert als perfekt leitend, sondern realistisch über einen komplexen Brechungsindex bzw. über ein äquivalentes Drude-Modell beschrieben. Nur so ist es möglich die Ausbreitung von surface plasmon polaritons zu modellieren und den Einfluß von Volmenplasmonen zu berücksichtigen. Untersucht wurden periodische und nicht periodische Strukturen aus dielektrischen und metallischen Materialien. Solche Systeme sind nützlich aufgrund der Bildung von photonischen Bandlücken (Dielektrika) und der Realisierung hoher Feldverstärkungen (Metalle). Die erste Eigenschaft kann für die Konstruktion von besonders effektiven Laser und die zweite im Rahmen der oberflächenverstärkten Raman-Streuung angewendet werden. Eine weiterer Schwerpunkt dieser Dissertation war die Analyse nahfeld-optischer Mikroskope (SNOM). Mit solchen Apparaturen kann eine Auflösung jenseits des Abbe-Limit erreicht werden. Untersucht wurden die Abbildungseigenschaften aperturloser Nahfeld-Mikroskope und die Ausbreitung von Femto-Sekunden Pulsen in einem konventionellen SNOM, welches mit einer metallbeschichteten Glasfaser ausgestattet ist. Die zweite Fragestellung ist relevant für die Kombination von hoher räumlicher mit hoher zeitlicher Auflösung.

Es werden nahfeldoptische Phänomene an Raster-Sonden-Geometrien untersucht. Hierzu zählen die Feldverstärkung an laserbestrahlten Rastersondenspitzen und die Feldverteilung in nahfeldoptischen Apertur- und Koaxialspitzen. Zur Berechnung der Feldverteilung werden verschiedene numerische Verfahren wie die Methode der Randelemente (BEM) und die Methode der Finiten Integration (FIT) angewendet. Die durchgeführten Berechnungen sagen eine erhebliche Feldüberhöhung im Nahfeld von laserbestrahlten Rastersondenspitzen voraus. Es wird diskutiert, ob diese Feldüberhöhung ursächlich für die Modifikation der Oberfläche unter der Sondenspitze verantwortlich ist. Zur weiteren Klärung der Strukturierungsmechanismen wird die Stromantwort des Tunnelübergangs auf Laserbestrahlung experimentell untersucht. Die Detektion des Stromsignals erfolgt mit einem neuen Vorverstärker, der die Eigenschaften einer hohen Eingangsimpedanz, einer hohen Transimpedanz und einer hohen Bandbreite miteinander vereinigt. Weiterhin wird die Feldverteilung in SNOM-Spitzen, insbesondere in den sogenannten Koaxialspitzen, bestimmt. Es wird untersucht wie sich koaxial Moden in koaxialen SNOM-Spitzen anregen lassen und wie sie sich ausbreiten. Ausgehend von bereits existierenden Entwürfen solcher Spitzen wird diskutiert, welche Strukturen sich zur Anregung von koaxialen Moden eignen.

This dissertation describes two possible methods for the calibration of an ultrasound tomography system developed at University of Manitoba's Electromagnetic Imaging Laboratory for imaging with the contrast-source inversion algorithm. The calibration techniques are adapted from existing procedures employed for microwave tomography. A theoretical model of these calibration principles is developed in order to provide a rationale for the effectiveness of the proposed procedures. The applicability of such an imaging algorithm and calibration methods in the context of ultrasound are discussed. Also presented are 2D and 3D finite-difference time-domain update equations for the simulation of acoustic wave propagation in inhomogeneous media. Details regarding the application of an absorbing boundary-condition, point-source modelling and the treatment of penetrable objects are included in this document.

This dissertation describes two possible methods for the calibration of an ultrasound tomography system developed at University of Manitoba's Electromagnetic Imaging Laboratory for imaging with the contrast-source inversion algorithm. The calibration techniques are adapted from existing procedures employed for microwave tomography. A theoretical model of these calibration principles is developed in order to provide a rationale for the effectiveness of the proposed procedures. The applicability of such an imaging algorithm and calibration methods in the context of ultrasound are discussed. Also presented are 2D and 3D finite-difference time-domain update equations for the simulation of acoustic wave propagation in inhomogeneous media. Details regarding the application of an absorbing boundary-condition, point-source modelling and the treatment of penetrable objects are included in this document.

La propagation acoustique en milieu externe est fortement influencée par l'environnement. Les effets liés à la géométrie, comme la topographie ou la présence d'obstacles, sont principalement les réflexions et les diffractions. Concernant l'effet de l'atmosphère, les gradients moyens génèrent des réfractions tandis que la turbulence provoque des fluctuations aléatoires et une perte de cohérence du signal. La plupart de ces effets sont généralement bien décrits de manière théorique, cependant dans les configurations réelles le cumul de tous ces effets rend l'utilisation des expressions analytiques très difficile. Les études expérimentales présentent également des limites liées à la difficulté de connaître l'environnement parfaitement et d'isoler un effet physique particulier. Dans cette perspective, la simulation numérique est une alternative pratique et complémentaire à la théorie et l'expérimentation. Parmi les modèles numériques de propagation existants, ceux basés sur une résolution par différences finies dans le domaine temporel (FDTD pour Finite-Difference Time-Domain) des équations d'Euler linéarisées sont récents et particulièrement prometteurs. Cependant comme pour tout modèle nouveau, il reste à montrer qu'effectivement l'ensemble des phénomènes physiques d'intérêt sont retranscrits.Dans le cadre de ses études sur la propagation acoustique extérieure, l'Institut franco-allemand de recherches de Saint-Louis (ISL) a implémenté un tel modèle de propagation. Cette implémentation est ci-après appelée ITM, pour ISL FDTD Model. L'objectif de cette thèse, proposée par l'ISL en collaboration avec le Laboratoire de Mécanique des Fluides et d'Acoustique (LMFA), est de poursuivre le développement et les validations de cette implémentation. Une part importante du travail consiste également à illustrer les potentialités du code ITM pour des applications de propagation de signaux acoustiques complexes dans un environnement complexe. [...]

This thesis reports on the application of ground penetrating radar (GPR) as a non-destructive technique for the monitoring of ring separation in brick masonry arch bridges. In addition, research is reported on the assessment of the clay capping layer often used in construction as a waterproof backing to arches. The thrust of the research is numerical modelling, verified by large laboratory experiments. Due to the heterogeneity of these structures, the resultant signals from the interaction between the GPR system and the bridge are often complex and hence, hard to interpret. This highlighted the need to create a GPR numerical model that would allow the study of the attributes of reflected signals from various targets within the structure of the bridge. The GPR numerical analysis was undertaken using the finite-difference time-domain (FDTD) method. Since micro regions in the bridge structure need to be modelled, the introduction of subgrids of supporting finer spatial resolution into the standard FDTD method was considered essential in order to economise on the required computational resources. In the main part of this thesis, it is demonstrated how realistic numerical modelling of GPR using the FDTD method could greatly benefit from the implementation of subgrids into the conventional FDTD mesh. This is particularly important when (a) parts of the computational domain need to be modelled in detail (i.e., ring separation between the mortar layers and the brick units, which is the case studied in this thesis); and also (b) when there are features or regions in the overall computational mesh with values of high relative permittivity supporting propagation of waves at very short wavelengths. A scheme is presented that simplifies the process of implementing these subgrids into the traditional FDTD method. This scheme is based on the combination of the standard FDTD method and the unconditionally stable alternating-direction implicit (ADI) FDTD technique. Given that ADI-FDTD is unconditionally stable, its time-step can be set to any value that facilitates the accurate calculation of the electromagnetic fields. By doing so, the two grids can efficiently communicate information across their boundary without requiring to use a time-interpolation scheme. The performance of ADI-FDTD subgrids when implemented into the traditional FDTD method is discussed herein. The developed algorithm can handle cases where the subgrid crosses dielectrically inhomogeneous and/or conductive media. In addition, results from the comparison between the proposed scheme and a commonly employed purely FDTD subgridding technique are presented. After determination of the optimum ADI-FDTD scheme, numerical experiments were conducted and calibrated using GPR laboratory experiments. Good correlations were obtained between the numerical experiments and the actual GPR experiments. It was shown both numerically and experimentally that significant mortar loss between the masonry arch rings can be detected. Dry hairline delaminations between the mortar and the brick masonry are difficult to detect using standard GPR procedures. However, hairline faults containing water produce distinct and detectable GPR responses. In addition, the clay layer was successfully identified and its thickness calculated to a satisfactory accuracy.

The objective of this research work is to develop an efficient methodology for chip-package cosimulation. In the traditional design flow, the integrated circuit (IC) is first designed followed by the package design. The disadvantage of the conventional sequential design flow is that if there are problems with signal and power integrity after the integration of the IC and the package, it is expensive and time consuming to go back and change the IC layout for a different input/output (IO) pad assignment. To overcome this limitation, a concurrent design flow, where both the IC and the package are designed together, has been recommended by researchers to obtain a fast design closure. The techniques from this research work will enable multiscale cosimulation of the chip and the package making the concurrent design flow paradigm possible. Traditional time-domain techniques, such as the finite-difference time-domain method, are limited by the Courant condition and are not suitable for chip-package cosimulation. The Courant condition gives an upper bound on the time step that can be used to obtain stable simulation results. The smaller the mesh dimension the smaller is the Courant time step. In the case of chip-package cosimulation the on-chip structures require a fine mesh, which can make the time step prohibitively small. An unconditionally stable scheme using Laguerre polynomials has been recommended for chip-package cosimulation. Prior limitations in this method have been overcome in this research work. The enhanced transient simulation scheme using Laguerre polynomials has been named SLeEC, which stands for simulation using Laguerre equivalent circuit. A full-wave EM simulator has been developed using the SLeEC methodology. A scheme for efficient use of full-wave solver for chip-package cosimulation has been proposed. Simulation of the entire chip-package structure using a full-wave solver could be a memory and time-intensive operation. A more efficient way is to separate the chip-package structure into the chip, the package signal-delivery network, and the package power-delivery network; use a full-wave solver to simulate each of these smaller subblocks and integrate them together in the following step, before a final simulation is done on the integrated network. Examples have been presented that illustrate the technique.

Numerical simulation is an important tool in understanding the electromagnetic field and how it interacts with the environment. Different topics for time-domain finite-difference (FDTD) and finite-element (FETD) methods for Maxwell's equations are treated in this thesis. Subcell models are of vital importance for the efficient modeling of small objects that are not resolved by the grid. A novel model for thin sheets using shell elements is proposed. This approach has the advantage of taking into account discontinuities in the normal component of the electric field, unlike previous models based on impedance boundary conditions (IBCs). Several results are presented to illustrate the capabilities of the shell element approach. Waveguides are of fundamental importance in many microwave applications, for example in antenna feeds. The key issues of excitation and truncation of waveguides are addressed. A complex frequency shifted form of the uniaxial perfectly matched layer (UPML) absorbing boundary condition (ABC) in FETD is developed. Prism elements are used to promote automatic grid generation and enhance the performance. Results are presented where reflection errors below -70dB are obtained for different types of waveguides, including inhomogeneous cases. Excitation and analysis via the scattering parameters are achieved using waveguide modes computed by a general frequency-domain mode solver for the vector Helmholtz equation. Huygens surfaces are used in both FDTD and FETD for excitation in waveguide ports. Inverse problems have received an increased interest due to the availability of powerful computers. An important application is non-destructive evaluation of material. A time-domain, minimization approach is presented where exact gradients are computed using the adjoint problem. The approach is applied to a general form of Maxwell's equations including dispersive media and UPML. Successful reconstruction examples are presented both using synthetic and experimental measurement data. Parameter reduction of complex geometries using simplified models is an interesting topic that leads to an inverse problem. Gradients for subcell parameters are derived and a successful reconstruction example is presented for a combined dielectric sheet and slot geometry.