In the field of antenna design and analysis, often the need arises to numerically extrapolate the far-zone performance of a radiating structure from its known (or assumed known) near-zone electromagnetic field. Mathematical processes developed to accomplish such a task are known in the literature as near-zone to far-zone transformations (NZ-FZTs) as well as near-field far-field (NF-FF) transformations. These processes make use of sampled near-zone field quantities along some virtual surface, viz., the transformation surface, that surrounds the radiating structure of interest. Depending upon the application, samples of the required near-zone field quantities are supplied via analytical, empirical, or computational means. Over the years, a number of NZ-FZT processes have been developed to meet the demands of many applications. In short, their differences include, but are not limited to, the following: (1) the size and shape of the transformation surface, (2) the required near-zone field quantities and how they are sampled, (3) the computational methodology used, and (4) the imbedding of various application-driven features. Each process has its pros and cons depending upon its specific application as well as the type of radiation structure under consideration. In this dissertation we put forth a new and original NZ-FZT process that allows the transformation surface along which the near-zone is sampled to be spheroidal in shape: namely a prolate or oblate spheroid. Naturally, there are benefits gained in doing so. Our approach uses a formulated eigenfunction expansion of spheroidal wave-harmonics to develop two distinct, yet closely related, NZ-FZT algorithms for each type of spheroidal transformation surface. The process only requires knowledge of the E-field along the transformation surface and does not need the corresponding H-field. Given is a systematic exposition of the formulation, implementation, and verification of the newly developed NZ-FZT process. Accordingly, computer software is developed to implement both NZ-FZT algorithms. In the validation process, analytical and empirical radiation structures serve as computational benchmarks. Numerical models of both benchmark structures are created by integrating the software with a field solver, viz., a finite-difference time-domain (FDTD) code. Results of these computer models are compared with theoretical and empirical data to provide additional validation. / Ph. D.

The discontinuous Galerkin method (DGM) is applied to a number of problems in computational electromagnetics. This is achieved by obtaining numerical solutions to Maxwell's equations using the DGM. The aim of these simulations is to highlight the strengths of the method while showing its resilience in handling problems other schemes may not be able to accurately model. Although no method will ever be the best choice for every problem in electromagnetics, the discontinuous Galerkin method is able to accurately approximate any problem, although the computational costs can make the scheme impractical for some. Like other time domain schemes, the DGM becomes inefficient on large domains where the solution contains small wavelengths. We demonstrate that all of the different types of boundary conditions in electromagnetic wave propagation can be implemented into the DGM. Reflection and transmission boundaries fit easily into the framework, whereas perfect absorption requires a more advanced technique known as the perfectly matched layer. We begin by simulating mirrors with several different geometries, and analyze how the DGM method performs, and how it offers a more complete evaluation of the behavior in this problem than some other methods. Since Maxwell's equations describe the macroscopic features of electromagnetics, our simulations are able to capture the wave features of electromagnetics, such as interference and diffraction. We demonstrate this by accurately modelling Young's double slit experiment, a classic experiment which features well understood interference and diffraction phenomena. We also extend the basic electromagnetic wave propagation simulations to include situations where the waves travel into new media. The formulation of the DGM for Maxwell's equations allows the numerical solutions to accurately resolve the features at the interface of two media as predicted by the Fresnel coefficients. This allows the DGM to model lenses and other sources of refraction. We predict that the DGM will become an increasingly valuable method for computational electromagnetics because of its wide range of applicability as well as the lack of undesirable features in the numerical solutions. Furthermore, the only limiting factor for applying DGM, its computational cost, will become less influential as computing power continues to increase, allowing us to apply the DGM to an increasing set of applications.

This thesis presents a full-vectorial mode matching method for whispering gallery microcavity analysis. With this technique, optical properties such as resonance wavelength, quality factor and electromagnetic field distribution of an arbitrarily shaped microcavity can be computed with high accuracy. To illustrate this, a mode matching analysis that involves a single propagating whispering gallery mode is performed on a microtoroid in the presence of individual nonplasmonic nanoparticle on its surface. This method is also extended to the analysis of cavity adsorbed by a plasmonic nanoparticle at a wavelength close to plasmon resonance where the resulting field distortion invalidates other approaches. The simulation demonstrates high efficiency and is in close agreement with experimental measurements reported in previous work. Furthermore, we extend our mode matching analysis to the case where multiple whispering gallery modes are involved in the course of light propagation. The new formalism is performed on a cavity-waveguide coupling system to investigate the light delivery from a tapered optical waveguide to a microcavity at high precision. A novel hybrid integration scheme to implement an ultra-high quality factor microcavity on a silicon-on-insulator platform is proposed based on the related modelling results. / Graduate / 0752 / 0544 / duxuanmax@gmail.com

