Publications

The chemical and structural changes within thin films of (3-glycidoxypropyl)trimethoxysilane (GPS) after exposure for various periods of time to air saturated with either D 2O or H 2O at 80°C were studied. The X-ray and neutron reflectivity (XR and NR), combined wuth attenuated total reflection infrared (ATR-IR) spectroscopy were used. The chemical degradation mechanism was identified by IR as hydrolysis of siloxane bonds. GPS films were prepared by dip-coating, which resulted in a greater and more variable thickness than for the spin-coated samples, for ATR-IR. The little changes in the reflectivity data was observed for films conditioned with D 2O at 80°C for one month.

Identifying the areas to be protected is an important part of the development of measures for physical protection against sabotage at complex nuclear facilities. In June 1999, the International Atomic Energy Agency published INFCIRC/225/Rev.4, 'The Physical Protection of Nuclear Material and Nuclear Facilities.' This guidance recommends that 'Safety specialists, in close cooperation with physical protection specialists, should evaluate the consequences of malevolent acts, considered in the context of the State's design basis threat, to identify nuclear material, or the minimum complement of equipment, systems or devices to be protected against sabotage.' This report presents a structured, transparent approach for identifying the areas that contain this minimum complement of equipment, systems, and devices to be protected against sabotage that is applicable to complex nuclear facilities. The method builds upon safety analyses to develop sabotage fault trees that reflect sabotage scenarios that could cause unacceptable radiological consequences. The sabotage actions represented in the fault trees are linked to the areas from which they can be accomplished. The fault tree is then transformed (by negation) into its dual, the protection location tree, which reflects the sabotage actions that must be prevented in order to prevent unacceptable radiological consequences. The minimum path sets of this fault tree dual yield, through the area linkage, sets of areas, each of which contains nuclear material, or a minimum complement of equipment, systems or devices that, if protected, will prevent sabotage. This method also provides guidance for the selection of the minimum path set that permits optimization of the trade-offs among physical protection effectiveness, safety impact, cost and operational impact.

Most attempts to quantify differential diffusion (DD) are based on the difference between different definitions of the mixture fraction. This paper presents a general method for evaluating differential diffusion in premixed or nonpremixed systems based on conservation equations for the elemental mass fractions. These measures form a basis for analysing differential diffusion. Casting these in terms of a mixture fraction gives particular insight into differential diffusion for nonpremixed systems, and provides a single DD measure. Furthermore, it allows direct evaluation of the validity of the traditional assumptions involved in writing a mixture fraction transport equation. Results are presented for one-dimensional opposed flow simulations of hydrogen and methane flames as well as direct numerical simulations (DNS) of CH4/H2-air and CO/H2-air flames. For a common definition of the mixture fraction, the DD measure can be approximated well by considering only the contribution of H2 and CH4 in methane-air flames. Differential diffusion is largely driven by production of H2 in the flame zone for hydrocarbon flames. Effects of strain rate and filter width on the relative importance of differential diffusion are examined. © 2005 Taylor & Francis Group Ltd.

Rapid microchip reversed-phase HPLC of peptides and proteins at pressure gradients of 12 bar/cm (180 psi/cm) has been performed using a microdevice that integrates subnanoliter on-chip injection and separation with a miniaturized fluorescence detector. Proteins and peptides were separated on a C18 side-chain porous polymer monolith defined by contact lithography, and injection was achieved via a pressure-switchable fluoropolymer valve defined using projection lithography. Preliminary separations of peptide standards and protein mixtures were performed in 40-200 s, and switching between samples with no detectible sample carryover has been performed. The injections and separations were reproducible; the relative standard deviation (RSD) for retention time was 0.03%, and peak area RSD was 3.8%. Sample volumes ranging from 220 to 800 pL could be linearly metered by controlling the pressure injection pulse duration with conventional timing and valving. The current prototype system shows the potential for rapid and autonomous HPLC separations with varying modalities and the potential for direct connection to mass spectrometers at nanospray flow rates. © 2005 American Chemical Society.

