Publications

An oxidation treatment, often termed "pre-oxidation", is performed on austenitic stainless steel prior to glass/metal joining to produce hermetic seals. The resulting thin oxide acts as a transitional layer and a source of Cr and other elements which diffuse into the glass during the subsequent bonding process. Pre-oxidation is performed in a low pO 2 atmosphere to avoid iron oxide formation and the final oxide is composed of Cr 2O 3, MnCr 2O 4 spinel, and SiO 2. Significant heat-to-heat variations in the oxidation behavior of 304L stainless steel have been observed, which result in inconsistent glass/metal seal behavior. The objectives of this work were to characterize the oxidation kinetics, the oxide morphology and composition, and the stainless steel attributes that lead to robust glass/metal seals. The oxidation kinetics were determined by thermogravimetric (TG) analysis and the oxide layers were characterized using metallography, SEM, focused ion beam (FIB) analysis, and image analysis. The results show that poor sealing behavior is associated with slower oxidation kinetics and a more continuous layer of SiO 2 at the metal/oxide interface. In addition, the effects of 304L heat composition on oxidation behavior will be discussed. Copyright © 2006 ASM International®.

The solderability of an immersion Ag finish was evaluated after the exposure of test specimens to a Battelle Class II environment, which accelerates the storage conditions of light industrial surroundings. The solderability metric was the contact angle, (θC), as determined by the meniscometer/wetting balance technique. Auger surface and depth profile analyses were utilized to identify changes in the coating chemistry. The solderability test results indicate that there was no appreciable loss in solderability when the immersion Ag coated coupons were packaged in vapor phase corrosion (VPC) inhibitor bags and/or inhibitor bags with VPC inhibitor paper and aged for 8 hours, 1 week or 2 weeks in the Battelle Class II environment. An increase in surface carbon concentration after aging did not appear to significantly affect solderability. Copyright © 2006 ASM International®.

The properties of energetic thin films considered for alternative braze[1] techniques are investigated. Vapor-deposited Ni/Ti multilayer foils having a net 1:1 stoichiometry exhibit self-propagating, high temperature combustion reactions. The rate of reaction depends on Ni/Ti multilayer design with steady-state propagation speeds of freestanding foils measured from 0.2 to 1.0m/s. Transmission electron microscopy and x-ray diffraction further show that NiTi films reacted in a self-propagating mode have a fine-grain, polycrystalline microstructure. All films are composed of cubic B2 and monoclinic B19' phases with some having NiTi2 or Ni3Ti precipitates. Copyright © 2006 ASM International®.

Simulations of a low-speed square cylinder wake and a supersonic axisymmetric base wake are performed using the detached eddy simulation model. A reduced-dissipation form of a shock-capturing flux scheme is employed to mitigate the effects of dissipative error in regions of smooth flow. The reduced-dissipation scheme is demonstrated on a two-dimensional square cylinder wake problem, showing a marked improvement in accuracy for a given grid resolution. The results for simulations on three grids of increasing resolution for the three-dimensional square cylinder wake are compared with experimental data and to other computational studies. The comparisons of mean flow and global flow quantities to experimental data are favorable, whereas the results for second order statistics hi the wake are mixed and do not always improve with increasing spatial resolution. Comparisons to large eddy simulation are also generally favorable, suggesting detached eddy simulation provides an adequate subgrid scale model. Predictions of base drag and centerline wake velocity for the supersonic wake are also good, given sufficient grid refinement. These cases add to the validation library for detached eddy simulation and support its use as an engineering analysis tool for accurate prediction of global flow quantities and mean flow properties.

