Publications

This paper presents an overview of an explicit message-passing paradigm for a Eulerian finite volume method for modeling solid dynamics problems involving shock wave propagation, multiple materials, and large deformations. Three-dimensional simulations of high-velocity impact were conducted on the IBM SP2, the SGI Power Challenge Array, and the SGI Origin 2000. The scalability of the message-passing code on distributed-memory and symmetric multiprocessor architectures is presented and compared to the ideal linear performance.

We discuss revolutionary performance advances in selectively oxidized vertical-cavity surface emitting lasers (VCSELs), which have enabled low operating power laser diodes appropriate for aerospace applications. Incorporating buried oxide layers converted from AlGaAs layers within the laser cavity produces enhanced optical and electrical confinement enabling superior laser performance, such as high efficiency and modulation bandwidth. VCSELs are also shown to be viable over varied environmental conditions such as ambient temperature and ionized radiation. The development of novel VCSEL technologies for advanced system applications is also described. Two-dimensional individually addressable VCSEL arrays exhibit uniform threshold and operating characteristics. Bottom emitting 850 nm VCSEL arrays fabricated using wafer fusion are also reported.

LixMn2O4 materials are of considerable interest in battery research and development. The crystal structure of this material can significantly affect the electrochemical performance. The ability to monitor the changes of the crystal structure during use, that is during electrochemical cycling, would prove useful to verify these types of structural changes. We report in-situ XRD measurements of LiMn2O4 cathodes with the use of an electrochemical cell designed for in-situ X-ray analysis. Cells prepared using this cell design allow investigation of the changes in the LiMn2O4 structure during charge and discharge. We describe the variation in lattice parameters along the voltage plateaus and consider the structural changes in terms of the electrochemical results on each cell. Kinetic effects of LiMn2O4 phase changes are also addressed. Applications of the in-situ cell to other compounds such as LiCoO2 cathodes and carbon anodes are presented as well.

This paper discusses how phase plane analysis can be used to describe the overall behaviour of single and multiple autonomous robotic vehicles with finite state machine rules. The importance of this result is that we can begin to design provably asymptotically stable group behaviours from a set of simple control laws and appropriate switching points with decentralized variable structure control. The ability to prove asymptotically stable group behaviour is especially important for applications such as locating military targets or land mines.

We describe a class of microscale heaters fabricated with CMOS processes on silicon wafers. These heaters were designed to produce localized high temperatures above 400°C for test and sensor applications. The temperature levels produced for various input powers and the thermal profiles surrounding the heater for packaged and wafer-level heater structures were studied to guide the placement of microelectronics integrated with the heater structures on the same die. To show the performance of the design, we present resistance sensor measurements, IR temperature profiles, and results from a 3D thermal model of the die. This effort demonstrates that it is possible to successfully operate both a microscale heater and microcircuits on the same die.

The US Department of Energy has developed the Waste Isolation Pilot Plant (WIPP) in the bedded salt of southeastern New Mexico to demonstrate the safe disposal of radioactive transuranic wastes. Four vertical shafts provide access to the underground workings located at a depth of about 660 meters. These shafts connect the underground facility to the surface and potentially provide communication between lithologic units, so they will be sealed to limit both the release of hazardous waste from and fluid flow into the repository. The seal design must consider the potential for fluid flow through a disturbed rock zone (DRZ) that develops in the salt near the shafts. The DRZ, which forms initially during excavation and then evolves with time, is expected to have higher permeability than the native salt. The closure of the shaft openings (i.e., through salt creep) will compress the seals, thereby inducing a compressive back-stress on the DRZ. This back-stress is expected to arrest the evolution of the DRZ, and with time will promote healing of damage. This paper presents laboratory data from tertiary creep and hydrostatic compression tests designed to characterize damage evolution and healing in WIPP salt. Healing is quantified in terms of permanent reduction in permeability, and the data are used to estimate healing times based on considerations of first-order kinetics.

