Publications

The void percolation threshold is calculated for a distribution of overlapping spheres with equal radii, and for a binary sized distribution of overlapping spheres, where half of the spheres have radii twice as large as the other half. Using systems much larger than previous work, the authors determine a much more precise value for the percolation thresholds and correlation length exponent. The values for the percolation thresholds are shown to be significantly different, in contrast with previous, less precise works that speculated that the threshold might be universal with respect to sphere size distribution.

Pressure studies have provided new insights into the physics of compositionally-disordered ABO{sub 3} oxide relaxors. Specifically, results will be presented and discussed on a pressure-induced ferroelectric-to-relaxer crossover phenomenon, the continuous evolution of the energetic and dynamics of the relaxation process, and the interplay between pressure and electric field in determining the dielectric response.

The microstructure of lateral composition modulation in InAs/AlAs superlattices grown by MBE on InP is examined. The use of x-ray diffraction, TEM, AFM, and STEM to characterize the modulations is discussed. Combining the information from these techniques gives increased insight into the phenomenon and how to manipulate it. Diffraction measures the intensity of modulation and its wavelength, and is used to identify growth conditions giving strong modulation. The TEM and STEM analyses indicate that local compositions are modulated by as much as 0.38 InAs mole fraction. Plan-view images show that modulated structures consists of short ({approx_lt}0.2 {micro}m) In-rich wires with a 2D organization in a (001) growth plane. However, growth on miscut substrates can produce a single modulation along the miscut direction with much longer wires ({approx_gt}0.4 {micro}m), as desired for potential applications. Photoluminescence studies demonstrate that the modulation has large effects on the bandgap energy of the superlattice.

The authors present a new technique for the design of approximation algorithms that can be viewed as a generalization of randomized rounding. They derive new or improved approximation guarantees for a class of generalized congestion problems such as multicast congestion, multiple TSP etc. Their main mathematical tool is a structural decomposition theorem related to the integrality gap of a relaxation.

We present a split Hopkinson pressure bar technique to obtain compressive stress-strain data for rock materials. This technique modifies the conventional split Hopkinson bar apparatus by placing a thin copper disk on the impact surface of the incident bar. When the copper disk is impacted by the striker bar, a nondispersive ramp pulse propagates in the incident bar and produces a nearly constant strain rate in a rock sample. Data from experiments with limestone show that the samples are in dynamic stress equilibrium and have constant strain rates over most of the duration of the tests. We also present analytical models that predict the time durations for sample equilibrium and constant strain rate. Model predictions are in good agreement with measurements.

The influence of the substrate on the translational and orientational ordering in sub-monolayer films of passivated multiply-twinned gold clusters has been investigated using high resolution and dark field transmission electron microscopy. Although clear differences were observed in the degree of translational ordering on amorphous carbon and etched silicon substrates, there was no corresponding variation in the crystallographic orientation of the nanocrystal cores. The results demonstrate that the orientation of passivated clusters with multiply-twinned cores is effectively random with respect to both the superlattice and the substrate.

This paper examines the impact of inserting Micro-Electro-Mechanical Systems (MEMS) into US defense applications. As specific examples, the impacts of micro Inertial Measurement Units (IMUs), radio frequency MEMS (RF MEMS), and Micro-Opto-Electro-Mechanical Systems (MOEMS) to provide integrated intelligence, communication, and control to the defense infrastructure with increased affordability, functionality, and performance are highlighted.

In this paper the authors present performance results from several parallel benchmarks and applications on a 400-node Linux cluster at Sandia National Laboratories. They compare the results on the Linux cluster to performance obtained on a traditional distributed-memory massively parallel processing machine, the Intel TeraFLOPS. They discuss the characteristics of these machines that influence the performance results and identify the key components of the system software that they feel are important to allow for scalability of commodity-based PC clusters to hundreds and possibly thousands of processors.

We report the growth and characterization of quaternary AlGaInN. A combination of photoluminescence (PL), high-resolution x-ray diffraction (XRD), and Rutherford backscattering spectrometry (RBS) characterizations enables us to explore the contours of constant PL peak energy and lattice parameter as functions of the quaternary compositions. The observation of room temperature PL emission at 351nm (with 20% Al and 5% In) renders initial evidence that the quaternary could be used to provide confinement for GaInN (and possibly GaN). AlGaInN/GrdnN MQW heterostructures have been grown; both XRD and PL measurements suggest the possibility of incorporating this quaternary into optoelectronic devices.

We describe parallel simulations of viscous, incompressible, free surface, Newtonian fluid flow problems that include dynamic contact lines. The Galerlin finite element method was used to discretize the fully-coupled governing conservation equations and a ''pseudo-solid'' mesh mapping approach was used to determine the shape of the free surface. In this approach, the finite element mesh is allowed to deform to satisfy quasi-static solid mechanics equations subject to geometric or kinematic constraints on the boundaries. As a result, nodal displacements must be included in the set of problem unknowns. Issues concerning the proper constraints along the solid-fluid dynamic contact line in three dimensions are discussed. Parallel computations are carried out for an example taken from the coating flow industry, flow in the vicinity of a slot coater edge. This is a three-dimensional free-surface problem possessing a contact line that advances at the web speed in one region but transitions to static behavior in another part of the flow domain. Discussion focuses on parallel speedups for fixed problem size, a class of problems of immediate practical importance.

The analysis of many complex multiphase fluid flow systems is based on a scale decoupling procedure. At the macroscale continuum models are used to perform large-scale simulations. At the mesoscale statistical homogenization theory is used to derive continuum models based on representative volume elements (RVEs). At the microscale small-scale features, such as interfacial properties, are analyzed to be incorporated into mesoscale simulations. In this research mesoscopic simulations of hard particles suspended in a Newtonian fluid undergoing nonlinear shear flow are performed using a boundary element method. To obtain an RVE at higher concentrations, several hundred particles are included in the simulations, putting considerable demands on the computational resources both in terms of CPU and memory. Parallel computing provides a viable platform to study these large multiphase systems. The implementation of a portable, parallel computer code based on the boundary element method using a block-block data distribution is discussed in this paper. The code employs updated direct-solver technologies that make use of dual-processor compute nodes.

Al{sub 2}Cu thin films ({approx} 382 nm) are fabricated by melting and resolidifying Al/Cu bilayers in the presence of a {micro} 3 nm Al{sub 2}O{sub 3} passivating layer. X-ray Photoelectron Spectroscopy (XPS) measures a 1.0 eV shift of the Cu2p{sub 3/2} peak and a 1.6 eV shift of the valence band relative to metallic Cu upon Al{sub 2}Cu formation. Scanning Electron microscopy (SEM) and Electron Back-Scattered Diffraction (EBSD) show that the Al{sub 2}Cu film is composed of 30-70 {micro}m wide and 10-25 mm long cellular grains with (110) orientation. The atomic composition of the film as estimated by Energy Dispersive Spectroscopy (EDS) is 67 {+-} 2% Al and 33 {+-} 2% Cu. XPS scans of Al{sub 2}O{sub 3}/Al{sub 2}Cu taken before and after air exposure indicate that the upper Al{sub 2}Cu layers undergo further oxidation to Al{sub 2}O{sub 3} even in the presence of {approx} 5 nm Al{sub 2}O{sub 3}. The majority of Cu produced from oxidation is believed to migrate below the Al{sub 2}O{sub 3} layers, based upon the lack of evidence for metallic Cu in the XPS scans. In contrast to Al/Cu passivated with Al{sub 2}O{sub 3}, melting/resolidifying the Al/Cu bilayer without Al{sub 2}O{sub 3} results in phase-segregated dendritic film growth.

For the first time, we show that redox-sensitive metals, which are highly soluble in the oxidized state can be reduced and precipitated from aqueous solution using tin protoporphyrin and light in the presence of an electron donor. Hg{sup 2+} and Cu{sup 2+} were reduced to the metallic state, and Ub{sup 6+} precipitated as oxide with very low volubility, suggesting that removal of these metals via reductive photoreduction and precipitation may be an innovative way for wastewater treatment. Ag{sup 2+} and Au{sup 2+} were reduced to the metallic state and precipitated as nanoparticles. Finally, using tin porphyrins and light for a variety of purposes involving reactions that require a low redox potential may be a good step toward energy conservation and environmentally benign processing.

We report on reduction and precipitation of Se(VI), Pb(II), CU(II), U(VI), Mo(VI), and Cr(VI) in water by cytochrome c{sub 3} isolated from Desulfomicrobium baczdatum [strain 9974]. The tetraheme protein cytochrome c{sub 3} was reduced by sodium dithionite. Redox reactions were monitored by UV-visible spectroscopy of cytochrome c{sub 3}. Analytical electron microscopy work showed that Se(VI), Pb(II), and CU(II) were reduced to the metallic state, U(W) and Mo(W) to U(IV) and Mo(IV), respectively, and Cr(VI) probably to Cr(III). U(IV) and Mo(W) precipitated as oxides and Cr(III) as an amorphous hydroxide. Cytochrome c{sub 3} was used repeatedly in the same solution without loosing its effectiveness. The results suggest usage of cytochrome c{sub 3} to develop innovative and environmentally benign methods to remove heavy metals from waste- and groundwater.

We have investigated the thermochromic transition of an ultrathin poly(diacetylene) film. The Langmuir film is composed of three layers of polymerized 10,12-pentacosadiynoic acid [CH{sub 3}(CH{sub 2}){sub 11}C{triple_bond}CC{triple_bond}C(CH{sub 2}){sub 8}COOH] (poly-PCDA) organized into crystalline domains on a silicon substrate. Spectroscopic ellipsometry and fluorescence intensity measurements are obtained with in-situ temperature control. Poly-PCDA films exhibit a reversible thermal transition between the initial blue form and an intermediate ''purple'' form that exists only at elevated temperature (between 303-333 K), followed by an irreversible transition to the red form after annealing above 320 K. We propose that the purple form is thermally distorted blue poly-PCDA, and may represent a transitional configuration in the irreversible conversion to red. This hypothesis is supported by the appearance of unique features in the absorption spectra for each form as derived from the ellipsometry measurements. Significant fluorescence emission occurs only with the red form, and is reduced at elevated temperatures while the absorption remains unchanged. Reduced emission is likely related to thermal fluctuations of the hydrocarbon side chains. Time-resolved fluorescence measurements of the irreversible transition have been performed. Using a first-order kinetic analysis of these measurements we deduce an energy barrier of 17.6 {+-} 1.1 kcal mol{sup -1} between the blue and red forms.