The discontinuous Galerkin method (DGM) is applied to a number of problems in computational electromagnetics. This is achieved by obtaining numerical solutions to Maxwell's equations using the DGM. The aim of these simulations is to highlight the strengths of the method while showing its resilience in handling problems other schemes may not be able to accurately model. Although no method will ever be the best choice for every problem in electromagnetics, the discontinuous Galerkin method is able to accurately approximate any problem, although the computational costs can make the scheme impractical for some. Like other time domain schemes, the DGM becomes inefficient on large domains where the solution contains small wavelengths. We demonstrate that all of the different types of boundary conditions in electromagnetic wave propagation can be implemented into the DGM. Reflection and transmission boundaries fit easily into the framework, whereas perfect absorption requires a more advanced technique known as the perfectly matched layer. We begin by simulating mirrors with several different geometries, and analyze how the DGM method performs, and how it offers a more complete evaluation of the behavior in this problem than some other methods. Since Maxwell's equations describe the macroscopic features of electromagnetics, our simulations are able to capture the wave features of electromagnetics, such as interference and diffraction. We demonstrate this by accurately modelling Young's double slit experiment, a classic experiment which features well understood interference and diffraction phenomena. We also extend the basic electromagnetic wave propagation simulations to include situations where the waves travel into new media. The formulation of the DGM for Maxwell's equations allows the numerical solutions to accurately resolve the features at the interface of two media as predicted by the Fresnel coefficients. This allows the DGM to model lenses and other sources of refraction. We predict that the DGM will become an increasingly valuable method for computational electromagnetics because of its wide range of applicability as well as the lack of undesirable features in the numerical solutions. Furthermore, the only limiting factor for applying DGM, its computational cost, will become less influential as computing power continues to increase, allowing us to apply the DGM to an increasing set of applications.

No description available.

In order to accurately describe the behavior of an antenna, one needs to understand the radiation mechanisms that govern its operation. One way to gain such an insight is to view the fields and currents present on a radiating structure. Unfortunately, in close proximity to an antenna empirical techniques fail because the measurement probe alters the operation of the radiating structure. Computational methods offer a solution to this problem. By simulating the operation of an antenna, one can obtain electromagnetic field data near (or even internal to) a radiating structure. However, these computationally intense techniques often generate extremely large data sets that cannot be adequately interpreted using traditional graphical approaches. A visualization capability is developed that allows an analysis of the above-mentioned data sets. With this technique, the data is viewed from a unique, global perspective. This format is well suited for analytical investigations as well as debugging during modeling and simulation. An illustrative example is provided in the context of a rectangular microstrip patch antenna. A comparison is performed between the visualized data and the theory of operation for the microstrip patch in order to demonstrate that radiation mechanisms can be obtained visually. An additional analysis tool is developed using Gabor filters and image-processing techniques. This tool allows one to detect and filter electromagnetic waves propagating with different velocities (both speed and direction). By doing so, each mode of an antenna can be analyzed independently. The fields of a multi-moded, open-ended rectangular waveguide are analyzed in order to demonstrate the effectiveness of these techniques. / Master of Science

The Transmission Line Matrix (TLM) method is an existing and established mathematical method for conducting computational electromagnetic (CEM) simulations. TLM models Maxwell s equations by discretising the contiguous nature of an environment and its contents into individual small-scale elements and it is a computationally intensive process. This thesis focusses on parallel processing optimisations to the TLM method when considering the opposing ends of the contemporary computing hardware spectrum, namely large-scale computing systems versus small-scale mobile computing devices. Theoretical aspects covered in this thesis are: The historical development and derivation of the TLM method. A discrete random variable (DRV) for rain-drop diameter,allowing generation of a rain-field with raindrops adhering to a Gaussian size distribution, as a case study for a 3-D TLM implementation. Investigations into parallel computing strategies for accelerating TLM on large and small-scale computing platforms. Implementation aspects covered in this thesis are: A script for modelling rain-fields using free-to-use modelling software. The first known implementation of 2-D TLM on mobile computing devices. A 3-D TLM implementation designed for simulating the effects of rain-fields on extremely high frequency (EHF) band signals. By optimising both TLM solver implementations for their respective platforms, new opportunities present themselves. Rain-field simulations containing individual rain-drop geometry can be simulated, which was previously impractical due to the lengthy computation times required. Also, computationally time-intensive methods such as TLM were previously impractical on mobile computing devices. Contemporary hardware features on these devices now provide the opportunity for CEM simulations at speeds that are acceptable to end users, as well as providing a new avenue for educating relevant user cohorts via dynamic presentations of EM phenomena.

