A 1 MV linear transformer driver (LTD) is being tested with a large area e-beam diode load at Sandia National Laboratories (SNL). The experiments will be utilized to determine the repeatability of the output pulse and the reliability of the components. The 1 MV accelerator is being used to determine the feasibility of designing a 6 MV LTD for radiography experiments. The peak voltage, risetime, and pulse width as well as the cavity timing jitter are analyzed to determine the repeatability of the output pulse.
Graph partitioning is often used for load balancing in parallel computing, but it is known that hypergraph partitioning has several advantages. First, hypergraphs more accurately model communication volume, and second, they are more expressive and can better represent nonsymmetric problems. Hypergraph partitioning is particularly suited to parallel sparse matrix-vector multiplication, a common kernel in scientific computing. We present a parallel software package for hypergraph (and sparse matrix) partitioning developed at Sandia National Labs. The algorithm is a variation on multilevel partitioning. Our parallel implementation is novel in that it uses a two-dimensional data distribution among processors. We present empirical results that show our parallel implementation achieves good speedup on several large problems (up to 33 million nonzeros) with up to 64 processors on a Linux cluster.
A methodology was developed for computing the probability that the sensor dart for the 'Near Real-Time Site Characterization for Assured HDBT Defeat' Grand-Challenge LDRD project will survive deployment over a forested region. The probability can be decomposed into three approximately independent probabilities that account for forest coverage, branch density and the physics of an impact between the dart and a tree branch. The probability that a dart survives an impact with a tree branch was determined from the deflection induced by the impact. If a dart that was deflected so that it impacted the ground at an angle of attack exceeding a user-specified, threshold value, the dart was assumed to not survive the impact with the branch; otherwise it was assumed to have survived. A computer code was developed for calculating dart angle of attack at impact with the ground and a Monte Carlo scheme was used to calculate the probability distribution of a sensor dart surviving an impact with a branch as a function of branch radius, length, and height from the ground. Both an early prototype design and the current dart design were used in these studies. As a general rule of thumb, it we observed that for reasonably generic trees and for a threshold angle of attack of 5{sup o} (which is conservative for dart survival), the probability of reaching the ground with an angle of attack less than the threshold is on the order of 30% for the prototype dart design and 60% for the current dart design, though these numbers should be treated with some caution.
The study described in this report used mathematical modeling to estimate health risks from exposure to depleted uranium (DU) during the 1991 Gulf War for both U.S. troops and nearby Iraqi civilians. The analysis found that the risks of DU-induced leukemia or birth defects are far too small to result in an observable increase in these health effects among exposed veterans or Iraqi civilians. Only a few veterans in vehicles accidentally struck by U.S. DU munitions are predicted to have inhaled sufficient quantities of DU particulate to incur any significant health risk (i.e., the possibility of temporary kidney damage from the chemical toxicity of uranium and about a 1% chance of fatal lung cancer). The health risk to all downwind civilians is predicted to be extremely small. Recommendations for monitoring are made for certain exposed groups. Although the study found fairly large calculational uncertainties, the models developed and used are generally valid. The analysis was also used to assess potential uranium health hazards for workers in the weapons complex. No illnesses are projected for uranium workers following standard guidelines; nonetheless, some research suggests that more conservative guidelines should be considered.
The annual program report provides detailed information about all aspects of the SNL/CA Pollution Prevention Program for a given calendar year. It functions as supporting documentation to the ''SNL/CA Environmental Management System Program Manual''. The 2005 program report describes the activities undertaken during the past year, and activities planned in future years to implement the Pollution Prevention Program, one of six programs that supports environmental management at SNL/CA.
Particle image velocimetry (PIV) data have been acquired using three different configurations in the far-field of the interaction of a transverse supersonic jet with a transonic crossflow. The configurations included two-dimensional PIV in the centerline streamwise plane at two overlapping stations, as well as stereoscopic PIV in both the same streamwise plane and the crossplane. The streamwise data show the downstream evolution of the interaction whereas the crossplane data directly reveal its vortex structure. The measurement planes intersect at a common line, allowing a comparison of those mean velocity components and turbulent stresses common to all configurations. All data from the streamwise plane agree to within their estimated uncertainties, but data from the crossplane exhibit reduced velocity and turbulent stress magnitudes by a small but significant degree. Additionally, the vertical positions of the peak velocities are slightly nearer the wall for the crossplane configuration. This comparison suggests that routine methods of uncertainty quantification for data used in the validation of computational models may not fully capture the error sources of an experiment.
Binders such as Estane, Teflon, Kel F and HTPB are typically used in heterogeneous explosives to bond polycrystalline constituents together as an energetic composite. Combined theoretical and experimental studies are underway to unravel the mechanical response of these materials when subjected to isentropic compression loading. Key to this effort is the determination of appropriate constitutive and EOS property data at extremely high stress-strain states as required for detailed mesoscale modeling. The Sandia Z accelerator and associated diagnostics provides new insights into mechanical response of these nonreactive constituents via isentropic ramp-wave compression loading. Several thicknesses of samples, varied from 0.3 to 1.2 mm, were subjected to a ramp load of {approx}42 Kbar over 500 ns duration using the Sandia Z-machine. Profiles of transmitted ramp waves were measured at window interfaces using conventional VISAR. Shock physics analysis is then used to determine the nonlinear material response of the binder materials. In this presentation we discuss experimental and modeling details of the ramp wave loading ICE experiments designed specifically for binder materials.
Exposure of optical materials to transient-ionizing-radiation fields can give rise to transient and/or permanent photodarkening effects. In laser materials, such as YAG, such induced optical loss can result in significant degradation of the lasing characteristic of the material, making its selection for optical device applications in radiation environments unfeasible. In the present study, the effects of ionizing radiation on the optical response of undoped and 1.1% Nd-doped single-crystal and polycrystalline YAG have been investigated. In the undoped materials it is seen that both laser materials exhibit significant loss at the 1.06 ?m lasing wavelength following exposure to a 40 krad, 30 nsec pulse of gamma radiation. In the undoped single-crystal samples, the transmission loss is initially large but exhibits a rapid recovery. By contrast, the undoped polycrystalline YAG experiences an initial 100% loss in transmission, becoming totally opaque at 1.06 ?m following the radiation pulse. This loss is slow to recover and a large residual permanent photodarkening effect is observed. Nd-doping improves the optical response of the materials in that the radiation-induced optical loss is substantially smaller in both the polycrystalline and single-crystal YAG samples. Preliminary results on the radiation response of elevated-temperature samples will also be reported.
Applications requiring injection of a high-power multimode laser into multiple fibers with equal energies, or specific energy ratios, provide unique design challenges. As with most all systems, engineering trades must balance competing requirements to obtain an optimal overall design. This is particularly true when fabrication issues are considered in the design process. A few of these competing design requirements are discussed in this conceptually simple system. This fiber injection system consists of three components; a refractive beam homogenizer, a diffractive beamsplitter, and a fiber array. We show the design process, starting with first-order design, for an example fiber injection system that couples a high-power YAG laser into seven fibers. Design goals include high efficiency, good beamsplitting uniformity, compact overall size, maximum mode filling of the fibers, and low cost of fabrication and assembly.
In support of the US Government and the International Atomic Energy Agency (IAEA) Nuclear Security Programmes, Sandia National Laboratories (SNL) has advocated and practiced a risk-based, systematic approach to nuclear security. The risk equation has been implemented as the basis for a performance methodology for the design and evaluation of Physical Protection Systems against a Design Basis Threat (DBT) for theft or sabotage of nuclear and/or radiological materials. Since integrated systems must include people as well as technology and the man-machine interface, a critical aspect of the human element is to train all stakeholders in nuclear security on the systems approach. Current training courses have been beneficial but are still limited in scope. SNL has developed two primary international courses and is completing development of three new courses that will be offered and presented in the near term. In the long-term, SNL envisions establishing a comprehensive nuclear security training curriculum that will be developed along with a series of forthcoming IAEA Nuclear Security Series guidance documents.
This article presents continuum simulations of viscous polymer flow during nanoimprint lithography (NIL) for embossing tools having irregular spacings and sizes. Simulations vary nonuniform embossing tool geometry to distinguish geometric quantities governing cavity filling order, polymer peak deformation, and global mold filling times. A characteristic NIL velocity predicts cavity filling order. In general, small cavities fill more quickly than large cavities, while cavity spacing modulates polymer deformation mode. Individual cavity size, not total filling volume, dominates replication time, with large differences in individual cavity size resulting in nonuniform, squeeze flow filling. High density features can be modeled as a solid indenter in squeeze flow to accurately predict polymer flow and allow for optimization of wafer-scale replication. The present simulations make it possible to design imprint templates capable of distributing pressure evenly across the mold surface and facilitating symmetric polymer flow over large areas to prevent mold deformation and nonuniform residual layer thickness.
This paper investigates the transient response of 50-nm gate length fully and partially depleted SOI and bulk devices to pulsed laser and heavy ion microbeam irradiations. The measured transient signals on 50-nm fully depleted devices are very short, and the collected charge is small compared to older 0.25-{micro}m generation SOI and bulk devices. We analyze in detail the influence of the SOI architecture (fully or partially depleted) on the pulse duration and the amount of bipolar amplification. For bulk devices, the doping engineering is shown to have large effects on the duration of the transient signals and on the charge collection efficiency.
Structural reorientations in metallic fcc nanowires are controlled by a combination of size, thermal energy, and the type of defects formed during inelastic deformation. By utilizing atomistic simulations, we show that certain fcc nanowires can exhibit both shape memory and pseudoelastic behavior. We also show that the formation of defect-free twins, a process related to the material stacking fault energy, nanometer size scale, and surface stresses is the mechanism that controls the ability of fcc nanowires of different materials to show a reversible transition between two crystal orientations during loading and thus shape memory and pseudoelasticity.
This multinational, multi-phase spent fuel sabotage test program is quantifying the aerosol particles produced when the products of a high energy density device (HEDD) interact with and explosively particulate test rodlets that contain pellets of either surrogate materials or actual spent fuel. This program has been underway for several years. This program provides data that are relevant to some sabotage scenarios in relation to spent fuel transport and storage casks, and associated risk assessments. The program also provides significant technical and political benefits in international cooperation. We are quantifying the Spent Fuel Ratio (SFR), the ratio of the aerosol particles released from HEDD-impacted actual spent fuel to the aerosol particles produced from surrogate materials, measured under closely matched test conditions, in a contained test chamber. In addition, we are measuring the amounts, nuclide content, size distribution of the released aerosol materials, and enhanced sorption of volatile fission product nuclides onto specific aerosol particle size fractions. These data are the input for follow-on modeling studies to quantify respirable hazards, associated radiological risk assessments, vulnerability assessments, and potential cask physical protection design modifications. This document includes an updated description of the test program and test components for all work and plans made, or revised, during FY 2004. It also serves as a program status report as of the end of FY 2004. All available test results, observations, and aerosol analyses plus interpretations--primarily for surrogate material Phase 2 tests, series 2/5A through 2/9B, using cerium oxide sintered ceramic pellets are included. Advanced plans and progress are described for upcoming tests with unirradiated, depleted uranium oxide and actual spent fuel test rodlets. This spent fuel sabotage--aerosol test program is coordinated with the international Working Group for Sabotage Concerns of Transport and Storage Casks (WGSTSC) and supported by both the U.S. Department of Energy and the Nuclear Regulatory Commission.
Remote spectral sensing offers an attractive means of mapping river water quality over wide spatial regions. While previous research has focused on development of spectral indices and models to predict river water quality based on remote images, little attention has been paid to subsequent validation of these predictions. To address this oversight, we describe a retrospective analysis of remote, multispectral Compact Airborne Spectrographic Imager (CASI) images of the Ohio River and its Licking River and Little Miami River tributaries. In conjunction with the CASI acquisitions, ground truth measurements of chlorophyll-a concentration and turbidity were made for a small set of locations in the Ohio River. Partial least squares regression models relating the remote river images to ground truth measurements of chlorophyll-a concentration and turbidity for the Ohio River were developed. Employing these multivariate models, chlorophyll-a concentrations and turbidity levels were predicted in river pixels lacking ground truth measurements, generating detailed estimated water quality maps. An important but often neglected step in the regression process is to validate prediction results using a spectral residual statistic. For both the chlorophyll-a and turbidity regression models, a spectral residual value was calculated for each river pixel and compared to the associated statistical confidence limit for the model. These spectral residual statistic results revealed that while the chlorophyll-a and turbidity models could validly be applied to a vast majority of Ohio River and Licking River pixels, application of these models to Little Miami River pixels was inappropriate due to an unmodeled source of spectral variation.
Slow, dense granular flows often exhibit thin, localized regions of particle motion, called shear bands, separating largely solid-like regions. Recent experiments using a split-bottom Couette cell found that the width of the shear zone grew as the pack height increased and the azimuthal velocities when rescaled fall on a universal curve regardless of the particle properties. Here we present large-scale Discrete Element simulations of a similar system for packs of varying height up to 180,000 monodisperse spheres. The onset and evolution of granular shear flow is investigated as a function of height. We find a transition in the nature of the shear as a characteristic height is exceeded. Below this height there is a central quasi-solid core; above this height we observe the onset of additional axial shear associated with a torsional failure mode of the inner core. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height while the axial width remains narrow and fixed.
This Corrective Measures Evaluation Report was prepared as directed by the Compliance Order on Consent issued by the New Mexico Environment Department to document the process of selecting the preferred remedial alternative for contaminated groundwater at Technical Area V. Supporting information includes background information about the site conditions and potential receptors and an overview of work performed during the Corrective Measures Evaluation. Evaluation of remedial alternatives included identification and description of four remedial alternatives, an overview of the evaluation criteria and approach, qualitative and quantitative evaluation of remedial alternatives, and selection of the preferred remedial alternative. As a result of the Corrective Measures Evaluation, it was determined that monitored natural attenuation of all contaminants of concern (trichloroethene, tetrachloroethene, and nitrate) was the preferred remedial alternative for implementation as the corrective measure to remediate contaminated groundwater at Technical Area V of Sandia National Laboratories/New Mexico. Finally, design criteria to meet cleanup goals and objectives and the corrective measures implementation schedule for the preferred remedial alternative are presented.
A commercial stereolithography (SL) machine was modified to integrate fluid dispensing or direct-write (DW) technology with SL in an integrated manufacturing environment for automated and efficient hybrid manufacturing of complex electrical devices, combining three-dimensional (3D) electrical circuitry with SL-manufactured parts. The modified SL system operates similarly to a commercially available machine, although build interrupts were used to stop and start the SL build while depositing fluid using the DW system. An additional linear encoder was attached to the SL platform z-stage and used to maintain accurate part registration during the SL and DW build processes. Individual STL files were required as part of the manufacturing process plan. The DW system employed a three-axis translation mechanism that was integrated with the commercial SL machine. Registration between the SL part, SL laser and the DW nozzle was maintained through the use of 0.025-inch diameter cylindrical reference holes manufactured in the part during SL. After depositing conductive ink using DW, the SL laser was commanded to trace the profile until the ink was cured. The current system allows for easy exchange between SL and DW in order to manufacture fully functional 3D electrical circuits and structures in a semi-automated environment. To demonstrate the manufacturing capabilities, the hybrid SL/DW setup was used to make a simple multi-layer SL part with embedded circuitry. This hybrid system is not intended to function as a commercial system, it is intended for experimental demonstration only. This hybrid SL/DW system has the potential for manufacturing fully functional electromechanical devices that are more compact, less expensive, and more reliable than their conventional predecessors, and work is ongoing in order to fully automate the current system.
Electron backscattered diffraction (EBSD) is a widely used technique for both identifying the crystallographic phase and for mapping the orientation of crystalline materials on the micron length scale. Often the operating conditions necessary for phase identification are not suitable for orientation mapping and vice versa. In an effort to optimize the speed involved in the mapping technique, pattern quality is sacrificed and the wealth of information present in an EBSD pattern is compressed to basically 4 values: a matched phase and three Euler angles. However, ab initio identification of phases from EBSD patterns requires high quality patterns and fairly intense computation. Spectrum imaging is an analytical approach that may offer some solutions to the aforementioned problems. Spectrum imaging consists of collecting a whole spectrum at each pixel in a mapping style measurement. This large set of data is then analyzed using multivariate statistical analysis (MSA) techniques such as principle components analysis, multivariate curve resolution, or other least squares based techniques. The result of these calculations is a set of component spectral shapes with corresponding abundances that allow the analyst to extract the greatest amount of physically relevant information from an otherwise enormous data set. Spectrum imaging has been used successfully in EDX microanalysis (both in the SEM and TEM), TOF-SIMS, WDS, and EELS. To examine the potential benefits of the spectrum imaging approach for EBSD data, a series of basic experiments and calculations were run. Test data sets (20 x 20 patterns in .jpeg format) on polycrystalline Al and on the directionally solidified eutectic oxide, CoO/ZrO{sub 2}(CaO), were collected using the HKL Channel 5 system with a Nordlys detector under normal mapping conditions. The data was collected on a FEI dual beam FIB (model DB235) and a Zeiss (Supra 55 VP) SEM at 20keV for Al and CoO/ZrO{sub 2}(CaO), respectively. The data sets were analyzed according to the schematic shown in Figure 1. Each EBSD pattern was hough transformed, unzipped into a 1-D vector of channels with intensities ranging from 0-255, and then added to an overall data matrix. A range of treatments (edge/no edge detection, spatial simplicity/spectral simplicity, etc.) were examined to determine the optimal way of treating the data. The multivariate analyses were performed using the AXSIA code developed at Sandia National Laboratories. The MSA techniques were able to correctly identify individual grains in the Al sample and individual phases in the CoO/ZrO{sub 2}(CaO) sample. For each component EBSD pattern identified from the Al data, a corresponding color map of abundance can be seen which clearly corresponds to a single grain (Figure 2). The success in the CoO/ZrO{sub 2}(CaO) sample is particularly notable due to both phases sharing the Fm-3m space group which would confuse most autoindexing routines. The range of analytical treatments identified two extremes in results: a minimal number of components (patterns) with only kikuchi line positions present or a larger number of components with full intensity information present. The further application of these results to phase mapping will be discussed.
The adsorption and desorption behavior of a planar microfabricated preconcentrator (PC) has been modeled and simulated using the computational fluid dynamics (CFD) package CFDRC-ACE+trade. By comparison with the results of a designed experiment, model parameters were determined. Assuming a first-order reaction for the adsorption of a light hydrocarbon chemical analyte onto the PC adsorbent and a unity-value sticking coefficient, a rate constant of 36,500 s{sup -1} was obtained. This compares favorably with the value of 25,300 s{sup -1} obtained by application of the Modified-Wheeler equation. The modeled rate constant depends on the concentration of adsorbent sites, estimated to be 6.94 ldr 10{sup -8} kmol/m{sup 2} for the Carboxen 1000 adsorbent used. Using the integral method, desorption was found to be first order with an Arrhenius temperature dependence and an activation energy of 30.1 kj/mol. Validation of this model is reported herein, including the use of Aris-Taylor dispersion to predict the influence of fluidics surrounding the PC. A maximum in desorption peak area with flow rate, predicted from a quadratic fit to the results of the designed experiment, was not observed in the 2-D simulation. Either approximations in the simulated model or the nonphysical nature of the quadratic fit are responsible. Despite the apparent simplicity of the model, the simulation is internally self consistent and capable of predicting performance of new device designs. To apply the method to other analytes and other adsorbent materials, only a limited number of comparisons to experiment are required to obtain the necessary rate constants.
The shock behavior of two varieties of the ceramic silicon carbide was investigated through a series of time-resolved plate impact experiments reaching stresses of over 140 GPa. The Hugoniot data obtained are consistent for the two varieties tested as well as with most data from the literature. Through the use of reshock and release configurations, reloading and unloading responses for the material were found. Analysis of these responses provides a measure of the ceramic's strength behavior as quantified by the shear stress and the strength in the Hugoniot state. While previous strength measurements were limited to stresses of 20-25 GPa, measurements were made to 105 GPa in the current study. The initial unloading response is found to be elastic to stresses as high as 105 GPa, the level at which a solid-to-solid phase transformation is observed. While the unloading response lies significantly below the Hugoniot, the reloading response essentially follows it. This differs significantly from previous results for B{sub 4}C and Al{sub 2}O{sub 3}. The strength of the material increases by about 50% at stresses of 50-75 GPa before falling off somewhat as the phase transformation is approached. Thus, the strength behavior of SiC in planar impact experiments could be characterized as metal-like in character. The previously reported phase transformation at {approx}105 GPa was readily detected by the reshock technique, but it initially eluded detection with traditional shock experiments. This illustrates the utility of the reshock technique for identifying phase transformations. The transformation in SiC was found to occur at about 104 GPa with an associated volume change of about 9%.
We examine the scaling to ignition of the energy deposition of laser generated electrons in compressed fast ignition cores. Relevant cores have densities of several hundred g/cm{sup 3}, with a few keV initial temperature. As the laser intensities increase approaching ignition systems, on the order of a few 10{sup 21}W/cm{sup 2}, the hot electron energies expected to approach 100MeV. Most certainly anomalous processes must play a role in the energy transfer, but the exact nature of these processes, as well as a practical way to model them, remain open issues. Traditional PIC explicit methods are limited to low densities on current and anticipated computing platforms, so the study of relevant parameter ranges has received so far little attention. We use LSP to examine a relativistic electron beam (presumed generated from a laser plasma interaction) of legislated energy and angular distribution is injected into a 3D block of compressed DT. Collective effects will determine the stopping, most likely driven by magnetic field filamentation. The scaling of the stopping as a function of block density and temperature, as well as hot electron current and laser intensity is presented. Sub-grid models may be profitably used and degenerate effects included in the solution of this problem.
Grain size and texture of Ni electrodeposited from sulfamate baths depend greatly on current density. Increasing grain size is observed with increasing current density and the deposit texture changes from (110) at current densities lower than 5 mA cm{sup -2} to (100) for higher current densities. Co-deposition of Mn modifies the deposit structure by favoring the growth of the (110) texture and decreasing the average grain size even as the current density increases. While the average Mn film content increases with increasing current density, local Mn concentrations are a more complex function of deposition parameters, as indicated by atom probe tomography measurements. In both direct-current plated and pulse plated films, large variations on a nanometer scale in local Mn concentration are observed.
The effects of Mg alloying on the temporal evolution of Al{sub 3}Sc (L1{sub 2} structure) nanoscale precipitates are investigated, focusing on the morphology and coarsening kinetics of Al{sub 3}Sc precipitates in an Al-2.2 Mg-0.12 Sc at.% alloy aged between 300 and 400 C. Approximately spheroidal precipitates are obtained after aging at 300 C and irregular morphologies are observed at 400 C. The coarsening behavior is studied using conventional and high-resolution transmission electron microscopies to obtain the temporal evolution of the precipitate radius, and atom-probe tomography is employed to measure the Sc concentration in the {alpha}-matrix. The coarsening kinetics are analyzed using a coarsening model developed by Kuehmann and Voorhees for ternary systems [Kuehmann CJ, Voorhees PW. Metall Mater Trans A 1996;27:937]. Values of the interfacial free energy and diffusion coefficient for Sc diffusion in this Al-Mg-Sc alloy at 300 C are independently calculated, and are in good agreement with the calculated value of interfacial free energy and the experimental diffusivity obtained for the Al-Sc system.
Linear algebra is a powerful and proven tool in web search. Techniques, such as the PageRank algorithm of Brin and Page and the HITS algorithm of Kleinberg, score web pages based on the principal eigenvector (or singular vector) of a particular non-negative matrix that captures the hyperlink structure of the web graph. We propose and test a new methodology that uses multilinear algebra to elicit more information from a higher-order representation of the hyperlink graph. We start by labeling the edges in our graph with the anchor text of the hyperlinks so that the associated linear algebra representation is a sparse, three-way tensor. The first two dimensions of the tensor represent the web pages while the third dimension adds the anchor text. We then use the rank-1 factors of a multilinear PARAFAC tensor decomposition, which are akin to singular vectors of the SVD, to automatically identify topics in the collection along with the associated authoritative web pages.
Controlling the distribution of chemical constituents within complex, structurally heterogeneous systems represents one of the fundamental challenges of alloy design. We demonstrate how the combination of recent developments in sophisticated experimental high resolution characterization techniques and ab initio theoretical methods provide the basis for a detailed level of understanding of the microscopic factors governing compositional distributions in metallic alloys. In a study of the partitioning of Mg in two-phase ternary Al-Sc-Mg alloys by atom-probe tomography, we identify a large Mg concentration enhancement at the coherent {alpha}-Al/Al{sub 3}Sc heterophase interface with a relative Gibbsian interfacial excess of Mg with respect to Al and Sc, {Lambda}{sub Mg}{sup rel}, equal to 1.9 {+-} 0.5 atom nm{sup -2}. The corresponding calculated value of {Lambda}{sub Mg}{sup rel} is -1.2 atom nm{sup -2}. Theoretical ab initio investigations establish an equilibrium driving force for Mg interfacial segregation that is primarily chemical in nature and reflects the strength of the Mg-Sc interactions in an Al-rich alloy.
Effective monitoring of large computational clusters demands the analysis of a vast amount of raw data from a large number of machines. The fundamental interactions of the system are not, however, well-defined, making it difficult to draw meaningful conclusions from this data, even if one were able to efficiently handle and process it. In this paper we show that computational clusters, because they are comprised of a large number of identical machines, behave in a statistically meaningful fashion. We therefore can employ normal statistical methods to derive information about individual systems and their environment and to detect problems sooner than with traditional mechanisms. We discuss design details necessary to use these methods on a large system in a timely and low-impact fashion.
Electrochromic (EC) materials are used in 'smart' windows that can be darkened by applying a voltage across an EC stack on the window. The associated change in refractive index (n) in the EC materials might allow their use in tunable or temperature-insensitive Fabry-Perot filters and transmissive-spatial-light-modulators (SLMs). The authors are conducting a preliminary evaluation of these materials in many applications, including target-in-the-loop systems. Data on tungsten oxide, WO{sub 3}, the workhorse EC material, indicate that it's possible to achieve modest changes in n with only slight increases in absorption between the visible and {approx}10 {micro}m. This might enable construction of a tunable Fabry-Perot filter consisting of an active EC layer (e.g. WO{sub 3}) and a proton conductor (e.g.Ta{sub 2}O{sub 5}) sandwiched between two gold electrodes. A SLM might be produced by replacing the gold with a transparent conductor (e.g. ITO). This SLM would allow broad-band operation like a micromirror array. Since it's a transmission element, simple optical designs like those in liquid-crystal systems would be possible. Our team has fabricated EC stacks and characterized their switching speed and optical properties (n, k). We plan to study the interplay between process parameters, film properties, and performance characteristics associated with the FP-filter and then extend what we learn to SLMs. Our goals are to understand whether the changes in absorption associated with changes in n are acceptable, and whether it's possible to design an EC-stack that's fast enough to be interesting. We'll present our preliminary findings regarding the potential viability of EC materials for target-in-the-loop applications.
There are many important biological processes involving lipid bilayers on times scales beyond that accessible by atomistic simulations. We have developed coarse-grained, bead-spring models of lipid molecules to treat membrane fusion, domain formation and the general physical characteristics of lipid bilayers. A key aspect of these coarse-grained models is that the liquid nature of a bilayer is explicitly present in the simulations; the lipids diffuse far beyond their neighbors in contrast to atomistic simulations. With these models self-assembly into a bilayer starting from a random configuration of lipids and solvent is readily simulated. We have performed extensive simulations to characterize these lipid models in single component lipid bilayers. For a variety of tail lengths, the area per lipid as a function of temperature has been calculated; the liquid-gel transition has been characterized. Models have been developed for a variety of lipids including double bonds in the lipid tails. Simulation results will be presented for fusion and domain formation.
The latency and throughput of MPI messages are critically important to a range of parallel scientific applications. In many modern networks, both of these performance characteristics are largely driven by the performance of a processor on the network interface. Because of the semantics of MPI, this embedded processor is forced to traverse a linked list of posted receives each time a message is received. As this list grows long, the latency of message reception grows and the throughput of MPI messages decreases. This paper presents a novel hardware feature to handle list management functions on a network interface. By moving functions such as list insertion, list traversal, and list deletion to the hardware unit, latencies are decreased by up to 20% in the zero length queue case with dramatic improvements in the presence of long queues. Similarly, the throughput is increased by up to 10% in the zero length queue case and by nearly 100% in the presence queues of 30 messages.
The conclusions of this report are: (1) 1D and 2D RMHD simulations indicate feasibility of producing high thermonuclear neutron yields in deuterium and DT gas-puff Z-pinches -- (a) Z 1.7 x 10{sup 13} DD neutrons at 70 kV, 13 MA (Z1384); (b) (3 to 6) x 10{sup 14} at 90 kV, 17 MA (Z1422); (c) Predicted for ZR 2 x 10{sup 15} DD and 6 x 10{sup 16} DT neutrons; (2) Theory and modeling issues -- collisionless ions, nonthermal ions; (3) Experimental data on the origin of the neutrons not yet conclusive, need more shots; and (4) Applications -- (a) Fusion 2.5 and 14 MeV neutron source; (b) Pulsed subcritical neutron source with uranium blanket for {approx}10x neutron and {approx}1000x energy multiplication (Smirnov, Feoktistov and Klimov); and (c) Fusion-assisted keV x-ray plasma radiation source.
New experimental diagnostics and computational modeling provide an unprecedented means for improving the understanding of energetic material behavior at the mesoscale (grain or crystal ensemble levels). This study focuses on the determination of appropriate constitutive and EOS property data of the constituents of an energetic composite at high stress and moderate strain-rate states. The Sandia Z accelerator is used to determine the mechanical response of energetic composites via isentropic ramp wave compression loading. In this paper we describe an energy source method in CTH that models ramp loading for the analysis of ICE experiments. This approach is applied to design experimental configurations to probe the constituent response of PBX 9501 subjected to {approx}40 Kbar ramp load over 300 ns duration. Multiple VISAR are used to determine the averaged response of the composite material in comparison to the individual constituents including the effects of anisotropy of HMX crystals and the interactions of fine crystallites with binder material.
Mechanical dynamics can be a determining factor for the switching speed of radio-frequency microelectromechanical systems (RF MEMS) switches. This paper presents the simulation of the mechanical motion of a microswitch under actuation. The switch has a plate suspended by springs. When an electrostatic actuation is applied, the plate moves toward the substrate and closes the switch. Simulations are calculated via a high-fidelity finite element model that couples solid dynamics with electrostatic actuation. It incorporates non-linear coupled dynamics and accommodates fabrication variations. Experimental modal analysis gives results in the frequency domain that verifies the natural frequencies and mode shapes predicted by the model. An effective 1D model is created and used to calculate an actuation voltage waveform that minimizes switch velocity at closure. In the experiment, the switch is actuated with this actuation voltage, and the displacements of the switch at various points are measured using a laser Doppler velocimeter through a microscope. The experiments are repeated on several switches from different batches. The experimental results verify the model.
Several mixture models are evaluated for their suitability in predicting the equivalent permittivity of dielectric particles in a dielectric medium for intermediate solid volume fractions (0.4 to 0.6). Predictions of the Maxwell, Rayleigh, Bottcher and Bruggeman models are compared to computational simulations of several arrangements of solid particles in a gas and to the experimentally determined permittivity of a static particle bed. The experiment uses spherical glass beads in air, so air and glass permittivity values (1 and 7, respectively) are used with all of the models and simulations. The experimental system used to measure the permittivity of the static particle bed and its calibration are described. The Rayleigh model is found to be suitable for predicting permittivity over the entire range of solid volume fractions (0-0.6).
Proliferation of degrees-of-freedom has plagued discontinuous Galerkin methodology from its inception over 30 years ago. This paper develops a new computational formulation that combines the advantages of discontinuous Galerkin methods with the data structure of their continuous Galerkin counterparts. The new method uses local, element-wise problems to project a continuous finite element space into a given discontinuous space, and then applies a discontinuous Galerkin formulation. The projection leads to parameterization of the discontinuous degrees-of-freedom by their continuous counterparts and has a variational multiscale interpretation. This significantly reduces the computational burden and, at the same time, little or no degradation of the solution occurs. In fact, the new method produces improved solutions compared with the traditional discontinuous Galerkin method in some situations.
The intense magnetic field generated in the 20 MA Z-machine is used to accelerate metallic flyer plates to high velocity for the purpose of generating strong shocks in equation of state experiments. We present results pertaining to experiments in which a 0.085 cm thick Al flyer plate is magnetically accelerated across a vacuum gap into a quartz target. Peak magnetic drive pressures up to 4.9 Mbar were produced, which yielded a record 34 km/s flyer velocity without destroying it by shock formation or Joule heating. Two-dimensional MHD simulation was used to optimize the magnetic drive pressure on the flyer surface, shape the current pulse to accelerate the flyer without shock formation (i.e., quasi-isentropically), and predict the flyer velocity. Shock pressures up to 11.5 Mbar were produced in quartz. Accurate measurements of the shock velocity indicate that a fraction of the flyer is at solid density when it arrives at the target. Comparison of measurements and simulation results yields a consistent picture of the flyer state at impact with the quartz target.
We present Density Functional Theory (DFT) calculations of water in a region of phase space of interest in shock experiments. The onset of electrical conductivity in shocked water is determined by ionic conductivity, with the electron contribution dominating at higher pressures. The ionic contribution to the conduction is calculated from proton diffusion (Green-Kubo formula) and the electronic contribution is calculated using the Kubo-Greenwood formula [1]. The calculations are performed with VASP, a plane-wave pseudopotential code. At 2000K and a density of 2.3 g/cc, we find a significant dissociation of water into H, OH, and H3O, not only intermittent formation of OH - H3O pairs as suggested earlier for 2000 K and 1.95 g/cc [2]. The calculated conductivity is compared to experimental data [3]. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy's National Nuclear Safety Administration under contract DE-AC04-94AL85000. This project was supported by the Sandia LDRD office. [1] M. P. Desjarlais, J. D. Kress, and L. A. Collins; Phys. Rev. B 66, 025401 (2002). [2] E. Schwegler, et al. Phys. Rev. Lett. 87, 265501 (2001). [3] P.M. Celliers, et. al. Physics of Plasmas 11, L41 (2004).
This paper describes the photonic component development, which exploits pioneering work and unique expertise at Sandia National Laboratories, ARDEC and the Army Research Laboratory by combining key optoelectronic technologies to design and demonstrate components for this fuzing application. The technologies under investigation for the optical fuze design covered in this paper are vertical cavity surface emitting lasers (VECSELs), integrated resonant cavity photodetectors (RCPD), and diffractive micro-optics. The culmination of this work will be low cost, robust, fully integrated, g-hardened components designed suitable for proximity fuzing applications. The use of advanced photonic components will enable replacement of costly assemblies that employ discrete lasers, photodetectors, and bulk optics. The integrated devices will be mass produced and impart huge savings for a variety of Army applications.
A suite of impact experiments was conducted to assess spatial and shot-to-shot variability in dynamic properties of tantalum. Samples had a uniform refined {approx}20 micron grain structure with a strong axisymmetric [111] crystallographic texture. Two experiments performed with sapphire windows (stresses of approximately 7 and 12 GPa) clearly showed elastic-plastic loading and slightly hysteretic unloading behavior. An HEL amplitude of 2.8 GPa (corresponding to Y 1.5 GPa) was observed. Free-surface spall experiments showed clear wave attenuation and spallation phenomena. Here, loading stresses were {approx} 12.5 GPa and various ratios of impactor to target thicknesses were used. Spatial and shot-to-shot variability of the spall strength was {+-} 20%, and of the HEL, {+-} 10%. Experiments conducted with smaller diameter flyer plates clearly showed edge effects in the line and point VISAR records, indicating lateral release speeds of roughly 5 km/s.
Spectroscopic investigations in the visible and near UV are underway to study plasmas present in X-ray radiography diodes during the time of the electron beam propagation. These studies are being performed on the RITS-3 accelerator (5.25 MV and 120 kA) at Sandia National Laboratories using several diode configurations. The proper characterization of the plasmas occurring during the time of the X-ray pulse can lead to a greater understanding of diode behavior and X-ray spot size evolution. By studying these plasmas along with the use of selective dopants, insights into such phenomena as impedance collapse, thermal and non-thermal species behavior, charge and current neutralization, anode and cathode plasma formation and propagation, and beam/foil interactions, can be obtained. Information from line and continuum emission and absorption can give key plasma parameters such as temperatures, densities, charge states, and expansion velocities. This information is important for proper modeling and future predictive capabilities for the design and improvement of flash X-ray radiography diodes. Diagnostics include a gated, intensified multichannel plate camera combined with a 1 meter Czerny-Turner monochromator with a multi-fiber spectral input, allowing for both temporal and spatial resolution. Recent results are presented.
Single layer transition metal sulfides (SLTMS) such as MoS{sub 2}, WS{sub 2}, and ReS{sub 2}, play an important role in catalytic processes such as the hydrofining of petroleum streams, and are involved in at least two of the slurry-catalyst hydroconversion processes that have been proposed for upgrading heavy petroleum feed and other sources of hydrocarbon fuels such as coal and shale oils. Additional promising catalytic applications of the SLTMS are on the horizon. The physical, chemical, and catalytic properties of these materials are reviewed in this report. Also discussed are areas for future research that promise to lead to advanced applications of the SLTMS.
The magnetically immersed (B{sub z}) diode is being investigated as a source for pulsed-power driven flash radiography. Experiments fielding this diode have revealed a limit on its achievable current density on target. Either a small spot produces a low dose, or a high dose is achieved with a large spot. It has been proposed that this limit is due to non-protonic ions liberated from the anode surface and subsequently ionizing to higher states. The three-dimensional particle-in-cell code LSP is used to investigate this proposal. Data from the recent immersed diode experiments conducted on the RITS-3 accelerator are compared to LSP models of the experimental configuration, including the B{sub z} field map. We report on how the non-protonic and protonic ion models compare to data, and proposals for future investigation.
Controlled seeding of perturbations is employed to study the evolution of wire array z-pinch implosion instabilities which strongly impact x-ray production when the 3D plasma stagnates on axis. Wires modulated in radius exhibit locally enhanced magnetic field and imploding bubble formation at discontinuities in wire radius due to the perturbed current path. Wires coated with localized spectroscopic dopants are used to track turbulent material flow. Experiments and MHD modeling offer insight into the behavior of z-pinch instabilities.
Experiments to study the implosion dynamics and radiation characteristics of copper z-pinches have been fielded at the 1 MA Zebra facility and the 20 MA Z facility. The impact of initial load mass, initial load diameter, and nesting of wire arrays on the precursor and the stagnated plasma has been evaluated through spectroscopy, shadowgraphy, and fluence measurements. Plasma parameters extracted from modeling of the time-integrated L-shell spectra indicate the presence of more than one plasma source contributing to the radiation, likely due to non-uniform hot spot x-ray emission or temporal gradients.
Over the last several years, rapid progress has been made evaluating the double-z-pinch indirect-drive, inertial confinement fusion (ICF) high-yield target concept (Hammer et al 1999 Phys. Plasmas 6 2129). We have demonstrated efficient coupling of radiation from two wire-array-driven primary hohlraums to a secondary hohlraum that is large enough to drive a high yield ICF capsule. The secondary hohlraum is irradiated from two sides by z-pinches to produce low odd-mode radiation asymmetry. This double-pinch source is driven from a single electrical power feed (Cuneo et al 2002 Phys. Rev. Lett. 88 215004) on the 20 MA Z accelerator. The double z-pinch has imploded ICF capsules with even-mode radiation symmetry of 3.1 {+-} 1.4% and to high capsule radial convergence ratios of 14-21 (Bennett et al 2002 Phys. Rev. Lett. 89 245002; Bennett et al 2003 Phys. Plasmas 10 3717; Vesey et al 2003 Phys. Plasmas 10 1854). Advances in wire-array physics at 20 MA are improving our understanding of z-pinch power scaling with increasing drive current. Techniques for shaping the z-pinch radiation pulse necessary for low adiabat capsule compression have also been demonstrated.
The evolution of granular shear flow is investigated as a function of height in a split-bottom Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations, we find a transition in the nature of the shear as a characteristic height H* is exceeded. Below H* there is a central stationary core; above H* we observe the onset of additional axial shear associated with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height, while the axial width remains narrow and fixed.
The use of Ion Mobility Spectrometry (IMS) in the Detection of Contraband Sandia researchers use ion mobility spectrometers for trace chemical detection and analysis in a variety of projects and applications. Products developed in recent years based on IMS-technology include explosives detection personnel portals, the Material Area Access (MAA) checkpoint of the future, an explosives detection vehicle portal, hand-held detection systems such as the Hound and Hound II (all 6400), micro-IMS sensors (1700), ordnance detection (2500), and Fourier Transform IMS technology (8700). The emphasis to date has been on explosives detection, but the detection of chemical agents has also been pursued (8100 and 6400). Combining Ion Mobility Spectrometry (IMS) with Mass Spectrometry (MS) is described. The IMS-MS combination overcomes several limitations present in simple IMS systems. Ion mobility alone is insufficient to identify an unknown chemical agent. Collision cross section, upon which mobility is based, is not sufficiently unique or predictable a priori to be able to make a confident peak assignment unless the compounds present are already identified. Molecular mass, on the other hand, is much more readily interpreted and related to compounds. For a given compound, the molecular mass can be determined using a pocket calculator (or in one's head) while a reasonable value of the cross-section might require hours of computation time. Thus a mass spectrum provides chemical specificity and identity not accessible in the mobility spectrum alone. In addition, several advanced mass spectrometric methods, such as tandem MS, have been extensively developed for the purpose of molecular identification. With an appropriate mass spectrometer connected to an ion mobility spectrometer, these advanced identification methods become available, providing greater characterization capability.
Sedimentary basins can increase the magnitude and extend the duration of seismic shaking. This potential for seismic amplification is investigated for Pahrump Valley, Nevada-California. The Pahrump Valley is located approximately 50 km northwest of Las Vegas and 75 km south of the Nevada Test Site. Gravity data suggest that the city of Pahrump sits atop a narrow, approximately 5 km deep sub-basin within the valley. The seismic amplification, or ''site effect'', was investigated using a combination of in situ velocity modeling and comparison of the waveforms and spectra of weak ground motion recorded in the city of Pahrump, Nevada, and those recorded in the nearby mountains. Resulting spectral ratios indicate seismic amplification factors of 3-6 over the deepest portion of Pahrump Valley. This amplification predominantly occurs at 2-2.5 Hz. Amplification over the deep sub-basin is lower than amplification at the sub-basin edge, location of the John Blume and Associates PAHA seismic station, which recorded many underground nuclear tests at the Nevada Test Site. A comprehensive analysis of basin amplification for the city of Pahrump should include 3-D basin modeling, due to the extreme basement topography of the Pahrump Valley.
This report describes results of tests using a laser system designed by Coherent Technologies, Inc., in conjunction with Sandia developed nonlinear optics technology. Test results are described for three different optical parametric oscillators built at Sandia. The report concludes with recommendations for future work.
During the past several years Sandia National Laboratories has carried out proof-of-concept experiments to demonstrate tunable, efficient, high-energy ultraviolet nanosecond light sources for satellite-based ozone DIAL. We designed our UV sources to generate pulse energies ≳ 200 mJ at 10 Hz in the range of 308-320 nm with optical-to-optical efficiency approaching 25%. We use sum-frequency generation to mix the 532 nm second harmonic of Nd:YAG with near-IR light derived from a self-injection-seeded image-rotating nonplanar-ring optical parametric oscillator. Laboratory configurations using extra- and intra-cavity sum-frequency generation were designed and tested, yielding 1064 nm to 320 nm conversion efficiencies of 21% and 23% respectively, with pulse energies of 190 mJ and 70 mJ. These energies and efficiencies require pump depletion in the parametric oscillator of at least 80% and SFG efficiency approaching 60%. While the results reported here fall slightly short of our original goals, we believe UV pulse energies exceeding 250 mJ are possible with additional refinements to our technology. Although the sources tested to date are laboratory prototypes with extensive diagnostics, the core components are compact and mechanically robust and can easily be packaged for satellite deployment.
Accelerated aging of Nylon 6.6 fibers used in parachutes has been conducted by following the tensile strength loss under both thermal-oxidative and 100% relative humidity conditions. Thermal-oxidative studies (air circulating ovens) were performed for time periods of weeks to years at temperatures ranging from 37 °C to 138 °C. Accelerated aging humidity experiments (100% RH) were performed under both an argon atmosphere to examine the 'pure' hydrolysis pathway, and under an oxygen atmosphere (oxygen partial pressure close to that occurring in air) to mimic true aging conditions. As expected the results indicated that degradation caused by humidity is much more important than thermal-oxidative degradation. Surprisingly when both oxygen and humidity were present the rate of degradation was dramatically enhanced relative to humidity aging in the absence of oxygen. This significant and previously unknown phenomena underscores the importance of careful accelerated aging that truly mimics real world storage conditions. Published by Elsevier Ltd.
This report is a review of open literature concerning threats including sabotage and theft related to fissile material transport in Japan. It is intended to aid Japanese officials in the development of a design basis threat. This threat includes the external threats of the terrorist, criminal, and extremist, and the insider threats of the disgruntled employee, the employee forced into cooperation via coercion, the psychotic employee, and the criminal employee. Examination of the external terrorist threat considers Japanese demographics, known terrorist groups in Japan, and the international relations of Japan. Demographically, Japan has a relatively homogenous population, both ethnically and religiously. Japan is a relatively peaceful nation, but its history illustrates that it is not immune to terrorism. It has a history of domestic terrorism and the open literature points to the Red Army, Aum Shinrikyo, Chukaku-Ha, and Seikijuku. Japan supports the United States in its war on terrorism and in Iraq, which may make Japan a target for both international and domestic terrorists. Crime appears to remain low in Japan; however sources note that the foreign crime rate is increasing as the number of foreign nationals in the country increases. Antinuclear groups' recent foci have been nuclear reprocessing technology, transportation of MOX fuel, and possible related nuclear proliferation issues. The insider threat is first defined by the threat of the disgruntled employee. This threat can be determined by studying the history of Japan's employment system, where Keiretsu have provided company stability and lifetime employment. Recent economic difficulties and an increase of corporate crime, due to sole reliability on the honor code, have begun to erode employee loyalty.
Oil caverns at the U.S. Strategic Petroleum Reserve (SPR) are subjected to geothermal heating from the surrounding domal salt. This process raises the temperature of the crude oil from around 75 F upon delivery to SPR to as high as 130 F after decades of storage. While this temperature regime is adequate for long-term storage, it poses challenges for offsite delivery, with warm oil evolving gases that pose handling and safety problems. SPR installed high-capacity oil coolers in the mid-1990's to mitigate the emissions problem by lowering the oil delivery temperature. These heat exchanger units use incoming raw water as the cooling fluid, and operate only during a drawdown event where incoming water displaces the outgoing oil. The design criteria for the heat exchangers are to deliver oil at 100 F or less under all drawdown conditions. Increasing crude oil vapor pressures due in part to methane intrusion in the caverns is threatening to produce sufficient emissions at or near 100 F to cause the cooled oil to violate delivery requirements. This impending problem has initiated discussion and analysis of alternative cooling methods to bring the oil temperature even lower than the original design basis of 100 F. For the study described in this report, two alternative cooling methods were explored: (1) cooling during a limited drawdown, and (2) cooling during a degas operation. Both methods employ the heat exchangers currently in place, and do not require extra equipment. An analysis was run using two heat transfer models, HEATEX, and CaveMan, both developed at Sandia National Laboratories. For cooling during a limited drawdown, the cooling water flowrate through the coolers was varied from 1:1 water:oil to about 3:1, with an increased cooling capacity of about 3-7 F for the test cavern Bryan Mound 108 depending upon seasonal temperature effects. For cooling in conjunction with a degas operation in the winter, cavern oil temperatures for the test cavern Big Hill 102 were cooled sufficiently that the cavern required about 9 years to return to the temperature prior to degas. Upon reviewing these results, the authors recommended to the U.S. Department of Energy that a broader study of the cooling during degas be pursued in order to examine the potential benefits of cooling on all caverns in the current degasification schedule.
In October of 2003 experts involved in various aspects of homeland security from the Pacific region met to engage in a free-wheeling discussion and brainstorming (a 'fest') on the role that technology could play in winning the war on terrorism in the Pacific region. The result was a concise and relatively thorough definition of the terrorism problem in the Pacific region, emphasizing the issues unique to Island nations in the Pacific setting, along with an action plan for developing working demonstrations of advanced technological solutions to these issues. Since PacFest 2003, the maritime dimensions of the international security environment have garnered increased attention and interest. To this end, PacFest 2004 sought to identify gaps and enabling technologies for maritime domain awareness and responsive decision-making in the Asia-Pacific region. The PacFest 2004 participants concluded that the technologies and basic information building blocks exist to create a system that would enable the Pacific region government and private organizations to effectively collaborate and share their capabilities and information concerning maritime security. The proposed solution summarized in this report integrates national environments in real time, thereby enabling effective prevention and first response to natural and terrorist induced disasters through better use of national and regional investments in people, infrastructure, systems, processes and standards.
The report presented below is to appear in ''Electrochemistry at the Nanoscale'', Patrik Schmuki, Ed. Springer-Verlag, (ca. 2005). The history of the LIGA process, used for fabricating dimensional precise structures for microsystem applications, is briefly reviewed, as are the basic elements of the technology. The principal focus however, is on the unique aspects of the electrochemistry of LIGA through-mask metal deposition and the generation of the fine and uniform microstructures necessary to ensure proper functionality of LIGA components. We draw from both previously published work by external researchers in the field as well as from published and unpublished studies from within Sandia.
Sandia National Laboratories, California (SNL/CA) is a government-owned/contractor-operated laboratory. Sandia Corporation, a Lockheed Martin Company, operates the laboratory for the Department of Energy's (DOE) National Nuclear Security Administration. The DOE Sandia Site Office oversees operations at the site, using Sandia Corporation as a management and operating contractor. This Site Environmental Report for 2004 was prepared in accordance with DOE Order 231.1A. The report provides a summary of environmental monitoring information and compliance activities that occurred at SNL/CA during calendar year 2004. General site and environmental program information is also included.
Model-based computer simulations have revolutionized product development in the last 10 to 15 years. Technologies that have existed for many decades or even centuries have been improved with the aid of computer simulations. Everything from low-tech consumer goods such as detergents, lubricants and light bulb filaments to the most advanced high-tech products such as airplane wings, wireless communication technologies and pharmaceuticals is engineered with the aid of computer simulations today. In this paper, we present a framework for describing computational tools and their application within the context of product engineering. We examine a few cases of product development that integrate numerical computer simulations into the development stage. We will discuss how the simulations were integrated into the development process, what features made the simulations useful, the level of knowledge and experience that was necessary to run meaningful simulations and other details of the process. Based on this discussion, recommendations for the incorporation of simulations and computational tools into product development will be made.
The objective of this subtask of the Unconventional Nuclear Warfare Design project was to demonstrate mitigation technologies for radiological material dispersal and to assist planners with incorporation of the technologies into a concept of operations. The High Consequence Assessment and Technology department at Sandia National Laboratories (SNL) has studied aqueous foam's ability to mitigate the effects of an explosively disseminated radiological dispersal device (RDD). These benefits include particle capture of respirable radiological particles, attenuation of blast overpressure, and reduction of plume buoyancy. To better convey the aqueous foam attributes, SNL conducted a study using the Explosive Release Atmospheric Dispersion model, comparing the effects of a mitigated and unmitigated explosive RDD release. Results from this study compared health effects and land contamination between the two scenarios in terms of distances of effect, population exposure, and remediation costs. Incorporating aqueous foam technology, SNL created a conceptual design for a stationary containment area to be located at a facility entrance with equipment that could minimize the effects from the detonation of a vehicle transported RDD. The containment design was evaluated against several criteria, including mitigation ability (both respirable and large fragment particle capture as well as blast overpressure suppression), speed of implementation, cost, simplicity, and required space. A mock-up of the conceptual idea was constructed at SNL's 9920 explosive test site to demonstrate the containment design.
The reaction zone of a diesel fuel jet stabilizes at a location downstream of the fuel injector once the initial autoignition phase is over. This distance is referred to as flame lift-off length. Recent investigations have examined the effects of a wide range of parameters (injection pressure, orifice diameter, and ambient gas temperature, density and oxygen concentration) on lift-off length under quiescent diesel conditions. Many of the experimental trends in lift-off length were in agreement with scaling laws developed for turbulent, premixed flame propagation in gas-jet lifted flames at atmospheric conditions. However, several effects did not correlate with the gas-jet scaling laws, suggesting that other mechanisms could be important to lift-off stabilization at diesel conditions. This paper shows experimental evidence that ignition processes affect diesel lift-off stabilization. Experiments were performed in the same optically-accessible combustion vessel as the previous lift-off research. The experimental results show that the ignition quality of a fuel affects lift-off. Fuels with shorter ignition delays generally produce shorter lift-off lengths. In addition, a cool flame is found upstream of, or near the same axial location as, the quasi-steady lift-off length, indicating that first-stage ignition processes affect lift-off. High-speed chemiluminescence imaging also shows that high-temperature self-ignition occasionally occurs in kernels that are upstream of, and detached from, the high-temperature reaction zone downstream, suggesting that the lift-off stabilization is not by flame propagation into upstream reactants in this instance. Finally, analysis of the previous lift-off length database shows that the time-scale for jet mixing from injector-tip orifice to lift-off length collapses to an Arrhenius-type expression, a common method for describing ignition delay in diesel sprays. This Arrhenius-based lift-off length correlation shows comparable accuracy as a previous power-law fit of the No.2 diesel lift-off length database.
The radiant heat test facility develops test sets providing well-characterized thermal environments, often representing fires. Many of the components and procedures have become standardized to such an extent that the development of a specialized design tool was appropriate. SPLASH (Single Panel Lamp and Shroud Helper) is that tool. SPLASH is implemented as a user-friendly program that allows a designer to describe a test setup in terms of parameters such as lamp number, power, position, and separation distance. Thermal radiation is the dominant mechanism of heat transfer and the SPLASH model solves a radiation enclosure problem to estimate temperature distributions in a shroud providing the boundary condition of interest. Irradiance distribution on a specified viewing plane is also estimated. This document provides the theoretical development for the underlying model. A series of tests were conducted to characterize SPLASH's ability to analyze lamp and shroud systems. The comparison suggests that SPLASH succeeds as a design tool. Simplifications made to keep the model tractable are demonstrated to result in estimates that are only approximately as uncertain as many of the properties and characteristics of the operating environment.
This report describes the methodology, analysis and conclusions of a comparison survey of recycling programs at ten Department of Energy sites including Sandia National Laboratories/New Mexico (SNL/NM). The goal of the survey was to compare SNL/NM's recycling performance with that of other federal facilities, and to identify activities and programs that could be implemented at SNL/NM to improve recycling performance.
The absence of agreed definitions and metrics for supercomputer RAS obscures meaningful discussion of the issues involved, hinders their solution, and increases total system cost. Seeking to foster a common basis for communication about supercomputer RAS, [1] proposed a general system state model, definitions, and measurements based on the SEMI-E10 specification [2] used in the semiconductor manufacturing industry. This document enumerates the platform-specific details necessary to apply that general framework to the Red Storm system at Sandia National Laboratories. Familiarity with [1] is a strong prerequisite for understanding of this document, as is familiarity with the Red Storm RAS subsystem (although to a much lesser degree). Given the current pre-production status of Red Storm, this document does not specify actual policy or practice, but rather proposes a framework by which to measure RAS performance on Red Storm.
The Advanced Simulation and Computing (ASC) Distance Computing (DisCom) Wide Area Network (WAN) is a high performance, long distance network environment that is based on the ubiquitous TCP/IP protocol set. However, the Transmission Control Protocol (TCP) and the algorithms that govern its operation were defined almost two decades ago for a network environment vastly different from the DisCom WAN. In this paper we explore and evaluate possible modifications to TCP that purport to improve TCP performance in environments like the DisCom WAN. We also examine a much newer protocol, SCTP (Stream Control Transmission Protocol) that claims to provide reliable network transport while also implementing multi-streaming, multi-homing capabilities that are appealing in the DisCom high performance network environment. We provide performance comparisons and recommendations for continued development that will lead to network communications protocol implementations capable of supporting the coming ASC Petaflop computing environments.
An experimental study was conducted to measure the mechanical properties of the Thermal Protection System (TPS) materials used for the Space Shuttle. Three types of TPS materials (LI-900, LI-2200, and FRCI-12) were tested in 'in-plane' and 'out-of-plane' orientations. Four types of quasi-static mechanical tests (uniaxial tension, uniaxial compression, uniaxial strain, and shear) were performed under low (10{sup -4} to 10{sup -3}/s) and intermediate (1 to 10/s) strain rate conditions. In addition, split Hopkinson pressure bar tests were conducted to obtain the strength of the materials under a relatively higher strain rate ({approx}10{sup 2} to 10{sup 3}/s) condition. In general, TPS materials have higher strength and higher Young's modulus when tested in 'in-plane' than in 'through-the-thickness' orientation under compressive (unconfined and confined) and tensile stress conditions. In both stress conditions, the strength of the material increases as the strain rate increases. The rate of increase in LI-900 is relatively small compared to those for the other two TPS materials tested in this study. But, the Young's modulus appears to be insensitive to the different strain rates applied. The FRCI-12 material, designed to replace the heavier LI-2200, showed higher strengths under tensile and shear stress conditions. But, under a compressive stress condition, LI-2200 showed higher strength than FRCI-12. As far as the modulus is concerned, LI-2200 has higher Young's modulus both in compression and in tension. The shear modulus of FRCI-12 and LI-2200 fell in the same range.
An unavoidable by-product of a metallic structure's use is the appearance of crack and corrosion flaws. Economic barriers to the replacement of these structures have created an aging infrastructure and placed even greater demands on efficient and safe repair methods. In the past decade, an advanced composite repair technology has made great strides in commercial aviation use. Extensive testing and analysis, through joint programs between the Sandia Labs FAA Airworthiness Assurance Center and the aviation industry, have proven that composite materials can be used to repair damaged aluminum structure. Successful pilot programs have produced flight performance history to establish the durability of bonded composite patches as a permanent repair on commercial aircraft structures. With this foundation in place, this effort is adapting bonded composite repair technology to civil structures. The use of bonded composite doublers has the potential to correct the difficulties associated with current repair techniques and the ability to be applied where there are no rehabilitation options. It promises to be cost-effective with minimal disruption to the users of the structure. This report concludes a study into the application of composite patches on thick steel structures typically used in mining operations. Extreme fatigue, temperature, erosive, and corrosive environments induce an array of equipment damage. The current weld repair techniques for these structures provide a fatigue life that is inferior to that of the original plate. Subsequent cracking must be revisited on a regular basis. The use of composite doublers, which do not have brittle fracture problems such as those inherent in welds, can help extend the structure's fatigue life and reduce the equipment downtime. Two of the main issues for adapting aircraft composite repairs to civil applications are developing an installation technique for carbon steel and accommodating large repairs on extremely thick structures. This study developed and proved an optimum field installation process using specific mechanical and chemical surface preparation techniques coupled with unique, in-situ heating methods. In addition, a comprehensive performance assessment of composite doubler repairs was completed to establish the viability of this technology for large, steel structures. The factors influencing the durability of composite patches in severe field environments were evaluated along with related laminate design issues.
A validation study has been conducted for a turbulence model used to close the temporally filtered Navier Stokes (TFNS) equations. A turbulence model was purposely built to support fire simulations under the Accelerated Strategic Computing (ASC) program. The model was developed so that fire transients could be simulated and it has been implemented in SIERRA/Fuego. The model is validated using helium plume data acquired for the Weapon System Certification Campaign (C6) program in the Fire Laboratory for Model Accreditation and Experiments (FLAME). The helium plume experiments were chosen as the first validation problem for SIERRA/Fuego because they embody the first pair-wise coupling of scalar and momentum fields found in fire plumes. The validation study includes solution verification through grid and time step refinement studies. A formal statistical comparison is used to assess the model uncertainty. The metric uses the centerline vertical velocity of the plume. The results indicate that the simple model is within the 95% confidence interval of the data for elevations greater than 0.4 meters and is never more than twice the confidence interval from the data. The model clearly captures the dominant puffing mode in the fire but under resolves the vorticity field. Grid dependency of the model is noted.
The use of surrogate models to approximate computationally expensive simulation models, e.g., large comprehensive finite element models, is widespread. Applications include surrogate models for design, sensitivity analysis, and/or uncertainty quantification. Typically, a surrogate model is defined by a postulated functional form; values for the surrogate model parameters are estimated using results from a limited number of solutions to the comprehensive model. In general, there may be multiple surrogate models, each defined by possibly a different functional form, consistent with the limited data from the comprehensive model. We refer to each as a candidate surrogate model. Methods are developed and applied to select the optimal surrogate model from the collection of candidate surrogate models. The classical approach is to select the surrogate model that best fits the data provided by the comprehensive model; this technique is independent of the model use and, therefore, may be inappropriate for some applications. The proposed approach applies techniques from decision theory, where postulated utility functions are used to quantify the model use. Two applications are presented to illustrate the methods. These include surrogate model selection for the purpose of: (1) estimating the minimum of a deterministic function, and (2) the design under uncertainty of a physical system.
The conclusions of this paper are: (1) Adsorption/desorption on bulk unmodified zeolites showed isoprene adsorbed by zeolite-L and n-pentane adsorbed by zeolite-Y and ZSM-5; (2) Bulk carbonization is used to passivate zeolite activity toward organic adsorption/decomposition; (3) Based on the bulk modified zeolite separation results, we have determined that the MFI type has the most potential for isoprene enrichment; (4) Modified MFI type membrane are jointly made by Sandia and the Univ. of Colorado. Separation experiments are performed by Goodyear Chemical; (5) Isoprene/n-pentane separations have been demonstrated by using both zeolite membranes and modified bulk zeolites at various temperatures on the Goodyear Pilot-scale unit; and (6) Target zeolite membrane separations values of 6.7% isoprene enrichment have been established by economic analysis calculations by Burns & McDonnell.
We consider four asynchronous parallel algorithms for Implicit Monte Carlo (IMC) thermal radiation transport on spatially decomposed meshes. Two of the algorithms are from the production codes KULL from Lawrence Livermore National Laboratory and Milagro from Los Alamos National Laboratory. Improved versions of each of the existing algorithms are also presented. All algorithms were analyzed in an implementation of the KULL IMC package in ALEGRA, a Sandia National Laboratory high energy density physics code. The improved Milagro algorithm performed the best by scaling almost linearly out to 244 processors for well load balanced problems.
The success of Lagrangian contact modeling leads one to believe that important aspects of this capability may be used for multi-material modeling when only a portion of the simulation can be represented in a Lagrangian frame. We review current experience with two dual mesh technologies where one of these meshes is a Lagrangian mesh and the other is an Arbitrary Lagrangian/Eulerian (ALE) mesh. These methods are cast in the framework of an operator-split ALE algorithm where a Lagrangian step is followed by a remesh/remap step. An interface-coupled methodology is considered first. This technique is applicable to problems involving contact between materials of dissimilar compliance. The technique models the more compliant (soft) material as ALE while the less compliant (hard) material and associated interface are modeled in a Lagrangian fashion. Loads are transferred between the hard and soft materials via explicit transient dynamics contact algorithms. The use of these contact algorithms remove the requirement of node-tonode matching at the soft-hard interface. In the context of the operator-split ALE algorithm, a single Lagrangian step is performed using a mesh to mesh contact algorithm. At the end of the Lagrangian step the meshes will be slightly offset at the interface but non-interpenetrating. The ALE mesh nodes at the interface are then remeshed to their initial location relative to the Lagrangian body faces and the ALE mesh is smoothed, translated and rotated to follow Lagrangian body. Robust remeshing in the ALE region is required for success of this algorithm, and we describe current work in this area. The second method is an overlapping grid methodology that requires mapping of information between a Lagrangian mesh and an ALE mesh. The Lagrangian mesh describes a relatively hard body that interacts with softer material contained in the ALE mesh. A predicted solution for the velocity field is performed independently on both meshes. Element-centered velocity and momentum are transferred between the meshes using the volume transfer capability implemented in contact algorithms. Data from the ALE mesh is mapped to a phantom mesh that surrounds the Lagrangian mesh, providing for the reaction to the predicted motion of the Lagrangian material. Data from the Lagrangian mesh is mapped directly to the ALE mesh. A momentum balance is performed on both meshes to adjust the velocity field to account for the interaction of the material from the other mesh. Subsequent, remeshing and remapping of the ALE mesh is performed to allow large deformation of the softer material. We overview current progress using this approach and discuss avenues for future research and development.
Study of the triggered plasma opening switch (TPOS) characteristics is in progress via an ion current collection diagnostic (ICCD), in addition to offline apparatus. This initial ion current collection diagnostic has been designed, fabricated, and tested on the TPOS in order to explore the opening profile of the main switch. The initial ion current collection device utilizes five collectors which are positioned perpendicularly to the main switch stage in order to collect radially traveling ions. It has been shown through analytical prowess that this specific geometry can be treated as a planar case of the Child-Langmuir law with only a 6% deviation from the cylindrical case. Additionally, magnetostatic simulations with self consistent space charge emitting surfaces of the main switch using the Trak code are under way. It is hoped that the simulations will provide evidence in support of both the analytical derivations and experimental data. Finally, an improved design of the ICCD (containing 12 collectors in the axial direction) is presently being implemented.
We provide a common framework for compatible discretizations using algebraic topology to guide our analysis. The main concept is the natural inner product on cochains, which induces a combinatorial Hodge theory. The framework comprises of mutually consistent operations of differentiation and integration, has a discrete Stokes theorem, and preserves the invariants of the DeRham cohomology groups. The latter allows for an elementary calculation of the kernel of the discrete Laplacian. Our framework provides an abstraction that includes examples of compatible finite element, finite volume and finite difference methods. We describe how these methods result from the choice of a reconstruction operator and when they are equivalent.