This paper reexamines orientations of shear bands (fault angles) predicted by a theory of shear localization as a bifurcation from homogeneous deformation. In contrast to the Coulomb prediction, which does not depend on deviatoric stress state, the angle between the band normal and the least (most compressive) principal stress increases as the deviatoric stress state varies from axisymmetric compression to axisymmetric extension. This variation is consistent with the data of Mogi (1967) on Dunham dolomite for axisymmetric compression, extension and biaxial compression, but the predicted angles are generally less than observed. This discrepancy may be due to anisotropy that develops due to crack growth in preferred orientations. Results from specialized constitutive relations for axisymmetric compression and plane strain that include this anisotropy indicate that it tends to increase the predicted angles. Measurements for a weak, porous sandstone (Castlegate) indicate that the band angle decreases with increasing inelastic compaction that accompanies increasing mean stress. This trend is consistent with the predictions of the theory but, for this rock, the observed angles are less than predicted.

The development of deep underground structures (e.g., shafts, mines, storage and disposal caverns) significantly alters the stress state in the rock near the structure or opening. The effect of such an opening is to concentrate the far-field stress near the free surface. For soft rock such as salt, the concentrating effect of the opening induces deviatoric stresses in the salt that may be large enough to initiate microcracks which then propagate with time. The volume of rock susceptible to damage by microfracturing is often referred to as the disturbed rock zone and, by its nature, is expected to exhibit high permeability relative to that of the native, far-field rock. This paper presents laboratory data that characterize microfracture-induced damage and the effect this damage has on permeability for bedded salt from the Waste Isolation Pilot Plant located in southeastern New Mexico. Damage is induced in the salt through a series of tertiary creep experiments and quantified in terms of dilatant volumetric strain. The permeability of damaged specimens is then measured using nitrogen gas as the permeant. The range in damage investigated included dilatant volumetric strains from less than 0.03 percent to nearly 4.0 percent. Permeability values corresponding to these damage levels ranged from 1 {times} 10{sup {minus}18} m{sup 2} to 1 {times} 10{sup {minus}12} m{sup 2}. Two simple models were fitted to the data for use in predicting permeability from dilatant volumetric strain.

This paper describes current research and development on a robotic visual servoing system for assembly of LIGA (lithography galvanoforming abforming) parts. The workcell consists of an AMTI robot, precision stage, long working distance microscope, and LIGA fabricated tweezers for picking up the parts. Fourier optics methods are used to generate synthetic microscope images from CAD drawings. These synthetic images are used off-line to test image processing routines under varying magnifications and depths of field. They also provide reference image features which are used to visually servo the part to the desired position.

An electrically injected coupled-resonator vertical-cavity laser (CRVCL) diode is described. The CRVCL consists of a lower 1-λ-thick active resonator containing three InGaAs quantum wells and a passive upper resonator composed of 1-λ-thick GaAs. Some of the characteristics arising from the cavity coupling, including methods for external modulation of the laser are demonstrated.

Microelectromechanical engines that convert the linear outputs from dual orthogonal electrostatic actuators to rotary motion were first developed in 1993. Referred to as microengines, these early devices demonstrated the potential of microelectromechanical technology, but, as expected from any first-of-its-kind device, were not yet optimized. Yield was relatively low, and the 10 micronewtons of force generated by the actuators was not always enough to ensure reliable operation. Since initial development, these engines have undergone a continuous series of significant improvements on three separate fronts: design, fabrication, and electrical activation. Although all three areas will be discussed, emphasis will be on aspects related to mechanical design and generation of the electrical waveforms used to drive these devices. Microtransmissions that dramatically increase torque will also be discussed. Electrostatically driven microengines can be operated at hundreds of thousands of revolutions per minute making large gear reduction ratios feasible; overall ratios of 3,000,000:1 have been successfully demonstrated. Today's microengines have evolved into high endurance (one test device has seen over 7,000,000,000 revolutions), high yield, robust devices that have become the primary actuation source for MicroElectroMechanical Systems (MEMS) at Sandia National Laboratories.

