Publications

This study models the incremental radiological risk of transporting NARM to the Hanford commercial LLW facility, both for incident-free transportation and for possible transportation accidents, compared with the radiological risk of transporting LLW to that facility. Transportation routes are modeled using HIGHWAY 3.1 and risks are modeled using RADTRAN 4. Both annual population doses and risks, and annual average individual doses and risks are reported. Three routes to the Hanford site were modeled from Albany, OR, from Coeur d'Alene, ID (called the Spokane route), and from Seattle, WA. Conservative estimates are used in the RADTRAN inputs, and RADTRAN itself is conservative.

A parametric analysis is presented which summarizes the amount of total ({gamma}/NbC + {gamma}/Laves) and individual {gamma}/NbC and {gamma}/Laves constituents which form during solidification of {gamma}{sub (Fe,Ni,Cr)} alloys with variations in nominal Nb and C contents. Calculated results are presented for Fe base alloys and Ni base alloys. The results provide a quantitative rationale for understanding the relation between alloy composition and solidification microstructures and should provide useful insight into commercial alloys of similar composition.

There are two sources of contamination in solder alloys. The first source is trace elements from the primary metals used in the as-manufactured product, be that product in ingot, wire, or powder form. Their levels in the primary metal are determined by the refining process. While some of these trace elements are naturally occurring materials, additional contamination can result from the refining and/or forming processes. Sources include: furnace pot liners, debris on the cutting edges of shears, rolling mill rollers, etc. The types and levels of contaminants per solder alloy are set by recognized industrial, federal, military, and international specifications. For example, the 63Sn-37Pb solder purchased to the ASTM B 32 standard can have maximum levels of contamination for the following metals: 0.08(wt.)%Cu, 0.001 %Cd, 0.005%Al, 0.25%Bi, 0.03%As, 0.02%Fe, and 0.005 %Zn. A second cause of contamination in solders, and solder baths in particular, is their actual use in soldering operations. Each time a workpiece is introduced into the bath, some dissolution of the joint base metal(s), protective or solderable coatings, and fixture metal takes place which adds to contamination levels in the solder. The potential impurities include Cu; Ni; Au or other noble metals used as protective finishes and Al; Fe; and Zn to name a few. Even dissolution of the pot wall or liner is a source of impurities, typically Fe.

This paper details the application of a standard physical security system design process to a US Border Port of Entry (PoE) for vehicle entry/exit. The physical security design methodology is described as well as the physical security similarities to facilities currently at a US Border PoE for vehicles. The physical security design process description includes the various elements that make up the methodologies well as the considerations that must be taken into account when dealing with system integration of those elements. The distinctions between preventing unlawful entry/exit of illegal contraband and personnel are described. The potential to enhance the functions of drug/contraband detection in the Pre-Primary Inspection area through the application of emerging technologies are also addressed.

We consider the problem of planning shortest paths for a tethered robot with a finite length tether in a 2D environment with polygonal obstacles. We present an algorithm that runs in time O((k{sub 1} + 1){sup 2}n{sup 4}) and finds the shortest path or correctly determines that none exists that obeys the constraints; here n is the number obstacle vertices, and k{sub 1} is the number loops in the initial configuration of the tether. The robot may cross its tether but nothing can cross obstacles, which cause the tether to bend. The algorithm applies as well for planning a shortest path for the free end of an anchored cable.

The authors have developed a new type of heat spreader based on the integration of heat pipes directly within a thin planar structure suitable for use as a heat spreader or as the base layer in a substrate. The process uses micromachining methods to produce micron scale patterns that act as a wick in these small scale heat pipes. By using silicon or a low expansion metal as the wall material of these spreaders, they achieve a good match to the thermal coefficient of expansion of the die. The match allows the use of a thin high performance die attachment even on large size die. The embedded heat pipes result in high effective thermal conductivity for the new spreader technology.

