Publications

We report a fully integrated high-Q factor micro-ring resonator using silicon nitride/dioxide on a silicon wafer. The micro-ring resonator is critically coupled to a low loss straight waveguide. An intrinsic quality factor of 2.4 x 10{sup 5} has been measured.

Abstract not provided.

A series of amorphous silicate materials with the general formula Na{sub x+2y}M{sub x}{sup 3+}Si{sub 1-x}O{sub 2+y}(M{sup 3+} = Al, Mn, Fe, Y) were studied. Samples were synthesized by a precipitation reaction at room temperature. The results indicate that the ion-exchange capacity (IEC) decreases as follows: Al > Fe > Mn > Y. Additionally, the IEC increases with increasing aluminum concentration. Structural studies show that the relative amount of octahedrally coordinated aluminum increases with increasing Al content, as does the total amount of AlO{sub 4} species increases. The data suggest that the IEC value of these amorphous aluminosilicates is dependent on the tetrahedrally coordinated aluminum. Regeneration of the Al-silicate with acetic acid does not decrease the IEC significantly.

Organic light-emitting diodes (OLEDs), with few exvceptions, are fabricated in the standard way of sequentially depositing active layers and elecrodes onto a substrate. The conventional devices have 'a detrimental layer' at the interface between the organic and the top metal electrode because evaporation results in metal in-diffusion and chemical disruption at the metal-organic interface, Here, a different approach is introduced to construct OLEDs: soft contact lamination (SCL) is based on thysical lamination of thin metal electrodes supported by an elastomeric layer against the electrolumnescent organic layer. Thei method produces spatially homogeneous, intimate contacts via van der Waals interaction between the metal and the organic, resulting in no chemical and physical damages to the organic. Devices fabricated by SCL are shown to have no detrimental layer and fewer luminescence-quenching channels than conventional devices that have evaporated top metal electrodes.

This report describes Phase III of a project entitled Innovative Applications of Energy Storage in a Restructured Electricity Marketplace. For this study, the authors assumed that it is feasible to operate an energy storage plant simultaneously for two primary applications: (1) energy arbitrage, i.e., buy-low-sell-high, and (2) to reduce peak loads in utility ''hot spots'' such that the utility can defer their need to upgrade transmission and distribution (T&D) equipment. The benefits from the arbitrage plus T&D deferral applications were estimated for five cases based on the specific requirements of two large utilities operating in the Eastern U.S. A number of parameters were estimated for the storage plant ratings required to serve the combined application: power output (capacity) and energy discharge duration (energy storage). In addition to estimating the various financial expenditures and the value of electricity that could be realized in the marketplace, technical characteristics required for grid-connected distributed energy storage used for capacity deferral were also explored.

This work is the result of a Sandia National Laboratories LDRD funded fellowship at the University of Michigan. Although, guidance and suggestions were offered by Sandia, the work contained here is primarily the work of Brian H. Stark, and his advisor, Professor Khalil Najafi. Junseok Chae, Andrew Kuo, and their coworkers at the University of Michigan helped to record some of the data. The following is an abstract of their work. We have developed a vacuum packaging technology using a thick nickel film to seal MEMS structures at the wafer level. The package is fabricated in a three-mask process by electroplating a 40 micro-meter thick nickel film over an 8 micro-meter sacrificial photoresist that is removed prior to package sealing. Implementation of electrical feedthroughs in this process requires no planarization. The large release channel enables an 800x800 micro-meter package to be released in less than three hours. Several mechanisms, based upon localized melting and lead/tin solder bumping, for sealing the release channel have been investigated. We have also developed Pirani gauges, integrated with this package, which can be used to establish the hermeticity of the different sealing technologies. They have measured a sealing pressure of approximately 1.5 Torr. Our work differs from previous Pirani gauges in that we utilize a novel doubly anchored structure that stiffens the structural membrane while not substantially degrading performance in order to measure fine leak rates.

We present experimental observation and theoretical analysis of looping carbon nanotubes connecting two electrodes. The measured conductance of the nanotubes is not strongly affected by the presence of these conformational defects, a result that is confirmed by quantum transport calculations. Our work indicates that solution-based fabrication methods for carbon nanotube devices can have high conformational defect tolerance, except for defects with 5-10 nanometer bending radius. © 2005 American Institute of Physics.

The objective of this LDRD was to develop microdevice strategies for dealing with samples to be examined in biological detection systems. This includes three sub-components: namely, microdevice fabrication, sample delivery to the microdevice, and sample processing within the microdevice. The first component of this work focused on utilizing Sandia's surface micromachining technology to fabricate small volume (nanoliter) fluidic systems for processing small quantities of biological samples. The next component was to develop interfaces for the surface-micromachined silicon devices. We partnered with Micronics, a commercial company, to produce fluidic manifolds for sample delivery to our silicon devices. Pressure testing was completed to examine the strength of the bond between the pressure-sensitive adhesive layer and the silicon chip. We are also pursuing several other methods, both in house and external, to develop polymer-based fluidic manifolds for packaging silicon-based microfluidic devices. The second component, sample processing, is divided into two sub-tasks: cell collection and cell lysis. Cell collection was achieved using dielectrophoresis, which employs AC fields to collect cells at energized microelectrodes, while rejecting non-cellular particles. Both live and dead Staph. aureus bacteria have been collected using RF frequency dielectrophoresis. Bacteria have been separated from polystyrene microspheres using frequency-shifting dielectrophoresis. Computational modeling was performed to optimize device separation performance, and to predict particle response to the dielectrophoretic traps. Cell lysis is continuing to be pursued using microactuators to mechanically disrupt cell membranes. Novel thermal actuators, which can generate larger forces than previously tested electrostatic actuators, have been incorporated with and tested with cell lysis devices. Significant cell membrane distortion has been observed, but more experiments need to be conducted to determine the effects of the observed distortion on membrane integrity and cell viability. Finally, we are using a commercial PCR DNA amplification system to determine the limits of detectable sample size, and to examine the amplification of DNA bound to microspheres. Our objective is to use microspheres as capture-and-carry chaperones for small molecules such as DNA and proteins, enabling the capture and concentration of the small molecules using dielectrophoresis. Current tests demonstrated amplification of DNA bound to micron-sized polystyrene microspheres using 20-50 microliter volume size reactions.

This report contains the summary of the 'Magnetophoretic Bead Trapping in a High-Flowrate Biological Detection System' LDRD project 74795. The objective of this project is to develop a novel biodetection system for high-throughput sample analysis. The chief application of this system is in detection of very low concentrations of target molecules from a complex liquid solution containing many different constituents--some of which may interfere with identification of the target molecule. The system is also designed to handle air sampling by using an aerosol system (for instance a WESP - Wet Electro-Static Precipitator, or an impact spray system) to get air sample constituents into the liquid volume. The system described herein automatically takes the raw liquid sample, whether air converted or initially liquid matrix, and mixes in magnetic detector beads that capture the targets of interest and then performs the sample cleanup function, allowing increased sensitivity and eliminating most false positives and false negatives at a downstream detector. The surfaces of the beads can be functionalized in a variety of ways in order to maximize the number of targets to be captured and concentrated. Bacteria and viruses are captured using antibodies to surface proteins on bacterial cell walls or viral particle coats. In combination with a cell lysis or PCR (Polymerase Chain Reaction), the beads can be used as a DNA or RNA probe to capture nucleic acid patterns of interest. The sample cleanup capability of this system would allow different raw biological samples, such as blood or saliva to be analyzed for the presence of different infectious agents (e.g. smallpox or SARS). For future studies, we envision functionalizing bead surfaces to bind to chemical weapons agents, radio-isotopes, and explosives. The two main objectives of this project were to explore methods for enhancing the mixing of the capture microspheres in the sample, and to develop a novel high-throughput magnetic microsphere trap. We have developed a novel technique using the magnetic capture microspheres as 'stirrer bars' in a fluid sample to enhance target binding to the microsphere surfaces. We have also made progress in developing a polymer-MEMS electromagnet for trapping magnetic spheres in a high-flowrate fluid format.

In order to develop the next generation of high peak intensity lasers, new grating technology providing higher damage thresholds and large apertures is required. The current assumption is that this technical innovation will be multilayer dielectric gratings, wherein the uppermost layer of a thin film mirror is etched to create the desired binary phase grating. A variant of this is explored with the upper grating layer being a lower density gelatin-based volume phase grating in either sol-gel or dichromated gelatin. One key benefit is the elimination of the etching step.

Networked systems of low-cost, small, integrable chemical sensors will enable monitoring of Nonproliferation and Materials Control targets and chemical weapons threats. Sandia-designed prototype chemical sensor systems are undergoing extended field testing supported by DOE and other government agencies. A required surety component will be verification of microanalytical system performance, which can be achieved by providing a programmable source of chemical signature(s) for autonomous calibration of analytical systems. In addition, such a controlled chemical source could be used to dispense microaliquots of derivatization reagents, extending the analysis capability of chemical sensors to a wider range of targets. We have developed a microfabricated system for controlled release of selected compounds (calibrants) into the analytical stream of microsensor systems. To minimize pumping and valve requirements of microfluidic systems, and to avoid degradation issues associated with storage of dilute solutions, we have utilized thermally labile organic salts as solid-phase reservoir materials. Reproducible deposition of tetrapropyl ammonium hydroxide onto arrays of microfabricated heating elements can provide a pair of calibration marker compounds (one fast and one slow-eluting compound) for GC analyses. The use of this microaliquot gas source array for hydrogen generation is currently under further development. The goal of the latter effort will be to provide a source of high-pressure, low viscosity GC carrier gas for Sandia's next-generation microfabricated gas-phase chemical analysis systems.

Thermal property measurements of 18650 cells for the Space Shuttle Orbiter Advanced Hydraulic Power System (AHPS, formerly known as EAPU) have been performed. Cells were measured from three commercial manufacturers: E-One MOLI (12 cells), Panasonic (5 cells) and Sanyo (5 cells). Thermal property measurements of the MOLI 18650 cells included measurements of specific heat, self discharge (microcalorimetry), overcharge response and thermal runaway by accelerating rate calorimetry (ARC). The Panasonic and Sanyo cells were measured only for thermal runaway response in the ARC at increasing states of charge (3.8V, 4.0V, 4.2V, 4.3V).

