Modeling human head tissues using fourth-order debye model in convolution-based three-dimensional finite-difference time-domain

A fourth order Debye model is derived using genetic algorithms to represent the dispersive properties of the 17 tissues that form the human head. The derived model gives accurate estimation of the electrical properties of those tissues across the frequency band from 0.1 GHz to 3 GHz that can be used in microwave systems for head imaging. A convolution-based three-dimensional finite-difference time-domain (3D-FDTD) formulation is implemented for modeling the electromagnetic wave propagation in the dispersive head tissues whose frequency dependent properties are represented by the derived fourth-order Debye model. The presented results show that the proposed 3D-FDTD and fourth-order Debye model can accurately show the electromagnetic interaction between a wide band radiation and head tissues with low computational overhead and more accurate results compared with using multi-pole Cole-Cole model

Transmission and reflection properties of layered left-handed materials

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, February 2007.Includes bibliographical references (p. 179-189).This thesis is concerned with the reflection and transmission properties of layered left-handed materials (LHM). In particular, the reflection properties of (LHM) slabs are studied for the Goos-Hanchen (GH) lateral shift phenomenon. We demonstrate a unique GH lateral shift phenomenon, which shows that both positive and negative shifts can be achieved using the same LHM slab configuration. This phenomenon is different from previously established cases where the GH lateral shift can be only negative or only positive when different LHM slab configurations are used. We also show that there exist two distinct cases with this unique phenomenon. One case has two regions of incident angles where the GH lateral shift directions are different, while another case has three regions with alternated GH shift directions. A generalized analytical formulation for analyzing the GH lateral shift direction is provided, which reveals that this unique phenomenon is related to the relative amplitudes of the growing and decaying evanescent waves inside the LHM slabs. The energy flux patterns within LHM slabs are further studied to show the influence of the evanescent waves on the GH shift direction change.(cont.) Furthermore, the transmission property of LHM slabs are studied on the finite slabs' maging capability. First, the development of the numerical simulation tool - the Finite-Difference Time-Domain method (FDTD) - investigates the ability of the method to model a perfect lens made of a slab of homogeneous LHM. It is shown that because of the frequency dispersive nature of the medium and the time discretization, an inherent mismatch in the constitutive parameters exists between the slab and its surrounding medium. This mismatch in the real part of the permittivity and permeability is found to have the same order of magnitude as the losses typically used in numerical simulations. Hence, when the LHM slab is lossless, this mismatch is shown to be the main factor contributing to the image resolution loss of the slab. In addition, finite-size LHM slabs are studied both analytically and numerically since they have practical importance in the actual experiments. The analytical method is based on Huygens' principles using truncated current sheets that cover only the apertures of the slabs. It is shown that the main effects on the images' spectra due to the size of the slabs can be predicted by the proposed analytical method, which can, therefore, be used as a fast alternative to numerical simulations.(cont.) Furthermore, the property of negative energy streams at the image plane is also investigated. This unique property is found to be due to the interactions between propagating and evanescent waves and can only occur with LHM slabs, of both finite-size and infinite size. The last part of the thesis deals with multi-layered media for the application to antenna isolations. The setup is with two horn antennas located beneath the ground plane with 10 A distance apart. In order to reduce the coupling between antennas, multi-layered media placed on top of the ground plane need to be designed to suppress the fields. After the problem is simplified to the dipole antenna coupling in infinite slabs, the method to evaluate the fields inside layered media is presented. This method obtains the spectral domain Green's function first and then transforms the fields to the spatial domain using the Sommerfeld-type integration. After the method is validated using right-handed materials (RHM) from references, it is extended to include media like LHM as well as p. negative material and : negative material . The validation with these materials are done by comparing the results with CST microwave studio simulations. The first configuration for the antenna isolation design if one layer slab backed by the grounded plane. Two different approaches are used to find the optimum slab parameters for the isolation.(cont.) One approach is to use Genetic Algorithm (GA) to optimize the slab's constitutive parameters and the thickness for a minimum coupling level. The other approach is to develop an analytic asymptotic expression for the field, and then used the expression to design the slab parameters for the best isolation. We conclude that both approaches yield the same design for the given configuration. The effectiveness of the design is also validated on a grounded finite slab, which is the representation of the actual application. Finally, multi-layered media for the antenna isolation is studied. GA method is applied with an optimization scheme tailed for a five layered structure. We show that GA converges very fast to the solution and the result yields satisfactory isolation between the antennas.by Jianbing James Chen.Ph.D

