Hypercube matrix computation task

The Hypercube Matrix Computation (Year 1986-1987) task investigated the applicability of a parallel computing architecture to the solution of large scale electromagnetic scattering problems. Two existing electromagnetic scattering codes were selected for conversion to the Mark III Hypercube concurrent computing environment. They were selected so that the underlying numerical algorithms utilized would be different thereby providing a more thorough evaluation of the appropriateness of the parallel environment for these types of problems. The first code was a frequency domain method of moments solution, NEC-2, developed at Lawrence Livermore National Laboratory. The second code was a time domain finite difference solution of Maxwell's equations to solve for the scattered fields. Once the codes were implemented on the hypercube and verified to obtain correct solutions by comparing the results with those from sequential runs, several measures were used to evaluate the performance of the two codes. First, a comparison was provided of the problem size possible on the hypercube with 128 megabytes of memory for a 32-node configuration with that available in a typical sequential user environment of 4 to 8 megabytes. Then, the performance of the codes was anlyzed for the computational speedup attained by the parallel architecture

Application of coupled-wave Wentzel-Kramers-Brillouin approximation to ground penetrating radar

This paper deals with bistatic subsurface probing of a horizontally layered dielectric half-space by means of ultra-wideband electromagnetic waves. In particular, the main objective of this work is to present a new method for the solution of the two-dimensional back-scattering problem arising when a pulsed electromagnetic signal impinges on a non-uniform dielectric half-space; this scenario is of interest for ground penetrating radar (GPR) applications. For the analytical description of the signal generated by the interaction of the emitted pulse with the environment, we developed and implemented a novel time-domain version of the coupled-wave Wentzel-Kramers-Brillouin approximation. We compared our solution with finite-difference time-domain (FDTD) results, achieving a very good agreement. We then applied the proposed technique to two case studies: in particular, our method was employed for the post-processing of experimental radargrams collected on Lake Chebarkul, in Russia, and for the simulation of GPR probing of the Moon surface, to detect smooth gradients of the dielectric permittivity in lunar regolith. The main conclusions resulting from our study are that our semi-analytical method is accurate, radically accelerates calculations compared to simpler mathematical formulations with a mostly numerical nature (such as the FDTD technique), and can be effectively used to aid the interpretation of GPR data. The method is capable to correctly predict the protracted return signals originated by smooth transition layers of the subsurface dielectric medium. The accuracy and numerical efficiency of our computational approach make promising its further development

A new hybrid implicit-explicit FDTD method for local subgridding in multiscale 2-D TE scattering problems

The conventional finite-difference time-domain (FDTD) method with staggered Yee scheme does not easily allow including thin material layers, especially so if these layers are highly conductive. This paper proposes a novel subgridding technique for 2-D problems, based on a hybrid implicit-explicit scheme, which efficiently copes with this problem. In the subgrid, the new method collocates field components such that the thin layer boundaries are defined unambiguously. Moreover, aspect ratios of more than a million do not impair the stability of the method and allow for very accurate predictions of the skin effect. The new method retains the Courant limit of the coarse Yee grid and is easily incorporated into existing FDTD codes. A number of illustrative examples, including scattering by a metal grating, demonstrate the accuracy and stability of the new method

Finite difference time domain calculations of antenna mutual coupling

The Finite Difference Time Domain (FDTD) technique was applied to a wide variety of electromagnetic analysis problems, including shielding and scattering. However, the method has not been exclusively applied to antennas. Here, calculations of self and mutual admittances between wire antennas are made using FDTD and compared with results obtained during the method of moments. The agreement is quite good, indicating the possibilities for FDTD application to antenna impedance and coupling

Plasmonic nanoparticle monomers and dimers: From nano-antennas to chiral metamaterials

We review the basic physics behind light interaction with plasmonic nanoparticles. The theoretical foundations of light scattering on one metallic particle (a plasmonic monomer) and two interacting particles (a plasmonic dimer) are systematically investigated. Expressions for effective particle susceptibility (polarizability) are derived, and applications of these results to plasmonic nanoantennas are outlined. In the long-wavelength limit, the effective macroscopic parameters of an array of plasmonic dimers are calculated. These parameters are attributable to an effective medium corresponding to a dilute arrangement of nanoparticles, i.e., a metamaterial where plasmonic monomers or dimers have the function of "meta-atoms". It is shown that planar dimers consisting of rod-like particles generally possess elliptical dichroism and function as atoms for planar chiral metamaterials. The fabricational simplicity of the proposed rod-dimer geometry can be used in the design of more cost-effective chiral metamaterials in the optical domain.Comment: submitted to Appl. Phys.

Hypercube matrix computation task

A major objective of the Hypercube Matrix Computation effort at the Jet Propulsion Laboratory (JPL) is to investigate the applicability of a parallel computing architecture to the solution of large-scale electromagnetic scattering problems. Three scattering analysis codes are being implemented and assessed on a JPL/California Institute of Technology (Caltech) Mark 3 Hypercube. The codes, which utilize different underlying algorithms, give a means of evaluating the general applicability of this parallel architecture. The three analysis codes being implemented are a frequency domain method of moments code, a time domain finite difference code, and a frequency domain finite elements code. These analysis capabilities are being integrated into an electromagnetics interactive analysis workstation which can serve as a design tool for the construction of antennas and other radiating or scattering structures. The first two years of work on the Hypercube Matrix Computation effort is summarized. It includes both new developments and results as well as work previously reported in the Hypercube Matrix Computation Task: Final Report for 1986 to 1987 (JPL Publication 87-18)

Comparison of FDTD Algorithms for Subcellular Modeling of Slots in Shielding Enclosures

Subcellular modeling of thin slots in the finite-difference time-domain (FDTD) method is investigated. Two subcellular algorithms for modeling thin slots with the FDTD method are compared for application to shielding end osures in electromagnetic compatibility (EMC). The stability of the algorithms is investigated, and comparisons between the two methods for slots in planes, and slots in loaded cavities are made. Results for scattering from a finite-length slot in an infinite plane employing one of the algorithms are shown to agree well with published experimental results, and power delivered to an enclosure with a slot agree well with results measured for this study

Direct and Inverse Computational Methods for Electromagnetic Scattering in Biological Diagnostics

Scattering theory has had a major roll in twentieth century mathematical physics. Mathematical modeling and algorithms of direct,- and inverse electromagnetic scattering formulation due to biological tissues are investigated. The algorithms are used for a model based illustration technique within the microwave range. A number of methods is given to solve the inverse electromagnetic scattering problem in which the nonlinear and ill-posed nature of the problem are acknowledged.Comment: 61 pages, 5 figure