Low residual stress silicon oxynitride thin films are investigated for use as a replacement for silicon dioxide (SiO2) as sacrificial layer in surface micromachined microelectrical-mechanical systems (MEMS). It is observed that the level of residual stress in oxynitrides is a function of the nitrogen content in the film. MEMS film stacks are prepared using both SiO2 and oxynitride sacrificial layers. Wafer bow measurements indicate that wafers processed with oxynitride release layers are significantly flatter. Polycrystalline Si (poly-Si) cantilevers fabricated under the same conditions are observed to be flatter when processed with oxynitride rather than SiO2 sacrificial layers. These results are attributed to the lower post-processing residual stress of oxynitride compared to SiO2.
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 this selective. We deposition process ensures the consistency necessary for process control. The tungsten is deposited after the removal of the sacrificial oxides to minimize stress and process integration problems. Tungsten coating adheres well and is hard and conducting, requirements for device performance. 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 stuck parts that are contacted over small areas such as dimples. The wear resistance of selectively coated W parts has been shown to be significantly improved on microengine test structures.
A new type of surface micromachined ratcheting actuation system has been developed at the Microelectronics Development Laboratory at Sandia National Laboratories. The actuator uses a torsional electrostatic comb drive that is coupled to an external ring gear through a ratcheting scheme. The actuator can be operated with a single square wave, has minimal rubbing surfaces, maximizes comb finger density, and can be used for open-loop position control. The prototypes function as intended with a minimum demonstrated operating voltage of 18V. The equations of motion are developed for the torsional electrostatic comb drive. The resonant frequency, voltage vs. displacement and force delivery characteristics are predicted and compared with the fabricated device's performance.
An increasing number of applications are requiring fiber transmission of high-intensity laser pulses. Our particular interests have led us to examine carefully the fiber transmission of Q-switched pulses from multimode Nd:YAG lasers at their fundamental wavelength. The maximum pulse energy that can be transmitted through a particular fiber is limited by the onset of laser-induced breakdown and damage mechanisms. Laser breakdown at the fiber entrance face is often the first limiting process to be encountered, but other mechanisms can result in catastrophic damage at either fiber face, within the initial `entry' segment of the fiber, and at other internal sites along the fiber path. In the course of our studies we have examined a number of factors that govern the relative importance of different mechanisms, including laser characteristics, the design and alignment of injection optics, fiber end-face preparation, and fiber routing. The present study emphasizes the important criteria for injection optics in high-intensity fiber transmission, and illustrates the opportunities that now exist for innovative designs of optics to meet these criteria. Our consideration of diffractive optics to achieve desired injection criteria began in 1993, and we have evaluated a progression of designs since that time. In the present study, two recent designs for injection optics are compared by testing a sufficient number of fibers with each design to establish statistics for the onset of laser-induced breakdown and damage. In this testing we attempted to hold constant other factors that can influence damage statistics. Both designs performed well, although one was less successful in meeting all injection criteria and consequently showed a susceptibility to a particular damage process.
A transient, multidimensional model has been developed to simulate proton exchange membrane fuel cells. The model accounts simultaneously for electrochemical kinetics, current distribution, hydrodynamics, and multicomponent transport. A single set of conservation equations valid for flow channels, gas-diffusion electrodes, catalyst layers, and the membrane region are developed and numerically solved using a finite-volume-based computational fluid dynamics technique. The numerical model is validated against published experimental data with good agreement. Subsequently, the model is applied to explore hydrogen dilution effects in the anode feed. The predicted polarization curves under hydrogen dilution conditions are in qualitative agreement with recent experiments reported in the literature. The detailed two-dimensional electrochemical and flow/transport simulations further reveal that in the presence of hydrogen dilution in the fuel stream, hydrogen is depleted at the reaction surface, resulting in substantial anode mass transport polarization and hence a lower current density that is limited by hydrogen transport from the fuel stream to the reaction site. Finally, a transient simulation of the cell current density response to a step change in cell voltage is reported.
Correctly identifying the possible alteration products and accurately predicting their occurrence in a repository-relevant environment are the key for source-term calculations in a repository performance assessment. Uraninite in uranium deposits has long been used as a natural analog to spent fuel in a repository because of their chemical and structural similarity. In this paper, a SEM/AEM investigation has been conducted on a partially alterated uraninite sample from a uranium ore deposit of Shinkolobwe of Congo. The mineral formation sequences were identified: uraninite→uranyl hydrates→uranyl silicates→Ca-uranyl silicates or uraninite→uranyl silicates→Ca-uranyl silicates. Reaction-path calculations were conducted for the oxidative dissolution of spent fuel in a representative Yucca Mountain groundwater. The predicted sequence is in general consistent with the SEM observations. The calculations also show that uranium carbonate minerals are unlikely to become major solubility-controlling mineral phases in a Yucca Mountain environment. Some discrepancies between model predictions and field observations are observed. Those discrepancies may result from poorly constrained thermodynamic data for uranyl silicate minerals.
Crystalline phases of pyrochlore (e.g., CaPuTi2O7, CaUTi2O7) have been proposed as a durable ceramic waste form for disposal of high level radioactive wastes including surplus weapons-usable plutonium. In this paper, we use a linear free energy relationship to predict the Gibbs free energies of formation of pyrochlore phases (CaMTi2O7). The Pu-pyrochlore phase is predicted to be stable with respect to PuO2, CaTiO3, and TiO2 at room temperatures. Pu-pyrochlore is expected to be stable in a geologic repository where silica and carbonate components are absent or limited. We suggest that a repository in a salt formation be an ideal environment for disposal of high level, pyrochlore-based ceramic wastes. In such environment, adding CaO as a backfill will make pyrochlore minerals thermodynamically stable and therefore effectively prevents actinide release from these mineral phases.
Interfaces play an important role in determining the effect of electric fields on the mechanism of the formation of spinel by solid-state reaction. The reaction occurs by the movement of phase boundaries but the rate of this movement can be affected by grain boundaries in the reactants or in the reaction product. Only by understanding these relationships will it be possible to engineer their behavior. As a particular example of such a study, MgIn2O4 can be formed by the reaction between single-crystal MgO substrate and a thin film of In2O3 with or without an applied electric field. High-resolution backscattered electron (BSE) imaging and electron backscattered diffraction (EBSD) in a scanning electron microscope (SEM) has been used to obtain complementary chemical and crystallographic information.
A new class of semiconductor lasers that can potentially produce much more short pulse energy is presented. This new laser is not limited in volume or aspect ratio by the depth of a p-n junction and are created from current filaments in semi-insulating GaAs. A current filament semiconductor lasers (CFSL) that have produced 75 nJ of 890 nm radiation in 1.5 ns were tested. A filaments as long as 3.4 cm and several hundred microns in diameter in high gain GaAs photoconductive switches were observed. Their smallest dimension can be more than 100 times the carrier diffusion length in GaAs. The spectral narrowing, lasing thresholds, beam divergence, temporal narrowing and energies which imply lasing for several configurations of CFSL are reported.
Failure analysis (FA) tools have been applied to analyze tungsten coated polysilicon microengines. These devices were stressed under accelerated conditions at ambient temperatures and pressure. Preliminary results illustrating the failure modes of microengines operated under variable humidity and ultra-high drive frequency will also be shown. Analysis of tungsten coated microengines revealed the absence of wear debris in microengines operated under ambient conditions. Plan view imaging of these microengines using scanning electron microscopy (SEM) revealed no accumulation of wear debris on the surface of the gears or ground plane on microengines operated under standard laboratory conditions. Friction bearing surfaces were exposed and analyzed using the focused ion beam (FIB). These cross sections revealed no accumulation of debris along friction bearing surfaces. By using transmission electron microscopy (TEM) in conjunction with electron energy loss spectroscopy (EELS), we were able to identify the thickness, elemental analysis, and crystallographic properties of tungsten coated MEMS devices. Atomic force microscopy was also utilized to analyze the surface roughness of friction bearing surfaces.
One method of providing the mode selectivity necessary to insure single mode operation in a large diameter VCSEL is to independently control the size of the gain region and that of the optical mode. Numerical simulations quantity this approach by predicting lateral mode discrimination for different sized gain apertures. Calculations are experimentally confirmed by the fabrication and testing of 850 nm VCSELs employing hybrid ion implantation/selective oxidation that produce a single-mode output of more than 5 mW.
In this paper we concern ourselves with modified versions of the traditional brachistochrone and tautochrone problems. In the modified version of each problem the constant gravity model is replaced with an attractive inverse square law, consequently we name these the 1/r2 brachistochrone and 1/r2 tautochrone problems. With regard to the 1/r2 brachistochrone problem, we show that the shape of the minimizing curve is formally constructed from an infinite series of elliptic integrals, and we use a numerical optimal control technique to generate the trajectories. The 1/r2 tautochrone problem is solved using fractional calculus techniques and we show that the solution satisfies Lagrange's rule for tautochronous curves.
Electrostatic discharge (ESD) and electrical overstress (EOS) damage of Micro-Electro-Mechanical Systems (MEMS) has been identified as a new failure mode. This failure mode has not been previously recognized or addressed primarily due to the mechanical nature and functionality of these systems, as well as the physical failure signature that resembles stiction. Because many MEMS devices function by electrostatic actuation, the possibility of these devices not only being susceptible to ESD or EOS damage but also having a high probability of suffering catastrophic failure due to ESD or EOS is very real. Results from previous experiments have shown stationary comb fingers adhered to the ground plane on MEMS devices tested in shock, vibration, and benign environments. Using Sandia polysilicon microengines, we have conducted tests to establish and explain the ESD/EOS failure mechanism of MEMS devices. These devices were electronically and optically inspected prior to and after ESD and EOS testing. This paper will address the issues surrounding MEMS susceptibility to ESD and EOS damage as well as describe the experimental method and results found from ESD and EOS testing. The tests were conducted using conventional IC failure analysis and reliability assessment characterization tools. In this paper we will also present a thermal model to accurately depict the heat exchange between an electrostatic comb finger and the ground plane during an ESD event.
In 1993, the Government Performance and Results Act (GPRA, PL 103-62) was enacted. GPRA, which applies to all federal programs, has three components: strategic plans, annual performance plans, and metrics to show how well annual plans are being followed. As part of meeting the GRPA requirement in FY2000, a 14-member external peer review panel (the Garwin Committee) was convened on May 17-19, 2000 to review Sandia National Laboratories' Pulsed Power Programs as a component of the Performance Appraisal Process negotiated with the Department of Energy (DOE). The scope of the review included activities in inertial confinement fission (ICF), weapon physics, development of radiation sources for weapons effects simulation, x-ray radiography, basic research in high energy density physics (HEDP), and pulsed power technology research and development. In his charge to the committee, Jeffrey Quintenz, Director of Pulsed Power Sciences (1600) asked that the review be based on four criteria (1) quality of science, technology, and engineering, (2) programmatic performance, management, and planning, (3) relevance to national needs and agency missions, and (4) performance in the operation and construction of major research facilities. In addition, specific programmatic questions were posed by the director and by the DOE-Defense Programs (DP). The accompanying report, produced as a SAND document, is the report of the committee's findings.
This paper identifies active research topics concerning human machine interfaces for intelligent machine systems. The paper was compiled by performing a series of literature searches and organizing the information according to the author's interest in better directing his own Human Machine Interface (HMI) research. Introductory literature from outside the HMI communities is also referenced to provide context.
This report summarizes the results of work conducted by the Semiconductor Manufacturing Technologies Program at Sandia National Laboratories (Sandia) during 1999. This work was performed by one working group: the Semiconductor Equipment Technology Center (SETEC). The group's projects included Numerical/Experimental Characterization of the Growth of Single-Crystal Calcium Fluoride (CaF{sub 2}); The Use of High-Resolution Transmission Electron Microscopy (HRTEM) Imaging for Certifying Critical-Dimension Reference Materials Fabricated with Silicon Micromachining; Assembly Test Chip for Flip Chip on Board; Plasma Mechanism Validation: Modeling and Experimentation; and Model-Based Reduction of Contamination in Gate-Quality Nitride Reactor. During 1999, all projects focused on meeting customer needs in a timely manner and ensuring that projects were aligned with the goals of the National Technology Roadmap for Semiconductors sponsored by the Semiconductor Industry Association and with Sandia's defense mission. This report also provides a short history of the Sandia/SEMATECH relationship and a brief on all projects completed during the seven years of the program.
Thermally-Induced Voltage Alteration (TIVA) is a relatively new technique for locating electrical defects in integrated circuits [1,2]. This paper describes a novel application of TIVA, to locate design anomalies. A newly designed integrated circuit with high and inconsistent Quiescent Power Supply Current (IDDQ) was initially diagnosed with limited success using various failsite isolation techniques. The TIVA technique was successful in accurately locating design anomalies. Results from TIVA identified a spurious ring oscillator in the design. Design modifications carried out using a focussed ion beam (FIB), verified the accuracy of the results from TIVA. This study clearly extends the use of TIVA beyond that of locating electrical defects and anomalies into the realm of design debugging.
LM111 voltage comparators exhibit a wide range of total-dose-induced degradation. Simulations show this variability may be a natural consequence of the low base doping of the substrate PNP (SPNP) input transistors. Low base doping increases the SPNPs collector to base breakdown voltage, current gain, and densities. The build-up of oxide trapped charge (N OT) and interface traps (N IT) is shown to be a function of pre-irradiation bakes. Experimental data indicate that, despite its structural similarities to the LM111, irradiated input transistors of the LM124 operational amplifier do not exhibit the same sensitivity to variations in pre-irradiation thermal cycles. Further disparities in LM111 and LM124 responses may result from a difference in the oxide defect build-up in the two part types. Variations in processing, packaging, and circuit effects are suggested as potential explanations.
The worst case bias during total dose irradiation of partially depleted SOI transistors from two technologies is correlated to the device architecture. Experiments and simulations are used to analyze SOI back transistor threshold voltage shift and charge trapping in the buried oxide.
This report describes the features of Aspen-EE (Electricity Enhancement), a new model for simulating the interdependent effects of market decisions and disruptions in the electric power system on other critical infrastructures in the US economy. Aspen-EE extends and modifies the capabilities of Aspen, an agent-based model previously developed by Sandia National Laboratories. Aspen-EE was tested on a series of scenarios in which the rules governing electric power trades were changed. Analysis of the scenario results indicates that the power generation company agents will adjust the quantity of power bid into each market as a function of the market rules. Results indicate that when two power markets are faced with identical economic circumstances, the traditionally higher-priced market sees its market clearing price decline, while the traditionally lower-priced market sees a relative increase in market clearing price. These results indicate that Aspen-EE is predicting power market trends that are consistent with expected economic behavior.
A record high fundamental-mode power of 5.1 mW was achieved from coupled-resonator vertical-cavity lasers (CRVCLs). In conventional VCSELs, the extent to which the gain volume may be increased is limited by the onset of multi-mode operation. Results indicate that this limitation is circumvented in a coupled-resonator device allowing high power fundamental-mode operation.
This report discusses an automatic target recognition (ATR) algorithm for synthetic aperture radar (SAR) imagery that is based on coherent change detection techniques. The algorithm relies on templates created from training data to identify targets. Objects are identified or rejected as targets by comparing their SAR signatures with templates using the same complex correlation scheme developed for coherent change detection. Preliminary results are presented in addition to future recommendations.
Current approaches to reliability are not adequate to keep pace with the need for faster, better and cheaper products and systems. This is especially true in high consequence of failure applications. The original proposal for the LDRD was to look at this challenge and see if there was a new paradigm that could make reliability predictions, along with a quantitative estimate of the risk in that prediction, in a way that was faster, better and cheaper. Such an approach would be based on the underlying science models that are the backbone of reliability predictions. The new paradigm would be implemented in two software tools: the Virtual Reliability Realization System (VRRS) and the Reliability Expert System (REX). The three-year LDRD was funded at a reduced level for the first year ($120K vs. $250K) and not renewed. Because of the reduced funding, we concentrated on the initial development of the expertise system. We developed an interactive semiconductor calculation tool needed for reliability analyses. We also were able to generate a basic functional system using Microsoft Siteserver Commerce Edition and Microsoft Sequel Server. The base system has the capability to store Office documents from multiple authors, and has the ability to track and charge for usage. The full outline of the knowledge model has been incorporated as well as examples of various types of content.
>This report describes the innovative modeling approach developed as a result of a 3-year Laboratory Directed Research and Development project. The overall goal of this project was to provide an effective suite of solvers for advanced production planning at facilities in the nuclear weapons complex (NWC). We focused our development activities on problems related to operations at the DOE's Pantex Plant. These types of scheduling problems appear in many contexts other than Pantex--both within the NWC (e.g., Neutron Generators) and in other commercial manufacturing settings. We successfully developed an innovative and effective solution strategy for these types of problems. We have tested this approach on actual data from Pantex, and from Org. 14000 (Neutron Generator production). This report focuses on the mathematical representation of the modeling approach and presents three representative studies using Pantex data. Results associated with the Neutron Generator facility will be published in a subsequent SAND report. The approach to task-based scheduling described here represents a significant addition to the literature for large-scale, realistic scheduling problems in a variety of production settings.
This report describes the design of PICO, a C++ framework for implementing general parallel branch-and-bound algorithms. The PICO framework provides a mechanism for the efficient implementation of a wide range of branch-and-bound methods on an equally wide range of parallel computing platforms. We first discuss the basic architecture of PICO, including the application class hierarchy and the package's serial and parallel layers. We next describe the design of the serial layer, and its central notion of manipulating subproblem states. Then, we discuss the design of the parallel layer, which includes flexible processor clustering and communication rates, various load balancing mechanisms, and a non-preemptive task scheduler running on each processor. We describe the application of the package to a branch-and-bound method for mixed integer programming, along with computational results on the ASCI Red massively parallel computer. Finally we describe the application of the branch-and-bound mixed-integer programming code to a resource constrained project scheduling problem for Pantex.
The Air Force's Electronic Systems Center has funded Sandia National Laboratories to develop an Automatic Target Recognition (ATR) System for the Air Force's Joint STARS platform using Mercury Computer systems hardware. This report provides general theory on the internal operations of the Real-Time ATR system and provides some basic techniques that can be used to reconfigure the system and monitor its runtime operation. In addition, general information on how to interface an image formation processor and a human machine interface to the ATR is provided. This report is not meant to be a tutorial on the ATR algorithms.
A pointing control system is developed and tested for a flying gimbaled telescope. The two-axis pointing system is capable of sub-microradian pointing stability and high accuracy in the presence of large host vehicle jitter. The telescope also has high agility--it is capable of a 50-degree retarget (in both axes simultaneously) in less than 2 seconds. To achieve the design specifications, high-accuracy, high-resolution, two-speed resolvers were used, resulting in gimbal-angle measurements stable to 1.5 microradians. In addition, on-axis inertial angle displacement sensors were mounted on the telescope to provide host-vehicle jitter cancellation. The inertial angle sensors are accurate to about 100 nanoradians, but do not measure low frequency displacements below 2 Hz. The gimbal command signal includes host-vehicle attitude information, which is band-limited. This provides jitter data below 20 Hz, but includes a variable latency between 15 and 25 milliseconds. One of the most challenging aspects of this design was to combine the inertial-angle-sensor data with the less perfect information in the command signal to achieve maximum jitter reduction. The optimum blending of these two signals, along with the feedback compensation were designed using Quantitative Feedback Theory.
The purpose of this research was to develop a science-based understanding of the early-time behavior of electric surface arcing in air at atmospheric pressure. As a first step towards accomplishing this, we used a kinetic approach to model an electron swarm as it evolved in a neutral gas under the influence of an applied electric field. A computer code was written in which pseudo-particles, each representing some number of electrons, were accelerated by an electric field. The electric field due to the charged particles was calculated efficiently using a tree algorithm. Collision of the electrons with the background gas led to the creation of new particles through the processes of ionization and photoionization. These processes were accounted for using measured cross-section data and Monte Carlo methods. A dielectric half-space was modeled by imaging the charges in its surface. Secondary electron emission from the surface, resulting in surface charging, was also calculated. Simulation results show the characteristics of a streamer in three dimensions. A numerical instability was encountered before the streamer matured to form branching.
Prediction of the evolution of microstructures in weapons systems is critical to meeting the objectives of stockpile stewardship in accordance with the Nuclear Weapons Test Ban Treaty. For example, accurate simulation of microstructural evolution in solder joints, cermets, PZT power generators, etc. is necessary for predicting the performance, aging, and reliability both of individual components and of entire weapons systems. A recently developed but promising approach called the ''Phase-Field Model'' (PFM) has the potential of allowing the accurate quantitative prediction of microstructural evolution, with all the spatial and thermodynamic complexity of a real microstructure. Simulating with the PFM requires solving a set of coupled nonlinear differential equations, one for each material variable (e.g., grain orientation, phase, composition, stresses, anisotropy, etc.). While the PFM is versatile and is able to incorporate the necessary complexity for modeling real material systems, it is very computationally intensive, and it has been a difficult and major challenge to formulate an efficient algorithmic implementation of the approach. We found that second order in space algorithm is more stable and leads to more accurate results. However, the computational requirements still remain high, so we have developed a single field algorithm to reduce the computations by 2 orders of magnitude. We have created a 3-D parallel version of the basic phase-field (PF model) and benchmarked it performance. Preliminary results indicate that we will be able to run very large problems effectively with the new parallel code. Microstructural evolution in a diffusion couple was simulated using PFM to simultaneously simulate grain growth, diffusion and phase transformation. Solute drag in a variable composition material, a process no other model can simulate, was successfully simulated using the phase-field model. The phase field model was used to study the evolution of fractal high curvature structures to show that these structures have very different morphological and kinetic behaviors than those of equi-axed structures.
Roll-isolated inertial measurement units are developed at Sandia for use in the instrumentation, guidance, and control of rapidly spinning vehicles. Roll-isolation is accomplished by supporting the inertial instrument cluster (gyros and accelerometers) on a single gimbal, the axis of which is parallel to the vehicle's spin axis. A rotary motor on the gimbal is driven by a servo loop to null the roll gyro output, thus inertially stabilizing the gimbal and instrument cluster while the vehicle spins around it. Roll-isolation prevents saturation of the roll gyro by the high vehicle spin rate, and vastly reduces measurement errors arising from gyro scale factor and alignment uncertainties. Nine versions of Sandia-developed roll-isolated inertial measurement units have been flown on a total of 27 flight tests since 1972.
In order to provide real-time data for validation of three dimensional numerical simulations of heterogeneous materials subjected to impact loading, an optically recording velocity interferometer system (ORVIS) has been adapted to a line-imaging instrument capable of generating precise mesoscopic scale measurements of spatially resolved velocity variations during dynamic deformation. Combining independently variable target magnification and interferometer fringe spacing, this instrument can probe a velocity field along line segments up to 15 mm in length. In high magnification operation, spatial resolution better than 10 {micro}m can be achieved. For events appropriate to short recording times, streak camera recording can provide temporal resolution better than 0.2 ns. A robust method for extracting spatially resolved velocity-time profiles from streak camera image data has been developed and incorporated into a computer program that utilizes a standard VISAR analysis platform. The use of line-imaging ORVIS to obtain measurements of the mesoscopic scale dynamic response of shocked samples has been demonstrated on several different classes of heterogeneous materials. Studies have focused on pressed, granular sugar as a simulant material for the widely used explosive HMX. For low-density (65% theoretical maximum density) pressings of sugar, material response has been investigated as a function of both impact velocity and changes in particle size distribution. The experimental results provide a consistent picture of the dispersive nature of the wave transmitted through these samples and reveal both transverse and longitudinal wave structures on mesoscopic scales. This observed behavior is consistent with the highly structured mesoscopic response predicted by 3-D simulations. Preliminary line-imaging ORVIS measurements on HMX as well as other heterogeneous materials such as foam and glass-reinforced polyester are also discussed.
VxInsight provides a visual mechanism for browsing, exploring and retrieving information from a database. The graphical display conveys information about the relationship between objects in several ways and on multiple scales. In this way, individual objects are always observed within a larger context. For example, consider a database consisting of a set of scientific papers. Imagine that the papers have been organized in a two dimensional geometry so that related papers are located close to each other. Now construct a landscape where the altitude reflects the local density of papers. Papers on physics will form a mountain range, and a different range will stand over the biological papers. In between will be research reports from biophysics and other bridging disciplines. Now, imagine exploring these mountains. If we zoom in closer, the physics mountains will resolve into a set of sub-disciplines. Eventually, by zooming in far enough, the individual papers become visible. By pointing and clicking you can learn more about papers of interest or retrieve their full text. Although physical proximity conveys a great deal of information about the relationship between documents, you can also see which papers reference which others, by drawing lines between the citing and cited papers. For even more information, you can choose to highlight papers by a particular researcher or a particular institution, or show the accumulation of papers through time, watching some disciplines explode and other stagnate. VxInsight is a general purpose tool, which enables this kind of interaction with wide variety of relational data: documents, patents, web pages, and financial transactions are just a few examples. The tool allows users to interactively browse, explore and retrieve information from the database in an intuitive way.
Proceedings of the 9th United Engineering Foundation Conference on Risk-Based Decisionmaking in Water Resources - Risk-Based Decisionmaking in Water Resources IX
Data mining involves the discovery and fusion of features from large databases to establish minimal probability of error (MPE) decision and estimation models. Our approach combines a weighted nearest neighbor (WNN) decision model for classification and estimation with genetic algorithms (GA) for feature discovery and model optimization. The WNN model is used to provide a mathematical framework for adaptively discovering and fusing features into near-MPE decision algorithms. The GA is used to discover weighted features and select decision points for the WNN decision model to achieve near-MPE decisions. The performance of the WNN fusion model is demonstrated on the first of two very different problems to demonstrate its robust and practical application to a wide variety of data-mining problems. The first problem involves the isolation of factors that cause hepatitis C virus (HCV) and requires the evaluation of large databases to establish the critical features that can detect with minimal error whether a person is at risk of having HCV. This requires discovering and extracting relevant information (features) from a questionnaire database and combining (fusing) them to achieve a minimal error decision rule. The primary objective of the research is to develop a practical basis for fusing information from questionnaires administered at hospitals to identify and verify features important to isolate risk factors for HCV. The basic problem involves creating a feature database from the questionnaire information, discovering features that provide sufficient information to reliably identify when a person is at risk under conditions with uncertainties caused by recording errors and evasive tactics of people answering the questionnaire. The results of this study demonstrate the WNN fusion algorithm ability to perform in supervised learning environments. The second phase of the research project is directed at the unsupervised learning environment. In this environment the feature data is presented without any classification. Clustering algorithms are developed to partition the feature data into clusters based upon similarity measure models. After the feature data is clustered and classified the supervised WNN fusion algorithms are used to classify the data based upon the minimal probability of error decision rule.
This report provides a summary of the LDRD project titled: An Electromagnetic Imaging System for Environmental Site Reconnaissance. The major initial challenge of this LDRD was to develop a ground penetrating radar (GPR) whose peak and average radiated power surpassed that of any other in existence. Goals were set to use such a system to detect the following: (1) disrupted soil layers where there is potential for buried waste, (2) buried objects such as 55-gallon drums at depths up to 3 m, and (3) detecting contaminated soil. Initial modeling of the problem suggested that for soil conditions similar to Puerto Rican clay loam, moisture content 10 percent (conductivity = 0.01 mhos at 350 MHz), a buried 55-gallon drum could be detected in a straightforward manner by an UWB GPR system at a depth of 3 meters. From the simulations, the highest attenuation ({minus}50 dB) was the result of scattering from a 3-m deep vertically orientated drum. A system loss of {minus}100 dB is a typical limit for all kinds of radar systems (either direct time-domain or swept frequency). The modeling work also determined that the waveshape of the pulse scattered off the buried drum would be relatively insensitive to drum orientation, and thus easier to detect with the GPR system.
Over the past several years, extensive databases have been developed for the S-N behavior of various mate&& used in wind turbine blades, primarily fiberglass composites. These data are typically presented both in their "raw" form and curve fit to define their average properties. For design, confidence limits must be placed on these descriptions. In particular, most designs call for the "95195" design values; namely, with a 95 percent level of confidence, thedesiguerisassuredthat95percentofthematerial will 'meet or exceed the design value. For such material properties as the ultimate streng& the procedures for estimating its value at a particular confidence level is wellffiedifthemeasured values follow a normal or a log-normal distribution. Namely, based upon the number of sample points and their standard deviation, a commonly-found table may be used to determine the survival percentagea t a particular confidencel evel with respect to its mean value. The same is true for fatigue data at a constaut stress level (the number of cycles to failure N at stress level SI). However, when the stress level is allowed to vary, as with a typical S-N fatigue curve, the proceduresf or determmingc onfidencel imits are not as well delked. This paper outlines techn.iques for determimng confklence limits of fatigue data Different approachesto estimating the 95195l evel are compared. Data from the MSUIDOE and the FACT fatigue databam are used to illustrate typical results.
A MEMS test structure capable of measuring friction between polysilicon surfaces under a variety of test conditions has been refined from previous designs. The device is applied here to measuring friction coefficients of polysilicon surfaces under different environmental, loading, and surface conditions. Two methods for qualitatively comparing friction coefficients (μ) using the device are presented. Samples that have been coated with a self-assembled monolayer of the lubricating film perfluorinated-decyltrichlorosilane (PFTS) have a coefficient of friction that is approximately one-half that of samples dried using super-critical CO2 (SCCO2) drying. Qualitative results indicate that μ is independent of normal pressure. Wear is shown to increase μ for both supercritically dried samples and PFTS coated samples, though the mechanisms appear to be different. Super critically dried surfaces appear to degrade continuously with increased wear cycles, while PFTS coated samples reach a steady state friction value after about 105 cycles.
Hydrogen atom-hydrogen atom scattering is a prototype for many of the fundamental principles of atomic collisions. In this work we present the formalism and the predictions of a time-dependent self-consistent-field description of the H + H system for scattering in the intermediate energy regime of 1-100 keV. Because of the unrestricted nature of the numerical orbital description, this method includes the effects of an unlimited basis set within each orbital. Electron exchange and a limited amount of electron correlation are included as well. We solve numerically coupled three-dimensional Schrodinger equations for the two-electron orbitals in singlet and triplet symmetries. Excitation and ionization cross sections are computed and compared with other theory and experiment. The results capture many features of the problem but illustrate a need for more quantitative experimental information concerning the H + H system in this energy range.
The cookoff of energetic materials involves the combined effects of several physical and chemical processes. These processes include heat transfer, chemical decomposition, and mechanical response. The interaction and coupling between these processes influence both the time-to-event and the violence of reaction. The prediction of the behavior of explosives during cookoff, particularly with respect to reaction violence, is a challenging task. To this end, a joint DoD/DOE program has been initiated to develop models for cookoff, and to perform experiments to validate those models. In this paper, a series of cookoff analyses are presented and compared with data from a number of experiments for the aluminized, RDX-based, Navy explosive PBXN-109. The traditional thermal-chemical analysis is used to calculate time-to-event and characterize the heat transfer and boundary conditions. A reaction mechanism based on Tarver and McGuire's work on RDX{sup 2} was adjusted to match the spherical one-dimensional time-to-explosion data. The predicted time-to-event using this reaction mechanism compares favorably with the validation tests. Coupled thermal-chemical-mechanical analysis is used to calculate the mechanical response of the confinement and the energetic material state prior to ignition. The predicted state of the material includes the temperature, stress-field, porosity, and extent of reaction. There is little experimental data for comparison to these calculations. The hoop strain in the confining steel tube gives an estimation of the radial stress in the explosive. The inferred pressure from the measured hoop strain and calculated radial stress agree qualitatively. However, validation of the mechanical response model and the chemical reaction mechanism requires more data. A post-ignition burn dynamics model was applied to calculate the confinement dynamics. The burn dynamics calculations suffer from a lack of characterization of the confinement for the flaw-dominated failure mode experienced in the tests. High-pressure burning rates are needed for more detailed post-ignition studies. Sub-models for chemistry, mechanical response and burn dynamics need to be validated against data from less complex experiments. The sub-models can then be used in integrated analysis for comparison with experimental data taken during integrated tests.
Since the discovery of surfactant-templated silica by Mobil scientists in 1992, mesostructured silica has been synthesized in various forms including thin films, powders, particles, and fibers. In general, mesostructured silica has potential applications, such as in separation, catalysis, sensors, and fluidic microsystems. In respect to these potential applications, mesostructured silica in the form of thin films is perhaps one of the most promising candidates. The preparation of mesostructured silica films through preferential solvent evaporation-induced self-assembly (EISA) has recently received much attention in the laboratories. However, no amphiphile/silica films with reverse mesophases have ever been made through this EISA procedure. Furthermore, templates employed to date have been either surfactants or poly(ethylene oxide)-b-poly(propylene oxide)-b-poly(ethylene oxide) triblock copolymers, such as pluronic P-123, both of which are water-soluble and alcohol-soluble. Due to their relatively low molecular weight, the templated silica films with mesoscopic order have been limited to relatively small characteristic length scales. In the present communication, the authors report a novel synthetic method to prepare mesostructured amphiphilic/silica films with regular and reverse mesophases of large characteristic length scales. This method involves evaporation-induced self-assembly (EISA) of amphiphilic polystyrene-block-poly(ethylene oxide) (PS-b-PEO) diblock copolymers. In the present study, the PS-b-PEO diblocks are denoted as, for example, PS(215)-b-PEO(100), showing that this particular sample contains 215 S repeat units and 100 EO repeat units. This PS(215)-b-PEO(100) diblock possesses high molecular weight and does not directly mix with water or alcohol. To the authors knowledge, no studies have reported the use of water-insoluble and alcohol-insoluble amphiphilic diblocks as structure-directing agents in the synthesis of mesostructured silica films through EISA. It is believed that the present system is the first to yield amphiphile/silica films with regular and reverse mesophases, as well as curved multi-bilayer mesostructures, through EISA. The ready formation of the diblock/silica films with multi-bilayer vesicular mesostructures is discussed.
Heat pipes are often proposed as cooling system components for small fission reactors. SAFE-300 and STAR-C are two reactor concepts that use heat pipes as an integral part of the cooling system. Heat pipes have been used in reactors to cool components within radiation tests (Deverall, 1973); however, no reactor has been built or tested that uses heat pipes solely as the primary cooling system. Heat pipe cooled reactors will likely require the development of a test reactor to determine the main differences in operational behavior from forced cooled reactors. The purpose of this paper is to describe the results of a systems code capable of modeling the coupling between the reactor kinetics and heat pipe controlled heat transport. Heat transport in heat pipe reactors is complex and highly system dependent. Nevertheless, in general terms it relies on heat flowing from the fuel pins through the heat pipe, to the heat exchanger, and then ultimately into the power conversion system and heat sink. A system model is described that is capable of modeling coupled reactor kinetics phenomena, heat transfer dynamics within the fuel pins, and the transient behavior of heat pipes (including the melting of the working fluid). The paper focuses primarily on the coupling effects caused by reactor feedback and compares the observations with forced cooled reactors. A number of reactor startup transients have been modeled, and issues such as power peaking, and power-to-flow mismatches, and loading transients were examined, including the possibility of heat flow from the heat exchanger back into the reactor. This system model is envisioned as a tool to be used for screening various heat pipe cooled reactor concepts, for designing and developing test facility requirements, for use in safety evaluations, and for developing test criteria for in-pile and out-of-pile test facilities.
An integrated approach, combining the continuum theory of sintering and Potts model based mesostructure evolution analysis, is used to solve the problem of bi-layered structure sintering. Two types of bi-layered structures are considered: layers of the same material with different initial porosity, and layers of two different materials. The effective sintering stress for the bi-layer powder sintering is derived, both at the meso- and the macroscopic levels. Macroscopic shape distortions and spatial distributions of porosity are determined as functions of the dimensionless specific time of sintering. The effect of the thickness of the layers on shrinkage, warpage, and pore-grain structure is studied. Ceramic ZnO powders are employed as a model experimental system to assess the model predictions.
This paper presents an analysis technique used to generate the structural response at locations not measured during the ejection of a captive-carried store. The ejection shock event is complicated by the fact that forces may be imparted to the store at eight distinct locations. The technique derives forcing functions by combining the initial field test data for a limited number of measurement locations with Frequency Response Functions (FRFs) measured using a traditional modal-type impact (tap) test at the same locations. The derived forcing functions were then used with tap test FRFs measured at additional locations of interest to produce the desired response data.
Operating in a potentially hostile city is every soldier's nightmare. The staggering complexity of the urban environment means that deadly threats--or non-combatants-may lurk behind every corner, doorway, or window. Urban operations present an almost unparalleled challenge to the modern professional military. The complexity of urban operations is further amplified by the diversity of missions that the military will be called upon to conduct in urban terrain. Peace-making and peace-keeping missions, urban raids to seize airports or WMD sites or to rescue hostages, and extended urban combat operations all present different sorts of challenges for planners and troops on the ground. Technology almost never serves as a magic bullet, and past predictions of technological miracles pile high on the ash heap of history. At the same time, it is a vital element of planning in the modern age to consider and, if possible, take advantage of emerging technologies. We believe that technologies can assist military operations in urbanized terrain (MOUT) in three primary areas, which are discussed.
This paper summarizes an evaluation of a treatment loop designed to upgrade the quality of spent rinse waters discharged from 10 wet benches located in the fab at Sandia's Microelectronics Development Laboratory (MDL). The goal of the treatment loop is to make these waters, presently being discharged to the fab's acid waste neutralization (AWN) station, suitable for recycling as feed water back into the fab's ultrapure water (UPW) plant. The MDL typically operates 2 shifts per day, 5 days per week. Without any treatment, the properties of the spent rinse waters now being collected have been shown to be compatible with recycling about 30% (50/168) of the time (weekends primarily, when the fab is idling) which corresponds to about 12% of the present water discharged from the fab to the AWN. The primary goal of adding a treatment loop is to increase the percentage of recyclable water from these 10 wet benches to near 100%, increasing the percentage of total recyclable water to near 40% of the total present fab discharge to the AWN. A second goal is to demonstrate compatibility with recycling this treated spent rinse water to the present R/O product water tank, reducing both the present volume of R/O reject water and the present load on the R/O. The approach taken to demonstrate achieving these goals is to compare all the common metrics of water quality for the treated spent rinse waters with those of the present R/O product water. Showing that the treated rinse water is equal or superior in quality to the water presently stored in the R/O tank by every metric all the time is assumed to be sufficient argument for proceeding with plans to incorporate recycling of these spent rinse waters back into MDL's R/O tank.
The use of plasma spray to deposit thin metal-sulfide cathode films is described in this paper. Conventional electroactive stack components in thermal batteries are constructed from pressed-powder parts that are difficult to fabricate in large diameters in thicknesses <0.010. Plasma-sprayed electrodes do not steer from this difficulty, allowing greater energy densities and specific energies to be realized. Various co-spraying agents have been found suitable for improving the mechanical as well as electrochemical properties of plasma-sprayed cathodes for thermal batteries. These electrodes generally show equal or improved performance over conventional pressed-powder electrodes. A number of areas for future growth and development of plasma-spray technology is discussed.
When embarking on an experimental program for purposes of discovery and understanding, it is only prudent to use appropriate analysis tools to aid in the discovery process. Due to the limited scope of experimental measurement analytical results can significantly complement the data after a reasonable validation process has occurred. In this manner the analytical results can help to explain certain measurements, suggest other measurements to take and point to possible modifications to the experimental apparatus. For these reasons it was decided to create a detailed nonlinear finite element model of the Sandia Microslip Experiment. This experiment was designed to investigate energy dissipation due to microslip in bolted joints and to identify the critical parameters involved. In an attempt to limit the microslip to a single interface a complicated system of rollers and cables was devised to clamp the two slipping members together with a prescribed normal load without using a bolt. An oscillatory tangential load is supplied via a shaker. The finite element model includes the clamping device in addition to the sequence of steps taken in setting up the experiment. The interface is modeled using Coulomb friction requiring a modest validation procedure for estimating the coefficient of friction. Analysis results have indicated misalignment problems in the experimental procedure, identified transducer locations for more accurate measurements, predicted complex interface motions including the potential for galling, identified regions where microslip occurs and during which parts of the loading cycle it occurs, all this in addition to the energy dissipated per cycle. A number of these predictions have been experimentally corroborated in varying degrees and are presented in the paper along with the details of the finite element model.
In this paper the authors investigate the relationship between glassy and ferromagnetic phases in disordered Ising ferromagnets in the presence of transverse magnetic fields, {Lambda}. Iterative mean field simulations probe the free energy landscape and suggest the existence of a glass transition line in the {Lambda}, temperature T plane well within the ferromagnetic phase. New experimental field-cooled and zero-field-cooled data on LiHo{sub x} Y{sub 1{minus}x}F{sub 4} provide support for our theoretical picture.
As a joint is loaded, the tangent stiffness of the joint reduces due to slip at interfaces. This stiffness reduction continues until the direction of the applied load is reversed or the total interface slips. Total interface slippage in joints is called macro-slip. For joints not undergoing macro-slip, when load reversal occurs the tangent stiffness immediately rebounds to its maximum value. This occurs due to stiction effects at the interface. Thus, for periodic loads, a softening and rebound hardening cycle is produced which defines a hysteretic, energy absorbing trajectory. For many jointed sub-structures, this hysteretic trajectory can be approximated using simple polynomial representations. This allows for complex joint substructures to be represented using simple non-linear models. In this paper a simple one dimensional model is discussed.
Unattended monitoring systems are being studied as a means of reducing both the cost and intrusiveness of present nuclear safeguards approaches. Such systems present the classic information overload problem to anyone trying to interpret the resulting data not only because of the sheer quantity of data but also because of the problems inherent in trying to correlate information from more than one source. As a consequence, analysis efforts to date have mostly concentrated on checking thresholds or diagnosing failures. Clearly more sophisticated analysis techniques are required to enable automated verification of expected activities level concepts in order to make automated judgments about safety, sensor system integrity, sensor data quality, diversion, and accountancy.
Multiple techniques have been developed to model the temporal evolution of infectious diseases. Some of these techniques have also been adapted to model the spatial evolution of the disease. This report examines the application of one such technique, the SEIR model, to the spatial and temporal evolution of disease. Applications of the SEIR model are reviewed briefly and an adaptation to the traditional SEIR model is presented. This adaptation allows for modeling the spatial evolution of the disease stages at the individual level. The transmission of the disease between individuals is modeled explicitly through the use of exposure likelihood functions rather than the global transmission rate applied to populations in the traditional implementation of the SEIR model. These adaptations allow for the consideration of spatially variable (heterogeneous) susceptibility and immunity within the population. The adaptations also allow for modeling both contagious and non-contagious diseases. The results of a number of numerical experiments to explore the effect of model parameters on the spread of an example disease are presented.
This Safety Analysis Report (SAR) is prepared in compliance with the requirements of DOE Order 5480.23, Nuclear Safety Analysis Reports, and has been written to the format and content guide of DOE-STD-3009-94 Preparation Guide for U. S. Department of Energy Nonreactor Nuclear Safety Analysis Reports. The Hot Cell Facility is a Hazard Category 2 nonreactor nuclear facility, and is operated by Sandia National Laboratories for the Department of Energy. This SAR provides a description of the HCF and its operations, an assessment of the hazards and potential accidents which may occur in the facility. The potential consequences and likelihood of these accidents are analyzed and described. Using the process and criteria described in DOE-STD-3009-94, safety-related structures, systems and components are identified, and the important safety functions of each SSC are described. Additionally, information which describes the safety management programs at SNL are described in ancillary chapters of the SAR.
The magnetic implosion of a high-Z quasi-spherical shell filled with DT fuel by the 20-MA Z accelerator can heat the fuel to near-ignition temperature. The attainable implosion velocity on Z, 13-cm/{micro}s, is fast enough that thermal losses from the fuel to the shell are small. The high-Z shell traps radiation losses from the fuel, and the fuel reaches a high enough density to reabsorb the trapped radiation. The implosion is then nearly adiabatic. In this case the temperature of the fuel increases as the square of the convergence. The initial temperature of the fuel is set by the heating of an ion acoustic wave to be about 200-eV after a convergence of 4. To reach the ignition temperature of 5-keV an additional convergence of 5 is required. The implosion dynamics of the quasi-spherical implosion is modeled with the 2-D radiation hydrodynamic code LASNEX. LASNEX shows an 8-mm diameter quasi-spherical tungsten shell on Z driving 6-atmospheres of DT fuel nearly to ignition at 3.5-keV with a convergence of 20. The convergence is limited by mass flow along the surface of the quasi-spherical shell. With a convergence of 20 the final spot size is 400-{micro}m in diameter.
A theory is developed for the response of moderately porous solids (no more than {approximately}20% void space) to high-strain-rate deformations. The model is consistent because each feature is incorporated in a manner that is mathematically compatible with the other features. Unlike simple p-{alpha} models, the onset of pore collapse depends on the amount of shear present. The user-specifiable yield function depends on pressure, effective shear stress, and porosity. The elastic part of the strain rate is linearly related to the stress rate, with nonlinear corrections from changes in the elastic moduli due to pore collapse. Plastically incompressible flow of the matrix material allows pore collapse and an associated macroscopic plastic volume change. The plastic strain rate due to pore collapse/growth is taken normal to the yield surface. If phase transformation and/or pore nucleation are simultaneously occurring, the inelastic strain rate will be non-normal to the yield surface. To permit hardening, the yield stress of matrix material is treated as an internal state variable. Changes in porosity and matrix yield stress naturally cause the yield surface to evolve. The stress, porosity, and all other state variables vary in a consistent manner so that the stress remains on the yield surface throughout any quasistatic interval of plastic deformation. Dynamic loading allows the stress to exceed the yield surface via an overstress ordinary differential equation that is solved in closed form for better numerical accuracy. The part of the stress rate that causes no plastic work (i.e-, the part that has a zero inner product with the stress deviator and the identity tensor) is given by the projection of the elastic stressrate orthogonal to the span of the stress deviator and the identity tensor.The model, which has been numerically implemented in MIG format, has been exercised under a wide array of extremal loading and unloading paths. As will be discussed in a companion sequel report, the CKP model is capable of closely matching plate impact measurements for porous materials.
Uncertainty distributions for specific parameters of the Cassini General Purpose Heat Source Radioisotope Thermoelectric Generator (GPHS-RTG) Final Safety Analysis Report consequence risk analysis were revised and updated. The revisions and updates were done for all consequence parameters for which relevant information exists from the joint project on Probabilistic Accident Consequence Uncertainty Analysis by the United States Nuclear Regulatory Commission and the Commission of European Communities.
Object Role Modeling (ORM) techniques produce a detailed domain model from the perspective of the business owner/customer. The typical process begins with a set of simple sentences reflecting facts about the business. The output of the process is a single model representing primarily the persistent information needs of the business. This type of model contains little, if any reference to a targeted computerized implementation. It is a model of business entities not of software classes. Through well-defined procedures, an ORM model can be transformed into a high quality objector relational schema.
This report presents research on public key, digital signature algorithms for cryptographic authentication in low-powered, low-computation environments. We assessed algorithms for suitability based on their signature size, and computation and storage requirements. We evaluated a variety of general purpose and special purpose computing platforms to address issues such as memory, voltage requirements, and special functionality for low-powered applications. In addition, we examined custom design platforms. We found that a custom design offers the most flexibility and can be optimized for specific algorithms. Furthermore, the entire platform can exist on a single Application Specific Integrated Circuit (ASIC) or can be integrated with commercially available components to produce the desired computing platform.
Polymer stresses around sharp corners and in constrained geometries of encapsulated components can generate cracks leading to system failures. Often, analysts use maximum stresses as a qualitative indicator for evaluating the strength of encapsulated component designs. Although this approach has been useful for making relative comparisons screening prospective design changes, it has not been tied quantitatively to failure. Accurate failure models are needed for analyses to predict whether encapsulated components meet life cycle requirements. With Sandia's recently developed nonlinear viscoelastic polymer models, it has been possible to examine more accurately the local stress-strain distributions in zones of likely failure initiation looking for physically based failure mechanisms and continuum metrics that correlate with the cohesive failure event. This study has identified significant differences between rubbery and glassy failure mechanisms that suggest reasonable alternatives for cohesive failure criteria and metrics. Rubbery failure seems best characterized by the mechanisms of finite extensibility and appears to correlate with maximum strain predictions. Glassy failure, however, seems driven by cavitation and correlates with the maximum hydrostatic tension. Using these metrics, two three-point bending geometries were tested and analyzed under variable loading rates, different temperatures and comparable mesh resolution (i.e., accuracy) to make quantitative failure predictions. The resulting predictions and observations agreed well suggesting the need for additional research. In a separate, additional study, the asymptotically singular stress state found at the tip of a rigid, square inclusion embedded within a thin, linear elastic disk was determined for uniform cooling. The singular stress field is characterized by a single stress intensity factor K{sub a} and the applicable K{sub a} calibration relationship has been determined for both fully bonded and unbended inclusions. A lack of interfacial bonding has a profound effect on inclusion-tip stress fields. A large radial compressive stress is generated in front of the inclusion-tip when the inclusion is well bonded, whereas a large tensile hoop stress is generated when the inclusion is unbended, and frictionless sliding is allowed. Consequently, an epoxy disk containing an unbended inclusion appears more likely to crack when cooled than a disk containing a fully bonded inclusion. A limited number of tests have been carried out to determine if encapsulant cracking can be induced by cooling a specimen fabricated by molding a square, steel insert within a thin, epoxy disk. Test results are in qualitative agreement with analysis. Cracks developed only in disks with mold-released inserts, and the tendency for cracking increased with inclusion size.
This paper looks at emerging technologies for converging voice and data networks and telephony transport over a data network using Internet Protocols. Considered are the benefits and drivers for this convergence. The paper describes these new technologies, how they are being used, and their application to Sandia.
The authors have successfully demonstrated an optical data interconnection from the output of a focal plane array to the downstream data acquisition electronics. The demonstrated approach included a continuous wave laser beam directed at a multiple quantum well reflectance modulator connected to the focal plane array analog output. The output waveform from the optical interconnect was observed on an oscilloscope to be a replica of the input signal. They fed the output of the optical data link to the same data acquisition system used to characterize focal plane array performance. Measurements of the signal to noise ratio at the input and output of the optical interconnection showed that the signal to noise ratio was reduced by a factor of 10 or more. Analysis of the noise and link gain showed that the primary contributors to the additional noise were laser intensity noise and photodetector receiver noise. Subsequent efforts should be able to reduce these noise sources considerably and should result in substantially improved signal to noise performance. They also observed significant photocurrent generation in the reflectance modulator that imposes a current load on the focal plane array output amplifier. This current loading is an issue with the demonstrated approach because it tends to negate the power saving feature of the reflectance modulator interconnection concept.
The HMX {beta}-{delta} solid-solid phase transition, which occurs as HMX is heated near 170 C, is linked to increased reactivity and sensitivity to initiation. Thermally damaged energetic materials (EMs) containing HMX therefore may present a safety concern. Information about the phase transition is vital to predictive safety models for HMX and HMX-containing EMs. We report work on monitoring the phase transition with real-time Raman spectroscopy aimed towards obtaining a better understanding of physical properties of HMX through the phase transition. HMX samples were confined in a cell of minimal free volume in a displacement-controlled or load-controlled arrangement. The cell was heated and then cooled at controlled rates while real-time Raman spectroscopic measurements were performed. Raman spectroscopy provides a clear distinction between the phases of HMX because the vibrational transitions of the molecule change with conformational changes associated with the phase transition. Temperature of phase transition versus load data are presented for both the heating and cooling cycles in the load-controlled apparatus, and general trends are discussed. A weak dependence of the temperature of phase transition on load was discovered during the heating cycle, with higher loads causing the phase transition to occur at a higher temperature. This was especially true in the temperature of completion of phase transition data as opposed to the temperature of onset of phase transition data. A stronger dependence on load was observed in the cooling cycle, with higher loads causing the reverse phase transitions to occur at a higher cooling temperature. Also, higher loads tended to cause the phase transition to occur over a longer period of time in the heating cycle and over a shorter period of time in the cooling cycle. All three of the pure HMX phases ({alpha}, {beta} and {delta}) were detected on cooling of the heated samples, either in pure form or as a mixture.
The goal of this LDRD Research project was to provide a preliminary examination of the use of infrared spectroscopy as a tool to detect the changes in cell cultures upon activation by an infectious agent. Due to a late arrival of funding, only 5 months were available to transfer and setup equipment at UTTM,develop cell culture lines, test methods of in-situ activation and collect kinetic data from activated cells. Using attenuated total reflectance (ATR) as a sampling method, live cell cultures were examined prior to and after activation. Spectroscopic data were collected from cells immediately after activation in situ and, in many cases for five successive hours. Additional data were collected from cells activated within a test tube (pre-activated), in both transmission mode as well as in ATR mode. Changes in the infrared data were apparent in the transmission data collected from the pre-activated cells as well in some of the pre-activated ATR data. Changes in the in-situ activated spectral data were only occasionally present due to (1) the limited time cells were studied and (2) incomplete activation. Comparison of preliminary data to infrared bands reported in the literature suggests the primary changes seen are due an increase in ribonucleic acid (RNA) production. This work will be continued as part of a 3 year DARPA grant.
The Rapid Terrain Visualization Advanced Concept Technology Demonstration (RTV-ACTD) is designed to demonstrate the technologies and infrastructure to meet the Army requirement for rapid generation of digital topographic data to support emerging crisis or contingencies. The primary sensor for this mission is an interferometric synthetic aperture radar (IFSAR) designed at Sandia National Laboratories. This paper will outline the design of the system and its performance, and show some recent flight test results. The RTV IFSAR will meet DTED level III and IV specifications by using a multiple-baseline design and high-accuracy differential and carrier-phase GPS navigation. It includes innovative near-real-time DEM production on-board the aircraft. The system is being flown on a deHavilland DHC-7 Army aircraft.
The Factsheets web application was conceived out of the requirement to create, update, publish, and maintain a web site with dynamic research and development (R and D) content. Before creating the site, a requirements discovery process was done in order to accurately capture the purpose and functionality of the site. One of the high priority requirements for the site would be that no specialized training in web page authoring would be necessary. All functions of uploading, creation, and editing of factsheets needed to be accomplished by entering data directly into web form screens generated by the application. Another important requirement of the site was to allow for access to the factsheet web pages and data via the internal Sandia Restricted Network and Sandia Open Network based on the status of the input data. Important to the owners of the web site would be to allow the published factsheets to be accessible to all personnel within the department whether or not the sheets had completed the formal Review and Approval (R and A) process. Once the factsheets had gone through the formal review and approval process, they could then be published both internally and externally based on their individual publication status. An extended requirement and feature of the site would be to provide a keyword search capability to search through the factsheets. Also, since the site currently resides on both the internal and external networks, it would need to be registered with the Sandia search engines in order to allow access to the content of the site by the search engines. To date, all of the above requirements and features have been created and implemented in the Factsheet web application. These have been accomplished by the use of flat text databases, which are discussed in greater detail later in this paper.
In a total-system performance assessment (TSPA), uncertainty in the performance measure (e.g., radiation dose) is estimated by first estimating the uncertain y in the input variables and then propagating that uncertain y through the model system by means of Monte Carlo simulation. This paper discusses uncertainty in surface infiltration, which is one of the input variables needed for performance assessments of the Yucca Mountain site. Infiltration has been represented in recent TSPA simulations by using three discrete infiltration maps (i.e., spatial distributions of infiltration) for each climate state in the calculation of unsaturated-zone flow and transport. A detailed uncertainty analysis of infiltration was carried out for two purposes: to better quantify the possible range of infiltration, and to determine what probability weights should be assigned to the three infiltration cases in a TSPA simulation. The remainder of this paper presents the approach and methodology for the uncertainty analysis, along with a discussion of the results.
In this paper we have shown how the direction cosine method of stripmap-mode IFSAR maybe modified for use in the spotlight-mode case. Spotlight-mode IFSAR geometry dictates a common aperture phase center, velocity vector, and baseline vector for every pixel in an image. Angle with respect to the velocity vector is the same for every pixel in a given column and can be computed from the column index, the Doppler of the motion compensation point and the Doppler column sample spacing used in image formation. With these modifications, the direction cosines and length of the line of sight vector to every scatterer in the scene may be computed directly from the raw radar measurements of range, Doppler, and interferometric phase.
One of the key elements of the Stochastic Finite Element Method, namely the polynomial chaos expansion, has been utilized in a nonlinear shock and vibration application. As a result, the computed response was expressed as a random process, which is an approximation to the true solution process, and can be thought of as a generalization to solutions given as statistics only. This approximation to the response process was then used to derive an analytically-based design specification for component shock response that guarantees a balanced level of marginal reliability. Hence, this analytically-based reference SRS might lead to an improvement over the somewhat ad hoc test-based reference in the sense that it will not exhibit regions of conservativeness. nor lead to overtesting of the design.
Experienced experimentalists have gone through the process of attempting to identify a final set of modal parameters from several different sets of extracted parameters. Usually, this is done by visually examining the mode shapes. With the advent of automated modal parameter extraction algorithms such as SMAC (Synthesize Modes and Correlate), very accurate extractions can be made to high frequencies. However, this process may generate several hundred modes that then must be consolidated into a final set of modal information. This has motivated the authors to generate a set of tools to speed the process of consolidating modal parameters by mathematical (instead of visual) means. These tools help quickly identify the best modal parameter extraction associated with several extractions of the same mode. The tools also indicate how many different modes have been extracted in a nominal frequency range and from which references. The mathematics are presented to achieve the best modal extraction of multiple modes at the same nominal frequency. Improvements in the SMAC graphical user interface and database are discussed that speed and improve the entire extraction process.
The solution-mediated synthesis and single crystal structure of (CN{sub 3}H{sub 6}){sub 2} {center_dot} Zn(HPO{sub 3}){sub 2} are reported. This phase is built up from a three-dimensional framework of vertex-linked ZnO{sub 4} and HPO{sub 3} building units encapsulating the extra-framework guanidinium cations. The structure is stabilized by template-to-framework hydrogen bonding. The inorganic framework shows a surprising similarity to those of some known zinc phosphates. Crystal data: (CN{sub 3}H{sub 6}){sub 2} {center_dot} Zn(HPO{sub 3}){sub 2}, AI,= 345.50, orthorhombic, space group Fdd2 (No. 43), a = 15.2109 (6) {angstrom}, b = 11.7281 (5) {angstrom}, c = 14.1821 (6) {angstrom}, V = 2530.0 (4){angstrom}{sup 3}, Z = 8, T = 298 (2)K, R(F) = 0.020, wR(F) = 0.025.
The thermal-hydrologic (TH) and coupled process models describe the evolution of a potential geologic repository as heat is released from emplaced waste. The evolution (thermal, hydrologic, chemical, and mechanical) of the engineered barrier and geologic systems is heavily dependent on the heat released by the waste packages and how the heat is transferred from the emplaced wastes through the drifts and through the repository host rock. The essential elements of this process are extracted (or abstracted) from the process-level models that incorporate the basic energy and mass conservation principles and applied to the total system models used to describe the overall performance of the potential repository. The process of total system performance assessment (TSPA) abstraction is the following. First is a description of the parameter inputs used in the process-level models. A brief description is given hereof past inputs for the viability assessment (e.g., for TSPA-VA) and current inputs for the site recommendation (TSPA-SR). This is followed by a highlight of the process-level models from which the abstractions are made. These include descriptions of TH, thermal-hydrologic-chemical (THC), and thermal-mechanical (TM) processes used to describe the performance of individual waste packages and waste emplacement drifts as well as the repository as a whole. Next is a description of what (and how) information is abstracted from the process-level models. This also includes an accounting of the features, events, and processes (FEPs) that are important to both the regulators and the international repository community in general. Finally, an identification of the TSPA model components that utilize the abstracted information to characterize the overall performance of a potential geologic repository is given.
In this paper an optimization-based method of drift prevention is presented for learning control of underdetermined linear and weakly nonlinear time-varying dynamic systems. By defining a fictitious cost function and the associated model-based sub-optimality conditions, a new set of equations results, whose solution is unique, thus preventing large drifts from the initial input. Moreover, in the limiting case where the modeling error approaches zero, the input that the proposed method converges to is the unique feasible (zero error) input that minimizes the fictitious cost function, in the linear case, and locally minimizes it in the (weakly) nonlinear case. Otherwise, under mild restrictions on the modeling error, the method converges to a feasible sub-optimal input.
The abstraction model used for seepage into emplacement drifts in recent TSPA simulations has been presented. This model contributes to the calculation of the quantity of water that might contact waste if it is emplaced at Yucca Mountain. Other important components of that calculation not discussed here include models for climate, infiltration, unsaturated-zone flow, and thermohydrology; drip-shield and waste-package degradation; and flow around and through the drip shield and waste package. The seepage abstraction model is stochastic because predictions of seepage are necessarily quite uncertain. The model provides uncertainty distributions for seepage fraction fraction of waste-package locations flow rate as functions of percolation flux. In addition, effects of intermediate-scale flow with seepage and seep channeling are included by means of a flow-focusing factor, which is also represented by an uncertainty distribution.
A numerical flow model is developed to simulate two-dimensional fluid flow past immersed, elastically supported tube arrays. This work is motivated by the objective of predicting forces and motion associated with both deep-water drilling and production risers in the oil industry. This work has other engineering applications including simulation of flow past tubular heat exchangers or submarine-towed sensor arrays and the flow about parachute ribbons. In the present work, a vortex method is used for solving the unsteady flow field. This method demonstrates inherent advantages over more conventional grid-based computational fluid dynamics. The vortex method is non-iterative, does not require artificial viscosity for stability, displays minimal numerical diffusion, can easily treat moving boundaries, and allows a greatly reduced computational domain since vorticity occupies only a small fraction of the fluid volume. A gridless approach is used in the flow sufficiently distant from surfaces. A Lagrangian remap scheme is used near surfaces to calculate diffusion and convection of vorticity. A fast multipole technique is utilized for efficient calculation of velocity from the vorticity field. The ability of the method to correctly predict lift and drag forces on simple stationary geometries over a broad range of Reynolds numbers is presented.
Societal needs related to demographics, resources, and human behavior will drive technological advances over the next 20 years. Nanotechnology is anticipated to be an important enabler of these advances, and thus maybe anticipated to have significant influence on new systems approaches to solving societal problems as well as on extending current science and technology-based applications. To examine the potential implications of nanotechnology a societal needs-driven approach is taken. Thus the methodology is to present the definition of the problem, and then examine system concepts, technology issues, and promising future directions. We approach the problem definition from a national and global security perspective and identify three key areas involving the condition of the planet, the human condition, and global security. In anticipating societal issues in the context of revolutionary technologies, such as maybe enabled by nanoscience, the importance of working on the entire life cycle of any technological solution is stressed.
The moduli used in RSA (see [5]) can be generated by many different sources. The generator of that modulus (assuming a single entity generates the modulus) knows its factorization. They would have the ability to forge signatures or break any system based on this moduli. If a moduli and the RSA parameters associated with it were generated by a reputable source, the system would have higher value than if the parameters were generated by an unknown entity. So for tracking, security, confidence and financial reasons it would be beneficial to know who the generator of the RSA modulus was. This is where digital marking comes in. An RSA modulus ia digitally marked, or digitally trade marked, if the generator and other identifying features of the modulus (such as its intended user, the version number, etc.) can be identified and possibly verified by the modulus itself. The basic concept of digitally marking an RSA modulus would be to fix the upper bits of the modulus to this tag. Thus anyone who sees the public modulus can tell who generated the modulus and who the generator believes the intended user/owner of the modulus is.
The free energy barriers of vapor tube formed in a metastable liquid confined between hydrophobic walls is investigated. Monte Carlo simulations, the transition state theory and constrained umbrella sampling techniques are performed to estimate the free energy barrier for vapor tube formation. Transmission coefficients calculated for the liquid layer show that capillary evaporation are also described from the size of a vapor pocket formed between the walls.
The capillary evaporation (cavitation) of water confined between two hydrophobic surfaces in close proximity is analyzed. The water is replaced by vapor due to the presence of bulk energetics and surface energetics. Monte Carlo simulations are performed to determine the effect of water confinement on the dynamics of surface-induced phase transitions. To relate the simulation rates to the experimental data, the mass-conserving Kawasaki algorithms are also performed.
Laser technology has advanced dramatically and is an integral part of today's healthcare delivery system. Lasers are used in the laboratory analysis of human blood samples and serve as surgical tools that kill, burn or cut tissue. Recent semiconductor microtechnology has reduced the size o f a laser to the size of a biological cell or even a virus particle. By integrating these ultra small lasers with biological systems, it is possible to create micro-electrical mechanical systems that may revolutionize health care delivery.
The rating and modeling of photovoltaic PW module performance has been of concern to manufacturers and system designers for over 20 years. Both the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) have developed methodologies to predict module and array performance under actual operating conditions. This paper compares the two methods of determining the performance of PV modules, The methods translate module performance to actual or reference conditions using slightly different approaches. The accuracy of both methods is compared for both hourly, daily, and annual energy production over a year of data recorded at NREL in Golden, CO. The comparison of the two methods will be presented for five different PV module technologies.
The Million Solar Roofs Initiative has motivated a renewed interest in the development of utility interconnected photovoltaic (UIPV) inverters. Government-sponsored programs (PVMaT, PVBONUS) and competition among utility interconnected inverter manufacturers have stimulated innovations and improved the performance of existing technologies. With this resurgence, Sandia National Laboratories (SNL) has developed a program to assist industry initiatives to overcome barriers to UIPV inverters. In accordance with newly adopted IEEE 929-2000, the utility interconnected PV inverters are required to cease energizing the utility grid when either a significant disturbance occurs or the utility experiences an interruption in service. Compliance with IEEE 929-2000 is being widely adopted by utilities as a minimum requirement for utility interconnection. This report summarizes work done at the SNL balance-of-systems laboratory to support the development of IEEE 929-2000 and to assist manufacturers in meeting its requirements.
For the last ten years, the Japanese High-Level Nuclear Waste (HLW) repository program has focused on assessing the feasibility of a basic repository concept, which resulted in the recently published H12 Report. As Japan enters the implementation phase, a new organization must identify, screen and choose potential repository sites. Thus, a rapid mechanism for determining the likelihood of site suitability is critical. The threshold approach, described here, is a simple mechanism for defining the likelihood that a site is suitable given estimates of several critical parameters. We rely on the results of a companion paper, which described a probabilistic performance assessment simulation of the HLW reference case in the H12 report. The most critical two or three input parameters are plotted against each other and treated as spatial variables. Geostatistics is used to interpret the spatial correlation, which in turn is used to simulate multiple realizations of the parameter value maps. By combining an array of realizations, we can look at the probability that a given site, as represented by estimates of this combination of parameters, would be good host for a repository site.
With the exponential growth in the power of computing hardware and software, modeling and simulation is becoming a key enabler for the rapid design of reliable Microsystems. One vision of the future microsystem design process would include the following primary software capabilities: (1) The development of 3D part design, through standard CAD packages, with automatic design rule checks that guarantee the manufacturability and performance of the microsystem. (2) Automatic mesh generation, for 3D parts as manufactured, that permits computational simulation of the process steps, and the performance and reliability analysis for the final microsystem. (3) Computer generated 2D layouts for process steps that utilize detailed process models to generate the layout and process parameter recipe required to achieve the desired 3D part. (4) Science-based computational tools that can simulate the process physics, and the coupled thermal, fluid, structural, solid mechanics, electromagnetic and material response governing the performance and reliability of the microsystem. (5) Visualization software that permits the rapid visualization of 3D parts including cross-sectional maps, performance and reliability analysis results, and process simulation results. In addition to these desired software capabilities, a desired computing infrastructure would include massively parallel computers that enable rapid high-fidelity analysis, coupled with networked compute servers that permit computing at a distance. We now discuss the individual computational components that are required to achieve this vision. There are three primary areas of focus: design capabilities, science-based capabilities and computing infrastructure. Within each of these areas, there are several key capability requirements.
Collective impact ionization has been used to explain lock-on in semi-insulating GaAs under high-voltage bias. We have used this theory to study some of the steady-state properties of lock-on current filaments. In steady state, the heat gained from the field is exactly compensated by the cooling due to phonon scattering. In the simplest approximation, the carrier distribution approaches a quasi-equilibrium Maxwell-Boltzmann distribution. In this report, we examine the validity of this approximation. We find that this approximation leads to a filament carrier density that is much lower than the high density needed to achieve a quasi-equilibrium distribution. Further work on this subject is in progress.
The nature of surfaces and the way they interact with each other during sliding contact can have a direct bearing on the performance of a microelectromechanical (MEMS) device. Therefore, a study was undertaken to characterize the surfaces of LIGA fabricated Ni and Cu components. Sidewall and planar surfaces were examined by scanning electron microscopy (SEM) and atomic force microscopy (AFM). Surface roughness was quantified using the AFM. Post-processing (e.g. lapping, removal of polymer film) can profoundly influence the morphology of LIGA components. Edge rounding and smearing of ductile materials during lapping can result in undesirable sidewall morphologies. By judicious selection of AFM scan sizes, the native roughness ({approximately}10 nm RMS) can be distinguished from that arising due to post processing, e.g. scratches, debris, polymer films. While certain processing effects on morphology such as those due to lapping or release etch can be controlled, the true side wall morphology appears to be governed by the morphology of the polymer mold or by the electroforming process itself, and may be much less amenable to modification.
Recent work at SNL has demonstrated unique capabilities to experimentally measure a variety of ion and neutral particle parameters inside surface features being etched, including ion energy, angular distributions, ion and neutral species measurements. This report details the construction of one recent laboratory tool designed to measure ion beam uniformity over the wafer surface in a reactive ion beam etch system, (RIBE). This information is critical to the development of accurate plasma processing computer models and simulation methods, and is essential for reducing the cost of introducing new processing technologies.
This paper describes a solution for the problem of authenticating the shapes of statistically variant gamma spectra while simultaneously concealing the shapes and magnitudes of the sensitive spectra. The shape of a spectrum is given by the relative magnitudes and positions of the individual spectral elements. Class-specific linear orthonormal transformations of the measured spectra are used to produce output that meet both the authentication and concealment requirements. For purposes of concealment, the n-dimensional gamma spectra are transformed into n-dimensional output spectra that are effectively indistinguishable from Gaussian white noise (independent of the class). In addition, the proposed transformations are such that statistical authentication metrics computed on the transformed spectra are identical to those computed on the original spectra.
Understanding high pressure behavior of homogeneous as well as heterogeneous materials is necessary in order to address the physical processes associated with hypervelocity impact events related to space science applications including orbital debris impact and impact lethality. At very high impact velocities, material properties will be subjugated to phase-changes, such as melting and vaporization. These phase states cannot be obtained through conventional gun technology. These processes need to be represented accurately in hydrodynamic codes to allow credible computational analysis of impact events resulting from hypervelocity impact. In this paper, techniques that are being developed and implemented to obtain the needed shock loading parameters (Hugoniot states) for material characterization studies, namely shock velocity and particle velocity, will be described at impact velocities up to 11 km/s. What is new in this report is that these techniques are being implemented for use at engagement velocities never before attained utilizing two-stage light-gas gun technology.
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 non-turbulent fluctuation growth based on that proposed by Warren and Hassan (Journal of Aircraft, Vol. 35, No. 5) is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittence function based on the work of Dhawan and Narasimha (Journal of Fluid Mechanics, Vol. 3, No. 4). 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 spatial accuracy of selected predictions is analyzed.