We have sputter-deposited 500-1200 Å thick WSi0.45 metallization onto n+ GaN (n≥1019 cm-3) doped either during MOCVD growth or by direct Si+ ion implantation (5×1015 cm-2, 100 keV) activated by RTA at 1100°C for 30 secs. In the epi samples Rc values of ∼10-14 ω cm2 were obtained, and were stable to ∼1000°C. The annealing treatments up to 600°C had little effect on the WSix/GaN interface, but the beta/-W2N phase formed between 700-800°C, concomitant with a strong reduction (approximately a factor of 2) in near-surface crystalline defects in the GaN. Spiking of the metallization down the threading and misfit dislocations was observed at 800°C, extending >5000 Å in some cases. This can create junction shorting in bipolar or thyristor devices, Rc values of <10-6 ωcm2 were obtained on the implanted samples for 950°C annealing, with values of after 1050°C anneals. The lower Rc values compared to epi samples appear to be a result of the higher peak doping achieved, ∼5×1020 cm-3. We observed wide spreads in Rc values over a wafer surface, with the values on 950°C annealed material ranging from 10-7 to 10-4 ω cm2. There appear to be highly nonuniform doping regions in the GaN, perhaps associated with the high defect density (1010 cm-2) in heteroepitaxial material, and this may contribute to the variations observed. We also believe that near-surface stoichiometry is variable in much of the GaN currently produced due to the relative ease of preferential N2 loss and the common use of HT containing growth (and cool-down) ambients. Finally the ohmic contact behavior of WSix on abrupt and graded composition InxAl1-xN layers has been studied as a function of growth temperature, InN mole fraction x=0.5-1) and post WSix deposition annealing treatment. Rc values in the range 10-3/-10sup-5/ ω cm2 are obtained for auto-doped n+ alloys, with the n-type background being little affected by growth conditions (n∼1020 cm-3). InN is the least temperature-stable alloy (les/700°C), and WSix contact morphology is found to depend strongly on the epi growth conditions.

The classification utility of a dual-antenna interferometric synthetic aperture radar (IFSAR) is explored by comparison of maximum likelihood classification results for synthetic aperture radar (SAR) intensity images and IFSAR intensity and coherence images. The addition of IFSAR coherence improves the overall classification accuracy for classes of trees, water, and fields. A threshold intensity-coherence classifier is also compared to the intensity-only classification results.

The constitutive model used to describe deformation of crushed salt is presented in this paper. Two mechanisms--dislocation creep and grain boundary diffusional pressure solutioning--are combined to form the basis for the constitutive model governing deformation of crushed salt. The constitutive model is generalized to represent three-dimensional states of stress. Recently completed creep consolidation tests are combined with an existing database that includes hydrostatic consolidation and shear consolidation tests conducted on Waste Isolation Pilot Plant (WIPP) and southeastern New Mexico salt to determine material parameters for the constitutive model. Nonlinear least-squares model fitting to data from shear consolidation tests and a combination of shear and hydrostatic tests produces two sets of material parameter values for the model. Changes in material parameter values from test group to test group indicate the empirical nature of the model but show significant improvement over earlier work. To demonstrate the predictive capability of the model, each parameter value set was used to predict each of the tests in the database. Based on fitting statistics and ability of the model to predict test data, the model appears to capture the creep consolidation behavior of crushed salt quite well.