This work presents an investigation into the optical scattering of heterogeneous cells with an application to two-photon imaging, optical scattering measurements and STED imaging. Using the finite difference time-domain (FDTD) method, the full-wave scattering by many cells containing multiple organelles with varying indices of refraction is computed. These simulations were previously limited to single cells for reasons of computational cost. A superposition approximation that uses the coherent linear superposition of FDTD-determined farfield scattering patterns of small numbers of cells to estimate the scattering from a larger tissue was developed and investigated. It was found that for the approximation to be accurate, the scattering sub-problems must at minimum extend along the incident field propagation axis for the full depth of the tissue, preserving the scattering that takes place in the direction of propagation. The FDTD method was used to study the scattering effects of multiple inhomogeneous cells on the propagation of a focused Gaussian beam with an application to two-photon imaging. It was found that scattering is mostly responsible for the reduction in two-photon fluorescence signal as depth is increased. It was also determined that for the chosen beam parameters and the cell and organelle configurations used, the nuclei are the dominant scatterers. FDTD was also utilized in an investigation of cellular scattering effects on the propagation of a common depletion beam used in STED microscopy and how scattering impacts the image obtained with a STED microscope. An axial doughnut beam was formulated and implemented in FDTD simulations, along with a corresponding focused Gaussian beam to simulate a fluorescence excitation beam. It was determined that the depletion beam will maintain a well-defined axial null in spite of scattering, although scattering will reduce the resulting fluorescence signal with focal depth. / text

A set of tools are proposed for the efficient modeling of several classes of problems related to periodic structures in microwave and optical regimes with Finite-Difference Time-Domain method. The first category of problems under study is the interaction of non-periodic sources and printed elements with infinitely periodic structures. Such problems would typically require a time-consuming simulation of a finite number of unit cells of the periodic structures, chosen to be large enough to achieve convergence. To alleviate computational cost, the sine-cosine method for the Finite-Difference Time-Domain based dispersion analysis of periodic structures is extended to incorporate the presence of non-periodic, wideband sources, enabling the fast modeling of driven periodic structures via a small number of low cost simulations. The proposed method is then modified for the accelerated simulation of microwave circuit geometries printed on periodic substrates. The scheme employs periodic boundary conditions applied at the substrate, to dramatically reduce the computational domain and hence, the cost of such simulations. Emphasis is also given on radiation pattern calculation, and the consequences of the truncated computational domain of the proposed method on the computation of the electric and magnetic surface currents invoked in the near-to-far-field transformation. It has been further demonstrated that from the mesh truncation point of view, the scheme, which has a unified form regardless dispersion and conductivity, serves as a much simpler but equally effective alternative to the Perfectly Matched Layer provided that the simulated domain is periodic in the direction of termination. The second category of problems focuses on the efficient characterization of nonlinear periodic structures. In Finite-Difference Time-Domain, the simulation of these problems is typically hindered by the fine spatial and time gridding. Originally proposed for linear structures, the Alternating-Direction Implicit Finite-Difference Time-Domain method, as well as a novel spatial filtering method, are extended to incorporate nonlinear media. Both methods are able to use time-step sizes beyond the conventional stability limit, offering significant savings in simulation time.

A set of tools are proposed for the efficient modeling of several classes of problems related to periodic structures in microwave and optical regimes with Finite-Difference Time-Domain method. The first category of problems under study is the interaction of non-periodic sources and printed elements with infinitely periodic structures. Such problems would typically require a time-consuming simulation of a finite number of unit cells of the periodic structures, chosen to be large enough to achieve convergence. To alleviate computational cost, the sine-cosine method for the Finite-Difference Time-Domain based dispersion analysis of periodic structures is extended to incorporate the presence of non-periodic, wideband sources, enabling the fast modeling of driven periodic structures via a small number of low cost simulations. The proposed method is then modified for the accelerated simulation of microwave circuit geometries printed on periodic substrates. The scheme employs periodic boundary conditions applied at the substrate, to dramatically reduce the computational domain and hence, the cost of such simulations. Emphasis is also given on radiation pattern calculation, and the consequences of the truncated computational domain of the proposed method on the computation of the electric and magnetic surface currents invoked in the near-to-far-field transformation. It has been further demonstrated that from the mesh truncation point of view, the scheme, which has a unified form regardless dispersion and conductivity, serves as a much simpler but equally effective alternative to the Perfectly Matched Layer provided that the simulated domain is periodic in the direction of termination. The second category of problems focuses on the efficient characterization of nonlinear periodic structures. In Finite-Difference Time-Domain, the simulation of these problems is typically hindered by the fine spatial and time gridding. Originally proposed for linear structures, the Alternating-Direction Implicit Finite-Difference Time-Domain method, as well as a novel spatial filtering method, are extended to incorporate nonlinear media. Both methods are able to use time-step sizes beyond the conventional stability limit, offering significant savings in simulation time.