We have determined limits on the cross section for both electronically nonadiabatic excitation and quenching in the Cl (Pj2) + D2 system. Our experiment incorporates crossed-molecular-beam scattering with state-selective Cl (P 12,32 2) detection and velocity-mapped ion imaging. By colliding atomic chlorine with D2, we address the propensity for collisions that result in a change of the spin-orbit level of atomic chlorine either through electronically nonadiabatic spin-orbit excitation Cl (P 32 2) + D2 → Cl* (P 12 2) + D2 or through electronically nonadiabatic spin-orbit quenching Cl* (P 12 2) + D2 →Cl (P 32 2) + D2. In the first part of this report, we estimate an upper limit for the electronically nonadiabatic spin-orbit excitation cross section at a collision energy of 5.3 kcalmol, which lies above the energy of the reaction barrier (4.9 kcalmol). Our analysis and simulation of the experimental data determine an upper limit for the excitation cross section as σNA ≤0.012 Å2. In the second part of this paper we investigate the propensity for electronically nonadiabatic spin-orbit quenching of Cl * following a collision with D2 or He. We perform these experiments at collision energies above and below the energy of the reaction barrier. By comparing the amount of scattered Cl* in our images to the amount of Cl* lost from the atomic beam we obtain the maximum cross section for electronically nonadiabatic quenching as σNA ≤ 15 -15 +44 Å2 for a collision energy of 7.6 kcalmol. Our experiments show the probability for electronically nonadiabatic quenching in Cl* + D2 to be indistinguishable to that for the kinematically identical system of Cl* +He. © 2005 American Institute of Physics.

This report is the latest in a continuing series that highlights the recent technical accomplishments associated with the work being performed within the Materials and Process Sciences Center. Our research and development activities primarily address the materials-engineering needs of Sandia's Nuclear-Weapons (NW) program. In addition, we have significant efforts that support programs managed by the other laboratory business units. Our wide range of activities occurs within six thematic areas: Materials Aging and Reliability, Scientifically Engineered Materials, Materials Processing, Materials Characterization, Materials for Microsystems, and Materials Modeling and Simulation. We believe these highlights collectively demonstrate the importance that a strong materials-science base has on the ultimate success of the NW program and the overall DOE technology portfolio.

Iron, the most ubiquitous of the transition metals and the fourth most plentiful element in the Earth's crust, is the structural backbone of our modern infrastructure. It is therefore ironic that as a nanoparticle, iron has been somewhat neglected in favor of its own oxides, as well as other metals such as cobalt, nickel, gold, and platinum. This is unfortunate, but understandable. Iron 's reactivity is important in macroscopic applications (particularly rusting), but is a dominant concern at the nanoscale. Finely divided iron has long been known to be pyrophoric, which is a major reason that iron nanoparticles have not been more fully studied to date. This extreme reactivity has traditionally made iron nanoparticles difficult to study and inconvenient for practical applications. Iron however has a great deal to offer at the nanoscale, including very potent magnetic and catalytic properties. Recent work has begun to take advantage of iron's potential, and work in this field appears to be blossoming. © 2005 Wiley-VCH Verlag GmbH & Co. KGaA.