This study describes the technical and economic benefits derived from adding an energy storage component to an existing building cooling, heating, and power system that uses microturbine generation to augment utility-provided power. Three different types of battery energy storage were evaluated: flooded lead-acid, valve-regulated lead-acid, and zinc/bromine. Additionally, the economic advantages of hybrid generation/storage systems were evaluated for a representative range of utility tariffs. The analysis was done using the Distributed Energy Technology Simulator developed for the Energy Storage Systems Program at Sandia National Laboratories by Energetics, Inc. The study was sponsored by the U.S. DOE Energy Storage Systems Program through Sandia National Laboratories and was performed in coordination with the University of Maryland's Center for Environmental Energy Engineering.

To date, all varistors used in ferroelectric neutron generators have been supplied from a single, proprietary source, General Electric Corporate Research and Development (GE CR&D). To protect against the vulnerability of a single source, Sandia initiated a program in the early 1980's to develop a second source for this material. A chemical preparation process for making homogeneous, high purity ZnO-based varistor powder was generated, scaled to production quantities, and transferred to external suppliers. In 1992, the chem-prep varistor program was suspended when it appeared there was sufficient inventory of GE CR&D material to supply ferroelectric neutron generator production for many years. In 1999, neutron generator production schedules increased substantially, resulting in a predicted exhaustion of the existing supply of varistor material within five years. The chem-prep program was restarted in January, 2000. The goals of the program were to (1) duplicate the chem-prep powder synthesis process that had been qualified for WR production, (2) demonstrate sintered billets from the chem-prep powder met requirements, (3) develop a process for rod fabrication and demonstrate that all component specifications could be met, and (4) optimize the process from powder synthesis through component fabrication for full-scale production. The first three of these goals have been met and are discussed in this report. A facility for the fabrication of production quantities of chem-prep powder has been established. All batches since the restart have met compositional requirements, but differences in sintering behavior between the original process and the restarted process were noted. Investigation into the equipment, precipitant stoichiometry, and powder processing procedures were not able to resolve the discrepancies. It was determined that the restarted process, which incorporated Na doping for electrical stability (a process that was not introduced until the end of the initial program and had not been investigated for processing effects), was responsible for the differences. Rod components fabricated since the restart have met requirements and have performed at a level comparable to chem-prep rods from the original program and GE CR&D rods currently in production.

This reports describes the successful extension of artificial immune systems from the domain of computer security to the domain of real time control systems for robotic vehicles. A biologically-inspired computer immune system was added to the control system of two different mobile robots. As an additional layer in a multi-layered approach, the immune system is complementary to traditional error detection and error handling techniques. This can be thought of as biologically-inspired defense in depth. We demonstrated an immune system can be added with very little application developer effort, resulting in little to no performance impact. The methods described here are extensible to any system that processes a sequence of data through a software interface.

Abstract not provided.

Abstract not provided.

The Unstructured Time-Domain ElectroMagnetics (UTDEM) portion of the EMPHASIS suite solves Maxwell's equations using finite-element techniques on unstructured meshes. This document provides user-specific information to facilitate the use of the code for applications of interest.

Semiconductor bridges (SCBs) are commonly used as initiators for explosive and pyrotechnic devices. Their advantages include reduced voltage and energy requirements and exceptional safety features. Moreover, the design of systems which implement SCBs can be expedited using electrical simulation software. Successful use of this software requires that certain parameters be correctly chosen. In this paper, we explain how these parameters can be identified using optimization. We describe the problem focusing on the application of a direct optimization method for its solution, and present some numerical results.

A commercial plasma enhanced chemical vapor deposition (PECVD) technique with planetary substrate rotation was used to apply a thin (200-400 nm thick) conformal diamond-like carbon (DLC) coating (known as a diamond-like nanocomposite (DLN)) on LIGA fabricated Ni-Mn alloy parts. The PECVD technique is known to overcome the drawbacks associated with the line-of-sight nature of physical vapor deposition (PVD) and substrate heating inherent with traditional chemical vapor deposition (CVD). The purpose of the present study is to characterize the coverage, adhesion, and tribological (friction and wear) behavior of DLN coatings applied to planar and sidewall surfaces of small featured LIGA Ni-Mn fabricated parts, e.g. 280 {micro}m thick sidewalls. Friction and wear tests were performed in dry nitrogen, dry air, and air with 50% RH at Hertzian contact pressures ranging from 0.3 to 0.6 GPa. The friction coefficient of bare Ni-Mn alloy was determined to be 0.9. In contrast, low friction coefficients ({approx}0.02 in dry nitrogen and {approx}0.2 in 50% RH air) and minimal amount of wear were exhibited for the DLN coated LIGA Ni-Mn alloy parts and test coupons. This behavior was due to the ability of the coating to transfer to the rubbing counterface providing low interfacial shear at the sliding contact; resultantly, coating one surface was adequate for low friction and wear. In addition, a 30 nm thick titanium bond layer was determined to be necessary for good adhesion of DLN coating to Ni-Mn alloy substrates. Raman spectroscopy and cross-sectional SEM with energy dispersive x-ray analysis revealed that the DLN coatings deposited by the PECVD with planetary substrate rotation covered both the planar and sidewall surfaces of LIGA fabricated parts, as well as narrow holes of 300 {micro}m (0.012 inch) diameter.

We have found that developing a computational framework for reconstructing error control codes for engineered data and ultimately for deciphering genetic regulatory coding sequences is a challenging and uncharted area that will require advances in computational technology for exact solutions. Although exact solutions are desired, computational approaches that yield plausible solutions would be considered sufficient as a proof of concept to the feasibility of reverse engineering error control codes and the possibility of developing a quantitative model for understanding and engineering genetic regulation. Such evidence would help move the idea of reconstructing error control codes for engineered and biological systems from the high risk high payoff realm into the highly probable high payoff domain. Additionally this work will impact biological sensor development and the ability to model and ultimately develop defense mechanisms against bioagents that can be engineered to cause catastrophic damage. Understanding how biological organisms are able to communicate their genetic message efficiently in the presence of noise can improve our current communication protocols, a continuing research interest. Towards this end, project goals include: (1) Develop parameter estimation methods for n for block codes and for n, k, and m for convolutional codes. Use methods to determine error control (EC) code parameters for gene regulatory sequence. (2) Develop an evolutionary computing computational framework for near-optimal solutions to the algebraic code reconstruction problem. Method will be tested on engineered and biological sequences.

The CABle ANAlysis (CABANA) portion of the EMPHASIS suite is designed specifically for the simulation of cable SGEMP. The code can be used to evaluate the response of a specific cable design to threat or to compare and minimize the relative response of difference designs. This document provides user-specific information to facilitate the application of the code to cables of interest.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The ability to detect Weapons of Mass Destruction biological agents rapidly and sensitively is vital to homeland security, spurring development of compact detection systems at Sandia and elsewhere. One such system is Sandia's microseparations-based pChemLab. Many bio-agents are serious health threats even at extremely low concentrations. Therefore, a universal challenge for detection systems is the efficient collection and selective transport of highly diffuse bio-agents against the enormous background of benign particles and species ever present in the ambient environment. We have investigated development of a ''front end'' system for the collection, preconcentration, and selective transport of aerosolized biological agents from dilute (1-10 active particles per liter of air) atmospheric samples, to ultimate concentrations of {approx}20 active particles per microliter of liquid, for interface with microfluidic-based analyses and detection systems. Our approach employs a Sandia-developed aerosol particle-focusing microseparator array to focus size-selected particles into a mating microimpinger array of open microfluidic transport channels. Upon collection (i.e., impingement, submergence, and liquid suspension), microfluidic dielectrophoretic particle concentrators and sorters can be employed to further concentrate and selectively transport bio-agent particles to the sample preparation stages of microfluidic analyses and detection systems. This report documents results in experimental testing, modeling and analysis, component design, and materials fabrication critical to establishing proof-of-principle for this collection ''front end''. Outstanding results have been achieved for the aerodynamic microseparator, and for the post-collection dielectrophoretic concentrator and sorter. Results have been obtained for the microimpinger, too, but issues of particle-trapping by surface tension in liquid surfaces have proven difficult. Subsequent particle submergence into liquid suspension for microfluidic transport has been demonstrated only inefficiently despite significant and varied effort. Importantly, the separate technologies whose development is described, (inertial microseparator, dielectrophoretic corduroy concentrator/sorter) should each, independently, prove greatly useful in a variety of additional applications.

Earth Power Resources, Inc. recently completed a combined rotary/core hole to a depth of 3,813 feet at it's Hot Sulphur Springs Tuscarora Geothermal Power Project Lease Area located 70-miles north of Elko, Nevada. Previous geothermal exploration data were combined with geologic mapping and newly acquired seismic-reflection data to identify a northerly tending horst-graben structure approximately 2,000 feet wide by at least 6,000 feet long with up to 1,700 feet of vertical offset. The well (HSS-2) was successfully drilled through a shallow thick sequence of altered Tertiary Volcanic where previous exploration wells had severe hole-caving problems. The ''tight-hole'' drilling problems were reduced using drilling fluids consisting of Polymer-based mud mixed with 2% Potassium Chloride (KCl) to reduce Smectite-type clay swelling problems. Core from the 330 F fractured geothermal reservoir system at depths of 2,950 feet indicated 30% Smectite type clays existed in a fault-gouge zone where total loss of circulation occurred during coring. Smectite-type clays are not typically expected at temperatures above 300 F. The fracture zone at 2,950 feet exhibited a skin-damage during injection testing suggesting that the drilling fluids may have caused clay swelling and subsequent geothermal reservoir formation damage. The recent well drilling experiences indicate that drilling problems in the shallow clays at Hot Sulphur Springs can be reduced. In addition, average penetration rates through the caprock system can be on the order of 25 to 35 feet per hour. This information has greatly reduced the original estimated well costs that were based on previous exploration drilling efforts. Successful production formation drilling will depend on finding drilling fluids that will not cause formation damage in the Smectite-rich fractured geothermal reservoir system. Information obtained at Hot Sulphur Springs may apply to other geothermal systems developed in volcanic settings.