Sampling during environmental drilling is essential to fully characterize the spatial distribution and migration of subsurface contaminants. However, analysis of the samples is expensive and time-consuming: off-site laboratory analysis can take weeks or months. Real-time information on environmental conditions, drill bit location and temperature during drilling is valuable in many environmental restoration operations. This type of information can be used to provide field screening data and improved efficiency of site characterization activities. The Environmental Measurement-While-Drilling (EMWD) System represents an innovative blending of new and existing technology in order to obtain real-time data during drilling. The system consists of two subsystems. The down-hole subsystem (at the drill bit) consists of sensors, a power supply, a signal conditioning and transmitter board, and a radio-frequency (RF) coaxial cable. The up-hole subsystem consists of a battery pack/coil, pickup coil, receiver, and personal computer. The system is compatible with fluid miser drill pipe, a directional drilling technique that uses minimal drilling fluids and generates little to no secondary waste. In EMWD, downhole sensors are located behind the drill bit and linked by a high-speed data transmission system to a computer at the surface. Sandia-developed Windows{trademark}-based software is used for data display and storage. As drilling is conducted, data is collected on the nature and extent of contamination, enabling on-the-spot decisions regarding drilling and sampling strategies. Initially, the downhole sensor consisted of a simple gamma radiation detector, a Geiger-Mueller tube (GMT). The design includes data assurance techniques to increase safety by reducing the probability of giving a safe indication when an unsafe condition exists. The EMWD system has been improved by the integration of a Gamma Ray Spectrometer (GRS) in place of the GMT. The GRS consists of a sodium iodide-thallium activated crystal coupled to a photomultiplier tube (PMT). The output of the PMT goes to a multichannel analyzer (MCA).The MCA data is transmitted to the surface via a signal conditioning and transmitter board similar to that used with the GMT. The EMWD system is described and the results of the GRS field tests and field demonstration are presented.

The Environmental Measurement-While-Drilling-Gamma Ray Spectrometer (EMWD-GRS) system represents an innovative blend of new and existing technology that provides real-time environmental and drill bit data during drilling operations. The EMWD-GRS technology was demonstrated at Savannah River Site (SRS) F-Area Retention Basin. The EMWD-GRS technology demonstration consisted of continuously monitoring for gamma-radiation-producing contamination while drilling two horizontal boreholes below the backfilled waste retention basin. These boreholes passed near previously sampled locations where concentrations of contaminant levels of cesium had been measured. Contaminant levels continuously recorded by the EMWD-GRS system during drilling were compared to contaminant levels previously determined through quantitative laboratory analysis of soil samples. The demonstration of the EMWD-GRS was a complete success. The results show general agreement between the soil sampling and EMWD-GRS techniques for CS-137. It was recognized that the EMWD-GRS tool would better satisfy our customers' needs if the instrument location could be continuously monitored. During the demonstration at SRS, an electromagnetic beacon with a walkover monitor (Subsite{reg_sign}) was used to measure bit location at depth. To use a beacon locator drilling must be stopped, thus it is normally only used when a new section of pipe was added. The location of contamination could only be estimated based on the position of the EMED-GRS package and the distance between locator beacon readings. A continuous location system that would allow us to know the location of each spectrum as it is obtained is needed.

The iterative advanced mean value algorithm (AMV+), introduced nearly ten years ago, is now widely used as a cost-effective probabilistic structural analysis tool when the use of sampling methods is cost prohibitive (Wu et al., 1990). The need to establish confidence bounds on calculated probabilities arises because of the presence of uncertainties in measured means and variances of input random variables. In this paper an algorithm is proposed that makes use of the AMV+ procedure and analytically derived probability sensitivities to determine confidence bounds on calculated probabilities.

Several ultrasonic inspection methods were developed at the Federal Aviation Administration's Airworthiness Assurance NDI Validation Center (AANC) to easily and rapidly detect hidden stress corrosion cracks in all vertical windshield posts on the US Coast Guard (USCG) HU-25 Guardian aircraft. The inspection procedure locates cracks as small as 2.0 millimeters emanating from internal fastener holes and determines their length. A test procedure was developed and a baseline assessment of the USCG fleet was conducted. Inspection results on twenty-five aircraft revealed a good correlation with results made during subsequent structural disassembly and visual inspection.

As III-V nitride device structures become more complicated and design rules shrink, well-controlled etch processes are necessary. Due to limited wet chemical etch results for the group-III nitrides, a significant amount of effort has been devoted to the development of dry etch processing. Dry etch development was initially focused on mesa structures where high etch rates, anisotropic profiles, smooth sidewalls, and equi-rate etching of dissimilar materials were required. For example, commercially available LEDs and laser facets for GaN-based laser diodes have been patterned using reactive ion etching (RIE). With the recent interest in high power, high temperature electronic devices, etch characteristics may also require smooth surface morphology, low plasma-induced damage, and selective etching of one layer over another. The principal criteria for any plasma etch process is its utility in the fabrication of a device. In this study, we will report plasma etch results for the group-III nitrides and their application to device structures.