13th Annual Review of Progress in Applied Computational Electromagnetics at the Naval Postgraduate School, Monterey, CA, March 17-21, 1997, Conference Proceedings Volumes I & II

Includes Volumes 1 &

Annual Review of Progress in Applied Computational Electromagnetics

Approved for public release; distribution is unlimited

Investigation of general-purpose computing on graphics processing units and its application to the finite element analysis of electromagnetic problems

In this dissertation, the hardware and API architectures of GPUs are investigated, and the corresponding acceleration techniques are applied on the traditional frequency domain finite element method (FEM), the element-level time-domain methods, and the nonlinear discontinuous Galerkin method. First, the assembly and the solution phases of the FEM are parallelized and mapped onto the granular GPU processors. Efficient parallelization strategies for the finite element matrix assembly on a single GPU and on multiple GPUs are proposed. The parallelization strategies for the finite element matrix solution, in conjunction with parallelizable preconditioners are investigated to reduce the total solution time. Second, the element-level dual-field domain decomposition (DFDD-ELD) method is parallelized on GPU. The element-level algorithms treat each finite element as a subdomain, where the elements march the fields in time by exchanging fields and fluxes on the element boundary interfaces with the neighboring elements. The proposed parallelization framework is readily applicable to similar element-level algorithms, where the application to the discontinuous Galerkin time-domain (DGTD) methods show good acceleration results. Third, the element-level parallelization framework is further adapted to the acceleration of nonlinear DGTD algorithm, which has potential applications in the field of optics. The proposed nonlinear DGTD algorithm describes the third-order instantaneous nonlinear effect between the electromagnetic field and the medium permittivity. The Newton-Raphson method is incorporated to reduce the number of nonlinear iterations through its quadratic convergence. Various nonlinear examples are presented to show the different Kerr effects observed through the third-order nonlinearity. With the acceleration using MPI+GPU under large cluster environments, the solution times for the various linear and nonlinear examples are significantly reduced

THE APPLICATION OF DISCONTINUOUS GALKERIN FINITE ELEMENT TIME-DOMAIN METHOD IN THE DESIGN, SIMULATION AND ANALYSIS OF MODERN RADIO FREQUENCY SYSTEMS

The discontinuous Galerkin finite element time-domain (DGFETD) method has been successfully applied to the solution of the coupled curl Maxwell’s equations. In this dissertation, important extensions to the DGFETD method are provided, including the ability to model lumped circuit elements and the ability to model thin-wire structures within a discrete DGFETD solution. To this end, a hybrid DGFETD/SPICE formulation is proposed for high-frequency circuit simulation, and a hybrid DGFETD/Thin-wire formulation is proposed for modeling thin-wire structures within a three-dimensional problem space. To aid in the efficient modeling of open-region structures, a Complex Frequency Shifted-Perfectly Matched Layer (CFS-PML) absorbing medium is applied to the DGFETD method for the first time. An efficient CFS-PML method that reduces the computational complexity and improves accuracy as compared to previous PML formulations is proposed. The methods have been successfully implemented, and a number of test cases are provided that validate the proposed methods. The proposed hybrid formulations and the new CFS-PML formulation dramatically enhances the ability of the DGFETD method to be efficiently applied to simulate complex, state of the art radio frequency systems

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