We have developed algorithms to automatically learn a detection map of a deployed sensor field for a virtual presence and extended defense (VPED) system without apriori knowledge of the local terrain. The VPED system is an unattended network of sensor pods, with each pod containing acoustic and seismic sensors. Each pod has the ability to detect and classify moving targets at a limited range. By using a network of pods we can form a virtual perimeter with each pod responsible for a certain section of the perimeter. The site's geography and soil conditions can affect the detection performance of the pods. Thus, a network in the field may not have the same performance as a network designed in the lab. To solve this problem we automatically estimate a network's detection performance as it is being installed at a site by a mobile deployment unit (MDU). The MDU will wear a GPS unit, so the system not only knows when it can detect the MDU, but also the MDU's location. In this paper, we demonstrate how to handle anisotropic sensor-configurations, geography, and soil conditions.
In this paper we will demonstrate that the computational effort of FWI can be reduced significantly by applying it to data formed by encoding and summing source gathers, if the encoding of the sources is changed between iterations. Changing the encoding between iterations changes the crosstalk noise caused by the summation of the sources. Thus, the source crosstalk-noise stacks out of the inverted earth model, allowing summation of a large number of encoded sources. We call this method encoded simultaneous-source FWI (ESSFWI).
The U.S. Department of Energy (DOE) provides scientific infrastructure and data archives to the international Arctic research community through a national user facility, the ARM Climate Research Facility, located on the North Slope of Alaska. The ARM sites at Barrow and Atqasuk, Alaska have been collecting and archiving atmospheric data for more than 10 years. These data have been used for scientific investigation as well as remote sensing validations. Funding from the Recovery Act (American Recovery and Reinvestment Act of 2009) will be used to install new instruments and upgrade existing instruments at the North Slope sites. These instruments include: scanning precipitation radar; scanning cloud radar; automatic balloon launcher; high spectral resolution lidar; eddy correlation flux systems; and upgraded ceilometer, AERI, micropulse lidar, and millimeter cloud radar. Information on these planned additions and upgrades will be provided in our poster. An update on activities planned at Oliktok Point will also be provided.
The objective of this Standard is the specification of a verification and validation approach that quantifies the degree of accuracy inferred from the comparison of solution and data for a specified variable at a specified validation point. The approach uses the concepts from experimental uncertainty analysis to consider the errors and uncertainties in both the solution and the data. The scope of this Standard is the quantification of the degree of accuracy of simulation of specified validation variables at a specified validation point for cases in which the conditions of the actual experiment are simulated. Consideration of solution accuracy at points within a domain other than the validation points, for example interpolation/extrapolation in a domain of validation, is a matter of engineering judgment specific to each family of problems and is beyond the scope of this Standard.
The purpose of this work is to improve detection methods that can reliably identify special nuclear material (SNM). One method that can be used to identify special nuclear material is neutron multiplicity analysis. This method detects multiple time-correlated neutrons released from a fission event in the SNM. This work investigates the ability of the software code MCNP-PoliMi to simulate neutron multiplicity measurements from a highly moderated SNM source. A measurement of a 4.5-kg alpha-phase metal plutonium sphere surrounded by up to 6 inches of polyethylene shells has recently been performed by Sandia National Laboratories personnel at the Nevada Test Site. A post-processing code was developed to account for dead-time effects within the detector and to determine the neutron multiplicity distributions for various time intervals. With the distributions calculated, the Feynman-Y can be determined. The Feynman-Y is a metric that measures the level of correlation present in a sample. At this time MCNP-PoliMi is able predict the Feynman-Y within 10% of the measured value.
This presentation discusses the following topics: (1) Red Sky Background; (2) 3D Torus Interconnect Concepts; (3) Difficulties of Torus in IB; (4) New Routing Code for IB a 3D Torus; (5) Red Sky 3D Torus Implementation; and (6) Managing a Large IB Machine. Computing at Sandia: (1) Capability Computing - Designed for scaling of single large runs, Usually proprietary for maximum performance, and Red Storm is Sandia's current capability machine; (2) Capacity Computing - Computing for the masses, 100s of jobs and 100s of users, Extreme reliability required, Flexibility for changing workload, Thunderbird will be decommissioned this quarter, Red Sky is our future capacity computing platform, and Red Mesa machine for National Renewable Energy Lab. Red Sky main themes are: (1) Cheaper - 5X capacity of Tbird at 2/3 the cost, Substantially cheaper per flop than our last large capacity machine purchase; (2) Leaner - Lower operational costs, Three security environments via modular fabric, Expandable, upgradeable, extensible, and Designed for 6yr. life cycle; and (3) Greener - 15% less power-1/6th power per flop, 40% less water-5M gallons saved annually, 10X better cooling efficiency, and 4x denser footprint.
The readout of a solid state qubit often relies on single charge sensitive electrometry. However the combination of fast and accurate measurements is non trivial due to large RC time constants due to the electrometers resistance and shunt capacitance from wires between the cold stage and room temperature. Currently fast sensitive measurements are accomplished through rf reflectrometry. I will present an alternative single charge readout technique based on cryogenic CMOS circuits in hopes to improve speed, signal-to-noise, power consumption and simplicity in implementation. The readout circuit is based on a current comparator where changes in current from an electrometer will trigger a digital output. These circuits were fabricated using Sandia's 0.35 {micro}m CMOS foundry process. Initial measurements of comparators with an addition a current amplifier have displayed current sensitivities of < 1nA at 4.2K, switching speeds up to {approx}120ns, while consuming {approx}10 {micro}W. I will also discuss an investigation of noise characterization of our CMOS process in hopes to obtain a better understanding of the ultimate limit in signal to noise performance.
Constructing high-fidelity control pulses that are robust to control and system/environment fluctuations is a crucial objective for quantum information processing (QIP). We combine dynamical decoupling (DD) with optimal control (OC) to identify control pulses that achieve this objective numerically. Previous DD work has shown that general errors up to (but not including) third order can be removed from {pi}- and {pi}/2-pulses without concatenation. By systematically integrating DD and OC, we are able to increase pulse fidelity beyond this limit. Our hybrid method of quantum control incorporates a newly-developed algorithm for robust OC, providing a nested DD-OC approach to generate robust controls. Motivated by solid-state QIP, we also incorporate relevant experimental constraints into this DD-OC formalism. To demonstrate the advantage of our approach, the resulting quantum controls are compared to previous DD results in open and uncertain model systems.
Growth of high quality graphene films on SiC is regarded as one of the more viable pathways toward graphene-based electronics. Graphitic films form on SiC at elevated temperature because of preferential sublimation of Si. Little is known, however, about the atomistic processes of interrelated SiC decomposition and graphene growth. We have observed the formation of graphene on SiC by Si sublimation in an Ar atmosphere using low energy electron microscopy, scanning tunneling microcopy and atomic force microscopy. This work reveals that the growth mechanism depends strongly on the initial surface morphology, and that carbon diffusion governs the spatial relationship between SiC decomposition and graphene growth. Isolated bilayer SiC steps generate narrow ribbons of graphene, whereas triple bilayer steps allow large graphene sheets to grow by step flow. We demonstrate how graphene quality can be improved by controlling the initial surface morphology specifically by avoiding the instabilities inherent in diffusion-limited growth.
Establishing the atomic structure and composition of interfaces in thermoelectric materials is important to understanding how these defects mediate thermal and electronic transport. Here, we discuss our experimental observations and theoretical calculations of the Bi{sub 2}Te{sub 3} (0001) basal twin in nanocrystalline Bi{sub 2}Te{sub 3}. This interface is important both because it is common in tetradymite-structured thermoelectric compounds and because it serves as a useful model system for more complex interfaces. Macroscopically, the (0001) twin corresponds to a 180 rotation of the crystal about the [0001] axis, which reverses the stacking of the basal planes. The basal planes of Bi{sub 2}Te{sub 3} are arranged in 5-plane groupings of alternating Bi and Te layers. Microscopically, one envisions three possible interface terminations: at the Te layer in the middle of the 5-layer packet, at a Bi layer, or at the Te-double layer at the junction of the 5-layer packet. Using aberration-corrected HAADF-STEM imaging, we have established that the twin boundary terminates at the Te-double layer. This result is consistent with ab initio calculations, which predict that the lowest energy for the three candidate structures is for this termination.
Electrons on the surface of superfluid helium have extremely high mobilities and long predicted spin coherence times, making them ideal mobile qubits. Previous work has shown that electrons localized in helium filled channels can be reliably transported between multiple underlying gates. Silicon chips have been designed, fabricated, and post processed by reactive ion etching to leverage the large scale integration capabilities of silicon technology. These chips, which serve as substrates for the electrons on helium research, utilize silicon CMOS for on-chip signal amplification and multiplexing and the uppermost metal layers for defining the helium channels and applying electrical potentials for moving the electrons. We will discuss experimental results for on-chip circuitry and clocked electron transport along etched channels.
A light beam induced current (LBIC) measurement is a non-destructive technique that produces a spatial graphical representation of current response in photovoltaic cells with respect to position when stimulated by a light beam. Generally, a laser beam is used for these measurements because the spot size can be made very small, on the order of microns, and very precise measurements can be made. Sandia National Laboratories Photovoltaic System Evaluation Laboratory (PSEL) uses its LBIC measurement technique to characterize single junction mono-crystalline and multi-crystalline solar cells ranging from miniature to conventional sizes. Sandia has modified the already valuable LBIC technique to enable multi-junction PV cells to be characterized.
Nanoporous metallic particles are of great interest for a range of applications including catalysis, gas storage, and electrical energy storage. In particular, recent work has shown that bulk powders of porous palladium can be synthesized in a scalable fashion. This material has pore sizes in the 2-5 nm range and has promise for use in hydrogen storage applications. However, because of the small pore size such materials are very susceptible to morphological evolution during aging, especially at elevated temperatures, leading to degradation of their storage properties. To better understand and predict the phenomena at work in nanoporous metal aging, we have developed a kinetic Monte Carlo (kMC) model for the simulation of atomic diffusion in a Pd lattice. The model is implemented in Sandia's parallelized kMC code SPPARKS. SPPARKS utilizes a spatial decomposition parallelization scheme, allowing large-scale simulations including millions of atoms. The diffusion model includes single-atom hops as well as Schwoebel barrier events that mimic concerted atom motions involving multiple lattice sites. Our simulations show that for statistically homogeneous nanoporous networks, coarsening at elevated temperature as measured by the surface area can be described by a scaling law that closely follows the L {approx} {sup 1/4} scaling predicted by continuum surface diffusion theory. This scaling holds despite the presence of surface faceting due to our simulations being run at temperatures below the roughening temperature of the material. Sensitivities of the rate of coarsening, the scaling exponent, and the amount of surface faceting to model parameters including temperature and event activation rates are explored. Because of the large spatial scales attainable in our computations, we are able to simulate nanoporous particle geometries similar to those synthesized in the laboratory, and compare directly to material aging experiments including porosimetry measurements and TEM images of particles.
JMP and design of experiments (DOE) have been successfully applied to security system technologies from sensors to communication and display systems. In all cases, the technologies have been complex enough to warrant the need for a statistical determination of significant factors and/or the generation of predictive models. For the sensors, it was the task of calibrating a fiber optic intrusion detection sensor (FOIDS) with 32 adjustable settings. In addition to the numerous settings, the FOIDS also had two software processors for detecting different types of alarms. The problem was made more complex when the different types of alarms occurred on the wrong processors, causing nuisance alarms. JMP's ability to optimize several predictive models simultaneously with JMP's Prediction Profiler flash files was an important factor in producing field solutions. For the Communications and Display testbed system, numerous hardware and software network components had been integrated to build a functional system. Although the components of the system had been tested individually, the system's performance could not be piecewise evaluated. Through the application of JMP's design of experiments and data mining capabilities, it was possible to test some of the factors affecting the system's performance and to differentiate between some of the software and hardware contributors. This paper will discuss design of experiments and the JMP tools applied to the solutions for both security systems.
Significant progress has been made over the last few years in understanding properties of matter subject to strong shocks and other extreme conditions. High-accuracy multi-Mbar experiments and first-principles theoretical studies together provide detailed insights into the physics and chemistry of high energy-density matter. While comprehensive advances have been made for pure elements like deuterium, helium, and carbon, progress has been slower for equally important, albeit more challenging, materials like molecular crystals, polymers, and foams. Hydrocarbon based polymer foams are common materials and in particular they are used in designing shock- and inertial confinement fusion experiments. Depending on their initial density, foams shock to relatively higher pressure and temperature compared to shocked dense polymers/plastics. As foams and polymers are shocked, they exhibit both structural and chemical transitions. We will present experimental and theoretical results for shocked polymers in the Mbar regime. By shock impact of magnetically launched flyer plates on poly(4-methyl-1-pentene) foams, we create multi-Mbar pressures in a dense plasma mixture of hydrogen, carbon, at temperatures of several eV. Concurrently with executing experiments, we analyze the system by multi-scale simulations, from density functional theory to continuum magneto-hydrodynamics simulations. In particular, density functional theory (DFT) molecular dynamics (MD) and classical MD simulations of the principal shock Hugoniot will be presented in detail for two hydrocarbon polymers: polyethylene (PE) and poly(4-methyl-1-pentene) (PMP).
The political dynamics associated with an election are typically a function of the interplay between political leaders and voters, as well as endogenous and exogenous factors that impact the perceptions and goals of the electorate. This paper describes an effort by Sandia National Laboratories to model the attitudes and behaviors of various political groups along with that population's primary influencers, such as government leaders. To accomplish this, Sandia National Laboratories is creating a hybrid system dynamics-cognitive model to simulate systems- and individual-level political dynamics in a hypothetical society. The model is based on well-established psychological theory, applied to both individuals and groups within the modeled society. Confidence management processes are being incorporated into the model design process to increase the utility of the tool and assess its performance. This project will enhance understanding of how political dynamics are determined in democratic society.
The Office of Nuclear Regulatory Research (RES) at the US Nuclear Regulatory Commission (USNRC) is sponsoring work in response to a Staff Requirements Memorandum (SRM) directing an effort to establish a single human reliability analysis (HRA) method for the agency or guidance for the use of multiple methods. As part of this effort an attempt to develop a comprehensive HRA qualitative approach is being pursued. This paper presents a draft of the method's middle layer, a part of the qualitative analysis phase that links failure mechanisms to performance shaping factors. Starting with a Crew Response Tree (CRT) that has identified human failure events, analysts identify potential failure mechanisms using the mid-layer model. The mid-layer model presented in this paper traces the identification of the failure mechanisms using the Information-Diagnosis/Decision-Action (IDA) model and cognitive models from the psychological literature. Each failure mechanism is grouped according to a phase of IDA. Under each phase of IDA, the cognitive models help identify the relevant performance shaping factors for the failure mechanism. The use of IDA and cognitive models can be traced through fault trees, which provide a detailed complement to the CRT.
A low-temperature upturn of the Coulomb drag resistivity {rho}{sub D} measured in undoped electron-hole bilayer devices, possibly manifesting from formation of a superfluid condensate or density modulated state, was recently observed. Here the effects of perpendicular and parallel magnetic fields on the drag upturn are examined. Measurements of {rho}{sub D} and drive layer resistivity {rho}{sub xx-e} as a function of temperature and magnetic field in two uEHBL devices are presented. In B{sub {perpendicular}}, the drag upturn was enhanced as the field increased up to roughly .2 T, beyond which oscillations in {rho}{sub D} and {rho}{sub xx-e}, reflecting Landau level formation, begin appearing. A small phase offset between those oscillations, which decreased at higher fields and temperatures, was also observed. In B{sub {parallel}}, the drag upturn magnitude diminished as the field increased. Above the upturn regime, both {rho}{sub D} and {rho}{sub xx-e} were enhanced by B{sub {parallel}}, the latter via decreased screening of the uniform background impurities.
Silicon is an ideal system for investigating single electron or isolated donor spins for quantum computation, due to long spin coherence times. Enhancement mode strained-silicon/silicon germanium (sSi/SiGe) devices would offer an as-yet untried path toward electron or electron/donor quantum dot systems. Thin, undoped SiGe dielectrics allow tight electrostatic confinement, as well as potential Lande g-factor engineered spin manipulation. In this talk we summarize recent progress toward sSi/SiGe enhancement mode devices on sSi on insulator, including characterization with X-ray diffraction and atomic force microscopy, as well as challenges faced and progress on integration of either top-down and bottom-up donor placement approaches in a sSi/SiGe enhancement mode structure.
The principal design drivers in the certification of wind turbine blades are ultimate strength, fatigue resistance, adequate tip-tower clearance, and buckling resistance. Buckling resistance is typically strongly correlated to both ultimate strength and fatigue resistance. A composite shell with spar caps forms the airfoil shape of a blade and reinforcing shear webs are placed inside the blade to stiffen the blade in the flap-wise direction. The spar caps are dimensioned and the shear webs are placed so as to add stiffness to unsupported panel regions and reduce their length. The panels are not the major flap-wise load carrying element of a blade; however, they must be designed carefully to avoid buckling while minimizing blade weight. Typically, buckling resistance is evaluated by consideration of the load-deflection behavior of a blade using finite element analysis (FEA) or full-scale static testing of blades under a simulated extreme loading condition. The focus of this paper is on the use of experimental modal analysis to measure localized resonances of the blade panels. It can be shown that the resonant behavior of these panels can also provide a means to evaluate buckling resistance by means of analytical or experimental modal analysis. Further, panel resonances have use in structural health monitoring by observing changes in modal parameters associated with panel resonances, and use in improving panel laminate model parameters by correlation with test data. In recent modal testing of wind turbine blades, a set of panel modes were measured. This paper will report on the findings of these tests and accompanying numerical and analytical modeling efforts aimed at investigating the potential uses of panel resonances for blade evaluation, health monitoring, and design.
When elastic-plastic materials, such as metals, are subjected to moderately high strain rates or dynamic loadings, the plastic stress wave trails behind the elastic wave because of its slower wave speed. Due to the inherent time-dependent nature of the plastic deformation, the elastic precursor generally loads the material to a metastable elastic state at a stres level that is higher than the static strength of the material. This metastable state gradually relaxes to the equilibrium state and the relaxation results in the so-called precursor decay. In a recent work by Asay et al. (J. Appl. Phys., 2009), the inelastic response of annealed and cold-rolled pure polycrystalline tantalum at intermediate strain rates ({approx} 106/sec) was experimentally characterized with ramp wave loading. It was found that the precursor of the annealed tantalum showed little decay over a propagation distance of 6 mm even though the peak precursor stress was well above the static strength of the mateiral. The precursor for the cold-rolled sample was more dispersive and did not exhibit the characteristics depicted by the annealed samples. In this study, a constitutive model based on the concept of dislocation motion and generation was developed to gain insights into this somewhat unusual precursor behavior, particularly for the annealed samples, and the possible underlying deformation mechanisms for tantalum. Despite its simplicity, the model worked quite well for both the annealed and cold-rolled materials. The tantalum studied here essentially exhibits strong rate sensitivity and this behavior is modeled through the low dislocation density and the strong stress dependence of the dislocation velocity. Both of these contributions may be related to the low mobility of the screw dislocations in bcc metals. This low mobility may result from its extended, three-dimensional core structure.
An overview of the study of the effects that electrical power and signal cables introduce on the dynamic response of precision structures is presented, along with a summary of lessons learned and most significant results. This was a three-year effort conducted at the Air Force Research Laboratory, Space Vehicles Directorate to discover a set of practical approaches for updating well defined dynamical models of cableless structures where knowledge of the cable type, position, and tie-down method are known. While cables can be found on many different types of structures, the focus of this effort was on precision, low-damping, and low-first modal frequency structures. Various obstacles, classified as tangents, rat holes, and dead ends, were encountered along the way. Rather than following a strictly technical flow, the paper presents the historical, experiential progression of the project. First, methods were developed to estimate cable properties. Problems were encountered because of the flexible, highly damped nature of cables. A simple beam was used as a test article to validate experimentally derived cable properties and to refine the assumptions regarding boundary conditions. A spacecraft bus-like panel with cables attached was designed, and finite element models were developed and validated through experiment. Various paths were investigated at each stage before a consistent test and analysis methodology was developed. These twists and turns are described.
The increased demand for Liquefied Natural Gas (LNG) as a fuel source in the U.S. has prompted a study to improve our capability to predict cascading damage to LNG tankers from cryogenic spills and subsequent fire. To support this large modeling and simulation effort, a suite of experiments were conducted on two tanker steels, ABS Grade A steel and ABS Grade EH steel. A thorough and complete understanding of the mechanical behavior of the tanker steels was developed that was heretofore unavailable for the span of temperatures of interest encompassing cryogenic to fire temperatures. This was accomplished by conducting several types of experiments, including tension, notched tension and Charpy impact tests at fourteen temperatures over the range of -191 C to 800 C. Several custom fixtures and special techniques were developed for testing at the various temperatures. The experimental techniques developed and the resulting data will be presented, along with a complete description of the material behavior over the temperature span.
We develop a switched feedback controller that optimizes the rate of convergence of the state trajectories to the origin for a class of second order LTI systems. Specifically, we derive an algorithm which optimizes the rate of convergence by employing a controller that switches between symmetric gains. As a byproduct of our investigation, we find that, in general, the controllers which optimize the rate of convergence switch between two linear subsystems, one of which is unstable. The algorithm we investigate will design optimal switching laws for the specific case of second order LTI plants of relative degree two.
Instrumented, fully coupled thermal-mechanical experiments were conducted to provide validation data for finite element simulations of failure in pressurized, high temperature systems. The design and implementation of the experimental methodology is described in another paper of this conference. Experimental coupling was accomplished on tubular 304L stainless steel specimens by mechanical loading imparted by internal pressurization and thermal loading by side radiant heating. Experimental parameters, including temperature and pressurization ramp rates, maximum temperature and pressure, phasing of the thermal and mechanical loading and specimen geometry details were studied. Experiments were conducted to increasing degrees of deformation, up to and including failure. Mechanical characterization experiments of the 304L stainless steel tube material was also completed for development of a thermal elastic-plastic material constitutive model used in the finite element simulations of the validation experiments. The material was characterized in tension at a strain rate of 0.001/s from room temperature to 800 C. The tensile behavior of the tube material was found to differ substantially from 304L bar stock material, with the plasticity characteristics and strain to failure differing at every test temperature.
Coupled thermal-mechanical experiments with well-defined, controlled boundary conditions were designed through an iterative process involving a team of experimentalists, material modelers and computational analysts. First the basic experimental premise was selected: an axisymmetric tubular specimen mechanically loaded by internal pressurization and thermally loaded asymmetrically by side radiant heating. Then several integrated experimental-analytical steps were taken to determine the experimental details. The boundary conditions were mostly thermally driven and were chosen so they could be modeled accurately; the experimental fixtures were designed to ensure that the boundary conditions were met. Preliminary, uncoupled analyses were used to size the specimen diameter, height and thickness with experimental consideration of maximum pressure loads and fixture design. Iterations of analyses and experiments were used to efficiently determine heating parameters including lamp and heating shroud design, set off distance between the lamps and shroud and between the shroud and specimen, obtainable ramp rates, and the number and spatial placement of thermocouples. The design process and the experimental implementation of the final coupled thermomechanical failure experiment design will be presented.
The blades of a modern wind turbine are critical components central to capturing and transmitting most of the load experienced by the system. They are complex structural items composed of many layers of fiber and resin composite material and typically, one or more shear webs. Large turbine blades being developed today are beyond the point of effective trial-and-error design of the past and design for reliability is always extremely important. Section analysis tools are used to reduce the three-dimensional continuum blade structure to a simpler beam representation for use in system response calculations to support full system design and certification. One model simplification approach is to analyze the two-dimensional blade cross sections to determine the properties for the beam. Another technique is to determine beam properties using static deflections of a full three-dimensional finite element model of a blade. This paper provides insight into discrepancies observed in outputs from each approach. Simple two-dimensional geometries and three-dimensional blade models are analyzed in this investigation. Finally, a subset of computational and experimental section properties for a full turbine blade are compared.
Resonant plasmonic detectors are potentially important for terahertz (THz) spectroscopic imaging. We have fabricated and characterized antenna coupled detectors that integrate a broad-band antenna, which improves coupling of THz radiation. The vertex of the antenna contains the tuning gates and the bolometric barrier gate. Incident THz radiation may excite 2D plasmons with wave-vectors defined by either a periodic grating gate or a plasmonic cavity determined by ohmic contacts and gate terminals. The latter approach of exciting plasmons in a cavity defined by a short micron-scale channel appears most promising. With this short-channel geometry, we have observed multiple harmonics of THz plasmons. At 20 K with detector bias optimized we report responsivity on resonance of 2.5 kV/W and an NEP of 5 x 10{sup -10} W/Hz{sup 1/2}.
Separation distances are used in hydrogen refueling stations to protect people, structures, and equipment from the consequences of accidental hydrogen releases. Specifically, hydrogen jet flames resulting from ignition of unintended releases can be extensive in length and pose significant radiation and impingement hazards. Depending on the leak diameter and source pressure, the resulting separation distances can be unacceptably large. One possible mitigation strategy to reduce exposure to hydrogen flames is to incorporate barriers around hydrogen storage, process piping, and delivery equipment. The effectiveness of barrier walls to reduce hazards at hydrogen facilities has been previously evaluated using experimental and modeling information developed at Sandia National Laboratories. The effect of barriers on the risk from different types of hazards including direct flame contact, radiation heat fluxes, and overpressures associated with delayed hydrogen ignition has subsequently been evaluated and used to identify potential reductions in separation distances in hydrogen facilities. Both the frequency and consequences used in this risk assessment and the risk results are described. The results of the barrier risk analysis can also be used to help establish risk-informed barrier design requirements for use in hydrogen codes and standards.
This study investigates a pathway to nanoporous structures created by hydrogen and helium implantation in aluminum. Previous experiments for fusion applications have indicated that hydrogen and helium ion implantations are capable of producing bicontinuous nanoporous structures in a variety of metals. This study focuses specifically on implantations of hydrogen and helium ions at 25 keV in aluminum. The hydrogen and helium systems result in remarkably different nanostructures of aluminum at the surface. Scanning electron microscopy, focused ion beam, and transmission electron microscopy show that both implantations result in porosity that persists approximately 200 nm deep. However, hydrogen implantations tend to produce larger and more irregular voids that preferentially reside at defects. Implantations of helium at a fluence of 10{sup 18} cm{sup -2} produce much smaller porosity on the order of 10 nm that is regular and creates a bicontinuous structure in the porous region. The primary difference driving the formation of the contrasting structures is likely the relatively high mobility of hydrogen and the ability of hydrogen to form alanes that are capable of desorbing and etching Al (111) faces.
There is an increasing need to assess the performance of high consequence systems using a modeling and simulation based approach. Central to this approach are the need to quantify the uncertainties present in the system and to compare the system response to an expected performance measure. At Sandia National Laboratories, this process is referred to as quantification of margins and uncertainties or QMU. Depending on the outcome of the assessment, there might be a need to increase the confidence in the predicted response of a system model; thus a need to understand where resources need to be allocated to increase this confidence. This paper examines the problem of resource allocation done within the context of QMU. An optimization based approach to solving the resource allocation is considered and sources of aleatoric and epistemic uncertainty are included in the calculations.
The relatively recent development of short (nsec) and ultra-short (fsec) pulsed laser systems has introduced process capabilities which are particularly suited for micro-manufacturing applications. Micrometer feature resolutions and minimal heat affected zones are commonly cited benefits, although unique material interactions also prove attractive for many applications. A background of short and ultra-short pulsed laser system capabilities and material interactions will be presented for micro-scale processing. Processing strengths and limitations will be discussed and demonstrated within the framework of applications related to micro-machining, material surface modifications, and fundamental material science research.
The amounts of charge collection by single-photon absorption to that by two-photon absorption laser testing techniques have been directly compared using specially made SOI diodes. Details of this comparison are discussed.
LAMMPS is a classical molecular dynamics code, and an acronym for Large-scale Atomic/Molecular Massively Parallel Simulator. LAMMPS has potentials for soft materials (biomolecules, polymers) and solid-state materials (metals, semiconductors) and coarse-grained or mesoscopic systems. It can be used to model atoms or, more generically, as a parallel particle simulator at the atomic, meso, or continuum scale. LAMMPS runs on single processors or in parallel using message-passing techniques and a spatial-decomposition of the simulation domain. The code is designed to be easy to modify or extend with new functionality.
We report reflectivity, design and laser damage comparisons of our AR coatings for use at 1054 nm and/or 527 nm, and at angles of incidence between 0 and 45 degrees.
Axial Ge/Si heterostructure nanowires (NWs) allow energy band-edge engineering along the axis of the NW, which is the charge transport direction, and the realization of asymmetric devices for novel device architectures. This work reports on two significant advances in the area of heterostructure NWs and tunnel FETs: (i) the realization of 100% compositionally modulated Si/Ge axial heterostructure NWs with lengths suitable for device fabrication and (ii) the design and implementation of Schottky barrier tunnel FETs on these NWs for high-on currents and suppressed ambipolar behavior. Initial prototype devices with 10 nm PECVD SiN{sub x} gate dielectric resulted in a very high current drive in excess of 100 {micro}A/{micro}m (I/{pi}D) and 10{sup 5} I{sub on}/I{sub off} ratios. Prior work on the synthesis of Ge/Si axial NW heterostructures through the VLS mechanism have resulted in axial Si/Si{sub 1-x}Ge{sub x} NW heterostructures with x{sub max} {approx} 0.3, and more recently 100% composition modulation was achieved with a solid growth catalyst. In this latter case, the thickness of the heterostructure cannot exceed few atomic layers due to the slow axial growth rate and concurrent radial deposition on the NW sidewalls leading to a mixture of axial and radial deposition, which imposes a big challenge for fabricating useful devices form these NWs in the near future. Here, we report the VLS growth of 100% doping and composition modulated axial Ge/Si heterostructure NWs with lengths appropriate for device fabrication by devising a growth procedure that eliminates Au diffusion on the NW sidewalls and minimizes random kinking in the heterostructure NWs as deduced from detailed microscopy analysis. Fig. 1 a shows a cross-sectional SEM image of epitaxial Ge/Si axial NW heterostructures grown on a Ge(111) surface. The interface abruptness in these Ge/Si heterostructure NWs is of the order of the NW diameter. Some of these NWs develop a crystallographic kink that is {approx}20{sup o} off the <111> axis at about 300 nm away from the Ge/Si interface. This provides a natural marker for placing the gate contact electrodes and gate metal at appropriate location for desired high-on current and reduced ambipolarity as shown in Fig. 2. The 1D heterostructures allow band-edge engineering in the transport direction, not easily accessible in planar devices, providing an additional degree of freedom for designing tunnel FETs (TFETs). For instance, a Ge tunnel source can be used for efficient electron/hole tunneling and a Si drain can be used for reduced back-tunneling and ambipolar behavior. Interface abruptness on the other hand (particularly for doping) imposes challenges in these structures and others for realizing high performance TFETs in p-i-n junctions. Since the metal-semiconductor contacts provide a sharp interface with band-edge control, we use properly designed Schottky contacts (aided by 3D Silvaco simulations) as the tunnel barriers both at the source and drain and utilize the asymmetry in the Ge/Si channel bandgap to reduce ambipolar transport behavior generally observed in TFETs. Fig. 3 shows the room-temperature transfer curves of a Ge/Si heterostructure TFET (H-TFET) for different V{sub DS} values showing a maximum on-current of {approx}7 {micro}A, {approx}170 mV/decade inverse subthreshold slope and 5 orders of magnitude I{sub on}/I{sub off} ratios for all V{sub DS} biases considered here. This high on-current value is {approx}1750 X higher than that obtained with Si p-i-n{sup +} NW TFETs and {approx}35 X higher than that obtained with CNT TFET. The I{sub on}/I{sub off} ratio and inverse subthreshold slope compare favorably to that of Si {approx} 10{sup 3} I{sub on}/I{sub off} and {approx} 800 mV/decade SS{sup -1} but lags behind those of CNT TFET due to poor PECVD nitride gate oxide quality ({var_epsilon}{sub r} {approx} 3-4). The asymmetry in the Schottky barrier heights used here eliminates the stringent requirements of abrupt doped interfaces used in p-i-n based TFETs, which is hard to achieve both in thin-film and in NW growth. These initial promising results are expected to be further improved by using a high-k gate dielectric.
Numerical simulations indicate that significant fusion yields (>100 kJ) may be obtained by pulsed-power-driven implosions of cylindrical metal liners onto magnetized and preheated deuterium-tritium fuel. The primary physics risk to this approach is the Magneto-Rayleigh-Taylor (MRT) instability, which operates during both the acceleration and deceleration phase of the liner implosion. We have designed and performed some experiments to study the MRT during the acceleration phase, where the light fluid is purely magnetic. Results from our first series of experiments and plans for future experiments will be presented. According to simulations, an initial axial magnetic field of 10 T is compressed to >100 MG within the liner during the implosion. The magnetic pressure becomes comparable to the plasma pressure during deceleration, which could significantly affect the growth of the MRT instability at the fuel/liner interface. The MRT instability is also important in some astronomical objects such as the Crab Nebula (NGC1962). In particular, the morphological structure of the observed filaments may be determined by the ratio of the magnetic to material pressure and alignment of the magnetic field with the direction of acceleration [Hester, ApJ, 456, 225 1996]. Potential experiments to study this MRT behavior using the Z facility will be presented.
3-D cubic unit cell arrays containing split ring resonators were fabricated and characterized. The unit cells are {approx}3 orders-of-magnitude smaller than microwave SRR-based metamaterials and exhibit both electrically and magnetically excited resonances for normally incident TEM waves in addition to showing improved isotropic response.
An AlN MEMS resonator technology has been developed, enabling massively parallel filter arrays on a single chip. Low-loss filter banks covering the 10 MHz--10-GHz frequency range have been demonstrated, as has monolithic integration with inductors and CMOS circuitry. The high level of integration enables miniature multi-bandm spectrally aware, and cognitive radios.
We describe a time-domain spectroscopy system in the thermal infrared used for complete transmission and reflection characterization of metamaterials in amplitude and phase. The system uses a triple-output near-infrared ultrafast fiber laser, phase-locked difference frequency generation and phase-matched electro-optic sampling. We will present measurements of several metamaterials designs.
Four approaches to modeling multi-junction concentrating photovoltaic system performance are assessed by comparing modeled performance to measured performance. Measured weather, irradiance, and system performance data were collected on two systems over a one month period. Residual analysis is used to assess the models and to identify opportunities for model improvement. Large photovoltaic systems are typically developed as projects which supply electricity to a utility and are owned by independent power producers. Obtaining financing at favorable rates and attracting investors requires confidence in the projected energy yield from the plant. In this paper, various performance models for projecting annual energy yield from Concentrating Photovoltaic (CPV) systems are assessed by comparing measured system output to model predictions based on measured weather and irradiance data. The results are statistically analyzed to identify systematic error sources.
Four approaches to modeling multi-junction concentrating photovoltaic system performance are assessed by comparing modeled performance to measured performance. Measured weather, irradiance, and system performance data were collected on two systems over a one month period. Residual analysis is used to assess the models and to identify opportunities for model improvement.
Problem Statement: (1) Uncertainties in PV system performance and reliability impact business decisions - Project cost and financing estimates, Pricing service contracts and guarantees, Developing deployment and O&M strategies; (2) Understanding and reducing these uncertainties will help make the PV industry more competitive (3) Performance has typically been estimated without much attention to reliability of components; and (4) Tools are needed to assess all inputs to the value proposition (e.g., LCOE, cash flow, reputation, etc.). Goals and objectives are: (1) Develop a stochastic simulation model (in GoldSim) that can represent PV system performance as a function of system design, weather, reliability, and O&M policies; (2) Evaluate performance for an example system to quantify sources of uncertainty and identify dominant parameters via a sensitivity study; and (3) Example System - 1 inverter, 225 kW DC Array latitude tilt (90 strings of 12 modules {l_brace}1080 modules{r_brace}), Weather from Tucumcari, NM (TMY2 with annual uncertainty).
The Fracture-Matrix Transport (FMT) code developed at Sandia National Laboratories solves chemical equilibrium problems using the Pitzer activity coefficient model with a database containing actinide species. The code is capable of predicting actinide solubilities at 25 C in various ionic-strength solutions from dilute groundwaters to high-ionic-strength brines. The code uses oxidation state analogies, i.e., Am(III) is used to predict solubilities of actinides in the +III oxidation state; Th(IV) is used to predict solubilities of actinides in the +IV state; Np(V) is utilized to predict solubilities of actinides in the +V state. This code has been qualified for predicting actinide solubilities for the Waste Isolation Pilot Plant (WIPP) Compliance Certification Application in 1996, and Compliance Re-Certification Applications in 2004 and 2009. We have established revised actinide-solubility uncertainty ranges and probability distributions for Performance Assessment (PA) by comparing actinide solubilities predicted by the FMT code with solubility data in various solutions from the open literature. The literature data used in this study include solubilities in simple solutions (NaCl, NaHCO{sub 3}, Na{sub 2}CO{sub 3}, NaClO{sub 4}, KCl, K{sub 2}CO{sub 3}, etc.), binary mixing solutions (NaCl+NaHCO{sub 3}, NaCl+Na{sub 2}CO{sub 3}, KCl+K{sub 2}CO{sub 3}, etc.), ternary mixing solutions (NaCl+Na{sub 2}CO{sub 3}+KCl, NaHCO{sub 3}+Na{sub 2}CO{sub 3}+NaClO{sub 4}, etc.), and multi-component synthetic brines relevant to the WIPP.
The need for the engineering analysis of systems in which the transport of thermal energy occurs primarily through a conduction process is a common situation. For all but the simplest geometries and boundary conditions, analytic solutions to heat conduction problems are unavailable, thus forcing the analyst to call upon some type of approximate numerical procedure. A wide variety of numerical packages currently exist for such applications, ranging in sophistication from the large, general purpose, commercial codes, such as COMSOL, COSMOSWorks, ABAQUS and TSS to codes written by individuals for specific problem applications. The original purpose for developing the finite element code described here, COYOTE, was to bridge the gap between the complex commercial codes and the more simplistic, individual application programs. COYOTE was designed to treat most of the standard conduction problems of interest with a user-oriented input structure and format that was easily learned and remembered. Because of its architecture, the code has also proved useful for research in numerical algorithms and development of thermal analysis capabilities. This general philosophy has been retained in the current version of the program, COYOTE, Version 5.0, though the capabilities of the code have been significantly expanded. A major change in the code is its availability on parallel computer architectures and the increase in problem complexity and size that this implies. The present document describes the theoretical and numerical background for the COYOTE program. This volume is intended as a background document for the user's manual. Potential users of COYOTE are encouraged to become familiar with the present report and the simple example analyses reported in before using the program. The theoretical and numerical background for the finite element computer program, COYOTE, is presented in detail. COYOTE is designed for the multi-dimensional analysis of nonlinear heat conduction problems. A general description of the boundary value problems treated by the program is presented. The finite element formulation and the associated numerical methods used in COYOTE are also outlined. Instructions for use of the code are documented in SAND2010-0714.
The viscosity of molten salts comprising ternary and quaternary mixtures of the nitrates of sodium, potassium, lithium and calcium was determined experimentally. Viscosity was measured over the temperature range from near the relatively low liquidus temperatures of he individual mixtures to 200C. Molten salt mixtures that do not contain calcium nitrate exhibited relatively low viscosity and an Arrhenius temperature dependence. Molten salt mixtures that contained calcium nitrate were relatively more viscous and viscosity increased as the roportion of calcium nitrate increased. The temperature dependence of viscosity of molten salts containing calcium nitrate displayed curvature, rather than linearity, when plotted in Arrhenius format. Viscosity data for these mixtures were correlated by the Vogel-Fulcher- ammann-Hesse equation.
In a multiyear research agreement with Tenix Investments Pty. Ltd., Sandia has been developing field deployable technologies for detection of biotoxins in water supply systems. The unattended water sensor or UWS employs microfluidic chip based gel electrophoresis for monitoring biological analytes in a small integrated sensor platform. This instrument collects, prepares, and analyzes water samples in an automated manner. Sample analysis is done using the {mu}ChemLab{trademark} analysis module. This report uses analysis results of two datasets collected using the UWS to estimate performance of the device. The first dataset is made up of samples containing ricin at varying concentrations and is used for assessing instrument response and detection probability. The second dataset is comprised of analyses of water samples collected at a water utility which are used to assess the false positive probability. The analyses of the two sets are used to estimate the Receiver Operating Characteristic or ROC curves for the device at one set of operational and detection algorithm parameters. For these parameters and based on a statistical estimate, the ricin probability of detection is about 0.9 at a concentration of 5 nM for a false positive probability of 1 x 10{sup -6}.
The Sierra Toolkit computational mesh is a software library intended to support massively parallel multi-physics computations on dynamically changing unstructured meshes. This domain of intended use is inherently complex due to distributed memory parallelism, parallel scalability, heterogeneity of physics, heterogeneous discretization of an unstructured mesh, and runtime adaptation of the mesh. Management of this inherent complexity begins with a conceptual analysis and modeling of this domain of intended use; i.e., development of a domain model. The Sierra Toolkit computational mesh software library is designed and implemented based upon this domain model. Software developers using, maintaining, or extending the Sierra Toolkit computational mesh library must be familiar with the concepts/domain model presented in this report.
Because digital image correlation (DIC) has become such an important and standard tool in the toolbox of experimental mechanicists, a complete uncertainty quantification of the method is needed. It should be remembered that each DIC setup and series of images will have a unique uncertainty based on the calibration quality and the image and speckle quality of the analyzed images. Any pretest work done with a calibrated DIC stereo-rig to quantify the errors using known shapes and translations, while useful, do not necessarily reveal the uncertainty of a later test. This is particularly true with high-speed applications where actual test images are often less than ideal. Work has previously been completed on the mathematical underpinnings of DIC uncertainty quantification and is already published, this paper will present corresponding experimental work used to check the validity of the uncertainty equations.
A mesoscale dimensional artifact based on silicon bulk micromachining fabrication has been developed and manufactured with the intention of evaluating the artifact both on a high precision coordinate measuring machine (CMM) and video-probe based measuring systems. This hybrid artifact has features that can be located by both a touch probe and a video probe system with a k=2 uncertainty of 0.4 {micro}m, more than twice as good as a glass reference artifact. We also present evidence that this uncertainty could be lowered to as little as 50 nm (k=2). While video-probe based systems are commonly used to inspect mesoscale mechanical components, a video-probe system's certified accuracy is generally much worse than its repeatability. To solve this problem, an artifact has been developed which can be calibrated using a commercially available high-accuracy tactile system and then be used to calibrate typical production vision-based measurement systems. This allows for error mapping to a higher degree of accuracy than is possible with a glass reference artifact. Details of the designed features and manufacturing process of the hybrid dimensional artifact are given and a comparison of the designed features to the measured features of the manufactured artifact is presented and discussed. Measurement results from vision and touch probe systems are compared and evaluated to determine the capability of the manufactured artifact to serve as a calibration tool for video-probe systems. An uncertainty analysis for calibration of the artifact using a CMM is presented.
Forward radiation transport is the problem of calculating the radiation field given a description of the radiation source and transport medium. In contrast, inverse transport is the problem of inferring the configuration of the radiation source and transport medium from measurements of the radiation field. As such, the identification and characterization of special nuclear materials (SNM) is a problem of inverse radiation transport, and numerous techniques to solve this problem have been previously developed. The authors have developed a solver based on nonlinear regression applied to deterministic coupled neutron-photon transport calculations. The subject of this paper is the experimental validation of that solver. This paper describes a series of experiments conducted with a 4.5-kg sphere of alpha-phase, weapons-grade plutonium. The source was measured in six different configurations: bare, and reflected by high-density polyethylene (HDPE) spherical shells with total thicknesses of 1.27, 2.54, 3.81, 7.62, and 15.24 cm. Neutron and photon emissions from the source were measured using three instruments: a gross neutron counter, a portable neutron multiplicity counter, and a high-resolution gamma spectrometer. These measurements were used as input to the inverse radiation transport solver to characterize the solver's ability to correctly infer the configuration of the source from its measured signatures.
We present Poblano v1.0, a Matlab toolbox for solving gradient-based unconstrained optimization problems. Poblano implements three optimization methods (nonlinear conjugate gradients, limited-memory BFGS, and truncated Newton) that require only first order derivative information. In this paper, we describe the Poblano methods, provide numerous examples on how to use Poblano, and present results of Poblano used in solving problems from a standard test collection of unconstrained optimization problems.
The difficulty of calculating the ambient properties of molecular crystals, such as the explosive PETN, has long hampered much needed computational investigations of these materials. One reason for the shortcomings is that the exchange-correlation functionals available for Density Functional Theory (DFT) based calculations do not correctly describe the weak intermolecular van der Waals' forces present in molecular crystals. However, this weak interaction also poses other challenges for the computational schemes used. We will discuss these issues in the context of calculations of lattice constants and structure of PETN with a number of different functionals, and also discuss if these limitations can be circumvented for studies at non-ambient conditions.
