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.     

Scale Up from Small Oven-Drying Tests of Mineral Concentrate to Pilot-Scale Drying with a Heated Pad

While Fickian diffusion models are commonly used in other applications, there are few reports of them being applied to the batch drying of a mineral concentrate. Diffusion coefficients estimated from small-scale oven-drying tests were used to predict the drying behavior of a concentrate sample 1 m × 1 m in area and 50 cm deep, with a heated bottom pad. These pilot-scale tests included both daily turning of the sample and turning every three days. The excellent quantitative agreement between the predicted and observed pilot-scale behavior gives a high level of confidence in the model predictions and suggests that a Fickian diffusion model is adequate to predict the behavior of mineral concentrates at the low moisture contents used here.     

Changes in the gaseous H2S barrier properties of damaged silicate coated nylon 6 films

The damage imposed on SiO _x deposited nylon 6 films as a result of abrasion with a cotton cloth and Gelboflex testing was examined by evaluating the rate at which copper plates, which were enveloped by the damaged films, were corroded by H₂S. Abrasion with a cotton cloth caused some micro-cracking of the SiO _x layer and the permeation rate of H₂S approached that of the uncoated nylon 6 film. Damage to the SiO _x layer by twisting and crushing progressed gradually with the number of Gelboflex test cycles and correspondingly the corrosion rate of the copper plates increased. Comparison of the corrosion rates of the copper plates kept in the pouches made of various commercial films with those obtained for the damaged SiO _x deposited nylon 6 films showed a clear relationship between the H₂ permeation rate of the films and the corrosion rate of the copper plates by H₂S.     

The rank and kernel of extended 1-perfect /spl Zopf//sub 4/-linear and additive non-/spl Zopf//sub 4/-linear codes

A binary extended 1-perfect code of length n + 1 = 2/sup t/ is additive if it is a subgroup of /spl Zopf//sub 2//sup /spl alpha// /spl times/ /spl Zopf//sub 4//sup /spl beta//. The punctured code by deleting a /spl Zopf//sub 2/ coordinate (if there is one) gives a perfect additive code. 1-perfect additive codes were completely characterized and by using that characterization we compute the possible parameters /spl alpha/, /spl beta/, rank, and dimension of the kernel for extended 1-perfect additive codes. A very special case is that of extended 1-perfect /spl Zopf//sub 4/-linear codes.     

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.     

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     

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     

Crosstalk between microstrip transmission lines

Methods for prediction of crosstalk between microstrip transmission lines are reviewed and simplified for the weak-coupling case. Classical coupled transmission line theory is used for uniform lines, and potential and induced EMF methods are used for crosstalk between nonuniform lines. It is shown that the potential method is equivalent to classical coupled transmission line theory for the case of uniform lines. An experiment was performed for uniform coupled microstrip lines for frequencies from 50 MHz to 5 GHz, and good agreement between theory and measurement was obtained for both near- and far-end crosstalk     

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