As the economic and environmental impact arguments for increasing the use of nuclear energy for electricity generation and hydrogen production strengthen, it becomes important to better understand human biases, critical thinking skills, and individual specific characteristics that influence decisions made during probabilistic safety assessments (PSAs), decisions regarding nuclear energy among the general public (e.g., trust of risk assessments, acceptance of new plants, etc.), and nuclear energy decisions made by high-level decision makers (e.g., energy policy makers & government regulators). To promote increased understanding and hopefully to improve decision making capacities, this paper provides four key elements. The foundation of these elements builds on decades of research and associated experimental data regarding risk perception and decision making. The first element is a unique taxonomy of twenty-six recognized biases. Examples of biases were generated by reviewing the relevant literature in nuclear safety, cognitive psychology, economics, science education, and neural science (to name a few) and customizing superficial elements of those examples to the nuclear energy domain. The second element is a listing of ten critical thinking skills (with precise definitions) applicable to risk perception and decision making. Third, three brief hypothetical decision making examples are presented and decomposed relative to the unique, decision making bias framework and critical thinking set. The fourth element is a briefly outlined strategy which may enable one to make better decisions in domains that demand careful reflection and strong adherence to the best available data (i.e., avoiding 'unhelpful biases' that conflict with proper interpretation of available data). The elements concisely summarized in this paper (and additional elements) are available in detail in an unclassified, unlimited release Sandia National Laboratories report (SAND2005-5730). The proposed taxonomy of biases contains the headings of normative knowledge, availability, and individual specific biases. Normative knowledge involves a person's skills in combinatorics, probability theory, and statistics. Research has shown that training and experience in these quantitative fields can improve one's ability to accurately determine event likelihoods. Those trained in statistics tend to seek appropriate data sources when assessing the frequency and severity of an event. The availability category of biases includes those which result from the structure of human cognitive machinery. Two examples of biases in the availability category include the anchoring bias and the retrievability bias. The anchoring bias causes a decision maker to bias subsequent values or items toward the first value or item presented to them. The retrievability bias refers to the bias that drives people to believe those values or items which are easier to retrieve from memory are more likely to occur. Individual specific biases include a particular person's values, personality, interests, group identity, and substantive knowledge (i.e., specific domain knowledge related to the decision to be made). Critical thinking skills are also offered as foundational for competent risk perception and decision making as they can mute the impact of undesirable biases, regulate the application of one's knowledge to a decision, and guide information gathering activities. The list of critical thinking skills presented here was originally articulated by the late Arnold B. Arons, a distinguished physicist and esteemed researcher of learning processes. Finally, in addition to borrowing insights from the literature domains mentioned above, the formal decision making approach supported in this paper incorporates methods used in multi-attribute utility theory. © 2006 by ASME.

The objectives of this study are to identify and evaluate alternative protective action recommendations (PARs) that could reduce dose to the public during a radiological emergency and to determine whether improvements or changes to the federal guidance would be beneficial. The emergency response strategies considered in this study are the following: (1) standard radial evacuation; (2) shelter-in-place followed by radial evacuation; (3) shelter-in-place followed by lateral evacuation; (4) preferential sheltering followed by radial evacuation; and (5) preferential sheltering followed by lateral evacuation. Radial evacuation is directly away from the plant; lateral evacuation is azimuthally (around the compass) away from the direction of the wind. Shelter-in-place is a protective action strategy in which individuals remain in their residence, place of work, or other facility at the time that a general warning is given. Preferential sheltering involves moving individuals to nearby, large buildings, e.g., high-school gymnasiums or courthouses, that afford greater protection than personal residences. This study shows that there are benefits to sheltering if followed by lateral evacuation. However, if the lateral evacuation strategy cannot be implemented, then early radial evacuation is often preferable. The most appropriate PAR depends on the evacuation time estimate (ETE) and, therefore, it is desirable to reduce the uncertainty associated with the ETE. Nuclear Regulatory Commission (NRC) guidance to commercial power plants currently allows for sheltering and/or evacuation as an emergency response to a serious nuclear accident. Frequently, however, licensees and states default to evacuation strategies and do not consider sheltering. Here we evaluate several alternative strategies to determine if standard radial evacuation is best or if other options could reduce the overall risk to the public. We consider two source terms based on the NUREG-1150 study. These involve a rapid release of radioactive material into the atmosphere and a more gradual release. Two variations in timing have also been investigated, but are not reported here. The evaluation was performed for a generic site, which uses a uniform population distribution and typical Midwest meteorological data (Moline, IL). Additional parameters that are varied in the study are the ETE (4-, 6-, 8-, and 10-hour ETEs are considered) and the duration of sheltering (2-, 4-, and 8-hr sheltering periods are considered). Additional sensitivity studies were performed to investigate nonuniform evacuation speed caused by traffic congestion, the time needed to reach a preferential shelter, and the effect of adverse weather conditions as opposed to favorable weather conditions. Adverse weather conditions are those for which precipitation occurs before the leading edge of the plume exits the 10-mile emergency planning zone (EPZ). Ultimately, emergency response strategies are ranked by their potential to reduce adverse health effects for residents within the EPZ. © 2006 by ASME.