This Final Report on the Monitoring Phase of the former Weeks Island Strategic Petroleum Reserve crude oil storage facility details the results of five years of monitoring of various surface accessible quantities at the decommissioned facility. The Weeks Island mine was authorized by the State of Louisiana as a Strategic Petroleum Reserve oil storage facility from 1979 until decommissioning of the facility in 1999. Discovery of a sinkhole over the facility in 1992 with freshwater inflow to the facility threatened the integrity of the oil storage and led to the decision to remove the oil, fill the chambers with brine, and decommission the facility. Thereafter, a monitoring phase, by agreement between the Department of Energy and the State, addressed facility stability and environmental concerns. Monitoring of the surface ground water and the brine of the underground chambers from the East Fill Hole produced no evidence of hydrocarbon contamination, which suggests that any unrecovered oil remaining in the underground chambers has been contained. Ever diminishing progression of the initial major sinkhole, and a subsequent minor sinkhole, with time was verification of the response of sinkholes to filling of the facility with brine. Brine filling of the facility ostensively eliminates any further growth or new formation from freshwater inflow. Continued monitoring of sinkhole response, together with continued surface surveillance for environmental problems, confirmed the intended results of brine pressurization. Surface subsidence measurements over the mine continued throughout the monitoring phase. And finally, the outward flow of brine was monitored as a measure of the creep closure of the mine chambers. Results of each of these monitoring activities are presented, with their correlation toward assuring the stability and environmental security of the decommissioned facility. The results suggest that the decommissioning was successful and no contamination of the surface environment by crude oil has been found.

Goal of this study was to develop and characterize novel polymeric materials as pseudostationary phases in electrokinetic chromatography. Fundamental studies have characterized the chromatographic selectivity of the materials as a function of chemical structure and molecular conformation. The selectivities of the polymers has been studied extensively, resulting in a large body of fundamental knowledge regarding the performance and selectivity of polymeric pseudostationary phases. Two polymers have also been used for amino acid and peptide separations, and with laser induced fluorescence detection. The polymers performed well for the separation of derivatized amino acids, and provided some significant differences in selectivity relative to a commonly used micellar pseudostationary phase. The polymers did not perform well for peptide separations. The polymers were compatible with laser induced fluorescence detection, indicating that they should also be compatible with chip-based separations.

The recent interests in developing multiscale model-based simulation procedures have brought about the challenging tasks of bridging different spatial and temporal scales within a unified framework. However, the research focus has been on the scale effect in the spatial domain with the loading rate being assumed to be quasi-static. Although material properties are rate-dependent in nature, little has been done in understanding combined loading rate and specimen size effects on the material properties at different scales. In addition, the length and time scales that can be probed by the molecular level simulations are still fairly limited due to the limitation of computational capability. Based on the experimental and computational capabilities available, therefore, an attempt is made in this report to formulate a hyper-surface in both spatial and temporal domains to predict combined size and rate effects on the mechanical properties of engineering materials. To demonstrate the features of the proposed hyper-surface, tungsten specimens of various sizes under various loading rates are considered with a focus on the uniaxial loading path. The mechanical responses of tungsten specimens under other loading paths are also explored to better understand the size effect. It appears from the preliminary results that the proposed procedure might provide an effective means to bridge different spatial and temporal scales in a unified multiscale modeling framework, and facilitate the application of nanoscale research results to engineering practice.