The AlGaInN material system is used for virtually all advanced solid state lighting and short wavelength optoelectronic devices. Although metal-organic chemical vapor deposition (MOCVD) has proven to be the workhorse deposition technique, several outstanding scientific and technical challenges remain, which hinder progress and keep RD&A costs high. The three most significant MOCVD challenges are: (1) Accurate temperature measurement; (2) Reliable and reproducible p-doping (Mg); and (3) Low dislocation density GaN material. To address challenge (1) we designed and tested (on reactor mockup) a multiwafer, dual wavelength, emissivity-correcting pyrometer (ECP) for AlGaInN MOCVD. This system simultaneously measures the reflectance (at 405 and 550 nm) and emissivity-corrected temperature for each individual wafer, with the platen signal entirely rejected. To address challenge (2) we measured the MgCp{sub 2} + NH{sub 3} adduct condensation phase diagram from 65-115 C, at typical MOCVD concentrations. Results indicate that it requires temperatures of 80-100 C in order to prevent MgCp{sub 2} + NH{sub 3} adduct condensation. Modification and testing of our research reactor will not be complete until FY2005. A new commercial Veeco reactor was installed in early FY2004, and after qualification growth experiments were conducted to improve the GaN quality using a delayed recovery technique, which addresses challenge (3). Using a delayed recovery technique, the dislocation densities determined from x-ray diffraction were reduced from 2 x 10{sup 9} cm{sup -2} to 4 x 10{sup 8} cm{sup -2}. We have also developed a model to simulate reflectance waveforms for GaN growth on sapphire.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Structured adaptive mesh refinement (SAMR) methods are being widely used for computer simulations of various physical phenomena. Parallel implementations potentially offer realistic simulations of complex, three-dimensional applications. But achieving good scalability for large-scale applications is non-trivial. Performance is limited by the partitioners ability to efficiently use the underlying computer's resources. The goal of our research project is to improve scalability for general SAMR applications executing on general parallel computers. We engineer the dynamically adaptive meta-partitioner, able to select and configure the most appropriate partitioning method at run-time, based on system and application state. This presentation gives an overview of our project, reports on recent achievements, and discusses the project's significance in a wider scientific context.

The Container Analysis Fire Environment (CAFE) computer code has been developed to model all relevant fire physics for predicting the thermal response of massive objects engulfed in large fires. It provides realistic fire thermal boundary conditions for use in design of radioactive material packages and in risk-based transportation studies. The CAFE code can be coupled to commercial finite-element codes such as MSC PATRAN/THERMAL and ANSYS. This coupled system of codes can be used to determine the internal thermal response of finite element models of packages to a range of fire environments. This document is a user manual describing how to use the three-dimensional version of CAFE, as well as a description of CAFE input and output parameters. Since this is a user manual, only a brief theoretical description of the equations and physical models is included.

The formulation, implementation and usage of a numerical solution verification code is described. This code uses the Richardson extrapolation procedure to estimate the order of accuracy and error of a computational program solution. It evaluates multiple solutions performed in numerical grid convergence studies to verify a numerical algorithm implementation. Analyses are performed on both structured and unstructured grid codes. Finite volume and finite element discretization programs are examined. Two and three-dimensional solutions are evaluated. Steady state and transient solution analysis capabilities are present in the verification code. Multiple input data bases are accepted. Benchmark options are included to allow for minimal solution validation capability as well as verification.

This report is a white paper summarizing the literature and different approaches to the problem of calibrating computer model parameters in the face of model uncertainty. Model calibration is often formulated as finding the parameters that minimize the squared difference between the model-computed data (the predicted data) and the actual experimental data. This approach does not allow for explicit treatment of uncertainty or error in the model itself: the model is considered the %22true%22 deterministic representation of reality. While this approach does have utility, it is far from an accurate mathematical treatment of the true model calibration problem in which both the computed data and experimental data have error bars. This year, we examined methods to perform calibration accounting for the error in both the computer model and the data, as well as improving our understanding of its meaning for model predictability. We call this approach Calibration under Uncertainty (CUU). This talk presents our current thinking on CUU. We outline some current approaches in the literature, and discuss the Bayesian approach to CUU in detail.

This report presents a modification of a previous model for the statistical distribution of linear antenna impedance. With this modification a simple formula is determined which yields accurate results for all ratios of modal spectral width to spacing. It is shown that the reactance formula approaches the known unit Lorentzian in the lossless limit.

Sandia National Laboratories performs many expensive tests using inertial measurement units (IMUs)--systems that use accelerometers, gyroscopes, and other sensors to measure flight dynamics in three dimensions. For the purpose of this report, the metrics used to evaluate an IMU are cost, size, performance, resolution, upgradeability and testing. The cost of a precision IMU is very high and can cost hundreds of thousands of dollars. Thus the goals and results of this project are as follows: (1) Examine the data flow in an IMU and determine a generic IMU design. (2) Discuss a high cost IMU implementation and its theoretically achievable results. (3) Discuss design modifications that would save money for suited applications. (4) Design and implement a low cost IMU and discuss its theoretically achievable results. (5) Test the low cost IMU and compare theoretical results with empirical results. (6) Construct a more streamlined printed circuit board design reducing noise, increasing capabilities, and constructing a self-contained unit. Using these results, we can compare a high cost IMU versus a low cost IMU using the metrics from above. Further, we can examine and suggest situations where a low cost IMU could be used instead of a high cost IMU for saving cost, size, or both.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The absence of agreed definitions and metrics for supercomputer RAS obscures meaningful discussion of the issues involved and hinders their solution. This paper seeks to foster a common basis for communication about supercomputer RAS, by proposing a system state model, definitions, and measurements. These are modeled after the SEMI-E10 specification which is widely used in the semiconductor manufacturing industry.

Abstract not provided.

Abstract not provided.

Relatively small motion measurement errors manifest themselves principally as a phase error in Synthetic Aperture Radar (SAR) complex data samples, and if large enough become observable as a smearing, blurring, or other degradation in the image. The phase error function can be measured and then deconvolved from the original data to compensate for the presumed motion error, ultimately resulting in a well-focused image. Techniques that do this are termed 'autofocus' algorithms. A very popular autofocus algorithm is the Phase Gradient Autofocus (PGA) algorithm. The nearly universal, and typically reasonable, assumption is that the motion errors are less than the range resolution of the radar, allowing solely a phase correction to suffice. Very large relative motion measurement errors manifest themselves as an unexpected additional shifting or migration of target locations beyond any deterministic migration during the course of the synthetic aperture. Degradation in images from data exhibiting errors of this magnitude are substantial, often rendering the image completely useless. When residual range migration due to either real or apparent motion errors exceeds the range resolution, conventional autofocus algorithms fail. Excessive residual migration is increasingly encountered as resolutions become finer, less expensive inertial sensors are used, and operating ranges become longer (due to atmospheric phenomena). A new migration-correction autofocus algorithm has been developed that estimates the excessive residual migration and applies phase and frequency corrections to properly focus the image. This overcomes the conventional constraint that motion errors not exceed the SAR range resolution.

Abstract not provided.

This report describes an integrated approach for designing communication, sensing, and control systems for mobile distributed systems. Graph theoretic methods are used to analyze the input/output reachability and structural controllability and observability of a decentralized system. Embedded in each network node, this analysis will automatically reconfigure an ad hoc communication network for the sensing and control task at hand. The graph analysis can also be used to create the optimal communication flow control based upon the spatial distribution of the network nodes. Edge coloring algorithms tell us that the minimum number of time slots in a planar network is equal to either the maximum number of adjacent nodes (or degree) of the undirected graph plus some small number. Therefore, the more spread out that the nodes are, the fewer number of time slots are needed for communication, and the smaller the latency between nodes. In a coupled system, this results in a more responsive sensor network and control system. Network protocols are developed to propagate this information, and distributed algorithms are developed to automatically adjust the number of time slots available for communication. These protocols and algorithms must be extremely efficient and only updated as network nodes move. In addition, queuing theory is used to analyze the delay characteristics of Carrier Sense Multiple Access (CSMA) networks. This report documents the analysis, simulation, and implementation of these algorithms performed under this Laboratory Directed Research and Development (LDRD) effort.

The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass-lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov-Galerkin streamline upwind schemes in the context of advection-diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection-dominated problems.

Electrical contact resistance testing was performed by hot-switching a simulated gold-platinum metal microelectromechanical systems contact. The experimental objective was to determine the sensitivity of the contact resistance degradation to current level and environment. The contact resistance increased sharply after 100 hot-switched cycles in air. Hot-switching at a reduced current and in nitrogen atmosphere curtailed contact resistance degradation by several orders of magnitude. The mechanism responsible for the resistance degradation was found to be arc-induced decomposition of adsorbed surface contaminants.

Intra-molecular cross-linking has been suggested as a method of obtaining distance constraints that would be useful in developing structural models of proteins. Recent work published on intra-molecular cross-linking for protein structural studies has employed commercially available primary amine selective reagents that can cross-link lysine residues to other lysine residues or the amino terminus. Previous work using these cross-linkers has shown that for several proteins of known structure, the number of cross-links that can be obtained experimentally may be small compared to what would be expected from the known structure, due to the relative reactivity, distribution, and solvent accessibility of the lysines in the protein sequence. To overcome these limitations we have investigated the use of cross-linking reagents that can react with other reactive sidechains in proteins. We used 1-Ethyl-3-(3-dimethylaminopropyl) carbodiimide hydrochloride (EDC) to activate the carboxylic acid containing residues, aspartic acid (D), glutamic acid (E), and the carboxy terminus (O), for cross-linking reactions. Once activated, the DEO sidechains can react to form 'zero-length' cross-links with nearby primary amine containing resides, lysines (K) and the amino terminus (X), via the formation of a new amide bond. We also show that the EDC-activated DEO sidechains can be cross-linked to each other using dihydrazides, two hydrazide moieties connected by an alkyl cross-linker ann of variable length. Using these reagents, we have found three new 'zero-length' cross-links in ubiquitin consistent with its known structure (M1-E16, M1-E18, and K63-E64). Using the dihydrazide cross-linkers, we have identified 2 new cross-links (D21-D32 and E24-D32) unambiguously. Using a library of dihydrazide cross-linkers with varying arm length, we have shown that there is a minimum arm length required for the DEO-DEO cross-links of 5.8 angstroms. These results show that additional structural information can be obtained by exploiting new cross-linker chemistry, increasing the probability that the protein target of choice will yield sufficient distance constraints to develop a structural model.

Membrane proteins make up a diverse and important subset of proteins for which structural information is limited. In this study, chemical cross-linking and mass spectrometry were used to explore the structure of the G-protein-coupled photoreceptor bovine rhodopsin in the dark-state conformation. All experiments were performed in rod outer segment membranes using amino acid 'handles' in the native protein sequence and thus minimizing perturbations to the native protein structure. Cysteine and lysine residues were covalently cross-linked using commercially available reagents with a range of linker arm lengths. Following chemical digestion of cross-linked protein, cross-linked peptides were identified by accurate mass measurement using liquid chromatography-fourier transform mass spectrometry and an automated data analysis pipeline. Assignments were confirmed and, if necessary, resolved, by tandem MS. The relative reactivity of lysine residues participating in cross-links was evaluated by labeling with NHS-esters. A distinct pattern of cross-link formation within the C-terminal domain, and between loop I and the C-terminal domain, emerged. Theoretical distances based on cross-linking were compared to inter-atomic distances determined from the energy-minimized X-ray crystal structure and Monte Carlo conformational search procedures. In general, the observed cross-links can be explained by re-positioning participating side-chains without significantly altering backbone structure. One exception, between C3 16 and K325, requires backbone motion to bring the reactive atoms into sufficient proximity for cross-linking. Evidence from other studies suggests that residues around K325 for a region of high backbone mobility. These findings show that cross-linking studies can provide insight into the structural dynamics of membrane proteins in their native environment.