In order to better understand how the US natural gas network might respond to disruptions, a model was created that represents the network on a regional basis. Natural gas storage for each region is represented as a stock. Transmission between each region is represented as a flow, as is natural gas production, importation, and consumption. Various disruption scenarios were run to test the robustness of the network. The system as modeled proved robust to a variety of disruption scenarios. However, a weakness of the system is that production shortfalls or interruptions cannot be replaced, and demand must therefore be reduced by the amount of the shortfall.

An oxidation treatment, often termed "pre-oxidation", is performed on austenitic stainless steel prior to joining to alkali barium silicate glass to produce hermetic seals. The resulting thin oxide acts as a transitional layer and a source of Cr and other elements which diffuse into the glass during the subsequent bonding process. Pre-oxidation is performed in a low pO2 atmosphere to avoid iron oxide formation and the final oxide is composed of Cr2O3, MnCr2O4 spinel, and SiO2. Significant heat-to-heat variations in the oxidation behavior of 304L stainless steel have been observed, which result in inconsistent glass-to-metal (GTM) seal behavior. The objectives of this work were to characterize the stainless steel pre-oxidized layer and the glass/oxide/304L interface region after glass sealing. The 304L oxidation kinetics were determined by thermogravimetric (TG) analysis and the glass/metal seals characteristics were studied using sessile drop tests, in which wetting angles were measured and glass adhesion was analyzed. The pre-oxidized layers and glass/metal interface regions were characterized using metallography, focused ion beam (FIB) sectioning, scanning and transmission electron microscopy, and electron probe microanalysis (EPMA). The results show that poor glass sealing behavior is associated with a more continuous layer of SiO 2 at the metal/oxide interface.

Radio Frequency Microelectromechanical Systems (RF MEMS) are advantageous for reconfigurable antennas providing the potential for steering, frequency agility, and tunable directivity. Until RF MEMS switches can consistently reach cycles into the billions (if not trillions), limited reliability outweighs the promised benefits, preventing the deployment of RF MEMS into systems. Understanding failure mechanisms is essential to improving reliability. This paper describes preliminary reliability results and tests conducted in a vacuum chamber to investigate and understand the impact of contamination related failure mechanisms. © 2006 IEEE.

Many systems can be approximated as linear with coefficients that vary periodically with time. For example, an anisotropic shaft rotating at constant speed on anisotropic bearings can be modeled as periodically time varying (PTV). Similar models can be obtained for wind turbines, some mechanisms, etc... However, the vast majority of modal analysis algorithms and techniques apply only to linear time invariant (LTI) systems. In this paper, two methods are demonstrated by which the free response of a periodically time varying system can be exactly parameterized by an LTI system. The parameters of the LTI representation can then be identified using standard techniques. The analysis techniques are demonstrated on a simple system, representing a rotor mounted on an anisotropic, flexible shaft, supported by anisotropic bearings. They are then applied to synthetic response data for a system with parameters that vary only weakly with time, as might be encountered when attempting to detect small cracks in a rotating shaft. These examples demonstrate the methods' ability to characterize the anisotropy of the shaft, even when both the shaft and supports are anisotropic.

