Publications

A major objective of the Nuclear Regulatory Commission`s (NRC) Individual Plant Examination (IPE) Insights Program is to identify the important determinants of core damage frequency (CDF) for the different reactor and containment types and plant designs as indicated in the IPEs. The human reliability analysis (HRA) is a critical component of the probabilistic risk assessments (PRAS) which were done for the IPES. The determination and selection of human actions for incorporation into the event and fault tree models and the quantification of their failure probabilities can have an important impact on the resulting estimates of CDF and risk. Therefore, two important goals of the NRCs IPE Insights Program are (1) to determine the extent to which human actions and their corresponding failure probabilities influenced the results of the IPEs and (2) to identify which factors played significant roles in determining the differences and similarities in the results of the HRA analyses across the different plants. To obtain the relevant information, the NRC`s IPE database, which contains information on plant design, CDF, and containment performance obtained from the IPES, was used in conjunction with a systematic examination of the HRA analyses and results from the IPES. Regarding the extent to which the results of the HRA analyses were significant contributors to the plants` CDFs, examinations of several different measures indicated that while individual human actions could have important influences on CDF for particular initiators, the HRA results did not appear to be the most significant driver of plant risk (CDF). Another finding was that while there were relatively wide variations in the calculated human error probabilities (HEPs) for similar events across plants, there was no evidence for any systematic variation as a function of the HRA methods used in the analyses.

SEMATECH and the Department of Energy have established a Contamination Free Manufacturing Research Center (CFMRC) located at Sandia National Laboratories. One of the programs underway at the CFMRC is directed towards defect reduction in semiconductor process reactors by the application of computational modeling. The goal is to use fluid, thermal, plasma, and particle transport models to identify process conditions and tool designs that reduce the deposition rate of particles on wafers. The program is directed toward defect reduction in specific manufacturing tools, although some model development is undertaken when needed. The need to produce quantifiable improvements in tool defect performance requires the close cooperation among Sandia, universities, SEMATECH, SEMATECH member companies, and equipment manufacturers. Currently, both plasma (e.g., etch, PECVD) and nonplasma tools (e.g., LPCVD, rinse tanks) are being worked on under this program. In this paper the authors summarize their recent efforts to reduce particle deposition on wafers during plasma-based semiconductor manufacturing.

Cold filling refers to flowing a fluid through piping or tubes that are at temperatures below the fluid’s freezing point. Since the piping and areas of the receiver in a molten-nitrate salt central-receiver solar power plant must be electrically heated to maintain their temperatures above the nitrate salt freezing point (430°F, 221 °C), considerable energy could be used to maintain such temperatures during nightly shutdown and bad weather. Experiments and analyses have been conducted to investigate cold filling receiver panels and piping as a way of reducing parasitic electrical power consumption and increasing the availability of the plant. The two major concerns with cold filling are (1) how far can the molten salt penetrate cold piping before freezing closed, and (2) what thermal stresses develop during the associated thermal shock. Cold fill experiments were conducted by flowing molten salt at 550°F (288°C) through cold panels, manifolds, and piping to determine the feasibility of cold filling the receiver and piping. The transient thermal responses were measured and heat transfer coefficients were calculated from the data. Nondimensional analysis is presented which quantifies the thermal stresses in a pipe or tube undergoing thermal shock. In addition, penetration distances were calculated to determine the distance salt could flow in cold pipes prior to freezing closed. © 1995 by ASME.

This is the final report on a series of experiments concerned with transient radiation-induced absorption in materials for a Cr,Nd:GSGG laser. Both the Sandia National Laboratories SPR III pulsed reactor and the Hermes III pulsed X-ray machine are used as radiation sources. The time dependence and the magnitude of the induced absorption in filter glasses and in doped and undoped LiNbO{sub 3} Q-switch materials have been measured. Gain has been observed in Cr,Nd:GSGG, the laser medium, when it is irradiated by X-rays.

Understanding the evolution of porous silicon (PS) layers at the early stages of growth is important for determining the mechanism of PS film growth and controlling the film properties. We have used X-ray reflectivity (XRR) to determine the evolution of layer thickness and interfacial roughness during the growth of thin PS layers (< 200 nm) prepared by electrochemical anodization. The porous layer grows at a constant rate for films as thin as 15 nm indicating a very short incubation period during which the surface may be electropolished before the PS structure begins to form. Interface roughness measurements indicate that the top surface of the film remains relatively smooth during growth while the roughness of the PS/silicon interface increases only slightly with film thickness. The XRR results are compared with results obtained from the same films by cross-sectional transmission electron microscopy (XTEM), atomic force microscopy (AFM) and gravimetry.

We have used {sup 13}C magic-angle spinning (MAS) nuclear magnetic resonance (NMR) to characterize the structure and rotational dynamics of C{sub 60} containing oxygen molecules located in the interstitial sites of the fcc lattice. Under normal conditions, a narrow peak at 143.7 ppm is observed for C{sub 60}. When exposed to oxygen at moderate pressures, several additional resonances appear in the {sup 13}C MAS NMR spectrum. These secondary resonances are shifted downfield from the main peak at 143.7 ppm and are due to the Fermi-contact interaction of the paramagnetic oxygen molecules with the {sup 13}C nuclear spins. The presence of oxygen depresses the orientational ordering transition by ca. 20 K as observed by DSC. The spin-lattice relaxation time (T{sub 1}) of each secondary peak shows a minimum near the ordering transition, indicating that this transition is not dependent on the number of oxygen molecules surrounding an individual C{sub 60} molecule. The T{sub 1}, due to paramagnetic relaxation, normalized by the number of surrounding oxygen molecules, is constant. This observation demonstrates that within a given sample, the dynamics of C{sub 60} molecules are independent of the number of surrounding oxygen molecules.

This paper describes a process of combining two state-of-the-art CFD tools, SPRINT and INCA, in a manner which extends the utility of both codes beyond what is possible from either code alone. The speed and efficiency of the PNS code, SPRINT, has been combined with the capability of a Navier-Stokes code to model fully elliptic, viscous separated regions on high performance, high speed flight systems. The coupled SPRINT/rNCA capability is applicable for design and evaluation of high speed flight vehicles in the supersonic to hypersonic speed regimes. This paper describes the codes involved, the interface process and a few selected test cases which illustrate the SPRINT/INCA coupling process. Results have shown that the combination of SPRINT and INCA produces correct results and can lead to improved computational analyses for complex, three-dimensional problems.

Sampling during environmental drilling is essential to fully characterize the spatial distribution and migration of near surface contaminants. However, the analysis of these samples is not only expensive, but can take weeks or months when sent to an off-site laboratory. In contrast, measurement-while-drilling (MWD) screening capability could save money and valuable time by quickly distinguishing between contaminated and uncontaminated areas. Real-time measurements provided by a MVM system would enable on-the-spot decisions to be made regarding sampling strategies, enhance worker safety, and provide the added flexibility of being able to ``steer`` the drill bit in or out hazardous zones. During measurement-while-drilling, down-hole sensors are located behind the drill bit and linked by a rapid data transmission system to a computer at the surface. As drilling proceeds, data are collected on the nature and extent of the subsurface contamination in real-time. The down-hole sensor is a Geiger-Mueller tube (GMT) gamma radiation detector. In addition to the GMT signal, the MWD system monitors these required down-hole voltages and two temperatures associated with the detector assembly. The Gamma Ray Detection System (GRDS) and electronics package are discussed in as well as the results of the field test. Finally, our conclusions and discussion of future work are presented.

Uncertainty and sensitivity analysis techniques based on Latin hypercube sampling, partial correlation analysis and stepwise regression analysis are used in an investigation with the MACCS model of the food pathways associated with a severe accident at a nuclear power station. The primary purpose of this study is to provide guidance on the variables to be considered in future review work to reduce the uncertainty in the important variables used in the calculation of reactor accident consequences. The effects of 87 imprecisely-known input variables on the following reactor accident consequences are studied: crop growing-season dose, crop long-term dose, milk growing-season dose, total food pathways dose, total ingestion pathways dose, total long-term pathways dose, area dependent cost, crop disposal cost, milk disposal cost, condemnation area, crop disposal area and milk disposal area. When the predicted variables are considered collectively, the following input variables were found to be the dominant contributors to uncertainty: fraction of cesium deposition on grain fields that is retained on plant surfaces and transferred directly to grain, maximum allowable ground concentrations of Cs-137 and Sr-90 for production of crops, ground concentrations of Cs-134, Cs-137 and I-131 at which the disposal of milk will be initiated due to accidents that occur during the growing season, ground concentrations of Cs-134, I-131 and Sr-90 at which the disposal of crops will be initiated due to accidents that occur during the growing season, rate of depletion of Cs-137 and Sr-90 from the root zone, transfer of Sr-90 from soil to legumes, transfer of Cs-137 from soil to pasture, transfer of cesium from animal feed to meat, and the transfer of cesium, iodine and strontium from animal feed to milk. © 1995.

Uncertainty and sensitivity analysis techniques based on Latin hypercube sampling, partial correlation analysis and stepwise regression analysis are used in an investigation with the MACCS model of the chronic exposure pathways associated with a severe accident at a nuclear power station. The primary purpose of this study is to provide guidance on the variables to be considered in future review work to reduce the uncertainty in the important variables used in the calculation of reactor accident consequences. The effects of 75 imprecisely known input variables on the following reactor accident consequences are studied: crop growing-season dose, crop long-term dose, water ingestion dose, milk growing-season dose, long-term groundshine dose, long-term inhalation dose, total food pathways dose, total ingestion pathways dose, total long-term pathways dose, total latent cancer fatalities, area-dependent cost, crop disposal cost, milk disposal cost, population-dependent cost, total economic cost, condemnation area, condemnation population, crop disposal area and milk disposal area. When the predicted variables are considered collectively, the following input variables were found to be the dominant contributors to uncertainty: dry deposition velocity, transfer of cesium from animal feed to milk, transfer of cesium from animal feed to meet, ground concentration of Cs-134 at which the disposal of milk products will be initiated, transfer of Sr-90 from soil to legumes, maximum allowable ground concentration of Sr-90 for production of crops, fraction of cesium entering surface water that is consumed in drinking water, groundshine shielding factor, scale factor defining resuspension, dose reduction associated with decontamination, and ground concentration of I-131 at which disposal of crops will be initiated due to accidents that occur during the growing season. Reducing the uncertainty in the preceding variables was found to substantially reduce the uncertainty in the predicted variables under consideration. For total number of latent cancer fatalities, the dominant variable was dry deposition velocity, with small effects indicated for a large number of additional variables. © 1995.

Uncertainty and sensitivity analysis techniques based on Latin hypercube sampling, partial correlation analysis and stepwise regression analysis are used in an investigation with the MACCS model of the early health effects associated with a severe accident at a nuclear power station. The primary purpose of this study is to provide guidance on the variables to be considered in future review work to reduce the uncertainty in the important variables used in the calculation of reactor accident consequences. The effects of 34 imprecisely-known input variables on the following reactor accident consequences are studied: number of early fatalities, number of cases prodromal vomiting, population dose within 10 mi of the reactor, population dose within 1000 mi of the reactor, individual early fatality probability within 1 mi of the reactor, and maximum early fatality distance. When the predicted variables are considered collectively, the following input variables were found to be the dominant contributors to uncertainty: scaling factor for horizontal dispersion, dry deposition velocity, inhalation protection factor for nonevacuees, groundshine shielding factor for nonevacuees, early fatality hazard function alpha value for bone marrow exposure, and scaling factor for vertical dispersion. © 1995.

The effect of the pump, signal, and idler wave phases on three-wave nonlinear parametric mixing is investigated in a series of single-pass- gain experiments. Measurements are made with two angle-tuned KTP crystals in a 532 nm pumped, walkoff-compensated, optical parametric amplifier that is seeded by an 800 nm cw diode laser. In one of the measurements the second crystal is orientated to have its effective nonlinearity deff. of opposite sign to that of the first crystal, so that all mixing that occurred in the first crystal is canceled by the second when the phase mismatch Δkcrystal 1 = Δkcrystal 2 = 0. Efficient two-crystal amplification is subsequently restored by selecting the correct phase relationship for the three waves entering the crystal by inserting a dispersive plate between the crystals. The experimental results are explained in a straightforward manner with diagrams involving the three input wave polarizations. These results demonstrate that walkoff-compensated geometries require phase correction to achieve efficient mixing in the second crystal whenever the nonlinear interaction involves two extraordinary waves (e-waves). One practical application of this work may be lower oscillation thresholds and enhanced performance in walkoff-compensated optical parametric oscillators which use two e-waves.

This report summarizes two approaches to time-optimal control of a nonlinear magnetically levitated platen. The system of interest is a candidate technology for next-generation photolithography machines used in the manufacture of integrated circuits. The dynamics and the variable peak control force of the electro-magnetic actuators preclude the direct application of classical time-optimal control methodologies for determining optimal rest-to-rest maneuver strategies. Therefore, this study explores alternate approaches using a previously developed computer simulation. In the first approach, conservative estimates of the available control forces are used to generate suboptimal switching curves. In the second approach, exact solutions are determined iteratively and used as a training set for an artificial neural network. The trained network provides optimal actuator switching times that incorporate the full nonlinearities of the magnetic levitation actuators. Sample problems illustrate the effectiveness of these techniques as compared to traditional proportional-derivative control.

The Critical Dynamics in Microgravity Experiment, DYNAMX, is under development for space flight at Sandia National Laboratories and the University of New Mexico with Dr. Rob Duncan as the Principal Investigator. This experiment will investigate the effects on the superfluid transition in 4He of currents generated by heat flow, measuring the thermal conductivity in the fluid as a function of applfed heat. DYNAMX will also take advantage of the weightless environment to measure the conductivity properties in the region of the interface between the two phases. Thus, DYNAMX represents an experiment that will explore a system driven far from equilibrium. The experiment development is sponsored by the Microgravity Science and Applications Division of NASA, with the Jet Propulsion Laboratory as the managing center. This paper will describe the science objectives, the current design of the experiment apparatus, the steps being taken to prepare this experiment for flight, and the results of ground-based feasibility demonstrations now underway.

We present magnetoluminescence data which provides a quantitative measure of the energy- band dispersion curves of novel compound semiconductor optoelectronic materials. Data for a n-type strained-layer InGaAs/GaAs (quantum-well width approximately 8 nm) and a n-type 4.5 nm-wide GaAs/AlGaAs lattice-matched single-quantum well are presented. We find that the conduction-bands are almost parabolic, with a mass of about 0.068m0 for the GaAs/AlGaAs structure. The valence-bands are nonparabolic with wave vector dependent in- plane valence-band masses varying from about 0.1m0 at zone center to about 0.3m0 for 20 meV energies.

Sol-gel processing methods are frequently used for the fabrication of lead zirconate titanate (PZT) thin films for many electronic applications. Our standard approach for film fabrication utilizes lead acetate and acetic acid modified metal alkoxides of zirconium and titanium in the preparation of our precursor solutions. This report highlights some of our recent results on the effects of the addition of a second chelating ligand, acetylacetone, to this process. We discuss the changes in film drying behavior, densification and ceramic microstructure which accompany acetylacetone additions to the precursor solution and relate the observed variations in processing behavior to differences in chemical precursor structure induced by the acetylacetone ligand. Improvements in thin film microstructure, ferroelectric and optical properties are observed when acetylacetone is added to the precursor solution.

Active control of structures has been under intensive development for the last ten years. Reference 2 reviews much of the identification and control technology for structural control developed during this time. The technology was initially focused on space structure and weapon applications; however, recently the technology is also being directed toward applications in manufacturing and transportation. Much of this technology focused on multiple-input/multiple-output (MIMO) identification and control methodology because many of the applications require a coordinated control involving multiple disturbances and control objectives where multiple actuators and sensors are necessary for high performance. There have been many optimal robust control methods developed for the design of MIMO robust control laws; however, there appears to be a significant gap between the theoretical development and experimental evaluation of control and identification methods to address structural control applications. Many methods have been developed for MIMO identification and control of structures, such as the Eigensystem Realization Algorithm (ERA), Q-Markov Covariance Equivalent Realization (Q-Markov COVER) for identification; and, Linear Quadratic Gaussian (LQG), Frequency Weighted LQG and H-/ii-synthesis methods for control. Upon implementation, many of the identification and control methods have shown limitations such as the excitation of unmodelled dynamics and sensitivity to system parameter variations. As a result, research on methods which address these problems have been conducted.

This paper describes numerical simulations that were performed to study laminar flow through a square duct with a 90° bend. The purpose of this work was two fold. First, an improved understanding was desired of the flow physics involved in the generation of secondary vortical flows in three-dimensions. Second, adaptive gridding techniques for structured grids in three-dimensions were investigated for die purpose of determining their utility in low Reynolds number, incompressible flows. It was also of interest to provide additional validation of the commercial computer code CFD-ACE. Velocity predictions for both non-adaptive and adaptive grids are compared with experimental data. Flow visualization was used to examine the characteristics of the flow though the curved duct in order to better understand the viscous flow physics of this problem. Generally fair agreement with the experimental data was found, but questions were raised concerning the accuracy of the experimental data. The adaptive grids did not significantly improve the accuracy of the results beyond the non-adaptive grid solution with a similar number of points.

Since both FeS2 and MoS2 are optically opaque and have bandgaps in the near IR, photochemical applications of these materials are non-existent. It is however demonstrated in the study that nanosize FeS2 and MoS2 have bandgaps that can be adjusted to the visible and even UV region of the spectrum by control of the cluster size. Discussed are measurements of size, structure, optical absorbance and photoemission for these new materials. Also, briefly outlines are synthetic procedures and novel methods for chemical processing of these nanoclusters.

We report a real-time, two-dimensional light scattering study of the evolution of structure of a two component nonionic micelle system undergoing phase separation. The micelles act like molecular slug-a-beds whose domain growth is pathetically lathargic (i.e. slower than the cube root of time prediction for simple binary fluids). In fact, the growth kinetics can be empirically described as a stretched exponential approach to a pinned domain size. Although the kinetics are not yet understood, this anomalous behavior may be due to the ability of the spherical micelles to reorganize into more complex structures.

A jib crane consists of a pendulum-like end line attached to a rotatable jib. Within this general category of cranes there exist devices with multiple degrees of freedom including variable load-line length and variable jib length. These cranes are commonly used for construction and transportation applications. Point-to-point payload maneuvers using jib cranes are performed so as not to excite the spherical pendulum modes of their cable and payload assemblies. Typically, these pendulum modes, although time-varying, exhibit low frequencies. The resulting maneuvers are therefore performed slowly, contributing to high construction and transportation costs. The crane considered here consists of a spherical pendulum attached to a rigid jib. The other end of the jib is attached to a direct drive motor of generating rotational motion. A general approach is presented for determining the open-loop trajectories for the jib rotation for accomplishing fixed-time, point-to-point, residual oscillation free, symmetric maneuvers. These residual oscillation free trajectories purposely excite the pendulum modes in such a way that at the end of the maneuver the oscillatory degrees of freedom are quiescent. Simulation results are presented with experimental verification.

First-principles density-functional calculations utilizing ab initio pseudopotentials and plane- wave expansions are used to determine lattice parameters, bulk moduli, and band structures for AlN, GaN, and InN. It is found that large numbers of plane waves are necessary to resolve the nitrogen 2p wave functions and that explicit treatment of the gallium 3d and indium 4d electrons is important for an accurate description of GaN and InN. Several properties of ternary zinc-blende alloys are determined including their bond-length and bond-angle relaxation and their energy-gap bowing parameters. The similarity of the calculated zinc- blende and wurtzite direct gaps also allows estimates to be made of the energy gap versus composition for wurtzite alloys.

A series of experiments was conducted to determine hydrogen combustion behavior under conditions of rapidly condensing steam caused by water sprays. Experiments were conducted in the Surtsey facility under conditions that were nearly prototypical of those that would be expected in a severe accident in the CE System 80+ containment. Mixtures were initially nonflammable owing to dilution by steam. The mixtures were ignited by thermal glow plugs when they became flammable after sufficient steam was removed by condensation caused by water sprays. No detonations or accelerated flame propagation was observed in the Surtsey facility. The combustion mode observed for prototypical mixtures was characterized by multiple deflagrations with relatively small pressure rises. The thermal glow plugs were effective in burning hydrogen safely by igniting the gases as the mixtures became marginally flammable.

This report identifies the key features noted as requirements in the diagnostic decision-making process of Single Photon Emission Computed Tomography (SPECT) cardiac imaging. The report discusses the critical issues that create the basic system framework for design of an automatic target recognizer (ATR) algorithm prototype to support diagnosis of coronary artery disease. Candidate feature discovery algorithms that may form the basis of future work include Adaptive Resonance Theory and Bayesian Decision Network. A framework for the practitioner-Human-System-Interface would include baseline patient history and demographic data; reference cardiac imagery history; and current overlay imagery to provide complementary information (i.e., coronary angiography, echocardiography, and SPECT images). The goal is to design a prototype that would represent a fused present and historical {open_quotes}whole{close_quotes} functional, structural, and physiologic cardiac patient model. This framework decision-assisting platform would be available to practitioner and student alike, with no {open_quotes}real-world{close_quotes} consequences.

This report describes the process used by Sandia National Laboratories, New Mexico to characterize weapon hardware for disposition. The report describes the following basic steps: (1) the drawing search process and primary hazard identification; (2) the development of Disassembly Procedures (DPs), including demilitarization and sanitization requirements; (3) the generation of a ``disposal tree``; (4) generating RCRA waste disposal information; and (5) documenting the information. Additional data gathered during the characterization process supporting hardware grouping and recycle efforts is also discussed.