Power sources capable of supplying tens of watts are needed for a wide variety of applications including portable electronics, sensors, micro aerial vehicles, and mini-robotics systems. The utility of these devices is often limited by the energy and power density capabilities of batteries. A small combustion engine using liquid hydrocarbon fuel could potentially increase both power and energy density by an order of magnitude or more. This report describes initial development work on a meso-scale external combustion engine based on the Stirling cycle. Although other engine designs perform better at macro-scales, we believe the Stirling engine cycle is better suited to small-scale applications. The ideal Stirling cycle requires efficient heat transfer. Consequently, unlike other thermodynamic cycles, the high heat transfer rates that are inherent with miniature devices are an advantage for the Stirling cycle. Furthermore, since the Stirling engine uses external combustion, the combustor and engine can be scaled and optimized semi-independently. Continuous combustion minimizes issues with flame initiation and propagation. It also allows consideration of a variety of techniques to promote combustion that would be difficult in a miniature internal combustion engine. The project included design and fabrication of both the engine and the combustor. Two engine designs were developed. The first used a cylindrical piston design fabricated with conventional machining processes. The second design, based on the Wankel rotor geometry, was fabricated by through-mold electroforming of nickel in SU8 and LIGA micromolds. These technologies provided the requisite precision and tight tolerances needed for efficient micro-engine operation. Electroformed nickel is ideal for micro-engine applications because of its high strength and ductility. A rotary geometry was chosen because its planar geometry was more compatible with the fabrication process. SU8 lithography provided rapid prototypes to verify the design. A final high precision engine was created via LIGA. The micro-combustor was based on an excess enthalpy concept. Development of a micro-combustor included both modeling and experiments. We developed a suite of simulation tools both in support of the design of the prototype combustors, and to investigate more fundamental aspects of combustion at small scales. Issues of heat management and integration with the micro-scale Stirling engine were pursued using CFD simulations. We found that by choice of the operating conditions and channel dimensions energy conversion occurs by catalysis-dominated or catalysis-then-homogeneous phase combustion. The purpose of the experimental effort in micro-combustion was to study the feasibility and explore the design parameters of excess enthalpy combustors. The efforts were guided by the necessity for a practical device that could be implemented in a miniature power generator, or as a stand-alone device used for heat generation. Several devices were fabricated and successfully tested using methane as the fuel.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We have synthesized defect-free aluminosilicate and silicalite zeolite thin films supported on commercially available alpha and gamma alumina disk substrates. We have also built a permeation unit that can test both pure and mixed gases from room temperature to 250 C. Results indicate fluxes on the order of 10{sup -6} to 10{sup -7} mole/(m{sup 2}Pa sec) and excellent separation values for H{sub 2} or CO{sub 2}. For the Al/Si membrane: H{sub 2}/N{sub 2} {ge} 61, H{sub 2}/CO{sub 2} {ge} 80, H{sub 2}/CH{sub 4} = 7, CH{sub 4}/CO{sub 2} {ge} 11; for the TPA/Si membrane: H{sub 2}/N{sub 2} {ge} 61, H{sub 2}/CO{sub 2} {ge} 80, H{sub 2}/CH{sub 4} = 7, CH{sub 4}/CO{sub 2} {ge} 11. Our data show that we can use the adsorption ability plus the effective pore diameter of the zeolite to 'tune' the selectivity of the membrane. Another avenue of research is into bulk novel molecular sieve materials, with the goal of 'tuning' pore sizes to molecular sieving needs. A novel crystalline 12-ring microporous gallophosphate material is described.

Abstract not provided.

This paper is about making reversible logic a reality for supercomputing. Reversible logic offers a way to exceed certain basic limits on the performance of computers, yet a powerful case will have to be made to justify its substantial development expense. This paper explores the limits of current, irreversible logic for supercomputers, thus forming a threshold above which reversible logic is the only solution. Problems above this threshold are discussed, with the science and mitigation of global warming being discussed in detail. To further develop the idea of using reversible logic in supercomputing, a design for a 1 Zettaflops supercomputer as required for addressing global climate warming is presented. However, to create such a design requires deviations from the mainstream of both the software for climate simulation and research directions of reversible logic. These deviations provide direction on how to make reversible logic practical.

Abstract not provided.

Binding of porphyrin to murine ferrochelatase, the terminal enzyme of the heme biosynthetic pathway, is investigated by employing a set of variants harboring mutations in a putative porphyrin-binding loop. Using resonance Raman (RR) spectroscopy, the structural properties of the ferrochelatase-bound porphyrins are examined, especially with respect to the porphyrin deformation occurring in the environment of the active site. This deformation is thought to be a key step in the enzymatic insertion of ferrous iron into the porphyrin ring to make heme. Our previous RR spectroscopic studies of binding of porphyrin to murine ferrochelatase led us to propose that the wild-type enzyme induces porphyrin distortion even in the absence of the metal ion substrate. Here, we broaden this view by presenting evidence that the degree of a specific nonplanar porphyrin deformation contributes to the catalytic efficiency of ferrochelatase and its variants. The results also suggest that the conserved Trp256 (murine ferrochelatase numbering) is partially responsible for the observed porphyrin deformation. Binding of porphyrin to the ferrochelatase variants causes a decrease in the intensity of RR out-of-plane vibrational mode {gamma}{sub 15}, a saddling-like mode that is strong in the wild-type enzyme. In particular, the variant with a catalytic efficiency 1 order of magnitude lower than that of the wild-type enzyme is estimated to produce less than 30% of the wild-type saddling deformation. These results suggest that specific conserved loop residues (especially Trp256) are directly involved in the saddling of the porphyrin substrate.