This paper presents a new method for the analysis of sinusoidal aquifer test data that extends current analytical techniques to include the analysis of pressure responses in multi-layer aquifers. The multi-layer sinusoidal solution is applied to data from a series of sinusoidal aquifer tests conducted in the two-layered Culebra dolomite at the Waste Isolation Pilot Plant near Carlsbad, NM. In addition to the multilayered solution, a summary of available sinusoidal solutions is provided along with tools to pre-process sinusoidal response data.

Gaussian processes are used as emulators for expensive computer simulations. Recently, Gaussian processes have also been used to model the "error field" or "code discrepancy" between a computer simulation code and experimental data, and the delta term between two levels of computer simulation (multi-fidelity codes). This work presents the use of Gaussian process models to approximate error or delta fields, and examines how one calculates the parameters governing the process. In multi-fidelity modeling, the delta term is used to correct a lower fidelity model to match or approximate a higher fidelity model. The terms governing the Gaussian process (e.g., the parameters of the covariance matrix) are updated using a Bayesian approach. We have found that use of Gaussian process models requires a good understanding of the method itself and an understanding of the problem in enough detail to identify reasonable covariance parameters. The methods are not "black-box" methods that can be used without some statistical understanding. However, Gaussian processes offer the ability to account for uncertainties in prediction. This approach can help reduce the number of high-fidelity function evaluations necessary in multi-fidelity optimization.

Forces generated by a static magnetic field interacting with eddy currents can provide a novel method of vibration damping. This paper discusses an experiment performed to validate modeling [3] for a case where a static magnetic field penetrates a thin sheet of conducting, non-magnetic material. When the thin sheet experiences motion, the penetrating magnetic field generates eddy currents within the sheet. These eddy currents then interact with the static field, creating magnetic forces that act on the sheet, providing damping to the sheet motion. In the presented experiment, the sheet was supported by cantilever springs attached to a frame, then excited with a vibratory shaker. The recorded motions of the sheet and the frame were used to characterize the effect of the eddy current damping.

Multiple references are often used to excite a structure in modal testing programs. This is necessary to excite all the modes and to extract accurate mode shapes when closely spaced roots are present. An algorithm known as SMAC (Synthesize Modes And Correlate), based on principles of modal filtering, has been in development for several years. This extraction technique calculates reciprocal modal vectors based on frequency response function (FRF) measurements. SMAC was developed to accurately extract modes from structures with moderately damped modes and/or high modal density. In the past SMAC has only worked with single reference data. This paper presents an extension of SMAC to work with multiple reference data. If roots are truly perfectly repeated, the mode shapes extracted by any method will be a linear combination of the "true" shapes. However, most closely spaced roots are not perfectly repeated but have some small difference in frequency and/or damping. SMAC exploits these very small differences. The multi-reference capability of SMAC begins with an evaluation of the MMIF (Multivariate Mode Indicator Function) or CMIF (Complex Mode Indicator Function) from the starting frequency list to determine which roots are likely repeated. Several seed roots are scattered in the region of the suspected multiple roots and convergence is obtained. Mode shapes are then created from each of the references individually. The final set of mode shapes are selected based on one of three different selection techniques. Each of these is presented in this paper. SMAC has long included synthesis of FRFs and MIFs from the roots and residues to check extraction quality against the original data, but the capability to include residual effects has been minimal. Its capabilities for including residual vectors to account for out-of-band modes have now been greatly enhanced. The ability to resynthesize FRFs and mode indicator functions from the final mode shapes and residual information has also been developed. Examples are provided utilizing the SMAC package on multi-reference experimental data from two different systems.

