Alternating-direction implicit (ADI) formulation of the finite-difference time-domain (FDTD) method: algorithm and material dispersion implementation

The alternating-direction implicit finite-difference time-domain (ADI-FDTD) technique is an unconditionally stable time-domain numerical scheme, allowing the /spl Delta/t time step to be increased beyond the Courant-Friedrichs-Lewy limit. Execution time of a simulation is inversely proportional to /spl Delta/t, and as such, increasing /spl Delta/t results in a decrease of execution time. The ADI-FDTD technique greatly increases the utility of the FDTD technique for electromagnetic compatibility problems. Once the basics of the ADI-FDTD technique are presented and the differences of the relative accuracy of ADI-FDTD and standard FDTD are discussed, the problems that benefit greatly from ADI-FDTD are described. A discussion is given on the true time savings of applying the ADI-FDTD technique. The feasibility of using higher order spatial and temporal techniques with ADI-FDTD is presented. The incorporation of frequency dependent material properties (material dispersion) into ADI-FDTD is also presented. The material dispersion scheme is implemented into a one-dimensional and three-dimensional problem space. The scheme is shown to be both accurate and unconditionally stable.     

Analytical multipoint moment matching reduction of distributed thermal networks

In this paper, the multipoint moment matching method for model order reduction of discretized linear thermal networks is extended to distributed linear thermal networks. As a result, from the analytical canonical forms of distributed linear thermal networks, reduced thermal networks are derived analytically. This direct construction of the reduced network, from the exact analytical solutions, avoids the inevitable inaccuracies inherent in conventional surface and volume meshing. It allows nearly exact reduced thermal network construction by domain decomposition for arbitrarily complicated structures.     

Failure analysis of a broken tooth

Several days after heart surgery, a patient discovered his upper right canine tooth had broken at the root. Such tooth damage, recognized post-operatively, is usually assumed to be caused by blunt mechanical force from an instrument used by the anesthesiologist during placement of a breathing tube at the start of surgery. In this case, the patient had saved the crown portion of the broken tooth, and it was possible to examine the root fracture characteristics. The curvature and direction of the crack path and natural tooth situation suggested that failure could be described through a cantilever beam model. This was confirmed when a whole extracted sample tooth was embedded and broken by a measured force in a manner consistent with the model. The resulting fracture surface matched that of the patient’s broken canine tooth. However, the high load and force direction necessary to fracture the root was inconsistent with forces applied during the anesthesia procedure. The failure analysis and further investigation indicated tooth clenching on the breathing tube during recovery was the likely cause of fracture. This paper presents an alternate explanation for intubation-related dental injury, demonstrates the practicality of fractographic analysis of biological materials, and introduces a methodology for simulating in vitro tooth settings for mechanical testing.     

Explosive chemistry: Simulating the chemistry of energetic materials at extreme conditions

In this article, we review recent atomistic computational techniques to study the electronic structure aspects and chemistry of energetic materials at high-pressure and/or high temperature. While several mechanisms have been proposed for the initial events of energetic materials at high-pressure, we explore the validity of a proposed shear-induced local metallization via molecular bond bending in the insensitive explosive TATB. We study the effect of high-stress (both uniform and uniaxial) on the electronic energy band-gap and the first chemical event of a prototypical energetic material, that of nitromethane. We also determine chemical reactions rate laws and decomposition mechanisms from a quantum-based molecular dynamics simulation of HMX, a widely used explosive material, at conditions of high density and temperature similar to that encounter under detonation. Finally, we review a new multi-scale computational tool recently developed to model the shock-induced chemistry of energetic materials at the atomistic level, and report its applicability to shocked solid nitromethane. This revised version was published online in June 2006 with corrections to the Cover Date.     

NANODEV: A Nanoelectronic-Device Simulation Software System

The concept and structure of the NANODEV simulation software are described. NANODEV deals with nanoelectronic devices that exploit single-electron tunneling, resonant tunneling, or quantum interference. It can use both simplified and sophisticated models and enables one to evaluate a wide variety of devices and configurations. The capabilities of NANODEV are illustrated by examples.     

Evanescently coupled photodiodes integrating a double-stage taper for 40-Gb/s applications-compared performance with side-illuminated photodiodes