A table top servohydraulic load frame equipped with a laser displacement measurement system was constructed for the mechanical characterization of LIGA fabricated electroforms. A drop-in tensile specimen geometry which includes a pattern to identify gauge length via laser scanning has proven to provide a convenient means to monitor and characterize mechanical property variations arising during processing. In addition to tensile properties, hardness and metallurgical data were obtained for nickel deposit specimens of current density varying between 20 and 80 mA/cm2 from a sulfamate based bath. Data from 80/20 nickel/iron deposits is also presented for comparison. As expected, substantial mechanical property differences from bulk metal properties are observed as well as a dependence of material strength on current density which is supported by grain size variation. While elastic modulus values of the nickel electrodeposit are near 160 GPa, yield stress values vary by over 60%. A strong orientation in the metal electrodeposits as well as variations in nucleating and growth morphology present a concern for anisotropic and geometry dependent mechanical properties within and between different LIGA components.

The first two truncation error terms resulting from finite differencing the convection terms in the two-dimensional Navier-Stokes equations are examined for the purpose of constructing two-dimensional grid generation schemes. These schemes are constructed such that the resulting grid distributions drive the error terms to zero. Two sets of equations result, one for each error term, that show promise in generating grids that provide more accurate flow solutions and possibly faster convergence. One set results in an algebraic scheme that drives the first truncation term to zero, and the other a hyperbolic scheme that drives the second term to zero. Also discussed is the possibility of using the schemes in sequentially constructing a grid in an iterative algorithm involving the flow solver. In essence, the process is envisioned to generate not only a flow field solution but the grid as well. Future work will include applications and extending the method to three dimensions.

Approximately 95% of the world's integrated chips are packaged using a hot, high pressure transfer molding process. The stress created by the flow of silica powder loaded epoxy can displace the fine bonding wires and can even distort the metalization patterns under the protective chip passivation layer [l, 2]. In this study we developed a technique to measure the mechanical stress over the surface of an integrated circuit during the molding process. A CMOS test chip with 25 diffused resistor stress sensors was applied to a commercial lead frame. Both compression and shear stresses were measured at all 25 locations on the surface of the chip every 50 milliseconds during molding. These measurements have a fine time and stress resolution which should allow comparison with computer simulation of the molding process, thus allowing optimization of both the manufacturing process and mold geometry.

Our research is focused on developing inorganic molecular sieve membranes for light gas separations such as hydrogen recovery and natural gas purification, and organic molecular separations, such as chiral enantiomers. We focus on zinc phosphates because of the ease in crystallization of new phases and the wide range of pore sizes and shapes obtained. With our hybrid systems of zinc phosphate crystalline phases templated by amine molecules, we are interested in better understanding the association of the template molecules to the inorganic phase, and how the organic transfers its size, shape, and (in some cases) chirality to the bulk. Furthermore, the new porous phases can also be synthesized as thin films on metal oxide substrates. These films allow us to make membranes from our organic/inorganic hybrid systems, suitable for diffusion experiments. Characterization techniques for both the bulk phases and the thin films include powder X-ray diffraction, TGA, Scanning Electron Micrograph (SEM) and Electron Dispersive Spectrometry (EDS).

The Southwest Surety Institute was formed in June, 1996 by Arizona State University (ASU), New Mexico Institute of Mining and Technology (NM Tech), New Mexico State University (NMSU), and Sandia National Laboratories (SNL) to provide new educational programs in Security Engineering. This is the first science-based program of its kind in the United States, directed at educating Security Engineers to help government and industry address their security needs. Current courses include security system design, evaluation, principles, and technology, the criminal justice system, and each member brings a unique educational capability to the institute. NMSU provides a security technology minor, merging programs in Criminal Justice and Electronics Technology. NM Tech has a formidable explosives testing and evaluation facility. ASU is developing a masters program in Security Engineering at their School of Technology located on a new campus in Mesa, Arizona. The Sandia National Laboratories security system design and evaluation process forms the basis for the security engineering curricula. In an effort to leverage the special capabilities of each university, distance education will be used to share courses among institute members and eventually with other sites across the country. The Institute will also pursue research and development funding in the areas of physical security information security, computer modeling and analysis, and counter-terrorist technology. Individual Institute members are currently working with sponsors from government and industry in areas such as counter-terrorism, microelectronics, banking, aviation, and sensor development.