A finite element (FE) model of a shell-payload structure is to be used to predict structural dynamic acceleration response to untestable blast environments. To understand the confidence level of these predictions, the model will be validated using test data from a blast tube experiment. The first step in validating the structural response is to validate the loading. A computational fluid dynamics (CFD) code, Saccara, was used to provide the blast tube pressure loading to the FE model. This paper describes the validation of the CFD pressure loading and its uncertainty quantification with respect to experimental pressure data obtained from geometrical mock-up structures instrumented with pressure gages in multiple nominal blast tube tests. A systematic validation approach was used from the uncertainty quantification group at Sandia National Labs. Significant effort was applied to distill the pressure loading to a small number of validation metrics important to obtaining valid final response which is in terms of acceleration shock response spectrum. Uncertainty in the pressure loading amplitude is quantified so that it can be applied to the validation blast tube test on the shell payload structure which has significant acceleration instrumentation but only a few pressure gages.

This paper presents a fictional scenario describing the effects of a category four hurricane on a metropolitan area, accompanied by a challenge: Describe, and eventually realize, a system able to carry out the necessary power system operations without human participation. We outline the capabilities of an automated system for managing electric power. The overarching task of the automated system is islanding: To separate the metro area's power system from the primary power grid and manage its operation during several hurricane-induced contingencies, with the power system operational throughout. The essential technology needed to support this automation is agents. We address the roles agents play in the transition to the islanded state and in power system operation within the island; features of an appropriate agent substrate; the way the agents are organized; and information exchange among agents, the power system, and human operators. ©2006 IEEE.

This work reports on the densification of iron nanoparticles by slow drying followed by pressureless sintering. In contrast, most previous work has used high heating rates to both dry and density the nanoparticle suspension in a single step. Laser heating has been required to achieve high densities by this approach. The slow drying/pressureless sintering approach is shown to be sensitive to reactions between the particles, the stabilizing ligands, the atmosphere, and the substrate. The sintering rate of iron nanoparticles and the final composition of the deposits are significantly impacted by these interactions. However, in both the cases studied, the nanoparticles densify under pressureless sintering. When the iron nanoparticle colloid is dried in a porous steel skeleton, it is shown to increase high-temperature strength and reduce the sintering shrinkage.

Oxygen-fuel fired glass melting furnaces have successfully reduced NO x and particulate emissions and improved the furnace energy efficiency relative to the more conventional air-fuel fired technology. However, full optimisation of the oxygen/fuel approach (particularly with respect to crown refractory corrosion) is unlikely to be achieved until there is improved understanding of the effects of furnace operating conditions on alkali vaporization, batch carryover, and the formation of gaseous air pollutants in operating furnaces. In this investigation, continuous online measurements of alkali concentration (by laser induced breakdown spectroscopy) were coupled with measurements of the flue gas composition in the exhaust of an oxygen/natural gas fired container glass furnace. The burner stoichiometry was purposefully varied while maintaining normal glass production. The data demonstrate that alkali vaporization and SO2 release increase as the oxygen concentration in the exhaust decreases. NOx emissions showed a direct correlation with the flow rate of infiltrated air into the combustion space. The extent of batch carryover was primarily affected by variations in the furnace differential pressure. The furnace temperature did not vary significantly during the measurement campaign, so no clear correlation could be obtained between the available measurements of furnace temperature and alkali vaporization.

Techniques to ensure shock data quality and to recognize bad data are discussed in this paper. For certain shock environments, acceleration response up to ten kHz is desired for structural model validation purposes. The validity and uncertainty associated with the experimental data need to be known in order to use it effectively in model validation. In some cases the frequency content of impulsive or pyrotechnic loading or metal to metal contact of joints in the structure may excite accelerometer resonances at hundreds of kHz. The piezoresistive accelerometers often used to measure such events can provide unreliable data depending on the level and frequency content of the shock. The filtered acceleration time history may not reveal that the data are unreliable. Some data validity considerations include accelerometer mounting systems, sampling rates, band-edge settings, peak acceleration specifications, signal conditioning bandwidth, accelerometer mounted resonance and signal processing checks. One approach for uncertainty quantification of the sensors, signal conditioning and data acquisition system is also explained.