The interfacial force microscope (IFM) was used to indent and image defect free Au(111) surfaces, providing atomic-scale observations of the onset of pileup and the excursion of material above the initial surface plane. Images and load-displacement measurements demonstrate that elastic accommodation of an indenter is followed by two stages of plasticity. The initial stage is identified by slight deviations of the load-displacement relationship from the predicted elastic response. Images acquired after indentations showing only this first stage indicate that these slight load relaxation events result in residual indentations 0.5 to 4 nm deep with no evidence of pileup or surface orientation dependence. The second stage of plasticity is marked by a series of dramatic load relaxation events and residual indentations tens of nanometers deep. Images acquired following this second stage document 0.25 nm pileup terraces which reflect the crystallography of the surface as well as the indenter geometry. Attempts to plastically displace the indenter 4-10 nanometers deep into the Au(111) surface were unsuccessful, demonstrating that the transition from stage I to stage H plasticity is associated with overcoming some sort of barrier. Stage I is consistent with previously reported models of dislocation nucleation. The dramatic load relaxations of stage II plasticity, and the pileup of material above the surface, require cross-slip and appear to reflect a dynamic process leading to dislocation intersection with the surface. The IFM measurements reported here offer new insights into the mechanisms underlying the very early stages of plasticity and the formation of pileup.

Gamma-densitometry tomography is applied to study the effect of sparger hole geometry, gas flow rate, column pressure, and phase properties on gas volume fraction profiles in bubble columns. Tests are conducted in a column 0.48 m in diameter, using air and mineral oil, superficial gas velocities ranging from 5 to 30 cm s{sup -1}, and absolute column pressures from 103 to 517 kPa. Reconstructed gas volume fraction profiles from two sparger geometries are presented. The development length of the gas volume fraction profile is found to increase with gas flow rate and column pressure. Increases in gas flow rate increase the local gas volume fraction preferentially on the column axis, whereas increases in column pressure produce a uniform rise in gas volume fraction across the column. A comparison of results from the two spargers indicates a significant change in development length with the number and size of sparger holes.

The failure of thermosetting polymer adhesives is an important problem which particularly lacks understanding from the molecular viewpoint. While linear elastic fracture mechanics works well for such polymers far from the crack tip, the method breaks down near the crack tip where large plastic deformation occurs and the molecular details become important [1]. Results of molecular dynamics simulations of highly crosslinked polymer networks bonded to a solid surface are presented here. Epoxies are used as the guide for modeling. The focus of the simulations is the network connectivity and the interfacial strength. In a random network, the bond stress is expected to vary, and the most stressed bonds will break first [2]. Crack initiation should occur where a cluster of highly constrained bonds exists. There is no reason to expect crack initiation to occur at the interface. The results to be presented show that the solid surface limits the interfacial bonding resulting in stressed interfacial bonds and interfacial fracture. The bonds in highly-crosslinked random networks do not become stressed as expected. The sequence of molecular structural deformations that lead to failure has been determined and found to be strongly dependent upon the network connectivity. The structure of these networks and its influence on the stress-strain behavior will be discussed in general. A set of ideal, ordered networks have been constructed to manipulate the deformation sequence to achieve different fracture modes (i.e. cohesive vs. adhesive).