We characterize in-situ the adhesion of surface micromachined polysilicon beams subject to controlled humidity ambients. Beams were freed by supercritical CO2 drying. Consistent adhesion results were obtained using a post-treatment in an oxygen plasma which rendered the microbeams uniformly hydrophilic. Individual beam deformations were measured by optical interferometry after equilibration at a given relative humidity (RH). Validation of each adhesion measurement was accomplished by comparing the deformations with elasticity theory. The data indicates that adhesion increases exponentially with RH from 30% to 95%, with values from 1 mJ/m2 to 50 mJ/m2. Using the Kelvin equation, we show that the data should be independent of RH if a smooth interface is considered. By modeling a rough interface consistent with atomic force microscopy (AFM) data, the exponential trend is satisfactorily explained.

Interferometric fringe maps are generated by accurately registering a pair of complex SAR images of the same scene imaged from two very similar geometries, and calculating the phase difference between the two images by averaging over a neighborhood of pixels at each spatial location. The phase difference (fringe) map resulting from this IFSAR operation is then unwrapped and used to calculate the height estimate of the imaged terrain. Although the method used to calculate interferometric fringe maps is well known, it is generally executed in a post-processing mode well after the image pairs have been collected. In that mode of operation, there is little concern about algorithm speed and the method is normally implemented on a single processor machine. This paper describes how the interferometric map generation is implemented on a distributed-memory parallel processing machine. This particular implementation is designed to operate on a 16 node Power-PC platform and to generate interferometric maps in near real-time. The implementation is able to accommodate large translational offsets, along with a slight amount of rotation which may exist between the interferometric pair of images. If the number of pixels in the IFSAR image is large enough, the implementation accomplishes nearly linear speed-up times with the addition of processors.

We describe two methods of combining two-pass RADARSAT interferometric phase maps with existing DTED (digital terrain elevation data) to produce improved terrain height estimates. The first is a least-squares estimation procedure that fits the unwrapped phase data to a phase map computed from the DTED. The second is a filtering technique that combines the interferometric height map with the DTED map based on spatial frequency content. Both methods preserve the high fidelity of the interferometric data.

In this paper we review the present status of z-pinches, and predict what the future holds. Although nobody can predict the future, the past 30 years have taught us some lessons that can be applied to the next 30 years.

Modern software development methods combined with key generalizations of standard computational algorithms enable the development of a new class of electromagnetic modeling tools. This paper describes current and anticipated capabilities of a frequency domain modeling code, EIGER, which has an extremely wide range of applicability. In addition, software implementation methods and high performance computing issues are discussed.

The Southwest Surety Institute was formed in June, 1996 by Arizona State University (ASU), New Mexico Institute of Mining and Technology (NM Tech), New Mexico State University (NMSU), and Sandia National Laboratories (SNL) to provide new educational programs in Security Engineering. This is the first science-based program of its kind in the United States, directed at educating Security Engineers to help government and industry address their security needs. Current courses include security system design, evaluation, principles, and technology, the criminal justice system, and each member brings a unique educational capability to the institute. NMSU provides a security technology minor, merging programs in Criminal Justice and Electronics Technology. NM Tech has a formidable explosives testing and evaluation facility. ASU is developing a masters program in Security Engineering at their School of Technology located on a new campus in Mesa, Arizona. The Sandia National Laboratories security system design and evaluation process forms the basis for the security engineering curricula. In an effort to leverage the special capabilities of each university, distance education will be used to share courses among institute members and eventually with other sites across the country. The Institute will also pursue research and development funding in the areas of physical security information security, computer modeling and analysis, and counter-terrorist technology. Individual Institute members are currently working with sponsors from government and industry in areas such as counter-terrorism, microelectronics, banking, aviation, and sensor development.