Clostridial neurotoxins, such as botulinum and tetanus, are generally thought to invade neural cells through a process of high affinity binding mediated by gangliosides, internalization via endosome formation, and subsequent membrane penetration of the catalytic domain activated by a pH drop in the endosome. This surface recognition and internalization process is still not well understood with regard to what specific membrane features the toxins target, the intermolecular interactions between bound toxins, and the molecular conformational changes that occur as a result of pH lowering. In an effort to elucidate the mechanism of tetanus toxin binding and permeation through the membrane a simple yet representative model was developed that consisted of the ganglioside G{sub tlb} incorporated in a bilayer of cholesterol and DPPC (dipalmitoylphosphatidyl choline). The bilayers were stable over time yet sensitive towards the binding and activity of whole toxin. A liposome leakage study at constant pH as well as with a pH gradient, to mimic the processes of the endosome, was used to elucidate the effect of pH on the toxin's membrane binding and permeation capability. Topographic imaging of the membrane surface, via in situ tapping mode AFM, provided nanoscale characterization of the toxin's binding location and pore formation activity.

This report describes a system level study on the use of a swarm of sea gliders to detect, confirm and kill littoral submarine threats. The report begins with a description of the problem and derives the probability of detecting a constant speed threat without networking. It was concluded that glider motion does little to improve this probability unless the speed of a glider is greater than the speed of the threat. Therefore, before detection, the optimal character for a swarm of gliders is simply to lie in wait for the detection of a threat. The report proceeds by describing the effect of noise on the localization of a threat once initial detection is achieved. This noise is estimated as a function of threat location relative to the glider and is temporally reduced through the use of an information or Kalman filtering. In the next section, the swarm probability of confirming and killing a threat is formulated. Results are compared to a collection of stationary sensors. These results show that once a glider has the ability to move faster than the threat, the performance of the swarm is equal to the performance of a stationary swarm of gliders with confirmation and kill ranges equal to detection range. Moreover, at glider speeds greater than the speed of the threat, swarm performance becomes a weak function of speed. At these speeds swarm performance is dominated by detection range. Therefore, to future enhance swarm performance or to reduce the number of gliders required for a given performance, detection range must be increased. Communications latency is also examined. It was found that relatively large communication delays did little to change swarm performance. Thus gliders may come to the surface and use SATCOMS to effectively communicate in this application.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Complex simulations (in particular, those involving multiple coupled physics) cannot be understood solely using geometry-based visualizations. Such visualizations are necessary in interpreting results and gaining insights into kinematics, however they are insufficient when striving to understand why or how something happened, or when investigating a simulation's dynamic evolution. For multiphysics simulations (e.g. those including solid dynamics with thermal conduction, magnetohydrodynamics, and radiation hydrodynamics) complex interactions between physics and material properties take place within the code which must be investigated in other ways. Drawing on the extensive previous work in view coordination, brushing and linking techniques, and powerful visualization libraries, we have developed Prism, an application targeted for a specific analytic need at Sandia National Laboratories. This multiview scientific visualization tool tightly integrates geometric and phase space views of simulation data and material models. Working closely with analysts, we have developed this production tool to promote understanding of complex, multiphysics simulations. We discuss the current implementation of Prism, along with specific examples of results obtained by using the tool.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Economists, systems analysts, engineers, regulatory specialists, and other experts were assembled from academia, the national laboratories, and the energy industry to discuss present restoration practices (many have already been defined to the level of operational protocols) in the sectors of the energy infrastructure as well as other infrastructures, to identify whether economics, a discipline concerned with the allocation of scarce resources, is explicitly or implicitly a part of restoration strategies, and if there are novel economic techniques and solution methods that could be used help encourage the restoration of energy services more quickly than present practices or to restore service more efficiently from an economic perspective. AcknowledgementsDevelopment of this work into a coherent product with a useful message has occurred thanks to the thoughtful support of several individuals:Kenneth Friedman, Department of Energy, Office of Energy Assurance, provided the impetus for the work, as well as several suggestions and reminders of direction along the way. Funding from DOE/OEA was critical to the completion of this effort.Arnold Baker, Chief Economist, Sandia National Laboratories, and James Peerenboom, Director, Infrastructure Assurance Center, Argonne National Laboratory, provided valuable contacts that helped to populate the authoring team with the proper mix of economists, engineers, and systems and regulatory specialists to meet the objectives of the work.Several individuals provided valuable review of the document at various stages of completion, and provided suggestions that were valuable to the editing process. This list of reviewers includes Jeffrey Roark, Economist, Tennessee Valley Authority; James R. Dalrymple, Manager of Transmission System Services and Transmission/Power Supply, Tennessee Valley Authority; William Mampre, Vice President, EN Engineering; Kevin Degenstein, EN Engineering; and Patrick Wilgang, Department of Energy, Office of Energy Assurance.With many authors, creating a document with a single voice is a difficult task. Louise Maffitt, Senior Research Associate, Institute for Engineering Research and Applications at New Mexico Institute of Mining & Technology (on contract to Sandia National Laboratories) served a vital role in the development of this document by taking the unedited material (in structured format) and refining the basic language so as to make the flow of the document as close to a single voice as one could hope for. Louise's work made the job of reducing the content to a readable length an easier process. Additional editorial suggestions from the authors themselves, particularly from Sam Flaim, Steve Folga, and Doug Gotham, expedited this process.

Visualization in scientific computing refers to the process of transforming data produced by a simulation into graphical representations that help scientific users interpret the results. While the back-end rendering phase of this work can be performed efficiently in graphics card hardware, the front-end 'post processing' portion of visualization is currently performed entirely in software. Field-Programmable Gate Arrays (FPGAs) are an attractive option for accelerating post-processing operations because they enable users to offload computations into reconfigurable hardware. A key challenge in utilizing FPGAs for this work is developing an infrastructure that allows FPGAs to be integrated into a distributed visualization system. We propose a networked approach, where each post-processing FPGA is equipped with specialized network interface (NI) hardware that is capable of transporting graphics commands across the network to existing rendering resources. In this paper we discuss a NI for FPGAs that is comprised of a Chromium OpenGL interface, a TCP Offload Engine, and a Gigabit Ethernet module. A prototype system has been tested for a distributed isosurfacing application.

Abstract not provided.

Abstract not provided.

The Sandia Z accelerator has become a unique platform to study matter at extreme conditions. The large currents (20 MA, 200-300 ns rise time) and magnetic fields (several MG) produced by Z generate magnetic compression in the multi-Mbar regime, enabling quasi-isentropic compression experiments (ICE) to several Mbar stresses. Thus, the Z platform is useful in determining high stress material isentropes, performing phase transition studies (including rapid solidification), obtaining constitutive property information, and estimating material strength at high stress. Furthermore, the magnetic pressure can also accelerate macroscopic flyer plates to velocities in excess of 30 km/s. Thus, impact experiments can be performed to ultra-high pressures. Furthermore, the adiabatic release response of materials can be investigated through shock and release experiments, allowing hot, dense liquid states to be probed. The Z platform allows a large expanse of the equation of state surface to be explored enabling new and exciting material dynamics experiments. Specific examples from each of the areas mentioned above will be discussed.

Concerns about acts of terrorism against critical infrastructures have been on the rise for several years. Critical infrastructures are those physical structures and information systems (including cyber) essential to the minimum operations of the economy and government. The President's Commission on Critical Infrastructure Protection (PCCIP) probed the security of the nation's critical infrastructures. The PCCIP determined the water infrastructure is highly vulnerable to a range of potential attacks. In October 1997, the PCCIP proposed a public/private partnership between the federal government and private industry to improve the protection of the nation's critical infrastructures. In early 2000, the EPA partnered with the Awwa Research Foundation (AwwaRF) and Sandia National Laboratories to create the Risk Assessment Methodology for Water Utilities (RAM-W{trademark}). Soon thereafter, they initiated an effort to create a template and minimum requirements for water utility Emergency Response Plans (ERP). All public water utilities in the US serving populations greater than 3,300 are required to undertaken both a vulnerability assessment and the development of an emergency response plan. This paper explains the initial steps of RAM-W{trademark} and then demonstrates how the security risk assessment is fundamental to the ERP. During the development of RAM-W{trademark}, Sandia performed several security risk assessments at large metropolitan water utilities. As part of the scope of that effort, ERPs at each utility were reviewed to determine how well they addressed significant vulnerabilities uncovered during the risk assessment. The ERP will contain responses to other events as well (e.g. natural disasters) but should address all major findings in the security risk assessment.

This report summarizes progress from the Laboratory Directed Research and Development (LDRD) program during fiscal year 2004. In addition to a programmatic and financial overview, the report includes progress reports from 352 individual R and D projects in 15 categories. The 15 categories are: (1) Advanced Concepts; (2) Advanced Manufacturing; (3) Biotechnology; (4) Chemical and Earth Sciences; (5) Computational and Information Sciences; (6) Differentiating Technologies; (7) Electronics and Photonics; (8) Emerging Threats; (9) Energy and Critical Infrastructures; (10) Engineering Sciences; (11) Grand Challenges; (12) Materials Science and Technology; (13) Nonproliferation and Materials Control; (14) Pulsed Power and High Energy Density Sciences; and (15) Corporate Objectives.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Experimental results [1] for the reflection coefficient of shock compressed xenon are compared with results from quantum molecular dynamics calculations with density functional theory (DFT). The real part of the optical conductivity is calculated within the Kubo-Greenwood formalism and the Kramers-Kroenig relations are used to generate the reflectivity and other optical properties. Improved agreement over non-ideal plasma theory [2] is found with the DFT calculations, but significant differences with the data remain. Since DFT in the various local density approximations tends to underestimate the band gap and overestimate the free electron population, we have used the ionizations from [2] to correct the DFT band gaps. This results in much improved agreement with the xenon reflectivity data and demonstrates a new approach to correcting DFT band gaps.