The generalized method of Backus and Gilbert (BG) is described and applied to the inverse problem of obtaining spectra from a 5-channel, filtered array of x-ray detectors (XRD's). This diagnostic is routinely fielded on the Z facility at Sandia National Laboratories to study soft x-ray photons ({le}2300 eV), emitted by high density Z-pinch plasmas. The BG method defines spectral resolution limits on the system of response functions that are in good agreement with the unfold method currently in use. The resolution so defined is independent of the source spectrum. For noise-free, simulated data the BG approximating function is also in reasonable agreement with the source spectrum (150 eV black-body) and the unfold. This function may be used as an initial trial function for iterative methods or a regularization model.

The dynamical behavior of Cl{sup {minus}} and H{sub 2}O molecules in the interlayer and on the (001) surface of the Ca-aluminate hydrate hydrocalumite (Friedel's salt) over a range of temperatures from {minus}100 to 300 C is studied using the technique of isothermal-isobaric molecular dynamics computer simulations. This phase is currently the best available model compound for other, typically more disordered, mixed-metal layered hydroxides. The computed crystallographic parameters and density are in good agreement with available X-ray diffraction data and the force field developed for these simulations preserves the structure and density to within less than 2% of their measured values. In contrast to the highly ordered arrangement of the interlayer water molecules interpreted from the X-ray data, the simulations reveal significant dynamic disorder in water orientations. At all simulated temperatures, the interlayer water molecules undergo rapid librations (hindered hopping rotations) around an axis essentially perpendicular to the layers. This results in breaking and reformation of hydrogen bonds with the neighboring Cl{sup {minus}} anions and in a time-averaged nearly uniaxial symmetry at Cl{sup {minus}}, in good agreement with recent {sup 35}Cl NMR measurements. Power spectra of translational, vibrational, and vibrational motions of interlayer and surface Cl{sup {minus}} and H{sub 2}O were calculated as Fourier transforms of the atomic velocity autocorrelation functions and compared with the corresponding spectra and dynamics for a bulk aqueous solution. The ordered interlayer space has significant effects on the motions. Strong electrostatic attraction between interlayer water molecules and Ca atoms in the principal layer makes the Ca{hor_ellipsis}OH{sub 2} bond direction the preferred axis for interlayer water librations. The calculated diffusion coefficient of Cl{sup {minus}} as an outer-sphere surface complex is almost three times that of inner-sphere Cl{sup {minus}}, but is still about an order of magnitude less than that of Cl{sup {minus}} in bulk aqueous solution at the same temperature.

This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity--transport equation for nonturbulent fluctuation growth based on that proposed by Warren and Hassan is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittency function based on the work of Dhawan and Narasimha. The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the grid-dependence of selected predictions is analyzed.

Many Navier-Stokes codes require that the governing equations be written in conservation form with a source term. The Spalart-Allmaras one-equation model was originally developed in substantial derivative form and when rewritten in conservation form, a density gradient term appears in the source term. This density gradient term causes numerical problems and has a small influence on the numerical predictions. Further work has been performed to understand and to justify the neglect of this term. The transition trip term has been included in the one-equation eddy viscosity model of Spalart-Allmaras. Several problems with this model have been discovered when applied to high-speed flows. For the Mach 8 flat plate boundary layer flow with the standard transition method, the Baldwin-Barth and both k-{omega} models gave transition at the specified location. The Spalart-Allmaras and low Reynolds number k-{var_epsilon} models required an increase in the freestream turbulence levels in order to give transition at the desired location. All models predicted the correct skin friction levels in both the laminar and turbulent flow regions. For Mach 8 flat plate case, the transition location could not be controlled with the trip terms as given in the Spalart-Allmaras model. Several other approaches have been investigated to allow the specification of the transition location. The approach that appears most appropriate is to vary the coefficient that multiplies the turbulent production term in the governing partial differential equation for the eddy viscosity (Method 2). When this coefficient is zero, the flow remains laminar. The coefficient is increased to its normal value over a specified distance to crudely model the transition region and obtain fully turbulent flow. While this approach provides a reasonable interim solution, a separate effort should be initiated to address the proper transition procedure associated with the turbulent production term. Also, the transition process might be better modeled with the Spalart-Allmaras turbulence model with modification of the damping function f{sub v1}. The damping function could be set to zero in the laminar flow region and then turned on through the transition flow region.

Certain applications for pulse power require narrow, high current pulses for their implementation. This work was performed to determine if MCTS (MOS Controlled Thyristors) could be used for these applications. The MCTS were tested as discharge switches in a low inductance circuit delivering 1 {micro}s pulses at currents between roughly 3 kA and 11 kA, single shot and repetitively at 1, 10 and 50 Hz. Although up to 9000 switching events could be obtained, all the devices failed at some combination of current and repetition rate. Failure was attributed to temperature increases caused by average power dissipated in the thyristor during the switching sequence. A simulation was performed to confirm that the temperature rise was sufficient to account for failure. Considerable heat sinking, and perhaps a better thermal package, would be required before the MCT could be considered for pulse power applications.

Inverter type thyristors were switched repetitively to failure with 1 {micro}s pulses at repetition rates of 10, 50 and 100 pps and at peak currents up to 12 kA. Millions of pulses could be obtained before failure if the peak current were held at around 6 kA.

The potential for cooperation between India and Pakistan is substantial. Topics as widely varying as national security, the environment and trade hold the potential for improved bilateral relations. This paper looks at a few areas in which monitoring technology could contribute to enhancing cooperative border agreements between the two nations. The goal of the paper is not to provide prescriptive solutions to regional problems, but to expand the number of options being considered for improving Indian-Pakistan relations. Many of the impediments to bilateral progress are a result of a history of conflict and mistrust. By utilizing technical monitoring and inspections, each side can begin to replace suspicion and doubt with knowledge and information useful in making informed political, economic and military decisions. At the same time, technical monitoring and inspections can build confidence through common interactions. India and Pakistan have pledged to resolve their disputes, including Kashmir, through dialogue. Implementation of that pledge is influenced by a number of factors, including changes in the political systems and the fortunes of the leadership. Events of the past year and a half have severely tested these two governments' ability to move forward along a constructive and positive path. Testing of new missile systems both preceded and followed testing of nuclear weapons in May 1998. Both countries disregarded subsequent international displeasure as they proceeded to openly declare their respective nuclear capability. Their brief engagement with each other in February 1999 and movement toward a rapprochement diluted international condemnation of their nuclear activity. Within a recent period of nine months however, progress in the dialogue has been stalled first by the Pakistani move in Kashmir in May 1999, then by the Indian election in the summer of 1999 and most recently by the military coup in Pakistan.

An assessment was made of the manufacturability of hybrid microcircuit test vehicles assembled using three Pb-free solder compositions 96.5Sn--3.5Ag (wt.%), 91.84Sn--3.33Ag--4.83Bi, and 86.85Sn--3.15Ag--5.0Bi--5.0Au. The test vehicle substrate was 96% alumina; the thick film conductor composition was 76Au--21Pt--3Pd. Excellent registration between the LCCC or chip capacitor packages and the thick film solder pads was observed. Reduced wetting of bare (Au-coated) LCCC castellations was eliminated by hot solder dipping the I/Os prior to assembly of the circuit card. The Pb-free solders were slightly more susceptible to void formation, but not to a degree that would significantly impact joint functionality. Microstructural damage, while noted in the Sn-Pb solder joints, was not observed in the Pb-free interconnects.

High-speed multiplication is frequently used in general-purpose and application-specific computer systems. These systems often support integer multiplication, where two n-bit integers are multiplied to produce a 2n-bit product. To prevent growth in word length, processors typically return the n least significant bits of the product and a flag that indicates whether or not overflow has occurred. Alternatively, some processors saturate results that overflow to the most positive or most negative representable number. This paper presents efficient methods for performing unsigned or two's complement integer multiplication with overflow detection or saturation. These methods have significantly less area and delay than conventional methods for integer multiplication with overflow detection and saturation.

The authors have demonstrated a functional MOCVD-grown AlGaAs/InGaAsN/GaAsPnP DHBT that is lattice matched to GaAs and has a peak current gain ({beta}) of 25. Because of the smaller bandgap (E{sub g}=1.20eV)of In{sub 0.03}Ga{sub 0.97}As{sub 0.99}N{sub 0.01} used for the base layer, this device has a low V{sub ON} of 0.79 V, 0.25 V lower than in a comparable Pnp AlGaAs/GaAs HBT. The BV{sub CEO} is 12 V, consistent with its GaAs collector thickness and doping level.

A high voltage GaAs HBT with an open-base collector breakdown voltage of 106 V and an open-emitter breakdown voltage of 134 V has been demonstrated. A high quality 9.0 {micro}m thick collector doped to 2.0{times}10{sup 15} cm{sup {minus}3} grown by MBE on a doped GaAs substrate is the key to achieving this breakdown. These results were achieved for HBTs with 4{times}40 {micro}m{sup 2} emitters. DC current gain of 38 at 6,000 A/cm{sup 2} was measured.

The Reproducing Kernel Particle Method (RKPM) is a discretization technique for partial differential equations that uses the method of weighted residuals, classical reproducing kernel theory and modified kernels to produce either ``mesh-free'' or ``mesh-full'' methods. Although RKPM has many appealing attributes, the method is new, and its numerical performance is just beginning to be quantified. In order to address the numerical performance of RKPM, von Neumann analysis is performed for semi-discretizations of three model one-dimensional PDEs. The von Neumann analyses results are used to examine the global and asymptotic behavior of the semi-discretizations. The model PDEs considered for this analysis include the parabolic and hyperbolic (first and second-order wave) equations. Numerical diffusivity for the former and phase speed for the later are presented over the range of discrete wavenumbers and in an asymptotic sense as the particle spacing tends to zero. Group speed is also presented for the hyperbolic problems. Excellent diffusive and dispersive characteristics are observed when a consistent mass matrix formulation is used with the proper choice of refinement parameter. In contrast, the row-sum lumped mass matrix formulation severely degraded performance. The asymptotic analysis indicates that very good rates of convergence are possible when the consistent mass matrix formulation is used with an appropriate choice of refinement parameter.

The authors compare the results of a microscopic laser theory with gain and recombination currents obtained from experimental spontaneous emission spectra. The calculated absorption spectrum is first matched to that measured on a laser, ensuring that the quasi-Fermi levels for the calculation and the experiment (spontaneous emission and gain) are directly related. This allows one to determine the inhomogeneous broadening in their experimental samples. The only other inputs to the theory are literature values of the bulk material parameter. The authors then estimate the non-radiative recombination current associated with the well and wave-guide core from a comparison of measured and calculated recombination currents.

Surface acoustic wave (SAW) sensors, which are sensitive to a variety of surface changes, have been widely used for chemical and physical sensing. The ability to control or compensate for the many surface forces has been instrumental in collecting valid data. In cases where it is not possible to neglect certain effects, such as frequency drift with temperature, methods such as the dual sensor technique have been utilized. This paper describes a novel use of a dual sensor technique, using two sensor materials, Quartz and GaAs, to separate out the contributions of mass and modulus of the frequency change during gas adsorption experiments. The large modulus change in the film calculated using this technique, and predicted by the Gassmann equation, provide a greater understanding of the challenges of SAW sensing.

The Nuclear Power Engineering Corporation of Japan and the US Nuclear Regulatory Commission Office of Nuclear Regulatory Research, are co-sponsoring and jointly funding a Cooperative Containment Research Program at Sandia National Laboratories, Albuquerque, New Mexico, USA. As a part of this program, a steel containment vessel model and contact structure assembly was tested to failure with over pressurization at Sandia on December 11--12, 1996. The steel containment vessel model was a mixed-scale model (1:10 in geometry and 1:4 in shell thickness) of a steel containment for an improved Mark-II Boiling Water Reactor plant in Japan. The contact structure, which is a thick, bell-shaped steel shell separated at a nominally uniform distance from the model, provides a simplified representation of features of the concrete reactor shield building in the actual plant. The objective of the internal pressurization test was to provide measurement data of the structural response of the model up to its failure in order to validate analytical modeling, to find its pressure capacity, and to observe the failure model and mechanisms.

The executive branch of the United States government has acknowledged and identified threats to the water supply infrastructure of the United States. These threats include contamination of the water supply, aging infrastructure components, and malicious attack. Government recognition of the importance of providing safe, secure, and reliable water supplies has a historical precedence in the water works of the ancient Romans, who recognized the same basic threats to their water supply infrastructure the United States acknowledges today. System surety is the philosophy of ''designing for threats, planning for failure, and managing for success'' in system design and implementation. System surety is an alternative to traditional compliance-based approaches to safety, security, and reliability. Four types of surety are recognized: reactive surety; proactive surety, preventative surety; and fundamental, inherent surety. The five steps of the system surety approach can be used to establish the type of surety needed for the water infrastructure and the methods used to realize a sure water infrastructure. The benefit to the water industry of using the system surety approach to infrastructure design and assessment is a proactive approach to safety, security, and reliability for water transmission, treatment, distribution, and wastewater collection and treatment.

Structural adhesive joints, when tested as made, typically fail cohesively through the centerline of the adhesive. However, in any study of adhesive joint durability, failure near the adhesive/substrate interface becomes an important consideration. In the current study, an interfacially debonding adhesive test, the notched coating adhesion (NCA) test, was applied to LaRC(trademark) PETI-5 adhesive bonded to chronic acid anodized (CAA) Ti-6Al-4V. Post-failure analysis of the interphase region included X-ray photoelectron spectroscopy (XPS), Auger electron spectroscopy (AES), field emission scanning electron microscopy (FE-SEM), and atomic force microscopy (AFM). Mechanical interlocking between an adhesive and a substrate occurs when the liquid adhesive flows into interstices of the substrate, solidifies, and becomes locked in place. Mechanical interlocking is believed to significantly contribute to the adhesion of substrates that exhibit microroughness, such as metal surfaces treated with chromic acid anodization or sodium hydroxide anodization. Filbey and Wightman found that an epoxy penetrated the pores of CAA Ti-6Al-4V, one of the limited number of pore penetration studies that have been reported. In the current study, the penetration of PETI-5 into the pores of CAA Ti-6Al-4V is investigated through analysis of adhesive/substrate failure surfaces.

A study of the errors introduced when one-dimensional inverse heat conduction techniques are applied to problems involving two-dimensional heat transfer effects was performed. The geometry used for the study was a cylinder with similar dimensions as a typical container used for the transportation of radioactive materials. The finite element analysis code MSC P/Thermal was used to generate synthetic test data that was then used as input for an inverse heat conduction code. Four different problems were considered including one with uniform flux around the outer surface of the cylinder and three with non-uniform flux applied over 360{degree}, 180{degree}, and 90{degree} sections of the outer surface of the cylinder. The Sandia One-Dimensional Direct and Inverse Thermal (SODDIT) code was used to estimate the surface heat flux of all four cases. The error analysis was performed by comparing the results from SODDIT and the heat flux calculated based on the temperature results obtained from P/Thermal. Results showed an increase in error of the surface heat flux estimates as the applied heat became more localized. For the uniform case, SODDIT provided heat flux estimates with a maximum error of 0.5% whereas for the non-uniform cases, the maximum errors were found to be about 3%, 7%, and 18% for the 360{degree}, 180{degree}, and 90{degree} cases, respectively.

The recent development of downhole tiltmeter arrays for monitoring hydraulic fractures has provided new information on fracture growth and geometry. These downhole arrays offer the significant advantages of being close to the fracture (large signal) and being unaffected by the free surface. As with surface tiltmeter data, analysis of these measurements requires the inversion of a crack or dislocation model. To supplement the dislocation models of Davis [1983], Okada [1992] and others, this work has extended several elastic crack solutions to provide tilt calculations. The solutions include constant-pressure 2D, penny-shaped, and 3D-elliptic cracks and a 2D-variable-pressure crack. Equations are developed for an arbitrary inclined fracture in an infinite elastic space. Effects of fracture height, fracture length, fracture dip, fracture azimuth, fracture width and monitoring distance on the tilt distribution are given, as well as comparisons with the dislocation model. The results show that the tilt measurements are very sensitive to the fracture dimensions, but also that it is difficult to separate the competing effects of the various parameters.

Efficient (1.2% yield) multikilovolt x-ray emission from Ba(L) (2.4--2.8{angstrom}) and Gd(L) (1.7--2.1{angstrom}) is produced by ultraviolet (248nm) laser-excited BaF{sub 2} and Gd solids. The high efficiency is attributed to an inner shell-selective collisional electron ejection. Much effort has been expended recently in attempts to develop an efficient coherent x-ray source suitable for high-resolution biological imaging. To this end, many experiments have been performed studying the x-ray emissions from high-Z materials under intense (>10{sup 18}W/cm{sup 2}) irradiation, with the most promising results coming from the irradiation of Xe clusters with a UV (248nm) laser at intensities of 10{sup 18}--10{sup 19}W/cm{sup 2}. In this paper the authors report the production of prompt x-rays with energies in excess of 5keV with efficiencies on the order of 1% as a result of intense irradiation of BaF{sub 2} and Gd targets with a terawatt 248nm laser. The efficiency is attributed to an inner shell-selective collisional electron ejection mechanism in which the previously photoionized electrons are ponderomotively driven into an ion while retaining a portion of their atomic phase and symmetry. This partial coherence of the laser-driven electrons has a pronounced effect on the collisional cross-section for the electron ion interaction.

UV filaments in air have been examined on the basis of the diameter and length of the filament, the generation of new spectral components, and the ionization by multiphoton processes. There have been numerous observations of filaments at 800 nm. The general perception is that, above a critical power, the beam focuses because nonlinear self-lensing overcomes diffraction. The self-focusing proceeds until an opposing higher order nonlinearity forms a stable balance.

The past several years rendered a resurgence of interest in phosphors for low-voltage flat panel displays utilizing cathodoluminescence (CL). A major selection criterion for these phosphors is CL efficiency. The objective is to maximize the efficiency at low voltages. This work focuses on understanding the materials properties that influence CL efficiency below 1 kV. Existing high-voltage CL efficiency models take into account intrinsic materials properties such as band-gap energy. Experimental data reveals that the CL efficiency also depends on physical properties such as particle and crystallite size. An update, predictive model of CL efficiency that includes the effects of crystallite size, radiative recombination probability, and electron accelerating potential was developed. The predicted efficiencies agree very well with experimental results. The experimental data were collected using a hot filament electron gun in a demountable high-vacuum chamber. To obtain measurement accuracy, secondary electrons were collected and the phosphor excited with a uniform beam profile. A CL characterization protocol for display phosphors was established at Sandia National Laboratories and made available to phosphor researchers.

Fine particle porous {alpha}-alumina has been prepared by a wet chemical method of combustion synthesis using an aqueous precursor containing aluminum nitrate (oxidizer) and carbohydrazide, an organic fuel as starting materials. The aluminum nitrate and carbohydrazide were reacted exothermically at 400--600 C. The synthesis of {alpha}-alumina ({alpha}-Al{sub 2}O{sub 3}) was used as a model for understanding the effects of processing parameters on physical properties such as surface area, average pore size, and residual carbon content. The porous powders were characterized using x-ray diffraction (XRD), scanning electron microscopy (SEM), BET surface area analysis and elemental analysis. The decomposition of the starting materials was investigated using differential thermal and thermogravimetric analyses (DTA/TGA). It has been shown that the furnace temperature, fuel/oxidizer ratio, and precursor water content can be tailored to produce powders with different physical properties.

This paper describes three-dimensional microstructures fabricated in a planar process and assembled in a single step. Multiple plates are constrained by hinges in such a way as to reduce the assembly process to a single degree-of-freedom of motion. Serial microassembly of these structures is simpler; moreover, self-assembly using hydrodynamic forces during release is much more feasible than with earlier, multiple degree-of-freedom hinged structures. A 250-{micro}m corner cube reflector, a 6-sided closed box, and a 3-D model of the Berkeley Campanile clock tower have been demonstrated in the 4-level polysilicon SUMMiT MEMS foundry.

The authors demonstrate, for the first time, a functional N-p-n heterojunction bipolar transistor using a novel material, InGaAsN, with a bandgap energy of 1.2eV as the p-type base layer. A 300{angstrom}-thick In{sub x}Ga{sub 1-x}As graded layer was introduced to reduce the conduction band offset at the p-type InGaAsN base and n-type GaAs collector junction. For an emitter size of 500 {mu}m{sup 2}, a peak current gain of 5.3 has been achieved.

The authors demonstrated a functional PnP double heterojunction bipolar transistor (DHBT) using AlGaAs, InGaAsN, and GaAs. The band alignment between InGaAsN and GaAs has a large {triangle}E{sub c} and negligible {triangle}E{sub v}, this unique characteristic is very suitable for PnP DHBT applications. The metalorganic vapor phase epitaxy (MOCVD) grown Al{sub 0.3}Ga{sub 0.7}As/In{sub 0.03}Ga{sub 0.97}As{sub 0.99}N{sub 0.01}/GaAs PnP DHBT is lattice matched to GaAs and has a peak current gain of 25. Because of the smaller bandgap (E{sub g}=1.20eV) of In{sub 0.03}Ga{sub 0.97}As{sub 0.99}N{sub 0.01} used for the base layer, this device has a low V{sub ON} of 0.79 V, which is 0.25 V lower than in a comparable Pnp AlGaAs/GaAs HBT. And because GaAs is used for the collector, its BV{sub CEO} is 12 V, consistent with BV{sub CEO} of AlGaAs/GaAs HBTs.

The natural space radiation environment presents a great challenge to present and future satellite systems with significant assets in space. Defining requirements for such systems demands knowledge about the space radiation environment and its effects on electronics and optoelectronics technologies, as well as suitable risk assessment of the uncertainties involved. For mission of high radiation levels, radiation-hardened integrated circuits will be required to preform critical mission functions. The most successful systems in space will be those that are best able to blend standard commercial electronics with custom radiation-hardened electronics in a mix that is suitable for the system of interest.

A new mode of operation for an ion mobility spectrometer (IMS) has been demonstrated that uses potential wells to trap and separate ions by their mobility. This mode of operation has been made feasible by the improvements in personal computers that now allow real-time control of the potentials on ring electrodes in the IMS drift tube. This mode of operation does not require a shutter grid and allows the accumulation of ions in the potential well to enhance the ion signal. Loss of ions from the potential well is controlled by the radial electric fields required by Gauss's law.

In their recent paper, Djikpesse and Tarantola (Geophysics 65 (4) pp. 1023-1035, hereinafter D and T) raise a central question about geophysical inversion: how accurately must the physics of seismic waves in the Earth be modeled in order that inversion succeed? Two general criteria for successful inversion appear in D and T's discussion: fit of predicted to observed data, and prediction of Earth structure. The hypothesis underlying inversion is that these criteria are unextricably linked, so that data fit should lead to accurate inference of subsurface features. The authors have also worked on the data discussed in D and T, using different modeling choices and inversion algorithms but also achieving quite successful inversions, in both senses. They feel that a brief comparison of methods and results might highlight the subtle relation between accuracy in modeling and success in inversion as well as raising questions about the appropriateness of D and T's modeling and inversion choices.

Phosphonate binding sites in guanidine and ammonium surface-functionalized silica xerogels were prepared via the molecular imprinting technique and characterized using solid state {sup 31}P MAS NMR. One-point, two-point, and non-specific host-guest interactions between phenylphosphonic acid (PPA) and the functionalized gels were distinguished by characteristic chemical shifts of the observed absorption peaks. Using solid state as well as solution phase NMR analyses, absorptions observed at 15.5 ppm and 6.5 ppm were identified as resulting from the 1:1 (one-point) and 2:1 (two-point) guanidine to phosphonate interactions, respectively. Similar absorptions were observed with the ammonium functionalized gels. By examining the host-guest interactions within the gels, the efficiency of the molecular imprinting procedure with regard to the functional monomer-to-template interaction could be readily assessed. Template removal followed by substrate adsorption studies conducted on the guanidine functionalized gels provided a method to evaluate the binding characteristics of the receptor sites to a phosphonate substrate. During these experiments, {sup 29}Si and {sup 31}P MAS NMR acted as diagnostic monitors to identify structural changes occurring in the gel matrix and at the receptor site from solvent mediated processes.

The adequacy for laboratory testing of four dolomite cores from the Culebra Dolomite of the Rustler Formation at the Waste Isolation Pilot Plant near Carlsbad, New Mexico, were evaluated using representative elementary volume (REV) theory. Gamma ray computerized tomography created three-dimensional grids of bulk density and macropore index over volumes from 1.4 x 10-7 to 1.6 L. Three different methods for both volume averaging and REV analysis were applied and compared. Both density and macropore index converged to single values with increasing volume, which meets the most common qualitative definition of a REV. Statistical test results for the relatively homogeneous samples indicate that volumes larger than 1 to 7 mL have constant properties. Contrarily, a highly varied sample required 250 and 373 mL to achieve invariant density and macropore characteristics, respectively. Prismatic volume averaging was found to be better than slice averaging, while a qualitative test for the REV provided similar results as a rigorous statistical method. All cores were larger than the REV but were significantly different from one another. This implies that multiple cores are necessary to determine the entire range of transport properties within the rock.

The effective-porosity approach and the dual-porosity approach are examined as two alternative conceptual models of radionuclide migration in fractured media of the saturated zone at Yucca Mountain. Numerical simulations of onedimensional radionuclide transport are performed for the domain relevant to repository performance assessment using the two alternative conceptual approaches. Dual-porosity solute transport modeling produces similar results to effective-porosity modeling for fracture spacing of less than approximately 1 m and greater than about 200 m, which corresponds to values of effective porosity equal to the matrix porosity and the fracture porosity, respectively. For intermediate values of fracture spacing, the dual-porosity approach results in concentration breakthrough curves that differ significantly from the effectiveporosity approach and are characterized by earlier first arrival, greater apparent dispersion, and lower concentrations at later times. The effective-porosity approach, as implemented in recent performance assessment analyses of saturated zone transport at Yucca Mountain, is conservative compared to the dual-porosity approach in terms of both radionuclide concentrations and, generally, travel times.

Injection of high power, multi-mode laser profiles into a fiber optic delivery system requires controlling a number of injection parameters to maximize throughput and minimize concerns for optical damage both at the entrance and exit faces of the fiber optic. A simple method for simultaneously achieving a compact fiber injection geometry and control of these injection parameters, independent of the input source characteristics, is provided by a refractive lenslet array and simple injection lens configuration. Design criteria together with analytical and experimental results for the refractive lenslet array and short focal length injection lens are presented. This arrangement provides a uniform spatial intensity distribution at the fiber injection plane to a large degree independent of the source mode structure, spatial profile, divergence, size, and/or alignment to the injection system. This technique has application to a number of laser systems where uniform illumination of a target or remote delivery of high peak power is desired.

In the present study, 10 impact tests were conducted on unpoled PZT 95/5, with 9% porosity and 2 at% Nb doping. These tests were instrumented to obtain time-resolved loading, unloading and span signatures. As well, PVDF gauges allowed shock timing to be established explicitly. The ferroelectric/antiferroelectric phases transition was manifested as a ramp to 0.4 GPa. The onset of crushup produced the most visible signature: a clear wave separation at 2.2 GPa followed by a highly dispersive wave. The end states also reflected crushup, and are consistent with earlier data and with related poled experiments. A span strength value of 0.17 GPa was measured for a shock stress of 0.5 GPa, this decreased to a very small value (no visible pullback signature) for a shock strength of 1.85 GPa.

SNLO is public domain software developed at Sandia Nat. Labs. It is intended to assist in the selection of the best nonlinear crystal for a particular application, and in predicting its performance. This paper briefly describes its functions and how to use them.

The variation of the value of the linewidth of an excitonic transition in InGaAsN alloys (1% and 2% nitrogen) as a function of hydrostatic pressure using photoluminescence spectroscopy is studied at 4 K. The excitonic linewidth increases as a function of pressure until about 100 kbar after which it tends to saturate. This pressure dependent excitonic linewidth is used to derive the pressure variation of the exciton reduced mass using a theoretical formalism based on the premise that the broadening of the excitonic transition is caused primarily by compositional fluctuations in a completely disordered alloy. The linewidth derived ambient pressure masses are compared and found to be in agreement with other mass measurements. The variation of this derived mass is compared with the results from a nearly first-principles approach in which calculations based on the local density approximation to the Kohn-Sham density functional theory are corrected using a small amount of experimental input.

In this paper, we overview several of the critical materials growth, design and performance issues for nitride-based UV (<400 nm) LEDs. The critical issue of optical efficiency is presented through temperature-dependent photoluminescence studies of various UV active regions. These studies demonstrate enhanced optical efficiencies for active regions with In-containing alloys (InGaN, AlInGaN). We discuss the trade-off between the challenging growth of high Al containing alloys (AlGaN, AlGaInN), and the need for sufficient carrier confinement in UV heterostructures. Carrier leakage for various composition AlGaN barriers is examined through a calculation of the total unconfined carrier density in the quantum well system. We compare the performance of two distinct UV LED structures: GaN/AlGaN quantum well LEDs for λ<360 nm emission, and InGaN/AlGaInN quantum well LEDs for 370 nm<λ<390 nm emission.

The following topics related to radionuclide and colloid transport in the Culebra Dolomite in the 1996 performance assessment for the Waste Isolation Pilot Plant (WIPP) are presented: (i) mathematical description of models: (ii) uncertainty and sensitivity analysis results arising from subjective (i,e. epistemic) uncertainty for individual releases; and (iii) construction of complementary cumulative distribution functions (CCDFs) arising from stochastic (i.e. aleatory) uncertainty. The presented results indicate that radionuclide and colloid transport in the Culebra Dolomite does not constitute a serious threat to the effectiveness of the WIPP as a disposal facility for transuranic waste. Even when the effects of uncertain analysis inputs are taken into account, no radionuclide transport to the boundary with the accessible environment was observed; thus, the associated CCDFs for comparison with the boundary line specified in the US Environmental Protection Agency's standard for the geologic disposal of radioactive waste (40 CFR 191, 40 CFR 194) are degenerate in the sense of having a probability of zero of exceeding a release of zero.

The following topics related to the treatment of direct brine releases to the surface environment in the 1996 performance assessment for the Waste Isolation Pilot Plant (WIPP) are presented: (i) mathematical description of models; (ii) uncertainty and sensitivity analysis results arising from subjective (i.e. epistemic) uncertainty for individual releases: (iii) construction of complementary cumulative distribution functions (CCDFs) arising from stochastic (i.e. aleatory) uncertainty; and (iv) uncertainty and sensitivity analysis results for CCDFs. The presented analyses indicate that direct brine releases do not constitute a serious threat to the effectiveness of the WIPP as a disposal facility for transuranic waste. Even when the effects of uncertain analysis inputs are taken into account, the CCDFs for direct brine releases fall substantially to the left of the boundary line specified in the US Environmental Protection Agency's standard for the geologic disposal of radioactive waste (40 CFR 191, 40 CFR 194).

We have designed and assembled two generations of integrated micro-optical systems that deliver pump light and detect broadband laser-induced fluorescence in micro-fluidic chemical separation systems employing electrochromatography. The goal is to maintain the sensitivity attainable with larger, tabletop machines while decreasing package size and increasing throughput (by decreasing the required chemical volume). One type of micro-optical system uses vertical-cavity surface-emitting lasers (VCSELs) as the excitation source. Light from the VCSELs is relayed with four-level surface relief diffractive optical elements (DOEs) and delivered to the chemical volume through substrate-mode propagation. Indirect fluorescence from dye-quenched chemical species is collected and collimated with a high numerical aperture DOE. A filter blocks the excitation wavelength, and the resulting signal is detected as the chemical separation proceeds. Variations of this original design include changing the combination of reflective and transmissive DOEs and optimizing the high numerical aperture DOE with a rotationally symmetric iterative discrete on-axis algorithm. We will discuss the results of these implemented optimizations.

Two programs have been written in C++ to greatly automate the process of computer simulation visualization inmost cases. These programs, rasterize.C and tracker.C, can be used to generate numerous images in order to create a video or still ties. In order to limit the amount of time and work involved in visualizing simulations, both of these programs have their own specific output formats. The first output format, from rasterize.C, is best suited for those who need only to visualize the actions of a single element, or elements that work on roughly the same time scale. The second format, from tracker.C, is best suited for simulations which involve multiple elements that work on different time scales and thus must be represented in a manner other than straight forward visualization.

The requirement to efficiently form images over a large range of geometries has a profound impact on the design of a Synthetic Aperture Radar (SAR) system. This article shows how a data set conducive to efficient processing might increase the total RF bandwidth. It also presents examples of how a fixed RF bandwidth might then limit SAR geometries.

Because of limited direct observation, understanding of the interior conditions of the massive storage caverns constructed in Gulf Coast salt domes is realizable only through predictions of salt response. Determination of the potential for formation of salt spalls, leading to eventual salt falls, is based on salt creep and fracture using the Multimechanism-Deformation Coupled Fracture (MCDF) model. This is a continuum model for creep, coupled to continuum damage evolution. The model has been successfully tested against underground results of damage around several test rooms at the Waste Isolation Pilot Plant (WIPP). Model simulations, here, evaluate observations made in the Strategic Petroleum Reserve (SPR), storage caverns, namely, the accumulation of material on cavern floors and evidence of salt falls. A simulation of a smooth cavern wall indicates damage is maximum at the surface but diminishes monotonically into the salt, which suggests the source of salt accumulation is surface sluffing. If a protuberance occurs on the wall, fracture damage can form beneath the protuberance, which will eventually cause fracture, and lead to a salt fall.

This paper will describe an algorithm for detecting and classifying seismic and acoustic signals for unattended ground sensors. The algorithm must be computationally efficient and continuously process a data stream in order to establish whether or not a desired signal has changed state (turned-on or off). The paper will focus on describing a Fourier-based technique that compares the running power spectral density estimate of the data to a predetermined signature in order to determine if the desired signal has changed state. How to establish the signature and the detection thresholds will be discussed as well as the theoretical statistics of the algorithm for the Gaussian noise case with results from simulated data. Actual seismic data results will also be discussed along with techniques used to reduce false alarms due to the inherent nonstationary noise environments found with actual data.

We study the low-temperature in-plane magnetoconductance of vertically coupled double quantum wires. Using a novel flip-chip technique, the wires are defined by two pairs of mutually aligned split gates on opposite sides of a≤1 micron thick AlGaAs/GaAs double quantum well heterostructure. We observe quantized conductance steps due to each quantum well and demonstrate independent control of each 1D wire. A broad dip in the magnetoconductance at approximately 6 T is observed when a magnetic field is applied perpendicular to both the current and growth directions. This conductance dip is observed only when 1D subbands are populated in both the top and bottom constrictions. This data is consistent with a counting model whereby the number of subbands crossing the Fermi level changes with field due to the formation of an anticrossing in each pair of 1D subbands.

A simple numerical procedure for predicting damage and permeability in the disturbed rock zone (DRZ) is given. The empirical procedure predicts damage based on a function of stress tensor invariants. For a wttte class of problems hydro logic/mechanical coupling is necessary for proper analysis. The RATDAMPER procedure incorporates dilatant volumetric strain and permeability. It has been implemented in a weakly coupled code, which combines a finite element structural code and a finite difference multi-phase fluid flow code. Using the development of inelastic volumetric strain, a value of permeability can be assigned. This flexibility allows empirical permeability functional relationships to be evaluated.

Designing products for easy assembly and disassembly during their entire life cycles for purposes including product assembly, product upgrade, product servicing and repair, and product disposal is a process that involves many disciplines. In addition, finding the best solution often involves considering the design as a whole and by considering its intended life cycle. Different goals and manufacturing plan selection criteria, as compared to initial assembly, require re-visiting significant fundamental assumptions and methods that underlie current assembly planning techniques. Previous work in this area has been limited to either academic studies of issues in assembly planning or to applied studies of life cycle assembly processes that give no attention to automatic planning. It is believed that merging these two areas will result in a much greater ability to design for, optimize, and analyze the cycle assembly processes. The study of assembly planning is at the very heart of manufacturing research facilities and academic engineering institutions; and, in recent years a number of significant advances in the field of assembly planning have been made. These advances have ranged from the development of automated assembly planning systems, such as Sandia's Automated Assembly Analysis System Archimedes 3.0{copyright}, to the startling revolution in microprocessors and computer-controlled production tools such as computer-aided design (CAD), computer-aided manufacturing (CAM), flexible manufacturing systems (EMS), and computer-integrated manufacturing (CIM). These results have kindled considerable interest in the study of algorithms for life cycle related assembly processes and have blossomed into a field of intense interest. The intent of this manuscript is to bring together the fundamental results in this area, so that the unifying principles and underlying concepts of algorithm design may more easily be implemented in practice.

Automatic assembly sequencing and visualization tools are valuable in determining the best assembly sequences, but without Human Factors and Figure Models (HFFMs) it is difficult to evaluate or visualize human interaction. In industry, accelerating technological advances and shorter market windows have forced companies to turn to an agile manufacturing paradigm. This trend has promoted computerized automation of product design and manufacturing processes, such as automated assembly planning. However, all automated assembly planning software tools assume that the individual components fly into their assembled configuration and generate what appear to be a perfectly valid operations, but in reality the operations cannot physically be carried out by a human. Similarly, human figure modeling algorithms may indicate that assembly operations are not feasible and consequently force design modifications; however, if they had the capability to quickly generate alternative assembly sequences, they might have identified a feasible solution. To solve this problem HFFMs must be integrated with automated assembly planning to allow engineers to verify that assembly operations are possible and to see ways to make the designs even better. Factories will very likely put humans and robots together in cooperative environments to meet the demands for customized products, for purposes including robotic and automated assembly. For robots to work harmoniously within an integrated environment with humans the robots must have cooperative operational skills. For example, in a human only environment, humans may tolerate collisions with one another if they did not cause much pain. This level of tolerance may or may not apply to robot-human environments. Humans expect that robots will be able to operate and navigate in their environments without collisions or interference. The ability to accomplish this is linked to the sensing capabilities available. Current work in the field of cooperative automation has shown the effectiveness of humans and machines directly interacting to perform tasks. To continue to advance this area of robotics, effective means need to be developed to allow natural ways for people to communicate and cooperate with robots just as they do with one another.

Experiments have been performed using a coaxial end-effector to combine a focused laser beam and a plasma arc. The device employs a hollow tungsten electrode, a focusing lens, and conventional plasma arc torch nozzles to co-locate the focused beam and arc on the workpiece. Plasma arc nozzles were selected to protect the electrode from laser generated metal vapor. The project goal is to develop an improved fusion welding process that exhibits both absorption robustness and deep penetration for small scale (<1.5 mm thickness) applications. On aluminum alloys 6061 and 6111, the hybrid process has been shown to eliminate hot cracking in the fusion zone. Fusion zone dimensions for both stainless steel and aluminum were found to be wider than characteristic laser welds, and deeper than characteristic plasma arc welds.

Due to extreme surface to volume ratios, adhesion and friction are critical properties for reliability of Microelectromechanical Systems (MEMS), but are not well understood. In this LDRD the authors established test structures, metrology and numerical modeling to conduct studies on adhesion and friction in MEMS. They then concentrated on measuring the effect of environment on MEMS adhesion. Polycrystalline silicon (polysilicon) is the primary material of interest in MEMS because of its integrated circuit process compatibility, low stress, high strength and conformal deposition nature. A plethora of useful micromachined device concepts have been demonstrated using Sandia National Laboratories' sophisticated in-house capabilities. One drawback to polysilicon is that in air the surface oxidizes, is high energy and is hydrophilic (i.e., it wets easily). This can lead to catastrophic failure because surface forces can cause MEMS parts that are brought into contact to adhere rather than perform their intended function. A fundamental concern is how environmental constituents such as water will affect adhesion energies in MEMS. The authors first demonstrated an accurate method to measure adhesion as reported in Chapter 1. In Chapter 2 through 5, they then studied the effect of water on adhesion depending on the surface condition (hydrophilic or hydrophobic). As described in Chapter 2, they find that adhesion energy of hydrophilic MEMS surfaces is high and increases exponentially with relative humidity (RH). Surface roughness is the controlling mechanism for this relationship. Adhesion can be reduced by several orders of magnitude by silane coupling agents applied via solution processing. They decrease the surface energy and render the surface hydrophobic (i.e. does not wet easily). However, only a molecular monolayer coats the surface. In Chapters 3-5 the authors map out the extent to which the monolayer reduces adhesion versus RH. They find that adhesion is independent of RH up to a threshold value, depending on the coating chemistry. The mechanism for the adhesion increase beyond this threshold value is that the coupling agent reconfigures from a surface to a bulk phase (Chapter 3). To investigate the details of how the adhesion increase occurs, the authors developed the mechanics for adhesion hysteresis measurements. These revealed that near-crack tip compression is the underlying cause of the adhesion increase (Chapter 4). A vacuum deposition chamber for silane coupling agent deposition was constructed. Results indicate that vapor deposited coatings are less susceptible to degradation at high RH (Chapter 5). To address issues relating to surfaces in relative motion, a new test structure to measure friction was developed. In contrast to other surface micromachined friction test structures, uniform apparent pressure is applied in the frictional contact zone (Chapter 6). The test structure will enable friction studies over a large pressure and dynamic range. In this LDRD project, the authors established an infrastructure for MEMS adhesion and friction metrology. They then characterized in detail the performance of hydrophilic and hydrophobic films under humid conditions, and determined mechanisms which limit this performance. These studies contribute to a fundamental understanding for MEMS reliability design rules. They also provide valuable data for MEMS packaging requirements.

Conventional methods of gathering forensic evidence at crime scenes are encumbered by difficulties that limit local law enforcement efforts to apprehend offenders and bring them to justice. Working with a local law-enforcement agency, Sandia National Laboratories has developed a prototype multispectral imaging system that can speed up the investigative search task and provide additional and more accurate evidence. The system, called the Criminalistics Light-imaging Unit (CLU), has demonstrated the capabilities of locating fluorescing evidence at crime scenes under normal lighting conditions and of imaging other types of evidence, such as untreated fingerprints, by direct white-light reflectance. CLU employs state of the art technology that provides for viewing and recording of the entire search process on videotape. This report describes the work performed by Sandia to design, build, evaluate, and commercialize CLU.

Z pinches, the oldest fusion concept, have recently been revisited in light of significant advances in the fields of plasma physics and pulsed power engineering. The possibility exists for z-pinch fusion to play a role in commercial energy applications. We report on work to develop z-pinch fusion concepts, the result of an extensive literature search, and the output for a congressionally-mandated workshop on fusion energy held in Snowmass, Co July 11-23,1999.

This paper compares the cost and effectiveness of several potential options that may be used to monitor silo-based ballistic missiles. Silo door monitoring can be used to verify that warheads removed to deactivate or download silo-based ballistic missiles have not been replaced. A precedent for monitoring warhead replacement using reentry vehicle on site inspections (RV-OSIs) and using satellites has been established by START-I and START-II. However, other monitoring options have the potential to be less expensive and more effective. Three options are the most promising if high verification confidence is desired: random monitoring using door sensors; random monitoring using manned or unmanned aircraft; and continuous remote monitoring using unattended door sensors.

The detailed understanding of the formation and evolution of plasma from rapidly heated metallic wires is a long-standing challenge in the field of plasma physics and in exploding wire engineering. This physical process is made even more complicated if the wire material is composed of a number of individual layers. The authors have successfully developed both optical and x-ray backlighting diagnostics. In particular, the x-ray backlighting technique has demonstrated the capability for quantitative determination of the plasma density over a wide range of densities. This diagnostic capability shows that the process of plasma formation is composed of two separate phases: first, current is passed through a cold wire and the wire is heated ohmically, and, second, the heated wire evolves gases that break down and forms a low-density plasma surrounding the wire.

Production of the molybdenum-99 isotope at the Annular Core Research Reactor requires highly enriched, uranium oxide loaded targets to be irradiated for several days in the high neutron-flux region of the core. This report presents the safety analysis for the irradiation of up to seven Cintichem-type targets in the central region of the core and compares the results to the Annular Core Research Reactor Safety Analysis Report. A 19 target grid configuration is presented that allows one to seven targets to be irradiated, with the remainder of the grid locations filled with aluminum ''void'' targets. Analyses of reactor, neutronic, thermal hydraulics, and heat transfer calculations are presented. Steady-state operation and accident scenarios are analyzed with the conclusion that the reactor can be operated safely with seven targets in the grid, and no additional risk to the public.

The bio-terrorism threat has become the ''poor man's'' nuclear weapon. The ease of manufacture and dissemination has allowed an organization with only rudimentary skills and equipment to pose a significant threat with high consequences. This report will analyze some of the most likely agents that would be used, the ease of manufacture, the ease of dissemination and what characteristics of the public health response that are particularly important to the successful characterization of a high consequence event to prevent excessive causalities.

The design of general-purpose dynamic load-balancing tools for parallel applications is more challenging than the design of static partitioning tools. Both algorithmic and software engineering issues arise. We have addressed many of these issues in the design of the Zoltan dynamic load-balancing library. Zoltan has an object-oriented interface that makes it easy to use and provides separation between the application and the load-balancing algorithms. It contains a suite of dynamic load-balancing algorithms, including both geometric and graph-based algorithms. Its design makes it valuable both as a partitioning tool for a variety of applications and as a research test-bed for new algorithmic development. In this paper, we describe Zoltan's design and demonstrate its use in an unstructured-mesh finite element application.

Two new failure analysis techniques have been developed for backside and front side localization of open and shorted interconnections on ICs. These scanning optical microscopy techniques take advantage of the interactions between IC defects and localized heating using a focused infrared laser ({lambda} = 1,340 nm). Images are produced by monitoring the voltage changes across a constant current supply used to power the IC as the laser beam is scanned across the sample. The methods utilize the Seebeck Effect to localize open interconnections and Thermally-Induced Voltage Alteration (TIVA) to detect shorts. Initial investigations demonstrated the feasibility of TIVA and Seebeck Effect Imaging (SEI). Subsequent improvements have greatly increased the sensitivity of the TIVA/SEI system, reducing the acquisition times by more than 20X and localizing previously unobserved defects. The interaction physics describing the signal generation process and several examples demonstrating the localization of opens and shorts are described. Operational guidelines and limitations are also discussed. The system improvements, non-linear response of IC defects to heating, modeling of laser heating and examples using the improved system for failure analysis are presented.

This report focuses on the relationship between the fundamental interactions acting across an interface and macroscopic engineering observable such as fracture toughness or fracture stress. The work encompasses experiment, theory, and simulation. The model experimental system is epoxy on polished silicon. The interfacial interactions between the substrate and the adhesive are varied continuously using self-assembling monolayer. Fracture is studied in two specimen geometries: a napkin-ring torsion geometry and a double cantilevered beam specimen. Analysis and modeling involves molecular dynamics simulations and continuum mechanics calculations. Further insight is gained from analysis of measurements in the literature of direct force measurements for various fundamental interactions. In the napkin-ring test, the data indicate a nonlinear relationship between interface strength and fracture stress. In particular, there is an abrupt transition in fracture stress which corresponds to an adhesive-to-cohesive transition. Such nonlinearity is not present in the MD simulations on the tens-of-nanometer scale, which suggests that the nonlinearity comes from bulk material deformation occurring on much larger length scales. We postulate that the transition occurs when the interface strength becomes comparable to the yield stress of the material. This postulate is supported by variation observed in the fracture stress curve with test temperature. Detailed modeling of the stress within the sample has not yet been attempted. In the DCB test, the relationship between interface strength and fracture toughness is also nonlinear, but the fracture mechanisms are quite different. The fracture does not transition from adhesive to cohesive, but remains adhesive over the entire range of interface strength. This specimen is modeled quantitatively by combining (i) continuum calculations relating fracture toughness to the stress at 90 {angstrom} from the crack tip, and (ii) a relationship from molecular simulations between fracture stress on a {approx} 90 {angstrom} scale and the fraction of surface sites which chemically bond. The resulting relationship between G{sub c} and fraction of bonding sites is then compared to the experimental data. This first order model captures the nonlinearity in the experimentally-determined relationship. A much more extensive comparison is needed (calculations extending to higher G{sub c} values, experimental data extending to lower G{sub c} values) to guide further model development.

Intelligent, integrated microsystems combine some or all of the functions of sensing, processing information, actuation, and communication within a single integrated package, and preferably upon a single silicon chip. As the elements of these highly integrated solutions interact strongly with each other, the microsystems can be neither designed nor fabricated piecemeal, in contrast to the more familiar assembled products. Driven by technological imperatives, microsystems will best be developed by multi-disciplinary teams, most likely within flatter, less hierarchical organizations. Standardization of design and process tools around a single, dominant technology will expedite economically viable operation under a common production infrastructure. The production base for intelligent, integrated microsystems has elements in common with the mathematical theory of chaos. Similar to chaos theory, the development of microsystems technology will be strongly dependent on, and optimized to, the initial product requirements that will drive standardization-thereby further rewarding early entrants to integrated microsystems technology.

Three reactive materials were evaluated at laboratory scale to identify the optimum treatment reagent for use in a Permeable Reactive Barrier Treatment System at Rocky Flats Environmental Technology Site (RFETS). The contaminants of concern (COCS) are uranium, TCE, PCE, carbon tetrachloride, americium, and vinyl chloride. The three reactive media evaluated included high carbon steel iron filings, an iron-silica alloy in the form of a foam aggregate, and a peculiar humic acid based sorbent (Humasorb from Arctech) mixed with sand. Each material was tested in the laboratory at column scale using simulated site water. All three materials showed promise for the 903 Mound Site however, the iron filings were determined to be the least expensive media. In order to validate the laboratory results, the iron filings were further tested at a pilot scale (field columns) using actual site water. Pilot test results were similar to laboratory results; consequently, the iron filings were chosen for the fill-scale demonstration of the reactive barrier technology. Additional design parameters including saturated hydraulic conductivity, treatment residence time, and head loss across the media were also determined and provided to the design team in support of the final design. The final design was completed by the Corps of Engineers in 1997 and the system was constructed in the summer of 1998. The treatment system began fill operation in December, 1998 and despite a few problems has been operational since. Results to date are consistent with the lab and pilot scale findings, i.e., complete removal of the contaminants of concern (COCs) prior to discharge to meet RFETS cleanup requirements. Furthermore, it is fair to say at this point in time that laboratory developed design parameters for the reactive barrier technology are sufficient for fuel scale design; however,the treatment system longevity and the long-term fate of the contaminants are questions that remain unanswered. This project along with others such as the Durango, CO and Monticello, UT reactive barriers will provide the data to determine the long-term effectiveness and return on investment (ROI) for this technology for comparison to the baseline pump and treat.

The burgeoning new technology of Micro-Electro-Mechanical Systems (MEMS) shows great promise in the weapons arena. We can now conceive of micro-gyros, micro-surety systems, and micro-navigators that are extremely small and inexpensive. Do we want to use this new technology in critical applications such as nuclear weapons? This question drove us to understand the reliability and failure mechanisms of silicon surface-micromachined MEMS. Development of a testing infrastructure was a crucial step to perform reliability experiments on MEMS devices and will be reported here. In addition, reliability test structures have been designed and characterized. Many experiments were performed to investigate failure modes and specifically those in different environments (humidity, temperature, shock, vibration, and storage). A predictive reliability model for wear of rubbing surfaces in microengines was developed. The root causes of failure for operating and non-operating MEMS are discussed. The major failure mechanism for operating MEMS was wear of the polysilicon rubbing surfaces. Reliability design rules for future MEMS devices are established.

Monte Carlo (MC) simulations of phosphate tetrahedron connectivity distributions in alkali and alkaline earth phosphate glasses are reported. By utilizing a discrete bond model, the distribution of next-nearest neighbor connectivities between phosphate polyhedron for random, alternating and clustering bonding scenarios was evaluated as a function of the relative bond energy difference. The simulated distributions are compared to experimentally observed connectivities reported for solid-state two-dimensional (2D) exchange and double-quantum (2Q) nuclear magnetic resonance (NMR) experiments of phosphate glasses. These MC simulations demonstrate that the polyhedron connectivity is best described by a random distribution in lithium phosphate and calcium phosphate glasses.

The authors have developed diode lasers for short pulse duration and high peak pulse power in the 0.01--100.0 m pulsewidth regime. A primary goal of the program was producing up to 10 W while maintaining good far-field beam quality and ease of manufacturability for low cost. High peak power, 17 W, picosecond pulses have been achieved by gain switching of flared geometry waveguide lasers and amplifiers. Such high powers area world record for this type of diode laser. The light emission pattern from diode lasers is of critical importance for sensing systems such as range finding and chemical detection. They have developed a new integrated optical beam transformer producing rib-waveguide diode lasers with a symmetric, low divergence, output beam and increased upper power limits for irreversible facet damage.

Junction field effect transistors (JFET) were fabricated on a GaN epitaxial structure grown by metal organic chemical vapor deposition. The DC and microwave characteristics, as well as the high temperature performance of the devices were studied. These devices exhibited excellent pinch-off and a breakdown voltage that agreed with theoretical predictions. An extrinsic transconductance (gm) of 48 mS/mm was obtained with a maximum drain current (ID) of 270 mA/mm. The microwave measurement showed an fT of 6 GHz and an fmax of 12 GHz. Both the ID and the gm were found to decrease with increasing temperature, possibly due to lower electron mobility at elevated temperatures. These JFETs exhibited a significant current reduction after a high drain bias was applied, which was attributed to a partially depleted channel caused by trapped electrons in the semi-insulating GaN buffer layer.

A variety of organic and hybrid organic-inorganic polymer systems were prepared and evaluated for their bulk response to optical, thermal and chemical environmental changes. These included modeling studies of polyene-bridged metal porphyrin systems, metal-mediated oligomerization of phosphaalkynes as heteroatomic analogues to polyacetylene monomers, investigations of chemically amplified degradation of acid- and base-sensitive polymers and thermally responsive thermoplastic thermosets based on Diels-Alder cycloaddition chemistry. The latter class of materials was utilized to initiate work to develop a new technique for rapidly building a library of systems with varying depolymerization temperatures.

An overview of the use of classical mechanical molecular simulations of porphyrins, hydroporphyrins and heme proteins is given. The topics cover molecular mechanics calculations of structures and conformer energies of porphyrins, energies of barriers for interconversion between stable conformers, molecular dynamics of porphyrins and heme proteins, and normal-coordinate structural analysis of experimental and calculated porphyrin structures. Molecular mechanics and dynamics are currently a fertile area of research on porphyrins. In the future, other computational methods such as Monte Carlo simulations, which have yet to be applied to porphyrins, will come into use and open new avenues of research into molecular simulations of porphyrins. Copyright (C) 2000 John Wiley and Sons, Ltd.

Using an interfacial force microscope, the measured elastic response of 100-nm-thick Au films was found to be strongly correlated with the films' stress state and thermal history. Large, reversible variations (2×) of indentation modulus were recorded as a function of applied stress. Low-temperature annealing caused permanent changes in the films' measured elastic properties. The measured elastic response was also found to vary in close proximity to grain boundaries in thin films and near surface steps on single-crystal surfaces. These results demonstrate a complex interdependence of stress state, defect structure, and elastic properties in thin metallic films.

The LDRD entitled ``Role of Defects in III-Nitride Based Devices'' is aimed to place Sandia National Laboratory at the forefront of the field of GaN materials and devices by establishing a scientific foundation in areas such as material growth, defect characterization/modeling, and processing (metalization and etching) chemistry. In this SAND report the authors summarize their studies such as (1) the MOCVD growth and doping of GaN and AlGaN, (2) the characterization and modeling of hydrogen in GaN, including its bonding, diffusion, and activation behaviors, (3) the calculation of energetic of various defects including planar stacking faults, threading dislocations, and point defects in GaN, and (4) dry etching (plasma etching) of GaN (n- and p-types) and AlGaN. The result of the first AlGaN/GaN heterojunction bipolar transistor is also presented.

QUICKSILVER is a 3-d electromagnetic particle-in-cell simulation code developed and used at Sandia to model relativistic charged particle transport. It models the time-response of electromagnetic fields and low-density-plasmas in a self-consistent manner: the fields push the plasma particles and the plasma current modifies the fields. Through an LDRD project a new parallel version of QUICKSILVER was created to enable large-scale plasma simulations to be run on massively-parallel distributed-memory supercomputers with thousands of processors, such as the Intel Tflops and DEC CPlant machines at Sandia. The new parallel code implements nearly all the features of the original serial QUICKSILVER and can be run on any platform which supports the message-passing interface (MPI) standard as well as on single-processor workstations. This report describes basic strategies useful for parallelizing and load-balancing particle-in-cell codes, outlines the parallel algorithms used in this implementation, and provides a summary of the modifications made to QUICKSILVER. It also highlights a series of benchmark simulations which have been run with the new code that illustrate its performance and parallel efficiency. These calculations have up to a billion grid cells and particles and were run on thousands of processors. This report also serves as a user manual for people wishing to run parallel QUICKSILVER.

A phenomenological model was developed for multicomponent transport of charged species with simultaneous electrochemical reactions in concentrated solutions, and was applied to model processes in a thermal battery cell. A new general framework was formulated and implemented in GOMA (a multidimensional, multiphysics, finite-element computer code developed and being enhanced at Sandia) for modeling multidimensional, multicomponent transport of neutral and charged species in concentrated solutions. The new framework utilizes the Stefan-Maxwell equations that describe multicomponent diffusion of interacting species using composition-insensitive binary diffusion coefficients. The new GOMA capability for modeling multicomponent transport of neutral species was verified and validated using the model problem of ternary gaseous diffusion in a Stefan tube. The new GOMA-based thermal battery computer model was verified using an idealized battery cell in which concentration gradients are absent; the full model was verified by comparing with that of Bernardi and Newman (1987) and validated using limited thermal battery discharge-performance data from the open literature (Dunning 1981) and from Sandia (Guidotti 1996). Moreover, a new Liquid Chemkin Software Package was developed, which allows the user to handle manly aspects of liquid-phase kinetics, thermodynamics, and transport (particularly in terms of computing properties). Lastly, a Lattice-Boltzmann-based capability was developed for modeling pore- or micro-scale phenomena involving convection, diffusion, and simplified chemistry; this capability was demonstrated by modeling phenomena in the cathode region of a thermal battery cell.

This project represented a coordinated LLNL-SNL collaboration to investigate the feasibility of developing radiation-hardened magnetic non-volatile memories using giant magnetoresistance (GMR) materials. The intent of this limited-duration study was to investigate whether giant magnetoresistance (GMR) materials similar to those used for magnetic tunnel junctions (MTJs) were process compatible with functioning CMOS circuits. Sandia's work on this project demonstrated that deposition of GMR materials did not affect the operation nor the radiation hardness of Sandia's rad-hard CMOS technology, nor did the integration of GMR materials and exposure to ionizing radiation affect the magnetic properties of the GMR films. Thus, following deposition of GMR films on rad-hard integrated circuits, both the circuits and the films survived ionizing radiation levels consistent with DOE mission requirements. Furthermore, Sandia developed techniques to pattern deposited GMR films without degrading the completed integrated circuits upon which they were deposited. The present feasibility study demonstrated all the necessary processing elements to allow fabrication of the non-volatile memory elements onto an existing CMOS chip, and even allow the use of embedded (on-chip) non-volatile memories for system-on-a-chip applications, even in demanding radiation environments. However, funding agencies DTRA, AIM, and DARPA did not have any funds available to support the required follow-on technology development projects that would have been required to develop functioning prototype circuits, nor were such funds available from LDRD nor from other DOE program funds.

This paper is an overview of the wafer and reticle positioning system of the Extreme Ultraviolet Lithography (EUVL) Engineering Test Stand (ETS). EUVL represents one of the most promising technologies for supporting the integrated circuit (IC) industry's lithography needs for critical features below 100 nm. EUVL research and development includes development of capabilities for demonstrating key EUV technologies. The ETS is under development at the EUV Virtual National Laboratory, to demonstrate EUV full-field imaging and provide data that supports production-tool development. The stages and their associated metrology operate in a vacuum environment and must meet stringent outgassing specifications. A tight tolerance is placed on the stage tracking performance to minimize image distortion and provide high position repeatability. The wafer must track the reticle with less than ±3 nm of position error and jitter must not exceed 10 nm rms. To meet these performance requirements, magnetically levitated positioning stages utilizing a system of sophisticated control electronics will be used. System modeling and experimentation have contributed to the development of the positioning system and results indicate that desired ETS performance is achievable.

Two-phase characteristic curves are necessary for the simulation of water and vapor flow in porous media. Existing functions such as van Genuchten [1980], Brooks and Corey [1966], and Luckner et al. [1989] have significant limitations in the dry region as the liquid saturation goes to zero. This region, which is important in a number of applications, including liquid and vapor flow and vapor-solid sorption, has been the subject of a number of previous investigations. Most previous studies extended standard capillary pressure curves into the adsorption region to zero water content and required a refitting of the revised curves to the data. In contrast, the present method provides for a simple extension of existing capillary pressure curves without the need to refit the experimental data. Therefore previous curve fits can be used, and the transition between the existing fit and the relationship in the adsorption region is easily calculated. The data-model comparison shows good agreement. This extension is a simple and convenient way to extend existing curves to the dry region.

A significant improvement to the classical least-squares (CLS) multivariate analysis method has been developed. The new method, called prediction-augmented classical least-squares (PACLS), removes the restriction for CLS that all interfering spectral species must be known and their concentrations included during the calibration. We demonstrate that PACLS can correct inadequate CLS models if spectral components left out of the calibration can be identified and if their 'spectral shapes' can be derived and added during a PACLS prediction step. The new PACLS method is demonstrated for a system of dilute aqueous solutions containing urea, creatinine, and NaCl analytes with and without temperature variations. We demonstrate that if CLS calibrations are performed with only a single analyte's concentrations, then there is little, if any, prediction ability. However, if pure-component spectra of analytes left out of the calibration are independently obtained and added during PACLS prediction, then the CLS prediction ability is corrected and predictions become comparable to that of a CLS calibration that contains all analyte concentrations. It is also demonstrated that constant-temperature CLS models can be used to predict variable-temperature data by employing the PACLS method augmented by the spectral shape of a temperature change of the water solvent. In this case, PACLS can also be used to predict sample temperature with a standard error of prediction of 0.07°C even though the calibration data did not contain temperature variations. The PACLS method is also shown to be capable of modeling system drift to maintain a calibration in the presence of spectrometer drift.

Nuclear power is an important and, the authors believe, essential component of a secure nuclear future. Although nuclear fuel cycles create materials that have some potential for use in nuclear weapons, with appropriate fuel cycles, nuclear power could reduce rather than increase real proliferation risk worldwide. Future fuel cycles could be designed to avoid plutonium production, generate minimal amounts of plutonium in proliferation-resistant amounts or configurations, and/or transparently and efficiently consume plutonium already created. Furthermore, a strong and viable US nuclear infrastructure, of which nuclear power is a large element, is essential if the US is to maintain a leadership or even participatory role in defining the global nuclear infrastructure and controlling the proliferation of nuclear weapons. By focusing on new fuel cycles and new reactor technologies, it is possible to advantageously burn and reduce nuclear materials that could be used for nuclear weapons rather than increase and/or dispose of these materials. Thus, the authors suggest that planners for a secure nuclear future use technology to design an ideal future. In this future, nuclear power creates large amounts of virtually atmospherically clean energy while significantly lowering the threat of proliferation through the thoughtful use, physical security, and agreed-upon transparency of nuclear materials. The authors must develop options for policy makers that bring them as close as practical to this ideal. Just as Atoms for Peace became the ideal for the first nuclear century, they see a potential nuclear future that contributes significantly to power for peace and prosperity.

MicroElectroMechanical Systems (MEMS) were subjected to a vibration environment that had a peak acceleration of 120 g and spanned frequencies from 20 to 2000 Hz. The device chosen for this test was a surface-micromachined microengine because it possesses many elements (springs, gears, rubbing surfaces) that may be susceptible to vibration. The microengines were unpowered during the test. We observed 2 vibration-related failures and 3 electrical failures out of 22 microengines tested. Surprisingly, the electrical failures also arose in four microengines in our control group indicating that they were not vibration related. Failure analysis revealed that the electrical failures were due to shorting of stationary comb fingers to the ground plane.

Dispersion of solutes in a variable aperture fracture results from a combination of molecular diffusion and velocity variations in both the plane of the fracture (macrodispersion) and across the fracture aperture (Taylor dispersion). We use a combination of physical experiments and computational simulations to test a theoretical model in which the effective longitudinal dispersion coefficient D(L) is expressed as a sum of the contributions of these three dispersive mechanisms. The combined influence of Taylor dispersion and macrodispersion results in a nonlinear dependence of D(L) on the Peclet number (Pe = V/D(m), where V is the mean solute velocity,is the mean aperture, and D(m) is the molecular diffusion coefficient). Three distinct dispersion regimes become evident: For small Pe (Pe << 1), molecular diffusion dominates resulting in D(L) proportional to Pe0; for intermediate Pe, macrodispersion dominates (D(L) proportional to Pe); and for large Pe, Taylor dispersion dominates (D(L) proportional to Pe2). The Pe range corresponding to these different regimes is controlled by the statistics of the aperture field. In particular, the upper limit of Pe corresponding to the macrodispersion regime increases as the macrodispersivity increases. Physical experiments in an analog, rough-walled fracture confirm the nonlinear Pe dependence of D(L) predicted by the theoretical model. However, the theoretical model underestimates the magnitude of D(L). Computational simulations, using a particle-tracking algorithm that incorporates all three dispersive mechanisms, agree very closely with the theoretical model predictions. The close agreement between the theoretical model and computational simulations is largely because, in both cases, the Reynolds equation describes the flow field in the fracture. The discrepancy between theoretical model predictions and D(L) estimated from the physical experiments appears to be largely, due to deviations from the local cubic law assumed by the Reynolds equation.

A long-standing goal of the equation of state (EOS) community has been the development of a loading capability for direct measurement of material properties along an isentrope. Previous efforts on smooth bore launchers have been somewhat successful, but quite difficult to accurately reproduce, had pressure limitations, or tended to be a series of small shocks as opposed to a smoothly increasing pressure load. A technique has recently been developed on the Sandia National Laboratories Z accelerator which makes use of the high current densities and magnetic fields available to produce nearly isentropic compression of samples that are approximately 1 mm in thickness over approximately 120 ns. Velocity interferometry is used to measure the rear surface motion of these samples. The resulting time resolved velocity profiles from multiple sample thicknesses provide information about mechanical response under isentropic loading conditions and phase transition kinetics. Feasibility experiments have been performed to pressures of approximately 130 kbar in copper and 300 kbar in iron with effects of the α–ε phase change kinetics in iron clearly observed. Work is in progress to achieve 1%–2% accuracy in [formula omitted] space along an isentrope, provide uniaxial strain, and to eliminate magnetic field and current diffusion within the sample of interest. © 2000, American Institute of Physics. All rights reserved.

Structural dynamic systems are often attached to a support structure to simulate proper boundary conditions during testing. In some cases, the support structure is fairly simple and can be modeled by discrete springs and dampers. In other cases, the desired test conditions necessitate the use of a support structure which introduces dynamics of its own. For such cases, a more complex structural dynamic model is required to simulate the response of the full combined system. In this paper, experimental frequency response functions, admittance function modeling concepts, and least squares reductions are used to develop a support structure model including both translational and rotational degrees of freedom at an attachment location. Subsequently, the modes of the support structure are estimated, and a NASTRAN model is created for attachment to the tested system.

Two major problems associated with Si-based MEMS (MicroElectroMechanical Systems) devices are stiction and wear. Surface modifications are needed to reduce both adhesion and friction in micromechanical structures to solve these problems. In this paper, we will present a CVD (Chemical Vapor Deposition) process that selectively coats MEMS devices with tungsten and significantly enhances device durability. Tungsten CVD is used in the integrated-circuit industry, which makes this approach manufacturable. This selective deposition process results in a very conformal coating and can potentially address both stiction and wear problems confronting MEMS processing. The selective deposition of tungsten is accomplished through the silicon reduction of WF6. The self-limiting nature of the process ensures consistent process control. The tungsten is deposited after the removal of the sacrificial oxides to minimize stress and process integration problems. The tungsten coating adheres well and is hard and conducting, which enhances performance for numerous devices. Furthermore, since the deposited tungsten infiltrates under adhered silicon parts and the volume of W deposited is less than the amount of Si consumed, it appears to be possible to release adhered parts that are contacted over small areas such as dimples. The wear resistance of tungsten coated parts has been shown to be significantly improved by microengine test structures.

Useful products generated from interferometric synthetic aperture radar (IFSAR) complex data include height measurement, coherent change detection, and classification. The IFSAR coherence is a spatial measure of complex correlation between two collects, a product of IFSAR signal processing. A tacit assumption in such IFSAR signal processing is that the terrain height is constant across an averaging box used in the process of correlating the two images. This paper presents simulations of IFSAR coherence if two targets with different heights exist in a given correlation cell, a condition in IFSAR collections produced by layover. It also includes airborne IFSAR data confirming the simulation results. The paper concludes by exploring the implications of the results on IFSAR height measurements and classification.