The experimentally determined creep responses of a number of domal salts have been reported in, the literature. Some of these creep results were obtained using standard (conventional) creep tests. However, more typically, the creep data have come from multistage creep tests, where the number of specimens available for testing was small. An incremental test uses abrupt changes in stress and temperature to produce several time increments (stages) of different creep conditions. Clearly, the ability to analyze these limited data and to correlate them with each other could be of considerable potential value in establishing the mechanical characteristics of salt domes, both generally and specifically. In any analysis, it is necessary to have a framework of rules to provide consistency. The basis for the framework is the Multimechanism-Deformation (M-D) constitutive model. This model utilizes considerable general knowledge of material creep deformation to supplement specific knowledge of the material response of salt. Because the creep of salt is controlled by just a few micromechanical mechanisms, regardless of the origin of the salt, certain of the material parameters are values that can be considered universal to salt. Actual data analysis utilizes the methodology developed for the Waste Isolation Pilot Plant (WIPP) program, and the response of a bedded pure WIPP salt as the baseline for comparison of the domal salts. Creep data from Weeks Island, Bryan Mound, West Hackberry, Bayou Choctaw, and Big Hill salt domes, which are all sites of Strategic Petroleum Reserve (SPR) storage caverns, were analyzed, as were data from the Avery Island, Moss Bluff, and Jennings salt domes. The analysis permits the parameter value sets for the domal salts to be determined in terms of the M-D model with various degrees of completeness. In turn this permits detailed numerical calculations simulating cavern response. Where the set is incomplete because of the sparse database, reasonable assumptions permit the set to be completed. From the analysis, two distinct response groups were evident, with the salts of one group measurably more creep resistant than the other group. Interestingly, these groups correspond well with the indirectly determined creep closure of the SPR storage caverns, a correlation that probably should be expected. Certainly, the results suggest a simple laboratory determination of the creep characteristics of a salt material from a dome site can indicate the relative behavior of any potential cavern placed within that dome.

Three-dimensional imaging techniques, numerical methods for simulating flow and transport, and emergent computational architectures are combined to enable fundamental studies of fluid flow at the pore scale. High resolution reconstructions of porous media obtained using laser scanning confocal microscopy reduce sampling artifacts to sub-micron features, and simultaneously capture multiple grain length scales. However, the volumetric image data sets are extremely large, and there are significant computational challenges in utilizing this information effectively. The principal problem lies in the complexity of the geometry and the retention of this structure in numerical analyses. Lattice Boltzmann (LB) methods provide a direct means to simulate transport processes in complex geometric domains due to the unique ability to treat accurately and efficiently the multitude of discrete boundary conditions. LB methods are numerically explicit as formulated, and this characteristic is exploited through a mapping of the numerical domain to distributed computing architectures. These techniques are applied to perform single phase flow simulations in 3D data sets obtained from cores of Berea sandstone using confocal microscopy. Simulations are performed using both a purpose-built distributed processor computer and a massively parallel processer (MPP) platform.

The authors designed and conducted experiments in a heterogeneous sand pack where gravity-destabilized nonwetting phase invasion (CO{sub 2} and TCE) could be recorded using high resolution light transmission methods. The heterogeneity structure was designed to be reminiscent of fluvial channel lag cut-and-fill architecture and contain a series of capillary barriers. As invasion progressed, nonwetting phase structure developed a series of fingers and pools; behind the growing front they found nonwetting phase saturation to pulsate in certain regions when viscous forces were low. Through a scale analysis, they derive a series of length scales that describe finger diameter, pool height and width, and regions where pulsation occurs within a heterogeneous porous medium. In all cases, they find that the intrinsic pore scale nature of the invasion process and resulting structure must be incorporated into the analysis to explain experimental results. The authors propose a simple macro-scale structural growth model that assembles length scales for sub-structures to delineate nonwetting phase migration from a source into a heterogeneous domain. For such a model applied at the field scale for DNAPL migration, they expect capillary and gravity forces within the complex subsurface lithology to play the primary roles with viscous forces forming a perturbation on the inviscid phase structure.