INVICE (INVerse analysis of Isentropic Compression Experiments) is a FORTRAN computer code that implements the inverse finite-difference method to analyze velocity data from isentropic compression experiments. This report gives a brief description of the methods used and the options available in the first beta version of the code, as well as instructions for using the code.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Agile ready-when-needed patterning of refractive index structures in photosensitive materials requires an understanding of the impact of local application environment on mechanisms contributing to the desired photoinduced index change. The present work examines the impact of atmosphere on the photosensitive response of poly(methylphenylsilane) (PMPS) thin films whose high photoinduced index'change under low incident optical fluence make them attractive candidates for such applications. Changes in optical absorption and refractive index are investigated after exposure to ultraviolet (UV) light resonant with the lowest energy transition exhibited by the Si-Si backbone structure in the material. A comparison between photoinduced absorption changes for thin films exposed in an air atmosphere versus those observed for samples subjected to a nitrogen environment during photoexposure is made for the first time. The study reveals that the anaerobic conditions of the nitrogen atmosphere significantly reduce the photosensitive response of the material to light. These results are discussed in terms of photooxidation processes within the polysilane structure and in the context of the need for predictable photosensitive refractive index change in varied photoimprinting environments. © 2004 Elsevier B.V. All rights reserved.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Metallic Phases in extraterrestrial materials are composed of Fe-Ni with minor amounts of Co, P, Si, Cr, etc. Electron microscopy techniques (SEM, TEM, EPMA, AEM) have been used for almost 50 years to study micron and submicron microscopic features in the metal phases (Fig. 1) such as clear taenite, cloudy zone, plessite, etc [1,2]. However lack of instrumentation to prepare TEM thin foils in specific sample locations and to obtain micro-scale crystallographic data have limited these investigations. New techniques such as the focused ion beam (FIB) and the electron backscatter electron diffraction (EBSD) techniques have overcome these limitations. The application of the FIB instrument has allowed us to prepare {approx}10 um long by {approx} 5um deep TEM thin sections of metal phases from specific regions of metal particles, in chondrites, irons and stony iron meteorites, identified by optical and SEM observation. Using a FEI dual beam FIB we were able to study very small metal particles in samples of CH chondrites [3] and zoneless plessite (ZP) in ordinary chondrites. Fig. 2 shows a SEM photomicrograph of a {approx}40 um ZP particle in Kernouve, a H6 chondrite. Fig. 3a,b shows a TEM photograph of a section of the FIB prepared TEM foil of the ZP particle and a Ni trace through a tetrataenite/kamacite region of the particle. It has been proposed that the Widmanstatten pattern in low P iron meteorites forms by martensite decomposition, via the reaction {gamma} {yields} {alpha}{sub 2} + {gamma} {yields} {alpha} + {gamma} in which {alpha}{sub 2}, martensite, decomposes to the equilibrium {alpha} and {gamma} phases during the cooling process [4]. In order to show if this mechanism for Widmanstatten pattern formation is correct, crystallographic information is needed from the {gamma} or taenite phases throughout a given meteorite. The EBSD technique was employed in this study to obtain the orientation of the taenite surrounding the initial martensite phase and the kamacite which forms as {alpha}{sub 2} or as Widmanstatten plates in a series of IVB irons. Fig. 4a,b shows EBSD orientation maps of taenite and kamacite from the Tawallah Valley IVB iron. We observe that the orientation of the taenite in the IVB meteorites is the same throughout the sample consistent with the orientation of the high temperature single phase taenite before formation of the Widmanstatten pattern.

Abstract not provided.

Abstract not provided.

The potential energy surface for the reaction between OH and acetylene has been calculated using the RQCISD(T) method and extrapolated to the complete basis-set limit. Rate coefficients were determined for a wide range of temperatures and pressures, based on this surface and the solution of the one-dimensional and two-dimensional master equations. With a small adjustment to the association energy barrier (1.1 kcal/mol), agreement with experiments is good, considering the discrepancies in such data. The rate coefficient for direct hydrogen abstraction is significantly smaller than that commonly used in combustion models. Also in contrast to previous models, ketene + H is found to be the main product at normal combustion conditions. At low temperatures and high pressures, stabilization of the C{sub 2}H{sub 2}OH adduct is the dominant process. Rate coefficient expressions for use in modeling are provided.

The National Spent Nuclear Fuel Program, located at the Idaho National Laboratory (INL), coordinates and integrates national efforts in management and disposal of US Department of Energy (DOE)-owned spent nuclear fuel. These management functions include development of standardised systems for long-term disposal in the proposed Yucca Mountain repository. Nuclear criticality control measures are needed in these systems to avoid restrictive fissile loading limits because of the enrichment and total quantity of fissile material in some types of the DOE spent nuclear fuel. This need is being addressed by development of corrosion-resistant, neutron-absorbing structural alloys for nuclear criticality control. This paper outlines results of a metallurgical development programme that is investigating the alloying of gadolinium into a nickel-chromium-molybdenum alloy matrix. Gadolinium has been chosen as the neutron absorption alloying element due to its high thermal neutron absorption cross section and low solubility in the expected repository environment. The nickel-chromium-molybdenum alloy family was chosen for its known corrosion performance, mechanical properties, and weldability. The workflow of this programme includes chemical composition definition, primary and secondary melting studies, ingot conversion processes, properties testing, and national consensus codes and standards work. The microstructural investigation of these alloys shows that the gadolinium addition is present in the alloy as a gadolinium-rich second phase. The mechanical strength values are similar to those expected for commercial Ni-Cr-Mo alloys. The alloys have been corrosion tested with acceptable results. The initial results of weldability tests have also been acceptable. Neutronic testing in a moderated critical array has generated favourable results. An American Society for Testing and Materials material specification has been issued for the alloy and a Code Case has been submitted to the American Society of Mechanical Engineers for code qualification.

Abstract not provided.

Wear is a critical factor in determining the durability of microelectromechanical systems (MEMS). While the reliability of polysilicon MEMS has received extensive attention, the mechanisms responsible for this failure mode at the microscale have yet to be conclusively determined. We have used on-chip polycrystalline silicon side-wall friction MEMS specimens to study active mechanisms during sliding wear in ambient air. Worn parts were examined by analytical scanning and transmission electron microscopy, while local temperature changes were monitored using advanced infrared microscopy. Observations show that small amorphous debris particles ({approx}50-100 nm) are removed by fracture through the silicon grains ({approx}500 nm) and are oxidized during this process. Agglomeration of such debris particles into larger clusters also occurs. Some of these debris particles/clusters create plowing tracks on the beam surface. A nano-crystalline surface layer ({approx}20-200 nm), with higher oxygen content, forms during wear at and below regions of the worn surface; its formation is likely aided by high local stresses. No evidence of dislocation plasticity or of extreme local temperature increases was found, ruling out the possibility of high temperature-assisted wear mechanisms.

Microanalysis is typically performed to analyze the near surface of materials. There are many instances where chemical information about the third spatial dimension is essential to the solution of materials analyses. The majority of 3D analyses however focus on limited spectral acquisition and/or analysis. For truly comprehensive 3D chemical characterization, 4D spectral images (a complete spectrum from each volume element of a region of a specimen) are needed. Furthermore, a robust statistical method is needed to extract the maximum amount of chemical information from that extremely large amount of data. In this paper, an example of the acquisition and multivariate statistical analysis of 4D (3-spatial and 1-spectral dimension) x-ray spectral images is described. The method of utilizing a single- or dual-beam FIB (w/o or w/SEM) to get at 3D chemistry has been described by others with respect to secondary-ion mass spectrometry. The basic methodology described in those works has been modified for comprehensive x-ray microanalysis in a dual-beam FIB/SEM (FEI Co. DB-235). In brief, the FIB is used to serially section a site-specific region of a sample and then the electron beam is rastered over the exposed surfaces with x-ray spectral images being acquired at each section. All this is performed without rotating or tilting the specimen between FIB cutting and SEM imaging/x-ray spectral image acquisition. The resultant 4D spectral image is then unfolded (number of volume elements by number of channels) and subjected to the same multivariate curve resolution (MCR) approach that has proven successful for the analysis of lower-dimension x-ray spectral images. The TSI data sets can be in excess of 4Gbytes. This problem has been overcome (for now) and images up to 6Gbytes have been analyzed in this work. The method for analyzing such large spectral images will be described in this presentation. A comprehensive 3D chemical analysis was performed on several corrosion specimens of Cu electroplated with various metals. Figure 1A shows the top view of the localized corrosion region prepared for FIB sectioning. The TSI region has been coated with Pt and a trench has been milled along the bottom edge of the region, exposing it to the electron beam as seen in Figure 1B. The TSI consisted of 25 sections and was approximately 6Gbytes. Figure 1C shows several of the components rendered in 3D: Green is Cu; blue is Pb; cyan represents one of the corrosion products that contains Cu, Zn, O, S, and C; and orange represents the other corrosion product with Zn, O, S and C. Figure 1 D shows all of the component spectral shapes from the analysis. There is severe pathological overlap of the spectra from Ni, Cu and Zn as well as Pb and S. in spite of this clean spectral shapes have been extracted from the TSI. This powerful TSI technique could be applied to other sectioning methods well.

Analyzing the performance of a complex System of Systems (SoS) requires a systems engineering approach. Many such SoS exist in the Military domain. Examples include the Army's next generation Future Combat Systems 'Unit of Action' or the Navy's Aircraft Carrier Battle Group. In the case of a Unit of Action, a system of combat vehicles, support vehicles and equipment are organized in an efficient configuration that minimizes logistics footprint while still maintaining the required performance characteristics (e.g., operational availability). In this context, systems engineering means developing a global model of the entire SoS and all component systems and interrelationships. This global model supports analyses that result in an understanding of the interdependencies and emergent behaviors of the SoS. Sandia National Laboratories will present a robust toolset that includes methodologies for developing a SoS model, defining state models and simulating a system of state models over time. This toolset is currently used to perform logistics supportability and performance assessments of the set of Future Combat Systems (FCS) for the U.S. Army's Program Manager Unit of Action.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We design a density-functional-theory (DFT) exchange-correlation functional that enables an accurate treatment of systems with electronic surfaces. Surface-specific approximations for both exchange and correlation energies are developed. A subsystem functional approach is then used: an interpolation index combines the surface functional with a functional for interior regions. When the local density approximation is used in the interior, the result is a straightforward functional for use in self-consistent DFT. The functional is validated for two metals (Al, Pt) and one semiconductor (Si) by calculations of (i) established bulk properties (lattice constants and bulk moduli) and (ii) a property where surface effects exist (the vacancy formation energy). Good and coherent results indicate that this functional may serve well as a universal first choice for solid-state systems and that yet improved functionals can be constructed by this approach.