Abstract not provided.

In contrast to traditional terascale simulations that have known, fixed data inputs, dynamic data-driven (DDD) applications are characterized by unknown data and informed by dynamic observations. DDD simulations give rise to inverse problems of determining unknown data from sparse observations. The main difficulty is that the optimality system is a boundary value problem in 4D space-time, even though the forward simulation is an initial value problem. We construct special-purpose parallel multigrid algorithms that exploit the spectral structure of the inverse operator. Experiments on problems of localizing airborne contaminant release from sparse observations in a regional atmospheric transport model demonstrate that 17-million-parameter inversion can be effected at a cost of just 18 forward simulations with high parallel efficiency. On 1024 Alphaserver EV68 processors, the turnaround time is just 29 minutes. Moreover, inverse problems with 135 million parameters - corresponding to 139 billion total space-time unknowns - are solved in less than 5 hours on the same number of processors. These results suggest that ultra-high resolution data-driven inversion can be carried out sufficiently rapidly for simulation-based 'real-time' hazard assessment.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Accurate and efficient turbulence simulation in complex geometries is a formidable chal-lenge. Traditional methods are often limited by low accuracy and/or restrictions to simplegeometries. We explore the merger of Discontinuous Galerkin (DG) spatial discretizationswith Variational Multi-Scale (VMS) modeling, termed Local VMS (LVMS), to overcomethese limitations. DG spatial discretizations support arbitrarily high-order accuracy on un-structured grids amenable for complex geometries. Furthermore, high-order, hierarchicalrepresentation within DG provides a natural framework fora prioriscale separation crucialfor VMS implementation. We show that the combined benefits of DG and VMS within theLVMS method leads to promising new approach to LES for use in complex geometries.The efficacy of LVMS for turbulence simulation is assessed by application to fully-developed turbulent channelflow. First, a detailed spatial resolution study is undertakento record the effects of the DG discretization on turbulence statistics. Here, the localhp[?]refinement capabilites of DG are exploited to obtain reliable low-order statistics effi-ciently. Likewise, resolution guidelines for simulating wall-bounded turbulence using DGare established. We also explore the influence of enforcing Dirichlet boundary conditionsindirectly through numericalfluxes in DG which allows the solution to jump (slip) at thechannel walls. These jumps are effective in simulating the influence of the wall commen-surate with the local resolution and this feature of DG is effective in mitigating near-wallresolution requirements. In particular, we show that by locally modifying the numericalviscousflux used at the wall, we are able to regulate the near-wall slip through a penaltythat leads to improved shear-stress predictions. This work, demonstrates the potential ofthe numerical viscousflux to act as a numerically consistent wall-model and this successwarrents future research.As in any high-order numerical method some mechanism is required to control aliasingeffects due to nonlinear interactions and to ensure nonlinear stability of the method. Inthis context, we evaluate the merits of two approaches to de-aliasing -- spectralfilteringand polynomial dealiasing. While both approaches are successful, polynomial-dealiasingis found to be better suited for use in large-eddy simulation. Finally, results using LVMSare reported and show good agreement with reference direct numerical simulation therebydemonstrating the effectiveness of LVMS for wall-bounded turbulence. This success pavesthe way for future applications of LVMS to more complex turbulentflows.3