Laser-induced damage mechanisms that can occur during high-intensity fiber transmission have been under study for a number of years. Our particular interest in laser initiation of explosives has led us to examine damage processes associated with the transmission of Q-switched, Nd:YAG pulses at 1.06 {micro}m through step-index, multimode, fused silica fiber. Laser breakdown at the fiber entrance face is often the first process to limit fiber transmission but catastrophic damage can also occur at either fiber end face, within the initial entry segment of the fiber, and at other internal sites along the fiber path. Past studies have examined how these various damage mechanisms depend upon fiber end-face preparation, fiber fixturing and routing, laser characteristics, and laser-to-fiber injection optics. In some applications of interest, however, a fiber transmission system may spend years in storage before it is used. Consequently, an important additional issue for these applications is whether or not there are aging processes that can result in lower damage thresholds over time. Fiber end-face contamination would certainly lower breakdown and damage thresholds at these surfaces, but careful design of hermetic seals in connectors and other end-face fixtures can minimize this possibility. A more subtle possibility would be a process for the slow growth of internal defects that could lead to lower thresholds for internal damage. In the current study, two approaches to stimulating the growth of internal defects were used in an attempt to produce observable changes in internal damage thresholds. In the first approach test fibers were subjected to a very high tensile stress for a time sufficient for some fraction to fail from static fatigue. In the second approach, test fibers were subjected to a combination of high tensile stress and large, cyclic temperature variations. Both of these approaches were rather arbitrary due to the lack of an established growth mechanism for internal defects. Damage characteristics obtained from fibers subjected to each of these aging environments were compared to results from fresh fibers tested under identical conditions. A surprising result was that internal damage was not observed in any of the tested fibers. Only breakdown at the fiber entrance face and catastrophic damage at both end faces were observed. Fiber end faces were not sealed during the accelerated aging environments, and thresholds at these faces were significantly lower in the aged fibers. However, most fibers transmitted relatively high pulse energies before damaging, and a large fraction never damaged before we reached the limits of our test laser. The absence of any observable affect on internal damage thresholds is encouraging, but the current results do not rule out the possibility that some other approach to accelerated aging could reveal a growth mechanism for internal defects.

The time-dependence of the wet oxidation of high-Al-content AlGaAs can be either linear, indicating reaction-rate limitation, or parabolic, indicating diffusion-limited rates. The transition from linear to parabolic time dependence can be explained by the increased rate of the formation of intermediate As{sub 2}O{sub 3} vs. its reduction to elemental As. A steadily increasing thickness of the As{sub 2}O{sub 3}-containing region at the oxidation front will shift the process from the linear to the parabolic regime. This shift from reaction-rate-limited (linear) to diffusion-limited (parabolic) time dependence is favored by increasing temperature or increasing Al mole fraction.

The work functions of polycrystalline metals are often used to systematize Schottky barrier height data for rectifying contacts to semiconductors. Rectifying contacts to silicon devices are predominantly formed using conductive metal silicides with work functions which are not as well characterized as metal work functions. The present work has two objectives. First, it classifies the transition metals using correlations between the metal work function and the atomic chemical potential. Second, the available data for metal silicides is collected and interpreted using an average charge transfer (ACT) model. The ACT model accounts for the electronic hardness of the component elements in addition to their chemical potentials. New trends in the behavior of silicide work functions are identified.

We discuss the role of localized high electric fields in the modification of Au surfaces with a W probe using the Interfacial Force Microscope. Upon bringing a probe close to a Au surface, we measure both the interfacial force and the field emission current as a function of separation with a constant potential of 100 V between tip and sample. The current initially increases exponentially as the separation decreases. However, at a distance of less than {approximately} 500{angstrom} the current rises sharply as the surface begins to distort and rapidly close the gap. Retraction of the tip before contact is made reveals the formation of a mound on the surface. We propose a simple model, in which the localized high electric field under the tip assists the production of mobile Au adatoms by detachment from surface steps, and a radial field gradient causes a net flux of atoms toward the tip by surface diffusion. These processes give rise to an unstable surface deformation which, if left unchecked, results in a destructive mechanical contact. We discuss our findings with respect to earlier work using voltage pulses in the STM as a means of nanofabrication.