In this paper, we explore the stability properties of time-domain numerical methods for multitime partial differential equations (MPDEs) in detail. We demonstrate that simple techniques for numerical discretization can lead easily to instability. By investigating the underlying eigenstructure of several discretization techniques along different artificial time scales, we show that not all combinations of techniques are stable. We identify choices of discretization method and step size, along fast and slow time scales, that lead to robust, stable time-domain integration methods for the MPDE. One of our results is that applying overstable methods along one time-scale can compensate for unstable discretization along others. Our novel integration schemes bring robustness to time-domain MPDE solution methods, as we demonstrate with examples.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Solid-state {sup 1}H magic angle spinning (MAS) NMR was used to investigate sulfonated Diels-Alder poly(phenlylene) polymer membranes. Under high spinning speed {sup 1}H MAS conditions, the proton environments of the sulfonic acid and phenylene polymer backbone are resolved. A double-quantum (DQ) filter using the rotor-synchronized back-to-back (BABA) NMR multiple-pulse sequence allowed the selective suppression of the sulfonic proton environment in the {sup 1}H MAS NMR spectra. This DQ filter in conjunction with a spin diffusion NMR experiment was then used to measure the domain size of the sulfonic acid component within the membrane. In addition, the temperature dependence of the sulfonic acid spin-spin relaxation time (T{sub 2}) was determined, providing an estimate of the activation energy for the proton dynamics of the dehydrated membrane.

Abstract not provided.

The fundamental chemical behavior of the AlCl{sub 3}/SO{sub 2}Cl{sub 2} catholyte system was investigated using {sup 27}Al NMR spectroscopy, Raman spectroscopy, and single-crystal X-ray diffraction. Three major Al-containing species were found to be present in this catholyte system, where the ratio of each was dependent upon aging time, concentration, and/or storage temperature. The first species was identified as [Cl{sub 2}Al({mu}-Cl)]{sub 2} in equilibrium with AlCl{sub 3}. The second species results from the decomposition of SO{sub 2}Cl{sub 2} which forms Cl{sub 2}(g) and SO{sub 2}(g). The SO{sub 2}(g) is readily consumed in the presence of AlCl{sub 3} to form the crystallographically characterized species [Cl{sub 2}Al({mu}-O{sub 2}SCl)]{sub 2} (1). For 1, each Al is tetrahedrally (T{sub d}) bound by two terminal Cl and two {mu}-O ligands whereas, the S is three-coordinated by two {mu}-O ligands and one terminal Cl. The third molecular species also has T{sub d}-coordinated Al metal centers but with increased oxygen coordination. Over time it was noted that a precipitate formed from the catholyte solutions. Raman spectroscopic studies show that this gel or precipitate has a component that was consistent with thionyl chloride. We have proposed a polymerization scheme that accounts for the precipitate formation. Further NMR studies indicate that the precipitate is in equilibrium with the solution.

We demonstrate direct diode-bar side pumping of a Yb-doped fiber laser using embedded-mirror side pumping (EMSP). In this method, the pump beam is launched by reflection from a micro-mirror embedded in a channel polished into the inner cladding of a double-clad fiber (DCF). The amplifier employed an unformatted, non-lensed, ten-emitter diode bar (20 W) and glass-clad, polarization-maintaining, large-mode-area fiber. Measurements with passive fiber showed that the coupling efficiency of the raw diode-bar output into the DCF (ten launch sites) was {approx}84%; for comparison, the net coupling efficiency using a conventional, formatted, fiber-coupled diode bar is typically 50-70%, i.e., EMSP results in a factor of 2-3 less wasted pump power. The slope efficiency of the side-pumped fiber laser was {approx}80% with respect to launched pump power and 24% with respect to electrical power consumption of the diode bar; at a fiber-laser output power of 7.5 W, the EMSP diode bar consumed 41 W of electrical power (18% electrical-to-optical efficiency). When end pumped using a formatted diode bar, the fiber laser consumed 96 W at 7.5 W output power, a factor of 2.3 less efficient, and the electrical-to-optical slope efficiency was lower by a factor of 2.0. Passive-fiber measurements showed that the EMSP alignment sensitivity is nearly identical for a single emitter as for the ten-emitter bar. EMSP is the only method capable of directly launching the unformatted output of a diode bar directly into DCF (including glass-clad DCF), enabling fabrication of low-cost, simple, and compact, diode-bar-pumped fiber lasers and amplifiers.

Abstract not provided.

Spreading of bacteria in a highly advective, disordered environment is examined. Predictions of super-diffusive spreading for a simplified reaction-diffusion equation are tested. Concentration profiles display anomalous growth and super-diffusive spreading. A perturbation analysis yields a crossover time between diffusive and super-diffusive behavior. The time's dependence on the convection velocity and disorder is tested. Like the simplified equation, the full linear reaction-diffusion equation displays super-diffusive spreading perpendicular to the convection. However, for mean positive growth rates the full nonlinear reaction-diffusion equation produces symmetric spreading with a Fisher wavefront, whereas net negative growth rates cause an asymmetry, with a slower wavefront velocity perpendicular to the convection.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Chemical crosslinking is an important tool for probing protein structure and protein-protein interactions. The approach usually involves crosslinking of specific amino acids within a folded protein or protein complex, enzymatic digestion of the crosslinked protein(s), and identification of the resulting crosslinked peptides by liquid chromatography/mass spectrometry (LC/MS). In this manner, distance constraints are obtained for residues that must be in close proximity to one another in the native structure or complex. As the complexity of the system under study increases, for example, a large multi-protein complex, simply measuring the mass of a crosslinked species will not always be sufficient to determine the identity of the crosslinked peptides. In such a case, tandem mass spectrometry (MS/MS) could provide the required information if the data can be properly interpreted. In MS/MS, a species of interest is isolated in the gas phase and allowed to undergo collision induced dissociation (CID). Because the gas-phase dissociation pathways of peptides have been well studied, methods are established for determining peptide sequence by MS/MS. However, although crosslinked peptides dissociate through some of the same pathways as isolated peptides, the additional dissociation pathways available to the former have not been studied in detail. Software such as MS2Assign has been written to assist in the interpretation of MS/MS from crosslinked peptide species, but it would be greatly enhanced by a more thorough understanding of how these species dissociate. We are thus systematically investigating the dissociation pathways open to crosslinked peptide species. A series of polyalanine and polyglycine model peptides have been synthesized containing one or two lysine residues to generate defined inter- and intra-molecular crosslinked species, respectively. Each peptide contains 11 total residues, and one arginine residue is present at the carboxy terminus to mimic species generated by tryptic digestion. The peptides have been allowed to react with a series of commonly used crosslinkers such as DSS, DSG, and DST. The tandem mass spectra acquired for these crosslinked species are being examined as a function of crosslinker identity, site(s) of crosslinking, and precursor charge state. Results from these model studies and observations from actual experimental systems are being incorporated into the MS2Assign software to enhance our ability to effectively use chemical crosslinking in protein complex determination.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Surface Evolver was used to compute the equilibrium microstructure of random soap foams with bidisperse cell-size distributions and to evaluate topological and geometric properties of the foams and individual cells. The simulations agree with the experimental data of Matzke and Nestler for the probability {rho}(F) of finding cells with F faces and its dependence on the fraction of large cells. The simulations also agree with the theory for isotropic Plateau polyhedra (IPP), which describes the F-dependence of cell geometric properties, such as surface area, edge length, and mean curvature (diffusive growth rate); this is consistent with results for polydisperse foams. Cell surface areas are about 10% greater than spheres of equal volume, which leads to a simple but accurate relation for the surface free energy density of foams. The Aboav-Weaire law is not valid for bidisperse foams.

Bulk migration of particles towards regions of lower shear occurs in suspensions of neutrally buoyant spheres in Newtonian fluids undergoing creeping flow in the annular region between two rotating, coaxial cylinders (a wide-gap Couette). For a monomodal suspension of spheres in a viscous fluid, dimensional analysis indicates that the rate of migration at a given concentration should scale with the square of the sphere radius. However, a previous experimental study showed that the rate of migration of spherical particles at 50% volume concentration actually scaled with the sphere radius to approximately the 2.9 power.

Three nested molybdenum wire arrays with initial outer diameters of 45, 50, and 55 mm were imploded by the - 20 MA, 90 ns rise-time current pulse of Sandia's Z accelerator. The implosions generated Mo plasmas with {approx} 10% of the array's initial mass reaching Ne-like and nearby ionization stages. These ions emitted 2-4 keV L-shell x rays with radiative powers approaching 10 TW. Mo L-shell spectra with axial and temporal resolution were captured and have been analyzed using a collisional-radiative model. The measured spectra indicate significant axial variation in the electron density, which increases from a few times 10{sup 20} cm{sup -3} at the cathode up to - 3 x 10{sup 21} cm{sup -3} near the middle of the 20 mm plasma column (8 mm from the anode). Time-resolved spectra indicate that the peak electron density is reached before the peak of the L-shell emission and decreases with time, while the electron temperature remains within 10% of 1.7 keV over the 20-30 ns L-shell radiation pulse. Finally, while the total yield, peak total power, and peak L-shell power all tended to decrease with increasing initial wire array diameters, the L-shell yield and the average plasma conditions varied little with the initial wire array diameter.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

ODRPAC (TOMS Algorithm 676) has provided a complete package for weighted orthogonal distance regression for many years. The code is complete with user selectable reporting facilities, numerical and analytic derivatives, derivative checking, and many more features. The foundation for the algorithm is a stable and efficient trust region Levenberg-Marquardt minimizer that exploits the structure of the orthogonal distance regression problem. ODRPAC95 is a modification of the original ODRPAC code that adds support for bound constraints, uses the newer Fortran 95 language, and simplifies the interface to the user called subroutine.

We consider the accuracy of predictions made by integer programming (IP) models of sensor placement for water security applications. We have recently shown that IP models can be used to find optimal sensor placements for a variety of different performance criteria (e.g. minimize health impacts and minimize time to detection). However, these models make a variety of simplifying assumptions that might bias the final solution. We show that our IP modeling assumptions are similar to models developed for other sensor placement methodologies, and thus IP models should give similar predictions. However, this discussion highlights that there are significant differences in how temporal effects are modeled for sensor placement. We describe how these modeling assumptions can impact sensor placements.

Abstract not provided.