The design, fabrication, and performance of double-stage taper photodiodes (DSTPs) are reported. The objective of this work is to develop devices compatible with 40-Gb/s applications. Such devices require high efficiency, ultrawide band, high optical power handling capability, and compatibility with low-cost module fabrication. The integration of mode size converters improves both the coupling efficiency and the responsivity with a large fiber mode diameter. Responsivity of 0.6 A/W and 0.45 A/W are achieved with a 6-/spl mu/m fiber mode diameter and cleaved fiber, respectively, providing relaxed alignment tolerances (/spl plusmn/1.6 /spl mu/m and /spl plusmn/2 /spl mu/m, respectively), compatible with cost-effective packaging techniques. DSTPs also offer a wide bandwidth greater than 40 GHz and transverse-electric/transverse-magnetic polarization dependence lower than 0.2 dB. Furthermore, a DSTP saturation current as high as 11 mA results in optical power handling greater than +10 dBm and a high output voltage of 0.8 V. These capabilities allow the photodiode to drive the decision circuit without the need of a broad-band electrical amplifier. The DSTP devices presented here demonstrate higher responsivities with large fiber mode diameter and better optical power handling capabilities and are compared with classical side-illuminated photodiodes.     

LiNbO3 optical waveguides formed in a new proton source

Proton-exchanged planar waveguides have been fabricated on Z-cut and X-cut lithium niobate crystals by using a new proton source formed by a mixture of benzoic and adipic acids. Waveguide index profiles and optical characteristics have been obtained at different values of the adipic-benzoic acid concentration ratio. The samples have been structurally characterized by Raman and infrared (IR) absorption spectroscopy and double-crystal X-ray diffraction. Good quality samples have been fabricated by using 30 mol% ratio dilution, showing very low scattering levels (<0.1 dB/cm), relatively high electrooptic coefficient (r₃₃=0.88 pm/V), and low relative percentage of interstitial protons (26%). The main factor limiting the waveguide optical properties is the substitutional-interstitial proton ratio, which can be easily controlled to produce good quality waveguides. A demonstration of the repeatability of the exchange process in the acid mixture is also provided     

Exact pairwise error probability of space-time codes

We describe a simple technique for the numerical calculation, within any desired degree of accuracy, of the pairwise error probability (PEP) of space-time codes over fading channels. This method applies also to the calculation of E[Q(√ξ)] for any nonnegative random variable ξ whose moment-generating function Φ_ξ(s)=E[exp(-sξ)] is known. Its application to the multiple antenna independent Rayleigh-fading channel and to the Rayleigh block fading channel is discussed, and illustrated by two simple examples     

The measurement of solar array discharges in a simulated spaceenvironment

A spacecraft in a plasma builds up charge on all the dielectric surfaces and interfaces. Once the net charge exceeds the dielectric breakdown of the material, a discharge occurs. One of the more susceptible pieces of equipment is the antenna/receiver system. The radiated E-field may be strong enough to create an ambiguous signal which may be misinterpreted by the system electronics and cause a system malfunction. A technique is developed to monitor the radiated E-field of materials discharging in an electron environment, using vacuum chambers for measuring the material discharges which are made of highly reflective materials. These chambers affect the radiated E-field due to multiple reflections from the walls. The technique developed defines a method for correcting the effects caused by the measurement facilities. The methodology is: monitor the radiated E-field with a broadband dipole antenna, and digitize the radiated signal as a function of time. Determine the frequency response of the radiated E-field using an FFT algorithm. Measure the transmission and reflection characteristics of the two-port network inside the measurement chamber, and determine the impedance network from the measured E-parameters across the frequency band of interest. Transform the measured E-field frequency response through the impedance network to obtain the frequency response of the actual radiated discharge current. Find the inverse FFT of this response to obtain the actual radiated discharge current response. This technique aids in the prediction of the E-field coupling into receive antennas on-board actual satellites     

High-frequency analysis of the Ellipticus antenna

Efficient numerical solution techniques have been developed and used to examine the electromagnetic fields that can be developed in the working volume of the CW Ellipticus antenna operated at frequencies from 100 kHz to 1 GHz. An exponentially tapered transition section is designed to obtain the desired illumination pattern in the working volume. The input transition section is needed for impedance matching and to drive efficiently the Ellipticus antenna. A parametric study is performed to ascertain the performance of the Ellipticus antenna for frequencies up to 1 GHz