A flexible, robust method for linking grids of locally refined ground-water flow models constructed with different numerical methods is needed to address a variety of hydrologic problems. This work outlines and tests a new ghost-node model-linking method for a refined "child" model that is contained within a larger and coarser "parent" model that is based on the iterative method of Steffen W. Mehl and Mary C. Hill (2002, Advances in Water Res., 25, p. 497-511; 2004, Advances in Water Res., 27, p. 899-912). The method is applicable to steady-state solutions for ground-water flow. Tests are presented for a homogeneous two-dimensional system that has matching grids (parent cells border an integer number of child cells) or nonmatching grids. The coupled grids are simulated by using the finite-difference and finite-element models MODFLOW and FEHM, respectively. The simulations require no alteration of the MODFLOW or FEHM models and are executed using a batch file on Windows operating systems. Results indicate that when the grids are matched spatially so that nodes and child-cell boundaries are aligned, the new coupling technique has error nearly equal to that when coupling two MODFLOW models. When the grids are nonmatching, model accuracy is slightly increased compared to that for matching-grid cases. Overall, results indicate that the ghost-node technique is a viable means to couple distinct models because the overall head and flow errors relative to the analytical solution are less than if only the regional coarse-grid model was used to simulate flow in the child model's domain.

A discussion on an active gas imager that can potentially improve system performance and reliability in Smart Leak Detection and Repair covers conventional single-wavelength imaging; differential imaging; methane detection; modification for detecting fugitive emissions relevant to refineries and chemical plants; and system description. This is an abstract of a paper presented at the AWMA's 99th Annual Conference and Exhibition (New Orleans, LA 6/20-23/2006).

The RADTRAN Loss of Shielding (LOS) Model was benchmarked using MicroShield 6.20®. This analysis considers an intact spent fuel truck cask as well as a set of damaged truck casks. Ratios of dose rates are calculated for casks with a loss of lead shielding to those of intact casks, and are then compared to ratios generated by the LOS model. LOS Model results were considered verified if two main constraints were satisfied. First, the dose rate profiles for both the LOS and MicroShield 6.20® calculations must have the same general shape and behavior. Additionally, the largest factor difference between any two points of the dose rate profiles may not exceed an order of magnitude. Reasonable agreement is shown for large-fraction LOS scenarios; however the differences in results are not satisfactory for cases with small fractions of slump.

A wide variety of MicroElectroMechanical Systems (MEMS) are fabricated using existing novel technologies. State-of-the art integrated circuit (IC) fabrication methods are used for the fabrication of these MEMS. The fabrication of these structures requires many process steps that include deposition, patterning, etching, and CMP. The use of CMP enables the fabrication of complex, multi-level MEMS. Similar to IC fabrication, there are concerns about non-uniformity, erosion and dishing after CMP, but because of the thickness of the materials, CMP processing issues are amplified. Unlike ICs, there is no transistor basic building block so processing must be technology specific and process development is driven by the device/system performance requirements, which are very specific to the application.

Abstract not provided.

We apply ab initio molecular dynamics (AIMD) to study the hydration structures and electronic properties of the formohydroxamate anion in liquid water. We consider the cis- nitrogen-deprotonated, cis- oxygen-deprotonated, and trans- oxygen-deprotonated formohydroxamate tautomers. They form an average of 6.3, 6.9, and 6.0 hydrogen bonds with water molecules, respectively. The predicted pair correlation functions and time dependence of the hydration numbers suggest that water is highly structured around the nominally negatively charged oxime oxygen in O-deprotonated tautomers but significantly less so around the nitrogen atom in the N-deprotonated species. Wannier function analysis suggests that, in the O-deprotonated anions, the negative charge is concentrated on the oxime oxygen, while in the N-deprotonated case, it is partially delocalized between the nitrogen and the adjoining oxime oxygen atom. © 2006 Elsevier B.V. All rights reserved.