A focused ion beam (FIB) is used to accurately sculpt predetermined micron-scale, curved shapes in a number of solids. Using a digitally scanned ion beam system, various features are sputtered including hemispheres and sine waves having dimensions from 1-50 {micro}m. Ion sculpting is accomplished by changing pixel dwell time within individual boustrophedonic scans. The pixel dwell times used to sculpt a given shape are determined prior to milling and account for the material-specific, angle-dependent sputter yield, Y({theta}), as well as the amount of beam overlap in adjacent pixels. A number of target materials, including C, Au and Si, are accurately sculpted using this method. For several target materials, the curved feature shape closely matches the intended shape with milled feature depths within 5% of intended values.

Abstract not provided.

Abstract not provided.

Ccaffeine is a Common Component Architecture (CCA) framework devoted to high-performance computing. In this note we give an overview of the system features of Ccaffeine and CCA that support component-based HPC application development. Object-oriented, single-threaded and lightweight, Ccaffeine is designed to get completely out of the way of the running application after it has been composed from components. Ccaffeine is one of the few frameworks, CCA or otherwise, that can compose and run applications on a parallel machine interactively and then automatically generate a static, possibly self-tuning, executable for production runs. Users can experiment with and debug applications interactively, improving their productivity. When the application is ready, a script is automatically generated, parsed and turned into a static executable for production runs. Within this static executable, dynamic replacement of components can be performed by self-tuning applications.

Abstract not provided.

In recent years, several integer programming models have been proposed to place sensors in municipal water networks in order to detect intentional or accidental contamination. Although these initial models assumed that it is equally costly to place a sensor at any place in the network, there clearly are practical cost constraints that would impact a sensor placement decision. Such constraints include not only labor costs but also the general accessibility of a sensor placement location. In this paper, we extend our integer program to explicitly model the cost of sensor placement. We partition network locations into groups of varying placement cost, and we consider the public health impacts of contamination events under varying budget constraints. Thus our models permit cost/benefit analyses for differing sensor placement designs. As a control for our optimization experiments, we compare the set of sensor locations selected by the optimization models to a set of manually-selected sensor locations.

Electrical operation of III-Nitride light emitting diodes (LEDs) with photonic crystal structures is demonstrated. Employing photonic crystal structures in III-Nitride LEDs is a method to increase light extraction efficiency and directionality. The photonic crystal is a triangular lattice formed by dry etching into the III-Nitride LED. A range of lattice constants is considered (a {approx} 270-340nm). The III-Nitride LED layers include a tunnel junction providing good lateral current spreading without a semi-absorbing metal current spreader as is typically done in conventional III-Nitride LEDs. These photonic crystal III-Nitride LED structures are unique because they allow for carrier recombination and light generation proximal to the photonic crystal (light extraction area) yet displaced from the absorbing metal contact. The photonic crystal Bragg scatters what would have otherwise been guided modes out of the LED, increasing the extraction efficiency. The far-field light radiation patterns are heavily modified compared to the typical III-Nitride LED's Lambertian output. The photonic crystal affects the light propagation out of the LED surface, and the radiation pattern changes with lattice size. LEDs with photonic crystals are compared to similar III-Nitride LEDs without the photonic crystal in terms of extraction, directionality, and emission spectra.

Abstract not provided.

Abstract not provided.

Experimental evidence suggests that the energy balance between processes in play during wire array implosions is not well understood. In fact the radiative yields can exceed by several times the implosion kinetic energy. A possible explanation is that the coupling from magnetic energy to kinetic energy as magnetohydrodynamic plasma instabilities develop provides additional energy. It is thus important to model the instabilities produced in the after implosion stage of the wire array in order to determine how the stored magnetic energy can be connected with the radiative yields. To this aim three-dimensional hybrid simulations have been performed. They are initialized with plasma radial density profiles, deduced in recent experiments [C. Deeney et al., Phys. Plasmas 6, 3576 (1999)] that exhibited large x-ray yields, together with the corresponding magnetic field profiles. Unlike previous work, these profiles do not satisfy pressure balance and differ substantially from those of a Bennett equilibrium. They result in faster growth with an associated transfer of magnetic energy to plasma motion and hence kinetic energy.

Pulsed power driven metallic wire-array Z pinches are the most powerful and efficient laboratory x-ray sources. Furthermore, under certain conditions the soft x-ray energy radiated in a 5 ns pulse at stagnation can exceed the estimated kinetic energy of the radial implosion phase by a factor of 3 to 4. A theoretical model is developed here to explain this, allowing the rapid conversion of magnetic energy to a very high ion temperature plasma through the generation of fine scale, fast-growing m=0 interchange MHD instabilities at stagnation. These saturate nonlinearly and provide associated ion viscous heating. Next the ion energy is transferred by equipartition to the electrons and thus to soft x-ray radiation. Recent time-resolved iron spectra at Sandia confirm an ion temperature T{sub i} of over 200 keV (2 x 10{sup 9} degrees), as predicted by theory. These are believed to be record temperatures for a magnetically confined plasma.

Abstract not provided.

Bayesian medical monitoring is a concept based on using real-time performance-related data to make statistical predictions about a patient's future health. The following paper discusses the fundamentals behind the medical monitoring concept and the application to monitoring the health of nuclear reactors. Necessary assumptions are discussed regarding distributions and failure-rate calculations. A simple example is performed to illustrate the effectiveness of the methods. The methods perform very well for the thirteen subjects in the example, with a clear failure sequence identified for eleven of the subjects.

A novel dual stage chemiluminescence detection system incorporating individually controlled hot stages has been developed and applied to probe for material interaction effects during polymer degradation. Utilization of this system has resulted in experimental confirmation for the first time that in an oxidizing environment a degrading polymer A (in this case polypropylene, PP) is capable of infecting a different polymer B (in this case polybutadiene, HTPB) over a relatively large distance. In the presence of the infectious degrading polymer A, the thermal degradation of polymer B is observed over a significantly shorter time period. Consistent with infectious volatiles from material A initiating the degradation process in material B it was demonstrated that traces (micrograms) of a thermally sensitive peroxide in the vicinity of PP could induce degradation remotely. This observation documents cross-infectious phenomena between different polymers and has major consequences for polymer interactions, understanding fundamental degradation processes and long-term aging effects under combined material exposures.

Photocatalytic porphyrins are used to reduce metal complexes from aqueous solution and, further, to control the deposition of metals onto porphyrin nanotubes and surfactant assembly templates to produce metal composite nanostructures and nanodevices. For example, surfactant templates lead to spherical platinum dendrites and foam-like nanomaterials composed of dendritic platinum nanosheets. Porphyrin nanotubes are reported for the first time, and photocatalytic porphyrin nanotubes are shown to reduce metal complexes and deposit the metal selectively onto the inner or outer surface of the tubes, leading to nanotube-metal composite structures that are capable of hydrogen evolution and other nanodevices.

Abstract not provided.

Computation of space-charge current-limiting effects across a vacuum cavity between parallel electrodes has previously been carried out only for thermionic emission spectra. In some applications, where the current arises from an injected electron beam or photo-Compton emission from electrode walls, the electron energy spectra may deviate significantly from Maxwellian. Considering the space charge as a collisionless plasma, we derive an implicit equation for the peak cavity potential assuming steady-state currents. For the examples of graphite, nickel, and gold electrodes exposed to x rays, we find that cavity photoemission currents are typically more severely space-charge limited than they would be with the assumption of a purely Maxwellian energy distribution.

Network administrators and security analysts often do not know what network services are being run in every corner of their networks. If they do have a vague grasp of the services running on their networks, they often do not know what specific versions of those services are running. Actively scanning for services and versions does not always yield complete results, and patch and service management, therefore, suffer. We present Net-State, a system for monitoring, storing, and reporting application and operating system version information for a network. NetState gives security and network administrators the ability to know what is running on their networks while allowing for user-managed machines and complex host configurations. Our architecture uses distributed modules to collect network information and a centralized server that stores and issues reports on that collected version information. We discuss some of the challenges to building and operating NetState as well as the legal issues surrounding the promiscuous capture of network data. We conclude that this tool can solve some key problems in network management and has a wide range of possibilities for future uses.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The stainless steel alloy 17-4PH contains a martensitic microstructure and second phase delta ({delta}) ferrite. Strengthening of 17-4PH is attributed to Cu-rich precipitates produced during age hardening treatments at 900-1150 F (H900-H1150). For wrought 17-4PH, the effects of heat treatment and microstructure on mechanical properties are well-documented [for example, Ref. 1]. Fewer studies are available on cast 17-4PH, although it has been a popular casting alloy for high strength applications where moderate corrosion resistance is needed. Microstructural features and defects particular to castings may have adverse effects on properties, especially when the alloy is heat treated to high strength. The objective of this work was to outline the effects of microstructural features specific to castings, such as shrinkage/solidification porosity, on the mechanical behavior of investment cast 17-4PH. Besides heat treatment effects, the results of metallography and SEM studies showed that the largest effect on mechanical properties is from shrinkage/solidification porosity. Figure 1a shows stress-strain curves obtained from samples machined from castings in the H925 condition. The strength levels were fairly similar but the ductility varied significantly. Figure 1b shows an example of porosity on a fracture surface from a room-temperature, quasi-static tensile test. The rounded features represent the surfaces of dendrites which did not fuse or only partially fused together during solidification. Some evidence of local areas of fracture is found on some dendrite surfaces. The shrinkage pores are due to inadequate backfilling of liquid metal and simultaneous solidification shrinkage during casting. A summary of percent elongation results is displayed in Figure 2a. It was found that higher amounts of porosity generally result in lower ductility. Note that the porosity content was measured on the fracture surfaces. The results are qualitatively similar to those found by Gokhale et al. and Surappa et al. in cast A356 Al and by Gokhale et al. for a cast Mg alloys. The quantitative fractography and metallography work by Gokhale et al. illustrated the strong preference for fracture in regions of porosity in cast material. That is, the fracture process is not correlated to the average microstructure in the material but is related to the extremes in microstructure (local regions of high void content). In the present study, image analysis on random cross-sections of several heats indicated an overall porosity content of 0.03%. In contrast, the area % porosity was as high as 16% when measured on fracture surfaces of tensile specimens using stereology techniques. The results confirm that the fracture properties of cast 17-4PH cannot be predicted based on the overall 'average' porosity content in the castings.

Abstract not provided.

Abstract not provided.