The self-diffusion of Ca and the tracer diffusion of Mg in calcite have been experimentally measured using isotopic tracers of {sup 25}Mg and {sup 44}Ca. Natural single crystals of calcite were coated with a thermally-sputtered oxide thin film and then annealed in a CO{sub 2} gas at one atmosphere total pressure and temperatures from 550 to 800 C. Diffusion coefficient values were derived from the depth profiles obtained by ion microprobe analysis. The resultant activation energies for Mg tracer diffusion and Ca self-diffusion are respectively: E{sub a}(Mg) = 284 {+-} 74 kJ/mol and E{sub a}(Ca) = 271 {+-} 80 kJ/mol. For the temperature ranges in these experiments, the diffusion of Mg is faster than Ca. The results are generally consistent in magnitude with divalent cation diffusion rates obtained in previous studies and provide a means of interpreting the thermal histories of carbonate minerals, the mechanism of dolomitization, and other diffusion-controlled processes. The results indicate that cation diffusion in calcite is relatively slow and cations are the rate-limiting diffusing species for the deformation of calcite and carbonate rocks. Application of the calcite-dolomite geothermometer to metamorphic assemblages will be constrained by cation diffusion and cooling rates. The direct measurement of Mg tracer diffusion in calcite indicates that dolomitization is unlikely to be accomplished by Mg diffusion in the solid state but by a recrystallization process.

Calorimetric studies by Chai and Navrotsky (1996) on dolomite-ankerite energetic have been extended by including two additional types of samples: a very disordered stoichiometric MgCa(CO{sub 3}){sub 2} prepared from low temperature aqueous solution and three largely ordered natural samples of intermediate iron content. Combining these data with previous work, three distinct trends of energetic can be seen: those for samples with nearly complete order, nearly complete disorder, and intermediate order. From these trends, the enthalpy of complete disordering is estimated to be 33 {+-} 6 kJ/mol for MgCa(CO{sub 3}){sub 2} and 18 {+-} 5 kJ/mol for FeCa(CO{sub 3}){sub 2}.

Evolutionary pattern search algorithms (EPSAs) are a class of evolutionary algorithms (EAs) that have convergence guarantees on a broad class of nonconvex continuous problems. In previous work we have analyzed the empirical performance of EPSAs. This paper revisits that analysis and extends it to a more general model of mutation. We experimentally evaluate how the choice of the set of mutation offsets affects optimization performance for EPSAs. Additionally, we compare EPSAs to self-adaptive EAs with respect to robustness and rate of optimization. All experiments employ a suite of test functions representing a range of modality and number of multiple minima.

Evolutionary programs (EPs) and evolutionary pattern search algorithms (EPSAS) are two general classes of evolutionary methods for optimizing on continuous domains. The relative performance of these methods has been evaluated on standard global optimization test functions, and these results suggest that EPSAs more robustly converge to near-optimal solutions than EPs. In this paper we evaluate the relative performance of EPSAs and EPs on a real-world application: flexible ligand binding in the Autodock docking software. We compare the performance of these methods on a suite of docking test problems. Our results confirm that EPSAs and EPs have comparable performance, and they suggest that EPSAs may be more robust on larger, more complex problems.

As part of a pollution prevention program, a study was conducted at Sandia National Laboratories and at the Amarillo, ''Pantex Plant'' to identify a suitable replacement solvent(s) for cleaning hardware during routine maintenance operations. Current cleaning is performed using solvents (e.g. acetone, toluene, MEK, alcohols) that are classified as Resource Conservation and Recovery Act (RCW) materials. The Environmental Protection Agency (EPA) has assigned four characteristics as the criteria for determining whether a material is identified as hazardous under RCRA: Ignitability, Corrosivity, Reactivity and Toxicity. Within the DOE and DoD sector, these solvents are used with hand wipes to clean surfaces prior to O-ring replacement, to remove decals for new labeling, to clean painted surfaces prior to reconditioning, and for other general maintenance purposes. In some cases, low level radioactive contamination during cleaning necessitates that the RCIL4 solvent-containing wipes be classified as mixed waste. To avoid using RCRA materials, cleaning candidates were sought that had a flashpoint greater than 140 F, a pH between 2.5 and 12.5, and did not fail the reactivity and toxicity criteria. Three brominated cleaners, two hydrofluoroether azeotropes and two aliphatic hydrocarbon cleaner formulations were studied as potential replacements. Cleaning efficacy, materials compatibility, corrosion and accelerated aging studies were conducted and used to screen potential candidates. Hypersolve NPB (an n-propyl bromide based formulation) consistently ranked high in removing typical contaminants for weapons applications.

Deep and shallow electron traps form in irradiated thermal SiO{sub 2} as a natural response to hole transport and trapping. The density and stability of these defects are discussed, as are their implications for total-dose modeling.

Techniques for modeling oil well sand production have been developed using the formulations for superquadric discrete elements and Darcy fluid flow. Discrete element models are generated using the new technique of particle cloning. Discrete element sources and sinks allow simulation of sand production from the initial state through the transition to an equilibrium state where particles are created and removed at the same rate.