Qualification Alternatives to the Sandia Pulsed Reactor (QASPR) and US DOE's Advanced Simulation and Computing (ASC).

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Policy-based network management (PBNM) uses policy-driven automation to manage complex enterprise and service provider networks. Such management is strongly supported by industry standards, state of the art technologies and vendor product offerings. We present a case for the use of PBNM and related technologies for end-to-end service delivery. We provide a definition of PBNM terms, a discussion of how such management should function and the current state of the industry. We include recommendations for continued work that would allow for PBNM to be put in place over the next five years in the unclassified environment.

ALEGRA is an arbitrary Lagrangian-Eulerian finite element code that emphasizes large distortion and shock propagation in inviscid fluids and solids. This document describes user options for modeling resistive magnetohydrodynamics, thermal conduction, and radiation transport effects, and two material temperature physics.

To establish mechanical properties and failure criteria of silicon carbide (SiC-N) ceramics, a series of quasi-static compression tests has been completed using a high-pressure vessel and a unique sample alignment jig. This report summarizes the test methods, set-up, relevant observations, and results from the constitutive experimental efforts. Results from the uniaxial and triaxial compression tests established the failure threshold for the SiC-N ceramics in terms of stress invariants (I{sub 1} and J{sub 2}) over the range 1246 < I{sub 1} < 2405. In this range, results are fitted to the following limit function (Fossum and Brannon, 2004) {radical}J{sub 2}(MPa) = a{sub 1} - a{sub 3}e -a{sub 2}(I{sub 1}/3) + a{sub 4} I{sub 1}/3, where a{sub 1} = 10181 MPa, a{sub 2} = 4.2 x 10{sup -4}, a{sub 3} = 11372 MPa, and a{sub 4} = 1.046. Combining these quasistatic triaxial compression strength measurements with existing data at higher pressures naturally results in different values for the least-squares fit to this function, appropriate over a broader pressure range. These triaxial compression tests are significant because they constitute the first successful measurements of SiC-N compressive strength under quasistatic conditions. Having an unconfined compressive strength of {approx}3800 MPa, SiC-N has been heretofore tested only under dynamic conditions to achieve a sufficiently large load to induce failure. Obtaining reliable quasi-static strength measurements has required design of a special alignment jig and load-spreader assembly, as well as redundant gages to ensure alignment. When considered in combination with existing dynamic strength measurements, these data significantly advance the characterization of pressure-dependence of strength, which is important for penetration simulations where failed regions are often at lower pressures than intact regions.

A laser hazard analysis and safety assessment was performed for the 3rd Tech model DeltaSphere-3000{reg_sign} Laser 3D Scene Digitizer, infrared laser scanner model based on the 2000 version of the American National Standard Institute's Standard Z136.1, for the Safe Use of Lasers. The portable scanner system is used in the Robotic Manufacturing Science and Engineering Laboratory (RMSEL). This scanning system had been proposed to be a demonstrator for a new application. The manufacture lists the Nominal Ocular Hazard Distance (NOHD) as less than 2 meters. It was necessary that SNL validate this NOHD prior to its use as a demonstrator involving the general public. A formal laser hazard analysis is presented for the typical mode of operation for the current configuration as well as a possible modified mode and alternative configuration.

Several SIERRA applications make use of third-party libraries to solve systems of linear and nonlinear equations, and to solve eigenproblems. The classes and interfaces in the SIERRA framework that provide linear system assembly services and access to solver libraries are collectively referred to as solver services. This paper provides an overview of SIERRA's solver services including the design goals that drove the development, and relationships and interactions among the various classes. The process of assembling and manipulating linear systems will be described, as well as access to solution methods and other operations.

The Finite Element Interface to Linear Solvers (FEI) is a linear system assembly library. Sparse systems of linear equations arise in many computational engineering applications, and the solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver package capable of solving all of the linear systems that arise. This motivates the need to switch an application from one solver library to another, depending on the problem being solved. The interfaces provided by various solver libraries for data assembly and problem solution differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application can be greatly reduced by having an abstraction layer that puts a 'common face' on various solver libraries. The FEI has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory. The original FEI offered several advantages over using linear algebra libraries directly, but also imposed significant limitations and disadvantages. A new set of interfaces has been added with the goal of removing the limitations of the original FEI while maintaining and extending its strengths.

This report summarizes the research and development performed from October 2002 to September 2004 at Sandia National Laboratories under the Laboratory-Directed Research and Development (LDRD) project ''Massively-Parallel Linear Programming''. We developed a linear programming (LP) solver designed to use a large number of processors. LP is the optimization of a linear objective function subject to linear constraints. Companies and universities have expended huge efforts over decades to produce fast, stable serial LP solvers. Previous parallel codes run on shared-memory systems and have little or no distribution of the constraint matrix. We have seen no reports of general LP solver runs on large numbers of processors. Our parallel LP code is based on an efficient serial implementation of Mehrotra's interior-point predictor-corrector algorithm (PCx). The computational core of this algorithm is the assembly and solution of a sparse linear system. We have substantially rewritten the PCx code and based it on Trilinos, the parallel linear algebra library developed at Sandia. Our interior-point method can use either direct or iterative solvers for the linear system. To achieve a good parallel data distribution of the constraint matrix, we use a (pre-release) version of a hypergraph partitioner from the Zoltan partitioning library. We describe the design and implementation of our new LP solver called parPCx and give preliminary computational results. We summarize a number of issues related to efficient parallel solution of LPs with interior-point methods including data distribution, numerical stability, and solving the core linear system using both direct and iterative methods. We describe a number of applications of LP specific to US Department of Energy mission areas and we summarize our efforts to integrate parPCx (and parallel LP solvers in general) into Sandia's massively-parallel integer programming solver PICO (Parallel Interger and Combinatorial Optimizer). We conclude with directions for long-term future algorithmic research and for near-term development that could improve the performance of parPCx.

The manipulation of physical interactions between structural moieties on the molecular scale is a fundamental hurdle in the realization and operation of nanostructured materials and high surface area microsystem architectures. These include such nano-interaction-based phenomena as self-assembly, fluid flow, and interfacial tribology. The proposed research utilizes photosensitive molecular structures to tune such interactions reversibly. This new material strategy provides optical actuation of nano-interactions impacting behavior on both the nano- and macroscales and with potential to impact directed nanostructure formation, microfluidic rheology, and tribological control.

As part of the U.S. Department of Energy (DOE) Office of Industrial Technologies (OIT) Industries of the Future (IOF) Forest Products research program, a collaborative investigation was conducted on the sources, characteristics, and deposition of particles intermediate in size between submicron fume and carryover in recovery boilers. Laboratory experiments on suspended-drop combustion of black liquor and on black liquor char bed combustion demonstrated that both processes generate intermediate size particles (ISP), amounting to 0.5-2% of the black liquor dry solids mass (BLS). Measurements in two U.S. recovery boilers show variable loadings of ISP in the upper furnace, typically between 0.6-3 g/Nm{sup 3}, or 0.3-1.5% of BLS. The measurements show that the ISP mass size distribution increases with size from 5-100 {micro}m, implying that a substantial amount of ISP inertially deposits on steam tubes. ISP particles are depleted in potassium, chlorine, and sulfur relative to the fuel composition. Comprehensive boiler modeling demonstrates that ISP concentrations are substantially overpredicted when using a previously developed algorithm for ISP generation. Equilibrium calculations suggest that alkali carbonate decomposition occurs at intermediate heights in the furnace and may lead to partial destruction of ISP particles formed lower in the furnace. ISP deposition is predicted to occur in the superheater sections, at temperatures greater than 750 C, when the particles are at least partially molten.

The thermal hazard posed by large hydrocarbon fires is dominated by the radiative emission from high temperature soot. Since the optical properties of soot, especially in the infrared region of the electromagnetic spectrum, as well as its morphological properties, are not well known, efforts are underway to characterize these properties. Measurements of these soot properties in large fires are important for heat transfer calculations, for interpretation of laser-based diagnostics, and for developing soot property models for fire field models. This research uses extractive measurement diagnostics to characterize soot optical properties, morphology, and composition in 2 m pool fires. For measurement of the extinction coefficient, soot extracted from the flame zone is transported to a transmission cell where measurements are made using both visible and infrared lasers. Soot morphological properties are obtained by analysis via transmission electron microscopy of soot samples obtained thermophoretically within the flame zone, in the overfire region, and in the transmission cell. Soot composition, including carbon-to-hydrogen ratio and polycyclic aromatic hydrocarbon concentration, is obtained by analysis of soot collected on filters. Average dimensionless extinction coefficients of 8.4 {+-} 1.2 at 635 nm and 8.7 {+-} 1.1 at 1310 nm agree well with recent measurements in the overfire region of JP-8 and other fuels in lab-scale burners and fires. Average soot primary particle diameters, radius of gyration, and fractal dimensions agree with these recent studies. Rayleigh-Debye-Gans theory of scattering applied to the measured fractal parameters shows qualitative agreement with the trends in measured dimensionless extinction coefficients. Results of the density and chemistry are detailed in the report.

The measurement of heat flux in hydrocarbon fuel fires (e.g., diesel or JP-8) is difficult due to high temperatures and the sooty environment. Un-cooled commercially available heat flux gages do not survive in long duration fires, and cooled gages often become covered with soot, thus changing the gage calibration. An alternate method that is rugged and relatively inexpensive is based on inverse heat conduction methods. Inverse heat-conduction methods estimate absorbed heat flux at specific material interfaces using temperature/time histories, boundary conditions, material properties, and usually an assumption of one-dimensional (1-D) heat flow. This method is commonly used at Sandia.s fire test facilities. In this report, an uncertainty analysis was performed for a specific example to quantify the effect of input parameter variations on the estimated heat flux when using the inverse heat conduction method. The approach used was to compare results from a number of cases using modified inputs to a base-case. The response of a 304 stainless-steel cylinder [about 30.5 cm (12-in.) in diameter and 0.32-cm-thick (1/8-in.)] filled with 2.5-cm-thick (1-in.) ceramic fiber insulation was examined. Input parameters of an inverse heat conduction program varied were steel-wall thickness, thermal conductivity, and volumetric heat capacity; insulation thickness, thermal conductivity, and volumetric heat capacity, temperature uncertainty, boundary conditions, temperature sampling period; and numerical inputs. One-dimensional heat transfer was assumed in all cases. Results of the analysis show that, at the maximum heat flux, the most important parameters were temperature uncertainty, steel thickness and steel volumetric heat capacity. The use of a constant thermal properties rather than temperature dependent values also made a significant difference in the resultant heat flux; therefore, temperature-dependent values should be used. As an example, several parameters were varied to estimate the uncertainty in heat flux. The result was 15-19% uncertainty to 95% confidence at the highest flux, neglecting multidimensional effects.

Abstract not provided.