Publications

Abstract not provided.

Abstract not provided.

An experiment was conducted comparing the effectiveness of individual versus group electronic brainstorming in addressing real-world "wickedly difficult" challenges. Previous laboratory research has engaged small groups of students in answering questions irrelevant to an industrial setting. The current experiment extended this research to larger, real-world employee groups engaged in addressing organizationrelevant challenges. Within the present experiment, the data demonstrated that individuals performed at least as well as groups in terms of number of ideas produced and significantly (p<.02) outperformed groups in terms of the quality of those ideas (as measured along the dimensions of originality, feasibility, and effectiveness).

Tuberculosis (TB), caused by the bacterium Mycobacterium tuberculosis (Mtb), is a growing international health crisis. Mtb is able to persist in host tissues in a nonreplicating persistent (NRP) or latent state. This presents a challenge in the treatment of TB. Latent TB can re-activate in 10% of individuals with normal immune systems, higher for those with compromised immune systems. A quantitative understanding of latency-associated virulence mechanisms may help researchers develop more effective methods to battle the spread and reduce TB associated fatalities. Leveraging BioXyce's ability to simulate whole-cell and multi-cellular systems we are developing a circuit-based framework to investigate the impact of pathogenicity-associated pathways on the latency/reactivation phase of tuberculosis infection. We discuss efforts to simulate metabolic pathways that potentially impact the ability of Mtb to persist within host immune cells. We demonstrate how simulation studies can provide insight regarding the efficacy of potential anti-TB agents on biological networks critical to Mtb pathogenicity using a systems chemical biology approach. © 2008 IEEE.

Abstract not provided.

We develop an approach for estimating model parameters which result in the "best distribution fit" between experimental and simulation data. Best distribution fit means matching moments of experimental data to those of a simulation (and possibly matching a full probability distribution). This approach extends typical nonlinear least squares methods which identify parameters maximizing agreement between experimental points and computational simulation results. Several analytic formulations for the distribution matching problem are provided, along with results for solving test problems and comparisons of this parameter estimation technique with a deterministic least squares approach. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

We describe a novel approach to fabricate high-aspect-ratio membranes in microchannels by direct laser scanning, and demonstrate >10-fold improvement in sample preconcentration speed by achieving lower fM detection of proteins within 5 minutes. The integrated device can be used for continuous sample preparation, injection, preconcentration, and biochemical binding/reaction applications. © 2008 CBMS.

The rich history of scalable computing research owes much to a rapid rise in computing platform scale in terms of size and speed. As platforms evolve, so must algorithms and the software expressions of those algorithms. Unbridled growth in scale inevitably leads to complexity. This special issue grapples with two facets of this complexity: scalable execution and scalable development. The former results from efficient programming of novel hardware with increasing numbers of processing units (e.g., cores, processors, threads or processes). The latter results from efficient development of robust, flexible software with increasing numbers of programming units (e.g., procedures, classes, components or developers). The progression in the above two parenthetical lists goes from the lowest levels of abstraction (hardware) to the highest (people). This issue's theme encompasses this entire spectrum. The lead author of each article resides in the Scalable Computing Research and Development Department at Sandia National Laboratories in Livermore, CA. Their co-authors hail from other parts of Sandia, other national laboratories and academia. Their research sponsors include several programs within the Department of Energy's Office of Advanced Scientific Computing Research and its National Nuclear Security Administration, along with Sandia's Laboratory Directed Research and Development program and the Office of Naval Research. The breadth of interests of these authors and their customers reflects in the breadth of applications this issue covers. This article demonstrates how to obtain scalable execution on the increasingly dominant high-performance computing platform: a Linux cluster with multicore chips. The authors describe how deep memory hierarchies necessitate reducing communication overhead by using threads to exploit shared register and cache memory. On a matrix-matrix multiplication problem, they achieve up to 96% parallel efficiency with a three-part strategy: intra-node multithreading, non-blocking inter-node message passing, and a dedicated communications thread to facilitate concurrent communications and computations. On a quantum chemistry problem, they spawn multiple computation threads and communication threads on each node and use one-sided communications between nodes to minimize wait times. They reduce software complexity by evolving a multi-threaded factory pattern in C++ from a working, message-passing program in C.

Reaction dynamics in exothermic Co/Al multilayer foils are studied with high speed digital photography. Unsteady, spin-like reaction propagation is observed in which the net synthesis of a foil is accomplished through advancing transverse bands that propagate perpendicular to the net reaction direction. This unsteady behavior is connected to the final reacted foil surface morphology that exhibits periodic structures. The evolution of the reaction front shape and corresponding surface morphology are discussed with respect to Co/Al foil characteristics. © 2009 Materials Research Society.

The monolithic nature of transportation technologies offers opportunities for significant improvements in efficiency of 25-50% through strategic technical investments in both advanced fuels and new low-temperature engine concepts. The application of direct numerical simulation (DNS) provides a way to study fundamental issues related to small-scale combustion processes in well-defined canonical configurations, whereas the application of large eddy simulation (LES) provides a formal treatment of the full range of time and length scales that exist in turbulent reacting flows, and thus provides a direct link to experimental studies of relevant combustion devices. In the present study, through DOE INCITE and Oak Ridge National Laboratory 250 Tflop Transition-to-Operations grants in 2008, DNS is performed to understand how a lifted autoignitive flame is stabilized, and LES is performed to understand the high-pressure injection and mixing processes in internal combustion engines. Understanding of these and other fundamental issues is needed to develop robust and reliable ignition and combustion models for the combustion regimes observed under low-temperature combustion engine environments using alternative fuels. © 2008 IOP Publishing Ltd.

A coupled Euler-Lagrange solution approach is used to model the response of a buried reinforced concrete structure subjected to a close-in detonation of a high explosive charge. The coupling algorithm is discussed along with a set of benchmark calculations involving detonations in clay and sand.

The disturbed rock zone constitutes an important geomechanical element of the Waste Isolation Pilot Plant. The science and engineering underpinning the disturbed rock zone provide the basis for evaluating ongoing operational issues and their impact on performance assessment. Contemporary treatment of the disturbed rock zone applied to the evaluation of the panel closure system and to a new mining horizon improves the level of detail and quantitative elements associated with a damaged zone surrounding the repository openings. Technical advancement has been realized by virtue of ongoing experimental investigations and international collaboration. Initial sections summarize and document theoretical and experimental results, which quantify characteristics of the disturbed rock zone as applied to nuclear waste repositories in salt. This information is then applied to operational issues pertaining to recertification of the repository. © 2008 ASCE.

A total system performance assessment (TSPA) model has been developed to analyze the ability of the natural and engineered barriers of the Yucca Mountain repository to isolate nuclear waste over the period following repository closure. The principal features of the engineered barrier system are emplacement tunnels (or "drifts") containing a two-layer waste package (WP) for waste containment and a titanium drip shield to protect the WP from seeping water and falling rock. The 25-mm-thick outer shell of the WP is composed of Alloy 22, a highly corrosion-resistant nickel-based alloy. There are five nominal degradation modes of the Alloy 22: general corrosion, microbially influenced corrosion, stress corrosion cracking, early failure due to manufacturing defects, and localized corrosion (LC). This paper specifically examines the incorporation of the Alloy 22 LC model into the Yucca Mountain TSPA model, particularly the abstraction and modeling methodology, as well as issues dealing with scaling, spatial variability, uncertainty, and coupling to other submodels that are part of the total system model, such as the submodel for seepage water chemistry.

Future energy systems based on gasification of coal or biomass for co-production of electrical power and fuels may require gas turbine operation on unusual gaseous fuel mixtures. In addition, global climate change concerns may dictate the generation of a CO2 product stream for end-use or sequestration, with potential impacts on the oxidizer used in the gas turbine. In this study the operation at atmospheric pressure of a small, optically accessible swirl-stabilized premixed combustor, burning fuels ranging from pure methane to conventional and H2-rich and H2-lean syngas mixtures is investigated. Both air and CO2-diluted oxygen are used as oxidizers. CO and NOx emissions for these flames have been determined from the lean blowout limit to slightly rich conditions (1.03). In practice, CO2-diluted oxygen systems will likely be operated close to stoichiometric conditions to minimize oxygen consumption while achieving acceptable NOx performance. The presence of hydrogen in the syngas fuel mixtures results in more compact, higher temperature flames, resulting in increased flame stability and higher NOx emissions. Consistent with previous experience, the stoichiometry of lean blowout decreases with increasing H2 content in the syngas. Similarly, the lean stoichiometry at which CO emissions become significant decreases with increasing H2 content. For the mixtures investigated, CO emissions near the stoichiometric point do not become significant until 0.95. At this stoichiometric limit, CO emissions rise more rapidly for combustion in O2-CO2 mixtures than for combustion in air.

Understanding the relationship between chemical structure and function is a ubiquitous problem within the fields of chemistry and biology. Simulation approaches attack the problem utilizing physics to understand a given process at the particle level. Unfortunately, these approaches are often too expensive for many problems of interest. Informatics approaches attack the problem with empirical analysis of descriptions of chemical structure. The issue in these methods is how to describe molecules in a manner that facilitates accurate and general calculation of molecular properties. Here, we present a novel approach that utilizes aspects of simulation and informatics in order to formulate structure-property relationships. We show how supervised learning can be utilized to overcome the sampling problem in simulation approaches. Likewise, we show how learning can be achieved based on molecular descriptions that are rooted in the physics of dynamic intermolecular forces. We apply the approach to three problems including the analysis of corticosteroid binding globulin ligand binding affinity, identification of formylpeptide receptor ligands, and identification of resveratrol analogues capable of inhibiting activation of transcription factor nuclear factor kappaB. © 2008 American Chemical Society.

Abstract not provided.

In part I of this study we introduced a 17-parameter model that can predict heart rate regulation during postural change from sitting to standing. In this subsequent study, we focus on the 17 model parameters needed to adequately represent the observed heart rate response. In part I and in previous work (Olufsen et al. 2006), we estimated the 17 model parameters by minimizing the least squares error between computed and measured values of the heart rate using the Nelder-Mead method (a simplex algorithm). In this study, we compare the Nelder-Mead optimization method to two sampling methods: the implicit filtering method and a genetic algorithm. We show that these off-the-shelf optimization methods can work in conjunction with the heart rate model and provide reasonable parameter estimates with little algorithm tuning. In addition, we make use of the thousands of points sampled by the optimizers in the course of the minimization to perform an overall analysis of the model itself. Our findings show that the resulting least-squares problem has multiple local minima and that the non-linear-least squares error can vary over two orders of magnitude due to the complex interaction between the model parameters, even when provided with reasonable bound constraints. © Springer Science+Business Media, LLC 2007.

A method to measure interfacial mechanical properties at high temperatures and in a controlled atmosphere has been developed to study anodized aluminum surface coatings at temperatures where the interior aluminum alloy is molten. This is the first time that the coating strength has been studied under these conditions. We have investigated the effects of ambient atmosphere, temperature, and surface finish on coating strength for samples of aluminum alloy 7075. Surprisingly, the effective Young's modulus or strength of the coating when tested in air was twice as high as when samples were tested in an inert nitrogen or argon atmosphere. Additionally, the effective Young's modulus of the anodized coating increased with temperature in an air atmosphere but was independent of temperature in an inert atmosphere. The effect of surface finish was also examined. Sandblasting the surface prior to anodization was found to increase the strength of the anodized coating with the greatest enhancement noted for a nitrogen atmosphere. Machining marks were not found to significantly affect the strength.

Current published models for predicting squeeze film damping (SFD), which are based on different assumptions, give widely different results in the free-molecule regime. The work presented here provides experimental data for validating SFD models in that regime. The test device was an almost rectangular micro plate supported by beam springs. The structure was base-excited. The rigid plate oscillated vertically while staying parallel to the substrate. The velocities of the plate and of the substrate were measured with a laser Doppler vibrometer and a microscope. The damping ratio was calculated by performing modal analysis of the frequency response functions. The test structures were contained in a vacuum chamber with air pressures controlled to provide a five-order-of-magnitude range of Knudsen numbers. The damping coefficients from the measurements were compared with predictions from various published models. The results show that the continuum-base Reynolds equation predicts squeeze-film damping accurately if used with correct boundary conditions. The accuracy of molecular-based models depends heavily on the assumptions used in developing the models. Copyright © 2007 by ASME.

The shape control of thin, flexible structures has been studied primarily for edge-supported thin-plates. For applications such as electromagnetic wave reflectors, corner-supported configurations may prove more applicable since they allow for greater flexibility and larger achievable deflections when compared to edge-supported geometries under similar actuation conditions. Models of such structures provide insight for effective, realizable designs, enable design optimization, and provide a means of active shape control. Models for small deformations of corner-supported, thin laminates actuated by integrated piezoelectric actuators have been developed. However, membrane deflections expected for nominal actuation exceed those stipulated by linear, small deflection theories. In addition, large deflection models have been developed for membranes; however these models are not formulated for shape control. This paper extends a previously-developed linear model for a corner-supported thin, rectangular laminate to a more general large deflection model for a clamped-corner laminate composed of moment actuators and an array of actuating electrodes. First, a nonlinear model determining the deflected shape of a laminate given a distribution of actuation voltages is derived. Second, a technique is employed to formulate the model as a map between input voltage and deflection alone, making it suitable for shape control. Finally, comparisons of simulated deflections with measured deflections of a fabricated active laminate are investigated.

Optical filters based on resonant gratings have spectral characteristics that are lithographically defined. Nanoimprint lithography is a relatively new method for producing large area gratings with sub-micron features. Computational modeling using rigorous coupled-wave analysis allows gratings to be designed to yield sharp reflectance maxima and minima. Combining these gratings with microfluidic channels and micromechanical actuators produced using micro electromechanical systems (MEMS) technology forms the basis for producing tunable filters and other wavelength selective elements. These devices achieve tunable optical characteristics by varying the index of refraction on the surface of the grating. Coating the grating surface with water creates a 33% change in the resonant wavelength whereas bringing a grating into contact with a quartz surface shifts the resonant wavelength from 558 nm to 879 nm, a fractional change of 58%. The reflectivity at a single wavelength can be varied by approximately a factor of three. Future applications of these devices may include tunable filters or optical modulators. © 2008 Society of Photo-Optical Instrumentation Engineers.

We describe a method of performing trilinear analysis on large data sets using a modification of the PARAFAC-ALS algorithm. Our method iteratively decomposes the data matrix into a core matrix and three loading matrices based on the Tuckerl model. The algorithm is particularly useful for data sets that are too large to upload into a computer's main memory. While the performance advantage in utilizing our algorithm is dependent on the number of data elements and dimensions of the data array, we have seen a significant performance improvement over operating PARAFAC-ALS on the full data set. In one case of data comprising hyperspectral images from a confocal microscope, our method of analysis was approximately 60 times faster than operating on the full data set, while obtaining essentially equivalent results. Published in 2008 by John Wiley & Sons, Ltd.

Chemical solution deposition has been used to fabricate continuous ultrathin lead lanthanum zirconate titanate (PLZT) films as thin as 20 nm. Further, multilayer capacitor structures with as many as 10 dielectric layers have been fabricated from these ultrathin PLZT films by alternating spin-coated dielectric layers with sputtered platinum electrodes. Integrating a photolithographically defined wet etch step to the fabrication process enabled the production of functional multilayer stacks with capacitance values exceeding 600 nF. Such ultrathin multilayer capacitors offer tremendous advantages for further miniaturization of integrated passive components.

Advances are reported in several aspects of clathrate hydrate desalination fundamentals necessary to develop an economical means to produce municipal quantities of potable water from seawater or brackish feedstock. These aspects include the following, (1) advances in defining the most promising systems design based on new types of hydrate guest molecules, (2) selection of optimal multi-phase reactors and separation arrangements, and, (3) applicability of an inert heat exchange fluid to moderate hydrate growth, control the morphology of the solid hydrate material formed, and facilitate separation of hydrate solids from concentrated brine. The rate of R141b hydrate formation was determined and found to depend only on the degree of supercooling. The rate of R141b hydrate formation in the presence of a heat exchange fluid depended on the degree of supercooling according to the same rate equation as pure R141b with secondary dependence on salinity. Experiments demonstrated that a perfluorocarbon heat exchange fluid assisted separation of R141b hydrates from brine. Preliminary experiments using the guest species, difluoromethane, showed that hydrate formation rates were substantial at temperatures up to at least 12 C and demonstrated partial separation of water from brine. We present a detailed molecular picture of the structure and dynamics of R141b guest molecules within water cages, obtained from ab initio calculations, molecular dynamics simulations, and Raman spectroscopy. Density functional theory calculations were used to provide an energetic and molecular orbital description of R141b stability in both large and small cages in a structure II hydrate. Additionally, the hydrate of an isomer, 1,2-dichloro-1-fluoroethane, does not form at ambient conditions because of extensive overlap of electron density between guest and host. Classical molecular dynamics simulations and laboratory trials support the results for the isomer hydrate. Molecular dynamics simulations show that R141b hydrate is stable at temperatures up to 265K, while the isomer hydrate is only stable up to 150K. Despite hydrogen bonding between guest and host, R141b molecules rotated freely within the water cage. The Raman spectrum of R141b in both the pure and hydrate phases was also compared with vibrational analysis from both computational methods. In particular, the frequency of the C-Cl stretch mode (585 cm{sup -1}) undergoes a shift to higher frequency in the hydrate phase. Raman spectra also indicate that this peak undergoes splitting and intensity variation as the temperature is decreased from 4 C to -4 C.

This work utilized advanced engineering in several fields to find solutions to the challenges presented by the integration of MEMS/NEMS with optoelectronics to realize a compact sensor system, comprised of a microfabricated sensor, VCSEL, and photodiode. By utilizing microfabrication techniques in the realization of the MEMS/NEMS component, the VCSEL and the photodiode, the system would be small in size and require less power than a macro-sized component. The work focused on two technologies, accelerometers and microphones, leveraged from other LDRD programs. The first technology was the nano-g accelerometer using a nanophotonic motion detection system (67023). This accelerometer had measured sensitivity of approximately 10 nano-g. The Integrated NEMS and optoelectronics LDRD supported the nano-g accelerometer LDRD by providing advanced designs for the accelerometers, packaging, and a detection scheme to encapsulate the accelerometer, furthering the testing capabilities beyond bench-top tests. A fully packaged and tested die was never realized, but significant packaging issues were addressed and many resolved. The second technology supported by this work was the ultrasensitive directional microphone arrays for military operations in urban terrain and future combat systems (93518). This application utilized a diffraction-based sensing technique with different optical component placement and a different detection scheme from the nano-g accelerometer. The Integrated NEMS LDRD supported the microphone array LDRD by providing custom designs, VCSELs, and measurement techniques to accelerometers that were fabricated from the same operational principles as the microphones, but contain proof masses for acceleration transduction. These devices were packaged at the end of the work.

We analyze and compare findings from identical national surveys of the US general public on nuclear security and terrorism administered by telephone and Internet in mid-2007. Key areas of investigation include assessments of threats to US security; valuations of US nuclear weapons and nuclear deterrence; perspectives on nuclear proliferation, including the specific cases of North Korea and Iran; and support for investments in nuclear weapons capabilities. Our analysis of public views on terrorism include assessments of the current threat, progress in the struggle against terrorism, preferences for responding to terrorist attacks at different levels of assumed casualties, and support for domestic policies intended to reduce the threat of terrorism. Also we report findings from an Internet survey conducted in mid 2007 that investigates public views of US energy security, to include: energy supplies and reliability; energy vulnerabilities and threats, and relationships among security, costs, energy dependence, alternative sources, and research and investment priorities. We analyze public assessments of nuclear energy risks and benefits, nuclear materials management issues, and preferences for the future of nuclear energy in the US. Additionally, we investigate environmental issues as they relate to energy security, to include expected implications of global climate change, and relationships among environmental issues and potential policy options.

We report Lithium-Ion batteries are being considered as a high-energy density replacement for Nickel Metal Hydride (NiMH) batteries in Hybrid Electric Vehicles (HEVs) and in the new Plug-In Hybrids (PHEVs). Although these cells can result in significant reduction in weight and volume, they have several safety related issues that still need to be addressed. We report here on the thermal response of Li-ion cells specifically assembled in our laboratory to test new materials, electrolytes and additives. Finally, improvements in the thermal abuse tolerance of cells are reported and discussed in terms of the need for overall battery system safety.

A previous LDRD studying radiation hardened optoelectronic components for space-based applications led to the result that increased neutron irradiation from a fast-burst reactor caused increased responsivity in GaAs photodiodes up to a total fluence of 4.4 x 10{sup 13} neutrons/cm{sup 2} (1 MeV Eq., Si). The silicon photodiodes experienced significant degradation. Scientific literature shows that neutrons can both cause defects as well as potentially remove defects in an annealing-like process in GaAs. Though there has been some modeling that suggests how fabrication and radiation-induced defects can migrate to surfaces and interfaces in GaAs and lead to an ordering effect, it is important to consider how these processes affect the performance of devices, such as the basic GaAs p-i-n photodiode. In this LDRD, we manufactured GaAs photodiodes at the MESA facility, irradiated them with electrons and neutrons at the White Sands Missile Range Linac and Fast Burst Reactor, and performed measurements to show the effect of irradiation on dark current, responsivity and high-speed bandwidth.

Abstract not provided.

An analytic model for electron flow in a system driving a fixed inductive load is described and evaluated with particle in cell simulations. The simple model allows determining the impedance profile for a magnetically insulated transmission line given the minimum gap desired, and the lumped inductance inside the transition to the minimum gap. The model allows specifying the relative electron flow along the power flow direction, including cases where the fractional electron flow decreases in the power flow direction. The electrons are able to return to the cathode because they gain energy from the temporally rising magnetic field. The simulations were done with small cell size to reduce numerical heating. An experiment to compare electron flow to the simulations was done. The measured electron flow is {approx}33% of the value from the simulations. The discrepancy is assumed to be due to a reversed electric field at the cathode because of the inductive load and falling electron drift velocity in the power flow direction. The simulations constrain the cathode electric field to zero, which gives the highest possible electron flow.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Understanding and controlling friction in micromachine interfaces is critical to the reliability and operational efficiency of microelectromechanical systems (MEMS). The relatively high adhesion forces and friction forces encountered in these devices often present major obstacles to the design of reliable MEMS devices. Using surface micromachining, arrays of microstructures are being designed and tested to examine the adhesion characteristics, static friction behavior, and dynamic friction response. Emphasis is also being given to the control and actuation of the test structures and the modeling of the dynamic response and contact mechanics at the interface. Specifically, the purpose of the research is to fabricate and test MEMS devices in order to obtain insight into the effect of surface topography, material properties, surface chemical state, environmental conditions, and contact load on the static and dynamic characteristics of the contact interface.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This paper presents an overview of algorithms for directing messages through networks of varying topology. These are commonly referred to as routing algorithms in the literature that is presented. In addition to providing background on networking terminology and router basics, the paper explains the issues of deadlock and livelock as they apply to routing. After this, there is a discussion of routing algorithms for both store-and-forward and wormhole-switched networks. The paper covers both algorithms that do and do not adapt to conditions in the network. Techniques targeting structured as well as irregular topologies are discussed. Following this, strategies for routing in the presence of faulty nodes and links in the network are described.

Emerging high-bandwidth, low-latency network technology has made network-based architectures both feasible and potentially desirable for use in satellite payload architectures. The selection of network topology is a critical component when developing these multi-node or multi-point architectures. This study examines network topologies and their effect on overall network performance. Numerous topologies were reviewed against a number of performance, reliability, and cost metrics. This document identifies a handful of good network topologies for satellite applications and the metrics used to justify them as such. Since often multiple topologies will meet the requirements of the satellite payload architecture under development, the choice of network topology is not easy, and in the end the choice of topology is influenced by both the design characteristics and requirements of the overall system and the experience of the developer.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Peridynamics is a nonlocal formulation of continuum mechanics. The discrete peridynamic model has the same computational structure as a molecular dynamic model. This document details the implementation of a discrete peridynamic model within the LAMMPS molecular dynamic code. This document provides a brief overview of the peridynamic model of a continuum, then discusses how the peridynamic model is discretized, and overviews the LAMMPS implementation. A nontrivial example problem is also included.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report documents work done for a late-start LDRD project, which was carried out during the last quarter of FY07. The objective of this project was to experimentally explore the feasibility of converting vegetable (e.g., soybean) oils to biodiesel by employing slit-channel reactors and solid catalysts. We first designed and fabricated several slit-channel reactors with varying channel depths, and employed them to investigate the improved performance of slit-channel reactors over traditional batch reactors using a NaOH liquid catalyst. We then evaluated the effectiveness of several solid catalysts, including CaO, ZnO, MgO, ZrO{sub 2}, calcium gluconate, and heteropolyacid or HPA (Cs{sub 2.5}H{sub 0.5}PW{sub 12}O{sub 40}), for catalyzing the soybean oil-to-biodiesel transesterification reaction. We found that the slit-channel reactor performance improves as channel depth decreases, as expected; and the conversion efficiency of a slit-channel reactor is significantly higher when its channel is very shallow. We further confirmed CaO as having the highest catalytic activity among the solid catalysts tested, and we demonstrated for the first time calcium gluconate as a promising solid catalyst for converting soybean oil to biodiesel, based on our preliminary batch-mode conversion experiments.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This LDRD began as a three year program to integrate through-wafer vias, micro-mirrors and control electronics with high-voltage capability to yield a 64 by 64 array of individually controllable micro-mirrors on 125 or 250 micron pitch with piston, tip and tilt movement. The effort was a mix of R&D and application. Care was taken to create SUMMiT{trademark} (Sandia's ultraplanar, multilevel MEMS technology) compatible via and mirror processes, and the ultimate goal was to mate this MEMS fabrication product to a complementary metal-oxide semiconductor (CMOS) electronics substrate. Significant progress was made on the via and mirror fabrication and design, the attach process development as well as the electronics high voltage (30 volt) and control designs. After approximately 22 months, the program was ready to proceed with fabrication and integration of the electronics, final mirror array, and through wafer vias to create a high resolution OMEMS array with individual mirror electronic control. At this point, however, mission alignment and budget constraints reduced the last year program funding and redirected the program to help support the through-silicon via work in the Hyper-Temporal Sensors (HTS) Grand Challenge (GC) LDRD. Several months of investigation and discussion with the HTS team resulted in a revised plan for the remaining 10 months of the program. We planned to build a capability in finer-pitched via fabrication on thinned substrates along with metallization schemes and bonding techniques for very large arrays of high density interconnects (up to 2000 x 2000 vias). Through this program, Sandia was able to build capability in several different conductive through wafer via processes using internal and external resources, MEMS mirror design and fabrication, various bonding techniques for arrayed substrates, and arrayed electronics control design with high voltage capability.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Vapor Phase Lubrication (VPL) of silicon surfaces with pentanol has been demonstrated. Two potential show stoppers with respect to application of this approach to real MEMS devices have been investigated. Water vapor was found to reduce the effectiveness of VPL with alcohol for a given alcohol concentration, but the basic reaction mechanism observed in water-free environments is still active, and devices operated much longer in mixed alcohol and water vapor environments than with chemisorbed monolayer lubricants alone. Complex MEMS gear trains were successfully lubricated with alcohol vapors, resulting in a factor of 104 improvement in operating life without failure. Complex devices could be made to fail if operated at much higher frequencies than previously used, and there is some evidence that the observed failure is due to accumulation of reaction products at deeply buried interfaces. However, if hypothetical reaction mechanisms involving heated surfaces are valid, then the failures observed at high frequency may not be relevant to operation at normal frequencies. Therefore, this work demonstrates that VPL is a viable approach for complex MEMS devices in conventional packages. Further study of the VPL reaction mechanisms are recommended so that the vapor composition may be optimized for low friction and for different substrate materials with potential application to conventionally fabricated, metal alloy parts in weapons systems. Reaction kinetics should be studied to define effective lubrication regimes as a function of the partial pressure of the vapor phase constituent, interfacial shear rate, substrate composition, and temperature.

This document provides general guidance for the design and analysis of bolted joint connections. An overview of the current methods used to analyze bolted joint connections is given. Several methods for the design and analysis of bolted joint connections are presented. Guidance is provided for general bolted joint design, computation of preload uncertainty and preload loss, and the calculation of the bolted joint factor of safety. Axial loads, shear loads, thermal loads, and thread tear out are used in factor of safety calculations. Additionally, limited guidance is provided for fatigue considerations. An overview of an associated Mathcad{copyright} Worksheet containing all bolted joint design formulae presented is also provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Synthetic Aperture Radar (SAR) performance testing and estimation is facilitated by observing the system response to known target scene elements. Trihedral corner reflectors and other canonical targets play an important role because their Radar Cross Section (RCS) can be calculated analytically. However, reflector orientation and the proximity of the ground and mounting structures can significantly impact the accuracy and precision with which measurements can be made. These issues are examined in this report.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The annual program report provides detailed information about all aspects of the Sandia National Laboratories, California (SNL/CA) Environmental Planning and Ecology Program for a given calendar year. It functions as supporting documentation to the SNL/CA Environmental Management System Program Manual. The program report describes the activities undertaken during the past year, and activities planned in future years to implement the Planning and Ecology Program, one of six programs that supports environmental management at SNL/CA.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The United States produces only about 1/3 of the more than 20 million barrels of petroleum that it consumes daily. Oil imports into the country are roughly equivalent to the amount consumed in the transportation sector. Hence the nation in general, and the transportation sector in particular, is vulnerable to supply disruptions and price shocks. The situation is anticipated to worsen as the competition for limited global supplies increases and oil-rich nations become increasingly willing to manipulate the markets for this resource as a means to achieve political ends. The goal of this project was the development and improvement of technologies and the knowledge base necessary to produce and qualify a universal fuel from diverse feedstocks readily available in North America and elsewhere (e.g. petroleum, natural gas, coal, biomass) as a prudent and positive step towards mitigating this vulnerability. Three major focus areas, feedstock transformation, fuel formulation, and fuel characterization, were identified and each was addressed. The specific activities summarized herein were identified in consultation with industry to set the stage for collaboration. Two activities were undertaken in the area of feedstock transformation. The first activity focused on understanding the chemistry and operation of autothermal reforming, with an emphasis on understanding, and therefore preventing, soot formation. The second activity was focused on improving the economics of oxygen production, particularly for smaller operations, by integrating membrane separations with pressure swing adsorption. In the fuel formulation area, the chemistry of converting small molecules readily produced from syngas directly to fuels was examined. Consistent with the advice from industry, this activity avoided working on improving known approaches, giving it an exploratory flavor. Finally, the fuel characterization task focused on providing a direct and quantifiable comparison of diesel fuel and JP-8.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The author will describe two-photon-resonant LIF detection of CO, O, and H. Application of these techniques in flames frequently suffers from significant photolytic interferences caused by the intense UV excitation pulses required to produce measurable signal. When compared to nanosecond excitation, the use of short pulse (picosecond) excitation can significantly reduce the effect of the photolytic interference. Results of recent atomic oxygen imaging experiments using picosecond- and nanosecond-duration laser pulses will be presented, and potential improvements to CO and H imaging will be discussed.

A Dish Stirling parabolic concentrator typically consists of a number of mirror facets that must be aligned to focus the concentrated sunlight on the engine receiver. An alignment strategy must be developed to deliver the energy uniformly to the receiver while maximizing system performance. Several criteria must be met in order to maximize the performance and lifetime of the system. The peak flux should be minimized at the receiver to extend life. This is accomplished by locally optimizing the mirror aimpoints, minimizing overlap of facet images. The energy delivered to each cylinder of a multi-cylinder engine should be balanced to maximize the power production capability of the engine. This is accomplished through globally optimizing the mirror aimpoints. Depending on dish geometry, both of these constraints will be met by moving the aimpoints of certain facets away from a single point at the center of the aperture. However, this often results in a larger aperture or more flux spillage. The larger aperture results in greater thermal and reflective losses from the receiver cavity. This paper proposes and demonstrates a novel approach to optimizing the alignment strategy while obeying these constraints. The method uses an approach similar to molecular dynamics to globally and locally distribute the power on the receiver, while imposing movement constraints at the aperture to limit the focal plane spot size. The method can also impose additional geometric constraints at the receiver plane to accommodate un-cooled surfaces. The method is explored and demonstrated on the Stirling Energy Systems 25kW dish Stirling system at Sandia National Laboratories. The approach provides a receiver flux distribution and power balance equal to the strategy developed by McDonnell Douglas in the early 1980's, but with an aperture size equal to that of the single aimpoint strategy. This should result in about a 1kW increase in power generated at rated conditions, with no additional cost, due to reduced thermal losses from the receiver. The method can be extended to other point-focus concentrating solar technologies. On a tower, the heliostat aiming strategy could be dynamically updated to accommodate flux profile needs, sun position, or maintenance in the field. Copyright © 2007 by US Government.

Within this paper, we provide a solution to the static frame validation challenge problem (see this issue) in a manner that is consistent with the guidelines provided by the Validation Challenge Workshop tasking document. The static frame problem is constructed such that variability in material properties is known to be the only source of uncertainty in the system description, but there is ignorance on the type of model that best describes this variability. Hence both types of uncertainty, aleatoric and epistemic, are present and must be addressed. Our approach is to consider a collection of competing probabilistic models for the material properties, and calibrate these models to the information provided; models of different levels of complexity and numerical efficiency are included in the analysis. A Bayesian formulation is used to select the optimal model from the collection, which is then used for the regulatory assessment. Lastly, bayesian credible intervals are used to provide a measure of confidence to our regulatory assessment.

This paper recalls the events leading up to the author's 1973 discovery of Deep Level Transient Spectroscopy (DLTS). It discusses the status of junction capacitance techniques in the late 1960s and points out why the typical capacitance instrumentation of that era would not have lead the author to the DLTS discovery. This discovery is discussed in the context of the novel NMR-inspired instrumentation used by the author to study fast capacitance transients of the ZnO center in GaP LEDs. Finally, the author makes some general comments about the innovation process. © 2007 Elsevier B.V. All rights reserved.

A standard approach to cross-language information retrieval (CLIR) uses Latent Semantic Analysis (LSA) in conjunction with a multilingual parallel aligned corpus. This approach has been shown to be successful in identifying similar documents across languages - or more precisely, retrieving the most similar document in one language to a query in another language. However, the approach has severe drawbacks when applied to a related task, that of clustering documents 'language independently', so that documents about similar topics end up closest to one another in the semantic space regardless of their language. The problem is that documents are generally more similar to other documents in the same language than they are to documents in a different language, but on the same topic. As a result, when using multilingual LSA, documents will in practice cluster by language, not by topic. We propose a novel application of PARAFAC2 (which is a variant of PARAFAC, a multi-way generalization of the singular value decomposition [SVD]) to overcome this problem. Instead of forming a single multilingual term-by-document matrix which, under LSA, is subjected to SVD, we form an irregular three-way array, each slice of which is a separate term-by-document matrix for a single language in the parallel corpus. The goal is to compute an SVD for each language such that V (the matrix of right singular vectors) is the same across all languages. Effectively, PARAFAC2 imposes the constraint, not present in standard LSA, that the 'concepts' in all documents in the parallel corpus are the same regardless of language. Intuitively, this constraint makes sense, since the whole purpose of using a parallel corpus is that exactly the same concepts are expressed in the translations. We tested this approach by comparing the performance of PARAFAC2 with standard LSA in solving a particular CLIR problem. From our results, we conclude that PARAFAC2 offers a very promising alternative to LSA not only for multilingual document clustering, but also for solving other problems in crosslanguage information retrieval. © 2007 ACM.

The rotational dynamics of chemically similar systems based on freely jointed and freely rotating chains are studied. The second Legendre polynomial of vectors along chain backbones is used to investigate the rotational dynamics at different length scales. In a previous study, it was demonstrated that the additional bond-angle constraint in the freely rotating case noticeably perturbs the character of the translational relaxation away from that of the freely jointed system. Here, it is shown that differences are also apparent in the two systems' rotational dynamics. The relaxation of the end-to-end vector is found to display a long time, single-exponential tail and a stretched exponential region at intermediate times. The stretching exponents Β are found to be 0.75±0.02 for the freely jointed case and 0.68±0.02 for the freely rotating case. For both system types, time-packing-fraction superposition is seen to hold on the end-to-end length scale. In addition, for both systems, the rotational relaxation times are shown to be proportional to the translational relaxation times, demonstrating that the Debye-Stokes-Einstein law holds. The second Legendre polynomial of the bond vector is used to probe relaxation behavior at short length scales. For the freely rotating case, the end-to-end relaxation times scale differently than the bond relaxation times, implying that the behavior is non-Stokes-Einstein, and that time-packing-fraction superposition does not hold across length scales for this system. For the freely jointed case, end-to-endrelaxation times do scale with bond relaxation times, and both Stokes-Einstein and time-packing-fraction-across-length-scales superposition are obeyed. © 2007 American Institute of Physics.

A standard approach to cross-language information retrieval (CLIR) uses Latent Semantic Analysis (LSA) in conjunction with a multilingual parallel aligned corpus. This approach has been shown to be successful in identifying similar documents across languages - or more precisely, retrieving the most similar document in one language to a query in another language. However, the approach has severe drawbacks when applied to a related task, that of clustering documents 'language independently', so that documents about similar topics end up closest to one another in the semantic space regardless of their language. The problem is that documents are generally more similar to other documents in the same language than they are to documents in a different language, but on the same topic. As a result, when using multilingual LSA, documents will in practice cluster by language, not by topic. We propose a novel application of PARAFAC2 (which is a variant of PARAFAC, a multi-way generalization of the singular value decomposition [SVD]) to overcome this problem. Instead of forming a single multilingual term-by-document matrix which, under LSA, is subjected to SVD, we form an irregular three-way array, each slice of which is a separate term-by-document matrix for a single language in the parallel corpus. The goal is to compute an SVD for each language such that V (the matrix of right singular vectors) is the same across all languages. Effectively, PARAFAC2 imposes the constraint, not present in standard LSA, that the 'concepts' in all documents in the parallel corpus are the same regardless of language. Intuitively, this constraint makes sense, since the whole purpose of using a parallel corpus is that exactly the same concepts are expressed in the translations. We tested this approach by comparing the performance of PARAFAC2 with standard LSA in solving a particular CLIR problem. From our results, we conclude that PARAFAC2 offers a very promising alternative to LSA not only for multilingual document clustering, but also for solving other problems in crosslanguage information retrieval. © 2007 ACM.

In this work, we present two methods for solving overdetermined systems of the Time of Arrival (TOA) geolocation equations that achieve the minimum possible variance in all cases, not just when the satellites are at large equal radii. One of these techniques gives two solutions, and the other gives four solutions.

We measure frictional properties of liquid-expanded and liquid-condensed phases of lipid Langmuir-Blodgett monolayers by interfacial force microscopy. We find that over a reasonably broad surface-density range, the friction shear strength of the lipid monolayer film is proportional to the surface area (42-74 Â2/ molecule) occupied by each molecule. The increase in frictional force (i.e., friction shear strength with molecular area can be attributed to the increased conformational freedom and the resulting increase in the number of available modes for energy dissipation. © 2007 American Chemical Society.

A three-component balance system has been developed and implemented to measure the forces and moments on a sub-scale missile fin model interacting with the wake and shed vortex from an upstream fin. Measurements were made from Mach 0.5 - 0.8 with both the upstream and downstream fins pitched between -5° and 10° angle of attack. The results show that the downstream fin's forces and moments are shifted from the baseline single fin values dependent on the angle of attack of the upstream fin. Mach Number had only a secondary effect and its influence was found to grow stronger as the angles of attack of the upstream and downstream fins diverged.

A model for the simulation of the electron energy distribution in nanoscale metal-oxide-semiconductor field-effect transistor (MOSFET) devices, using a kinetic simulation technique, is implemented. The convective scheme (CS), a method of characteristics, is an accurate method of solving the Boltzmann transport equation, a nonlinear integrodifferential equation, for the distribution of electrons in a MOSFET device. The method is used to find probabilities for use in an iterative scheme which iterates to find collision rates in cells. The CS is also a novel approach to 2D semiconductor device simulation. The CS has been extended to handle boundary conditions in 2D as well as to calculation of polygon overlap for polygons of more than three sides. Electron energy distributions in the channel of a MOSFET are presented. © 2007 Elsevier Inc. All rights reserved.

The immersed- Bz diode is being developed as a high-brightness, flash x-ray radiography source at Sandia National Laboratories. This diode is a foil-less electron-beam diode with a long, thin, needlelike cathode that is inserted into the bore of a solenoid. The solenoidal magnetic field guides the electron beam emitted from the cathode to the anode while maintaining a small beam radius. The electron beam strikes a thin, high-atomic-number anode and produces forward-directed bremsstrahlung. In addition, electron beam heating of the anode produces surface plasmas allowing ion emission. Two different operating regimes for this diode have been identified: a nominal operating regime where the total diode current is characterized as classically bipolar and an anomalous operating regime characterized by a dramatic impedance collapse where the total diode current greatly exceeds the bipolar limit. Data from a comprehensive series of experiments fielded at 4 and 5 MV, where the diode operates in the nominal or stable impedance regime, with beam currents ranging from 20-40 kA on target are presented. In this mode, both the measured diode current and experimental radiation production are consistent with physics based models including two-dimensional particle-in-cell simulations. The analysis indicates that intermediate mass ions (e.g., 12-18 amu) control the nominal impedance evolution rather than expected lighter mass ions such as hydrogen. © 2007 American Institute of Physics.

A preliminary study was conducted which considered capturing carbon dioxide from fossil-fired power plants and combining it with nuclear hydrogen in order to produce alternative liquid fuels for transportation. Among the alternative liquid hydrocarbons which can be used as fuel in internal combustion engines, the two that are most promising are methanol and ethanol. We choose these two because they are relatively simple compounds and can be used with only minor changes to the fuel systems of most automobiles today. In fact, there are some vehicles today which can operate with any combination of conventional gasoline, ethanol, or methanol. We estimated the quantity of carbon dioxide that would be emitted by fossil-fired power plants in the future. We then use this information to determine how much ethanol or methanol can be created if enough hydrogen is made available. Using the quantity of hydrogen required and the thermodynamics of the reactions involved, we estimate the nuclear power that would be needed to produce the liquid fuel. This amount of liquid fuel is then used to estimate the effect of such a program on conventional gasoline usage, need for foreign oil, and decrease in CO 2 emissions.

Study of the Triggered Plasma Opening Switch (TPOS) characteristics is in progress via an ion current collection diagnostic (ICCD), in addition to offline apparatus. This initial ion current collection diagnostic has been designed, fabricated, and tested on the TPOS in order to explore the opening profile of the main switch. The initial ion current collection device utilizes five collectors which are positioned perpendicularly to the main switch stage in order to collect radially traveling ions. It has been shown through analytical prowess that this specific geometry can be treated as a planar case of the Child-Langmuir law with only a 6% deviation from the cylindrical case. Additionally, magnetostatic simulations with self consistent space charge emitting surfaces of the main switch using the Trak code are under way. It is hoped that the simulations will provide evidence in support of both the analytical derivations and experimental data. Finally, an improved design of the ICCD (containing 12 collectors in the axial direction) is presently being implemented. © 2005 IEEE.

The BROOM system was developed to collect, manage and analyze information from bioterrorist attacks on strategic buildings. GIS features help decision-makers and analysts rapidly assess the current status of contaminated facilities and develop optimized cleanup strategies. BROOM consists of networked server, desktop and PDA components. PDAs are deployed to collect samples of suspected bioagents, such as anthrax. Novel geostatistical methods are used to generate contaminant maps and define optimum locations for subsequent sampling. Efficiency and accuracy gains witnessed in field tests show that GIS technology can play a vital role in visualizing, managing and analyzing data from bioterrorism incidents. © 2007 The Haworth Press, Inc. All rights reserved.

The ZR accelerator is a refurbishment of Sandia National Laboratories Z accelerator [1]. The ZR accelerator components were designed using electrostatic and circuit modeling tools. Transient electromagnetic modeling has played a complementary role in the analysis of ZR components [2]. In this paper we describe a 3D transient electromagnetic analysis of the ZR water convolute and stack using edge-based finite element techniques. © 2005 IEEE.

The Z driver at Sandia National Laboratories delivers one to two megajoules of electromagnetic energy inside its ∼10 cm radius final feed in 100 ns. The high current (∼20 MA) at small diameter produces magnetic pressures well above yield strengths for metals. The metal conductors stay in place due to inertia long enough to deliver current to the load. Within milliseconds however, fragments of metal escape the load region at high velocity. Much of the hardware and diagnostics inside the vacuum chamber is protected from this debris by blast shields with small view ports, and fast-closing valves. The water-vacuum insulator requires different protection because the transmission line debris shield should not significantly raise the inductance or perturb the self-magnetically insulated electron flow. This report shows calculations and results from a design intended to protect the insulator assembly. © 2005 IEEE.

The Triggered Plasma Opening Switch (TPOS) at SNL is a unique device that exploits the high conductivity and low mass properties of plasma. The TPOS's objective is to take the initial ∼0.8MA (∼250ns rise time) storage inductor current and deliver ∼0.5MA at ∼2.4MV (∼10ns rise time) to a load of ∼5-10Ω. Configuration advantages include low current jitter and resistive voltage drop, power gain, and minimization of trigger input power as the result of using two stages in series. This two-stage design is novel and is the first to demonstrate operation of magnetically triggered stages. Study of TPOS characteristics is in progress via an offline interferometer diagnostic; specifically, a laser interferometer will be used to make density measurements of the source plasma. It is thought that the gross plasma source density is ∼1014 cm-3, but details of the spatial structure and temporal evolution have not previously been studied. In order to better understand switch operation, these details are essential. Presently two interferometer systems are planned for testing: a temporary 1 μm system for initial plasma characterization, and a 10.6 mu;m laser system for routine use. We will start with a single chord measurement then upgrade to a multi-chord system. Future plans involve varying plasma source parameters, such as magnetic field strength and plasma fill time, in order to understand the density dependence on these parameters. Improved knowledge of the plasma source density behavior should allow for improved switch operation. © 2005 IEEE.

Although there is much written in regards to voltage breakdown of polymeric insulators under AC and DC conditions, much less is written involving Rexolite® (1422) [1], non-uniform field geometries, and impulse conditions. Yet, in order to design optimized pulsed power systems with some desired degree of reliability, understanding the behavior of this type of insulating system is needed. Specifically, Sandia National Laboratory's ZR project, which will use anode plugs in the vacuum stack (thus increasing the electrical stress in the Rexolite insulators), needs to be able to estimate the reliability of these vacuum stack insulators [2]. In an effort to estimate the insulator's lifetime small scale testing is in progress. Nine samples have been tested so far and at least ten more will be tested. Results from the current testing suggest that the Rexolite "ages" from pulse to pulse, that there is some volume dependence on breakdown strength, and that the electrode-vacuum-insulator interface has an affect on the insulator lifetime. ©2005 IEEE.

The Refurbished Z machine, ZR, has 3 distinct missions. The first is to increase the capability for the user community by providing higher peak currents. The second is to increase the precision by increasing pulse repeatability and pulse shape flexibility. The third is to increase the capacity by providing operational turnaround time consistent with conducting a shot per shin. The pursuit of the third mission relies heavily on the reliability of the components and well defined maintenance cycles. To test the performance of the ZR design, a system assessment module has been tested and found to meet the capability mission. The system assessment module will be used to test the reliability of the components comprising the ZR pulsed power modules. To assess the program goal of 1 failure in 50 ZR shots we plan to perform 7200 shots of the system assessment module, Z20. At a typical shot rate, this task would take approximately three years. To minimize the impact on the facility while obtaining the required reliability data, the system will be configured to fire and reset autonomously with the ultimate goal of unmanned operation. A systems approach was developed using National Instruments LabVIEW ® software and Field Point I/O hardware. All communications between subsystems are provided via ethernet using fiber optic media converters. The major subsystems for operating the module which will be described are the gas pressure and purge, high voltage power supplies, oil diverter operation, triggering, precision monitoring of Marx system, personnel access control and remote operation of the laser trigger system. © 2005 IEEE.

The ZR refurbishment project [1] at Sandia National Laboratories (SNL) required a set of diverter switches to protect the Marx generators and intermediate storage (IS) capacitors from Marx pre-fire and/or laser triggered output switch (LTS) no-fire. Thirty-six such diverters, one for each Marx-IS set, will need to operate reliably over the full range of Marx charge voltages and LTS anticipated closure times. Operating voltage is up to 6 MV. A self-closing oil switch diverter was selected and design work began in late 2002. The first diverter (Phase I or just P1) was delivered in the summer of 2003 and tested on SNL's Z20 test-bed. Based on test results, operational experience and overall project budgetary concerns, it was decided to re-design the diverter, resulting a simpler, less costly switch. This new self-closing oil switch (Phase II or P2) was fielded at SNL on the Z20 test-bed in late 2004. Both designs include adjustable electrodes to control the closure time. Also incorporated is a mechanical clamp that minimizes or shorts the oil gap until Marx charge is complete. Both diverters feature liquid resistors sized to safely absorb the energy stored in the Marx or IS. This paper describes the design and test results from these diverters. © 2005 IEEE.

The six-cell RITS-6 accelerator is an upgrade of the existing RITS-3 accelerator and is next in the sequence of Sandia IVA accelerators built to investigate/validate critical accelerator and radiographic diode issues for scaling to the Radiographic Integrated Test Stand (RITS) (nominally 16 MV, 156 kA, and 70 ns). In the RITS-6 upgrade to RITS-3 the number of cells/cavities, PFLs, laser triggered gas switches and intermediate stores is being doubled. A rebuilt single 61-nF Marx generator will charge the two intermediate storage capacitors. The RITS-3 experiments have demonstrated a MITL configuration matched to the PFL/induction cell impedance and a higher impedance MITL. RITS-6 is designed to utilize the higher impedance MITL providing a 10.5-MV, 123-kA output. The three years of pulsed power performance data from RITS-3 will be summarized and the design improvements being incorporated into RITS-6 will be outlined. The predicted output voltage and current for RITS-6 as a function of diode impedance will be shown. Particle-in-cell simulations of the vacuum power flow from the cell to the load for a range of diode impedances from matched to ∼40 Ohms will be shown and compared with the re-trapped parapotential flow predictions. The status of the component fabrication and system integration will be given. Another potential upgrade under consideration is RITS-62. In this case the RITS-6 Marx, intermediate stores, gas switches, and PFLs would be duplicated and a tee would replace the elbow that now connects a single PFL to a cell thereby allowing two PFLs to be connected to one cell. The output of RITS-62 matched to the cell/PFL impedance would then be 8 MV, 312 kA or 25.6 ohms. The predicted operating curves for RITS-62 with other non-matched MITLs will be shown. The power delivered to a radiographic diode can be maximized by the correct choice of MITL impedance given the cell/PFL and radiographic diode impedances. If the radiated output for a given diode has a stronger than linear voltage dependence this dependence can also be included in the correct choice of MITL impedance. The optimizations and trade-offs will be shown for RITS-6 and RITS-62 for diode impedances characteristic of radiographic diodes. © 2005 IEEE.

The University of Missouri Terawatt Test Stand (MUTTS) has conducted many untriggered experiments on a Rimfire gas switch scaled to 2.5 MV. The focus of these experiments was to evaluate what methods may be used to control the distribution of cascade arcs. The untriggered data indicates that the rise time of switch current does not statistically improve, as expected, as the number of cascade arcs per gap increased beyond two channels. For the same data, the number of arcs in the cascade section more dramatically affects the output current period. This indicates that in late time increased multichanneling has a more pronounced effect than in early time. The switch is triggered with a frequency quadrupled Nd:YAG laser at 30 mJ with a 3-5 ns pulse width. Since the focused laser does not ionize the full length of the trigger section, there is little effect on current rise time when compared to untriggered data, but more channels form in the cascade section for an air filled switch. The cascade section was shorted and data are presented describing the contribution of the single channeling trigger section to overall switch impedance. The electrical effects of multichanneling using a laser trigger, the formation of arc channels in the cascade section, and the implications the results have on the future design of fast gas switches are discussed. © 2005 IEEE.

The ZR gas switch, located between the Intermediate Store capacitor (I-Store) and the Pulsed Forming Line (PFL), requires a laser pulse for its triggering. There are several routes for the beam to reach the gas switch but all of them cross over the high voltage regions. The Z laser tube crosses over the outer to inner PFL electrodes with a voltage difference no larger than 3.5 MV. The ZR gas switch was designed to be in oil, given the higher operational voltages, as a consequence the laser tube is in the oil side of the PFL interface. The ZR laser tube is required to hold in excess of 5 MV across it using high pressure SF6 gas, the ID is 2.5″ to accommodate the laser beam, mechanically should tolerate the non-axial shock loading during the water switches firing. After a couple of iterations it was decided to use Polyurethane, it provided most of the desired mechanical properties, except that it outgases ether and ether based compounds. The effect of just a few ppm of ether on SF6 is a significant reduction on the HV hold off especially surface tracking or flashover. As a consequence the final design is such that the electric field distribution on the tube is as conservative as it was possible due to space constrains. We present the basic design, the field distribution, its relationship with available SF6 breakdown data and the present performance. © 2005 IEEE.

The Z Refurbishment project is designed to increase the peak current to the load on Z to ∼26 MA in a 100-ns wide power pulse. This current is achieved by summing the current from 36 independent pulse-power modules. To meet these requirements, we have designed and constructed an SF6-insulated gas switch that can hold off 5.5 MV and conduct a peak current of 600 kA for over a hundred shots. The gas switch is charged by a Marx generator in ∼1 microsecond and transfers about 200-kilojoules of energy and 0.25 Coulombs of charge to a pulse-forming line in a ∼150-ns-wide power pulse peaking at 2.5 TW. The gas switch oonsists of a laser-triggered section holding off 15% of the voltage followed by 25 self-breakdown gaps. The self-breaking gaps are designed to provide multiple breakdown arcs in order to lower the overall inductance of the switch. The gas switch is submerged in transformer oil during operation. In this work, we show how simulation and experiment have worked together, first to verify proper operation of the switch, and then to solve problems with the switch design that arose during testing. © 2005 IEEE.

Sandia National Laboratories is investigating and developing high-dose, high-brightness flash radiographic sources. The immersed-Bz diode employs large-bore, high-field solenoid magnets to help guide and confine an intense electron beam from a needle-like cathode "immersed" in the axial field of the magnet. The electron beam is focused onto a high-atomic-number target/anode to generate an intense source of bremsstrahlung X-rays. Historically, these diodes have been unable to achieve high dose (> 500 rad @ m) from a small spot (< 3 mm diameter). It is believed that this limitation is due in part to undesirable effects associated with the interaction of the electron beam with plasmas formed at either the anode or the cathode. Previous research concentrated on characterizing the behavior of diodes, which used untreated, room temperature (RT) anodes. Research is now focused on improving the diode performance by modifying the diode behavior by using cryogenic anodes that are coated in-situ with frozen gases. The objective of these cryogenically treated anodes is to control and limit the ion species of the anode plasma formed and hence the species of the counter-streaming ions that can interact with the electron beam. Recent progress in the development, testing and fielding of the cryogenically cooled immersed diodes at Sandia is described. ©2005 IEEE.

This paper develops a novel control system design methodology that uniquely combines: concepts from thermodynamic exergy and entropy; Hamiltonian systems; Lyapunov's direct method and Lyapunov optimal analysis; electric AC power concepts; and power flow analysis. Relationships are derived between exergy/entropy and Lyapunov optimal functions for Hamiltonian systems. The methodology is demonstrated with a couple of fundamental numerical simulation examples: 1) a PID regulator control for a linear mass-spring-damper system and 2) a Duffing oscillator/Coulomb friction nonlinear model that employs PID regulator control. The control system performances are partitioned and evaluated based on exergy generation and exergy dissipation terms. This novel nonlinear control methodology results in both necessary and sufficient conditions for stability of nonlinear systems. © 2006 IEEE.

The dimensionless extinction coefficient, Ke, was measured for soot produced in 2m JP-8 pool fires. Light extinction and gravimetric sampling measurements were performed simultaneously at 635 and 1310nm wavelengths at three heights in the flame zone and in the overfire region. Measured average Ke values of 8.41.2 at 635nm and 8.71.1 at 1310nm in the overfire region agree well with values from 8-10 recently reported for different fuels and flame conditions. The overfire Ke values are also relatively independent of wavelength, in agreement with recent findings for JP-8 soot in smaller flames. Ke was nearly constant at 635nm for all sampling locations in the large fires. However, at 1310nm, the overfire Ke was higher than in the flame zone. Chemical analysis of physically sampled soot shows variations in carbon-to-hydrogen (C/H) ratio and polycyclic aromatic hydrocarbon (PAH) concentration that may account for the smaller Ke values measured in the flame zone. Rayleigh-Debye-Gans theory of scattering for polydisperse fractal aggregate (RDG-PFA) was applied to measured aggregate fractal dimensions and found to under-predict the extinction coefficient by 17-30% at 635nm using commonly accepted refractive indices of soot, and agreed well with the experiments using the more recently published refractive index of 1.99-0.89i. This study represents the first measurements of soot chemistry, morphology, and optical properties in the flame zone of large, fully-turbulent pool fires, and emphasizes the importance of accurate measurements of optical properties both in the flame zone and overfire regions for models of radiative transport and interpretation of laser-based diagnostics of soot volume fraction and temperature.

Cities without an early warning system of indwelling sensors can consider monitoring their networks manually, especially during times of heightened security levels. We consider the problem of calculating an optimal schedule for manual sampling in a municipal water network. Preliminary computations with a small-scale example indicate that during normal times, manual sampling can provide some benefit, but it is far inferior to an indwelling sensor network. However, given information that significantly constrains the nature of an imminent threat, manual sampling can perform as well as a small sensor network designed to handle normal threats. Copyright ASCE 2006.

We present a theory of impurity states in quantum wells, where the confining potential of the heterostructure and the random impurity potential are treated in a unified theory. After diagonalization of the 3D Hamiltonian we calculate the infrared absorption spectrum. We discuss the nature of impurity states that are confined in the quantum wells and their influence on the absorption spectra. We then calculate the absorption spectra for a quadruple quantum well, revealing impurity transitions as well as intersubband transitions. The results are compared to existing experimental data and show a remarkable agreement. © 2007 American Institute of Physics.

Accurate material models are fundamental to predictive structural finite element models. Because potting foams are routinely used to mitigate shock and vibration of encapsulated components in electro/mechanical systems, accurate material models of foams are needed. A linear-viscoelastic foam constitutive model has been developed to represent the foam's stiffness and damping throughout an application space defined by temperature, strain rate or frequency and strain level. Validation of this linear-viscoelastic model, which is integrated into the Salinas structural dynamics code, is being achieved by modeling and testing a series of structural geometries of increasing complexity that have been designed to ensure sensitivity to material parameters. Both experimental and analytical uncertainties are being quantified to ensure the fair assessment of model validity. Quantitative model validation metrics are being developed to provide a means of comparison for analytical model predictions to observations made in the experiments. This paper is one of several recent papers documenting the validation process for simple to complex structures with foam encapsulated components. This paper specifically focuses on model validation over a wide temperature range and using a simple dumbbell structure for modal testing and simulation. Material variations of density and modulus have been included. A double blind validation process is described that brings together test data with model predictions.

In this paper we present the results of a study to quantify uncertainty in experimental modal parameters due to test set-up uncertainty, measurement uncertainty, and data analysis uncertainty. Uncertainty quantification is required to accomplish a number of tasks including model updating, model validation, and assessment of unit-tounit variation. We consider uncertainty in the modal parameters due to a number of sources including force input location/direction, force amplitude, instrumentation bias, support conditions, and the analysis method (algorithmic variation). We compute the total uncertainty due to all of these sources, and discuss the importance of proper characterization of bias errors on the total uncertainty. This uncertainty quantification was applied to modal tests designed to assess modeling capabilities for emerging designs of wind turbine blades. In an example, we show that unit-to-unit variation of the modal parameters of two nominally identical wind turbine blades is successfully assessed by performing uncertainty quantification. This study aims to demonstrate the importance of the proper pre-test design and analysis for understanding the uncertainty in modal parameters, in particular uncertainty due to bias error.

In order to predict blast damage on structures, it is current industry practice to decouple shock calculations from computational structural dynamics calculations. Pressure-time histories from experimental tests were used to assess computational models developed using a shock physics code (CTH) and a structural dynamics code (PRONTO3D). CTH was shown to be able to reproduce three independent characteristics of a blast wave: arrival time, peak overpressure, and decay time. Excellent agreement was achieved for early times, where the rigid wall assumptions used in the model analysis were valid. A one-way coupling was performed for this blast-structure interaction problem by taking the pressure-time history from the shock physics simulation and applying it to the structure at the corresponding locations in the PRONTO3D simulation to capture the structural deformation. In general, the one-way coupling was shown to be a cost-effective means of predicting the structural response when the time duration of the load was less than the response time of the structure. Therefore, the computational models were successfully evaluated for the internal blast problems studied herein.

The water dynamics in a series of Sandia octahedral molecular sieves (SOMS) were investigated using high speed 1H magic angle spinning (MAS) NMR spectroscopy. For these materials both the 20% Ti-substituted material, Na 2Nb1.6Ti0.4(OH)0.4O 5.6·H2O and the 0% exchanged end member, Na 2Nb2O6·H2O were studied. By combining direct one dimensional (1D) MAS NMR experiments with double quantum (DQ) filtered MAS NMR experiments different water environments within the materials were identified based on differences in mobility. Two dimensional (2D) DQ correlation experiments were used to extract the DQ spinning sideband patterns allowing the residual 1H-1H homonuclear dipolar coupling to be measured. From these DQ experiments the effective order parameters for the different water environments were calculated. The water environments in the two different SOMS compositions investigated revealed very large differences in the water mobility. © 2007 Materials Research Society.

Both the NASA and DOE have programs that are investigating advanced power conversion cycles for planetary surface power on the moon or Mars, and for next generation nuclear power plants on earth. The gas Brayton cycle offers many practical solutions for space nuclear power systems and was selected as the nuclear power system of choice for the NASA Prometheus project. An alternative Brayton cycle that offers high efficiency at a lower reactor coolant outlet temperature is the supercritical Brayton cycle (SCBC). The supercritical cycle is a true Brayton cycle because it uses a single phase fluid with a compressor inlet temperature that is just above the critical point of the fluid. This paper describes the use of a supercritical Brayton cycle that achieves a cycle efficiency of 26.6% with a peak coolant temperature of 750 K and for a compressor inlet temperature of 390 K. The working fluid uses a clear odorless, nontoxic refrigerant C318 perflurocarbon (C4F8) that always operates in the gas phase. This coolant was selected because it has a critical temperature and pressure of 388.38 K and 2.777 MPa. The relatively high critical temperature allows for efficient thermal radiation that keeps the radiator mass small. The SCBC achieves high efficiency because the loop design takes advantage of the non-ideal nature of the coolant equation of state just above the critical point. The lower coolant temperature means that metal fuels, uranium oxide fuels, and uranium zirconium hydride fuels with stainless steel, ferretic steel, or superalloy cladding can be used with little mass penalty or reduction in cycle efficiency. The reactor can use liquid-metal coolants and no high temperature heat exchangers need to be developed. Indirect gas cooling or perhaps even direct gas cooling can be used if the C4F8 coolant is found to be sufficiently radiation tolerant. Other fluids can also be used in the supercritical Brayton cycle including Propane (C3H 8, Tcritical = 369 K) and Hexane (C6H 14, Tcritical = 506.1 K) provided they have adequate chemical compatibility and stability. Overall the use of supercritical Brayton cycles may offer "break through" operating capabilities for space nuclear power plants because high efficiencies can be achieved a very low reactor operating temperatures which in turn allows for the use of available fuels, cladding, and structural materials. © 2007 American Institute of Physics.

This paper investigates methods for coupling analytical dynamic models of subcomponents with experimentally derived models in order to predict the response of the combined system, focusing on modal substructuring or Component Mode Synthesis (CMS), the experimental analog to the ubiquitous Craig-Bampton method. While the basic methods for combining experimental and analytical models have been around for many years, it appears that these are not often applied successfully. The CMS theory is presented along with a new strategy, dubbed the Maximum Rank Coordinate Choice (MRCC), that ensures that the constrained degrees of freedom can be found from the unconstrained without encountering numerical ill conditioning. The experimental modal substructuring approach is also compared with frequency response function coupling, sometimes called admittance or impedance coupling. These methods are used both to analytically remove models of a test fixture (required to include rotational degrees of freedom) and to predict the response of the coupled beams. Both rigid and elastic models for the fixture are considered. Similar results are obtained using either method although the modal substructuring method yields a more compact database and allows one to more easily interrogate the resulting system model to assure that physically meaningful results have been obtained. A method for coupling the fixture model to experimental measurements, dubbed the Modal Constraint for Fixture and Subsystem (MCFS) is presented that greatly improves the result and robustness when an elastic fixture model is used.

The detection of anomalous water quality events has become an increased priority for distribution systems, both for quality of service and security reasons. Because of the high cost associated with false detections, both missed events and false alarms, algorithms which aim to provide event detection aid need to be evaluated and configured properly. CANARY has been developed to provide both real-time, and off-line analysis tools to aid in the development of these algorithms, allowing algorithm developers to focus on the algorithms themselves, rather than on how to read in data and drive the algorithms. Among the features to be discussed and demonstrated are: 1) use of a standard data exchange format for input and output of water quality and operations data streams; 2) the ability to "plug in" various water quality change detection algorithms, both in MATLAB® and compiled library formats for testing and evaluation by using a well defined interface; 3) an "operations mode" to simulate what a utility operator will receive; 4) side-by-side comparison tools for different evaluation metrics, including ROC curves, time to detect, and false alarm rates. Results will be shown using three algorithms previously developed (Klise and McKenna, 2006; McKenna, et al., 2006) using test and real-life data sets. © 2007 ASCE.

Like other interfaces, equilibrium grain boundaries are smooth at low temperature and rough at high temperature; however, little attention has been paid to roughening except for faceting boundaries. Using molecular dynamics simulations of face-centered cubic Ni, we studied two closely related grain boundaries with different boundary planes. In spite of their similarity, their boundary roughening temperatures differ by several hundred degrees, and boundary mobility is much larger above the roughening temperature. This has important implications for microstructural development during metallurgical processes.

Abstract not provided.

In order for electric power generating capacity to be supplanted to a meaningful extent by sources smaller than 200 kW, an automated means of managing systems of small sources must be found or their sheer number will overwhelm the power production community. Microgrids - power systems comprising multiple small interconnected generators - are a promising response to this need, but an automated microgrid management system has not been demonstrated. This paper describes the energy management task and its execution in a standardized grid services environment. ©2007 IEEE.

The evolution of Guinier-Preston zones in an Al-2.7 at.% Ag alloy was studied using atom probe tomography. The composition and morphology of the GP zones are time dependent, explaining discrepancies in previous work. This result requires the metastable irascibility gap for GP zones to be reevaluated, highlighting the importance of the temporal evolution of the GP zones. Preliminary results on the composition of γ′ and γ plates are also presented. © MICROSCOPY SOCIETY OF AMERICA 2007.

We have discussed the key areas of the IR process that should not be circumvented if an organization is to achieve a high level of assurance in high-dollar, high-risk cost estimates; lessons learned; and possible solutions to improve the process. In summary, the best practices described are to do the following. Develop a corporate policy for review of cost estimates based on TPC and potential financial and reputation risk; Develop a database of qualified, experienced personnel, who can perform well as IR team members; Spell out the process for approval of review team members, including the executive approval process; Address review team availability by developing review team member alternates; Increase lead-time notice on high-dollar, high risk estimates by developing an advanced notice system with internal organizations; Improve coordination of the estimating team's responses to the review team's questions and concerns; and Develop alternatives such as representatives and electronic briefings to alleviate challenges in scheduling executives for cost estimate briefings. Each organization has its own needs, culture, and level of maturity. If you have an IR process that works, great! If not, we hope that we have sparked your interest in developing a process that works for your company. The goal is to continuously improve and further refine the process to meet the needs of both external and internal customers. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company for the US Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Accurate material models are fundamental to predictive structural finite element models. Because potting foams are routinely used to mitigate shock and vibration of encapsulated components in electro/mechanical systems, accurate material models of foams are needed. A linear-viscoelastic foam constitutive model has been developed to represent the foam's stiffness and damping throughout an application space defined by temperature, strain rate or frequency and strain level. Validation of this linear-viscoelastic model, which is integrated into the Salinas structural dynamics code, is being achieved by modeling and testing a series of structural geometries of increasing complexity that have been designed to ensure sensitivity to material parameters. Both experimental and analytical uncertainties are being quantified to ensure the fair assessment of model validity. Quantitative model validation metrics are being developed to provide a means of comparison for analytical model predictions to observations made in the experiments. This paper is one of several recent papers documenting the validation process for simple to complex structures with foam encapsulated components. This paper specifically focuses on model validation over a wide temperature range and using a simple dumbbell structure for modal testing and simulation. Material variations of density and modulus have been included. A double blind validation process is described that brings together test data with model predictions.

Even within a well-established systems engineering organization, formalizing the practice of systems engineering can be an arduous task. This paper describes one organization's effort to start this formalization by defining and documenting the very foundation of its practice: its systems engineering principles. Topics include the process for developing principles (augmenting organizational guidance with industry best practices), the final form of the principles (justifying terminology), and the connection to broader formalization efforts (associating the principles diagram to the systems engineering logo). The principles are divided into four categories that tie mainstream systems engineering definitions together: Stakeholder, Systems Engineer, Problem, and Solution. © 2007 by Sandia National Laboratories.

The Z machine at Sandia National Laboratories drives 20 MA in 100 ns through a cylindrical array of fine wires which implodes due to the strong j × B force, generating up to 250 TW of soft x-ray radiation when the z-pinch plasma stagnates on axis. The copious broadband self-emission makes the dynamics of the implosion well suited to diagnosis with soft x-ray imaging and spectroscopy. A monochromatic self-emission imaging instrument has recently been developed on Z which reflects pinhole images from a multilayer mirror onto a 1 ns gated microchannel plate detector. The multilayer can be designed to provide narrowband (∼10 eV) reflection in the 100-700 eV photon energy range, allowing observation of the soft emission from accreting mass as it assembles into a hot, dense plasma column on the array axis. In the present instrument configuration, data at 277 eV photon energy have been obtained for plasmas ranging from Al to W, and the z-pinch implosion and stagnation will be discussed along with > 1 keV self-emission imaging and spectroscopy. Collisional-radiative simulations are currently being pursued in order to link the imaged emissivity to plasma temperature and density profiles and address the role of opacity in interpreting the data. © 2007 American Institute of Physics.

For most orbital maneuvers, small satellites in the sub-10 kg range require thrusters capable of spanning the micro-Newton to milli-Newton force range. At this scale, electrokinetic (EK) pumping offers precise metering of monergolic or hypergolic liquid propellants under purely electrical control at pressures and flow rates well-suited to microthruster applications. We have demonstrated direct and indirect EK pumping for delivery of anhydrous hydrazine and hydrogen peroxide monopropellants, respectively, into capillary-based microthrusters with integrated in-line catalyst beds. Catalytic decomposition generates gases which accelerate through a plasma-formed converging-diverging nozzle, producing thrust. Specific impulses up to 190 s have been shown for hydrazine in non-optimized nozzles. ©2007 IEEE.

Real time protein binding interactions between biotinylated bovine serum albumen (BSA) and streptavidin (SA) are detected using chemoreceptive neuron MOS (CvMOS) transistors with extended floating gate structures. This enables protein interaction events to be monitored by a sensing gate area that is only capacitively coupled to the sensing circuits and far removed from the active area, reducing the invasiveness without loss of sensitivity. The use of both sensing and control gates to control the floating gate potential eliminates the need of an analyte reference electrode. © 2007 IEEE.

Abstract not provided.

Abstract not provided.

The U.S. Department of Energy (DOE) has a rich history of significant contributions to geospatial science spanning the past four decades. In the early years, work focused on basic research, such as development of algorithms for processing geographic data and early use of LANDSAT imagery. The emphasis shifted in the mid-1970s to development of geographic information system (GIS) applications to support programs such as the National Uranium Resource Evaluation (NURE), and later to issue-oriented GIS applications supporting programs such as environmental restoration and management (mid-1980s through present). Throughout this period, the DOE national laboratories represented a strong chorus of voices advocating the importance of geospatial science and technology in the decades to come. The establishment of a Geospatial Science Program by the DOE Office of the Chief Information Officer in 2005 reflects the continued potential of geospatial science to enhance DOE's science, projects, and operations, as is well demonstrated by historical analysis. © 2007 The Haworth Press, Inc. All rights reserved.

When a system design approach is applied to wind turbine blades, manufacturing and structural requirements are included along with aerodynamic considerations in the design optimization. The resulting system-driven design includes several innovative structural features such as flat-back airfoils, a constant thickness carbon spar-cap, and a thin, large diameter root. Subscale blades were manufactured to evaluate the as-built integrated performance. The design resulted in a 22% reduction in mass, but withstood over 300% of its design load during testing. Compressive strains of nearly 0.9% were measured in the carbon spar-cap. The test results from this and an earlier design are compared, as are finite element models of each design. Included in the analysis is a review of the acoustic emission events that were detected through the use of surface mounted microphones.

The US Strategic Petroleum Reserve (SPR) stores oil in large underground salt caverns. This oil has compositional and thermal gradients induced by geothermal heating from both the bottom surface and side walls. Temperature layering has been recorded in SPR oil caverns, which is hypothesized to be predominantly due to double-diffusive layering that occurs when a stable compositional gradient is heated from below. Initial results of a laboratory experimental program aimed at studying dynamics of such double-diffusive layers in the context of SPR are described in this paper. Of particular interest are the thickness of converting layers, layer evolution (migration/merging) and conditions for the formation/non-formation of double-diffusive layers. Copyright © 2007 by ASME.

When measuring the structural dynamic response of test objects, the desired data is sometimes combined with some type of undesired periodic data. This can occur due to N-per-revolution excitation in systems with rotating components or when dither excitation is used. The response due to these (typically unmeasured) periodic excitations causes spikes in system frequency response functions (FRFs) and poor coherence. This paper describes a technique to remove these periodic components from the measured data. The data must be measured as a continuous time history which is initially processed as a single, long record. Given an initial guess for the periodic signal's fundamental frequency, an automated search will identify the actual fundamental frequency to very high accuracy. Then the fundamental and a user-specified number of harmonics are removed from the acquired data to create new time histories. These resulting time histories can then be processed using standard signal processing techniques. An example of this technique will be presented from a test where a vehicle is dithered with a fixed-frequency, sinusoidal force to linearize the behavior of the shock absorbers, while measuring the acceleration responses due to a random force applied elsewhere on the vehicle.

We present an LES-type variational multiscale theory of turbulence. Our approach derives completely from the incompressible Navier-Stokes equations and does not employ any ad hoc devices, such as eddy viscosities. We tested the formulation on forced homogeneous isotropic turbulence and turbulent channel flows. In the calculations, we employed linear, quadratic and cubic NURBS. A dispersion analysis of simple model problems revealed NURBS elements to be superior to classical finite elements in approximating advective and diffusive processes, which play a significant role in turbulence computations. The numerical results are very good and confirm the viability of the theoretical framework. © 2007 Elsevier B.V. All rights reserved.

Partitioned global address space (PGAS) programming models have been identified as one of the few viable approaches for dealing with emerging many-core systems. These models tend to generate many small messages, which requires specific support from the network interface hardware to enable efficient execution. In the past, Cray included E-registers on the Cray T3E to support the SHMEM API; however, with the advent of multi-core processors, the balance of computation to communication capabilities has shifted toward computation. This paper explores the message rates that are achievable with multi-core processors and simplified PGAS support on a more conventional network interface. For message rate tests, we find that simple network interface hardware is more than sufficient. We also find that even typical data distributions, such as cyclic or block-cyclic, do not need specialized hardware support. Finally, we assess the impact of such support on the well known RandomAccess benchmark. (c) 2007 ACM.

Implicit large eddy simulation (ILES) has provided many computer simulations with an efficient and effective model for turbulence. The capacity for ILES has been shown to arise from a broad class of numerical methods with specific properties producing nonoscillatory solutions using limiters that provide these methods with nonlinear stability. The use of modified equation has allowed us to understand the mechanisms behind the efficacy of ILES as a model. Much of the understanding of the ILES modeling has proceeded in the realm of incompressible flows. Here, we extend this analysis to compressible flows. While the general conclusions are consistent with our previous findings, the compressible case has several important distinctions. Like the incompressible analysis, the ILES of compressible flow is dominated by an effective self-similarity model (Bardina, J., Ferziger, J. H., and Reynolds, W. C., 1980, quot;Improved Subgrid Scale Models for Large Eddy Simulations, quot; AIAA Paper No. 80-1357; Borue, V., and Orszag, S. A., 1998, quot;Local Energy Flux and Subgrid-Scale Statistics in Three Dimensional Turbulence,quot; J. Fluid Mech., 366, pp. 1-31; Meneveau, C., and Katz, J., 2000, quot;Scale-Invariance and Turbulence Models for Large-Eddy Simulations,quot; Annu. Rev. Fluid. Mech., 32, pp. 1-32). Here, we focus on one of these issues, the form of the effective subgrid model for the conservation of mass equations. In the mass equation, the leading order model is a self-similarity model acting on the joint gradients of density and velocity. The dissipative ILES model results from the limiter and upwind differencing resulting in effects proportional to the acoustic modes in the flow as well as the convective effects. We examine the model in several limits including the incompressible limit. This equation differs from the standard form found in the classical Navier-Stokes equations, but generally follows the form suggested by Brenner (2005, quot;Navier-Stokes Revisited,quot; Physica A, 349(1-2), pp. 60-133) in a modification of Navier-Stokes necessary to successfully reproduce some experimentally measured phenomena. The implications of these developments are discussed in relation to the usual turbulence modeling approaches. Copyright ©2007 by ASME.

In this paper, we present an optimal method for calculating turning maneuvers for an unmanned aerial vehicle (UAV) developed for ecological research. The algorithm calculates several possible solutions using vectors represented in complex notation, and selects the shortest turning path given constraints determined by the aircraft. This algorithm considers the UAV's turning capabilities, generating a two-dimensional path that is feasible for the UAV to fly. We generate a test flight path and show that the UAV is capable of following the turn maneuvers.

The goal of this study is to model the electrical response of gold plated copper electrical contacts exposed to a mixed flowing gas stream consisting of air containing 10 ppb H2S at 30 °C and a relative humidity of 70%. This environment accelerates the attack normally observed in a light industrial environment (essentially a simplified version of the Battelle class 2 environment). Corrosion rates were quantified by measuring the corrosion site density, size distribution, and the macroscopic electrical resistance of the aged surface as a function of exposure time. A pore corrosion numerical model was used to predict both the growth of copper sulfide corrosion product which blooms through defects in the gold layer and the resulting electrical contact resistance of the aged surface. Assumptions about the distribution of defects in the noble metal plating and the mechanism for how corrosion blooms affect electrical contact resistance were needed to close the numerical model. Comparisons are made to the experimentally observed number density of corrosion sites, the size distribution of corrosion product blooms, and the cumulative probability distribution of the electrical contact resistance. Experimentally, the bloom site density increases as a function of time, whereas the bloom size distribution remains relatively independent of time. These two effects are included in the numerical model by adding a corrosion initiation probability proportional to the surface area along with a probability for bloom-growth extinction proportional to the corrosion product bloom volume. The cumulative probability distribution of electrical resistance becomes skewed as exposure time increases. While the electrical contact resistance increases as a function of time for a fraction of the bloom population, the median value remains relatively unchanged. In order to model this behavior, the resistance calculated for large blooms has been weighted more heavily. © 2007 IEEE.

A sub-scale experiment has been constructed using fins mounted on one wall of a transonic wind tunnel to investigate the influence of fin tip vortices upon downstream control surfaces. Data are collected using a fin balance mounted on the downstream fin to measure the aerodynamic forces of the interaction, combined with stereoscopic Particle Image Velocimetry to measure vortex properties. The fin balance data show that the response of the downstream fin essentially is shifted from the baseline single-fin data dependent upon the angle of attack of the upstream fin. Freestream Mach number and the spacing between fins have secondary effects. The velocimetry shows that the vortex strength increases markedly with upstream fin angle of attack, though even an uncanted fin generates a noticeable wake. No Mach number effect can be discerned in the normalized data, but measurements taken progressively further from the fin trailing edge show the decay in vortex strength with downstream distance. Correlations between the force data and the velocimetry suggest that the interaction is fundamentally a result of an angle of attack induced upon the downstream fin by the vortex shed from the upstream fin tip. © 2007 IEEE.

Parallel adaptive mesh refinement methods potentially lead to realistic modeling of complex three-dimensional physical phenomena. However, they also present significant challenges in data partitioning and load balancing. As the mesh adapts to the solution, the partitioning requirements change. By explicitly considering these dynamic conditions, the scalability for large, realistic simulations could possibly be significantly improved. Our hypothesis is that adaptive partitioning, meaning dynamic and automatic switching of partitioning techniques, based on the current run-time state, can be beneficial for these simulations. However, switching partitioners can be expensive due to differences in the algorithms' native mapping of data onto processors. We suggest forcing a uniform starting point for all included partitioners. We present a penalty-based method for determining whether switching is beneficial. We study the effects on data migration, as well as on overall cost, of using the uniform starting point and the switching-penalties to select the best partitioning algorithm, among a set of graph-based and geometric partitioning algorithms, for each adaptive time-step for four different adaptive scientific applications. The results show that data migration can be significantly reduced and that adaptive partitioning indeed can be effective for unstructured adaptive applications.

This paper describes the electromagnetic analysis that has been completed using the OPERA-3d product to characterize the folces on the ITER shield modules as part of the conceptual design. These forces exist due to the interaction of the eddy currents induced in the shield modules and the large magnetic fields present in the tokamak. ©2007 IEEE.

Isentropic compression of a solid to 100's of GPa by a ramped, planar compression wave allows measurement of material properties at high strain and at modest temperature. Introduction of a measurement plane disturbs the flow, requiring special analysis techniques. If the measurement interface is windowed, the unsteady nature of the wave in the window requires special treatment. When the flow is hyperbolic the equations of motion can be integrated backward in space in the sample to a region undisturbed by the interface interactions, fully accounting for the untoward interactions. For more complex materials like hysteretic elastic/plastic solids or phase changing material, hybrid analysis techniques are required. © 2007 American Institute of Physics.

At Sandia National Laboratories, we have coined the term "microenergetics" to describe sub-millimeter energetic material studies aimed at gaining knowledge of combustion and detonation behavior at the mesoscale. Our approach is to apply technologies developed by the microelectronics industry to fabricate test samples with well-defined geometries. Substrates have been fabricated from materials such as silicon and ceramics, with channels to contain the energetic material. Energetic materials have been loaded into the channels, either as powders, femtosecond laser-micromachined pellets, or as vapor-deposited films. Ignition of the samples has been achieved by simple hotwires, integrated semiconductor bridges, and also by lasers. Additionally, grain-scale patterning has been performed on explosive films using both oxygen plasma etching and femtosecond laser micromachining. We have demonstrated simple work functions in microenergetic devices, such as piston motion, which is also a relevant diagnostic to examine combustion properties. Detonation has been achieved in deposited explosive films, recorded by high-speed photography. © 2007 American Institute of Physics.

Requirements Engineering viewpoints are converging more between the Software Engineering Community and the Systems Engineering Community, but there remains driving viewpoints that linger and constrain our complete understanding of a problem and of requirements for a solution. These driving ideas need to be reassessed in order to represent and assess the problem or opportunity as well as the requirements for systems that are introduced. In particular, in order to perform the modeling and acquire understanding of an Enterprise it is crucial that we reevaluate these driving viewpoints. © 2007 by R.M. Griego.

An accurate method for controlling strain rates in dynamic compressions studies involves using the non-linear elastic property of fused silica to transform an initial shock into a ramp wave of known amplitude and duration. Fused silica when placed between a dry Indiana limestone specimen and a projectile produces strain rates in the range of 104/s. Ramp-loading strain rates are higher than what can be produced on Hopkinson bars and lower than what shock experiments attain. The strength determined at the elastic limit under ramp loading compared to Hopkinson bar measurements shows a significant strength increase with increasing strain rate. © 2007 American Institute of Physics.

Optically powered devices are typically irradiated by high intensity lasers and rely on the temperature excursion generated by the laser for operation. While numerical modeling can estimate the temperature profile of the irradiated devices, only direct measurements can determine the actual device temperatures. Available surface thermometry techniques, such as infrared imaging, scanning thermal microscopy and thermoreflectance are generally incompatible with an optical powering scheme, the micron-scale layer thicknesses of microsystem devices, or both. In this paper we discuss the use of micro-Raman thermometry to obtain the first spatiallyresolved temperature measurements of various polycrystalline silicon (polysilicon) surfaces heated with an 808 nm continuous wave (CW) laser at a 60° angle of incidence. The micron-scale resolution of the micro-Raman technique permitted mapping of the surface temperature in the vicinity of the heating laser spot and throughout the device. In addition to discussing the requirements for accurate data collection, the implications of optical interference on the heated structures are also considered. Copyright © 2007 by ASME.

Even within a well-established systems engineering organization, formalizing the practice of systems engineering can be an arduous task. This paper describes one organization's effort to start this formalization by defining and documenting the very foundation of its practice: its systems engineering principles. Topics include the process for developing principles (augmenting organizational guidance with industry best practices), the final form of the principles (justifying terminology), and the connection to broader formalization efforts (associating the principles diagram to the systems engineering logo). The principles are divided into four categories that tie mainstream systems engineering definitions together: Stakeholder, Systems Engineer, Problem, and Solution. © 2007 by Sandia National Laboratories.

Potential events involving biological or chemical contamination of buildings are of major concern in the area of homeland security. Tools are needed to provide rapid, onsite predictions of contaminant levels given only approximate measurements in limited locations throughout a building. In principal, such tools could use calculations based on physical process models to provide accurate predictions. In practice, however, physical process models are too complex and computationally costly to be used in a real-time scenario. In this paper, we investigate the feasibility of using machine learning to provide easily computed but approximate models that would be applicable in the field. We develop a machine learning method based on Support Vector Machine regression and classification. We apply our method to problems of estimating contamination levels and contaminant source location. © 2007 IEEE.

In this paper, we introduce EXACT, the EXperimental Algorithmics Computational Toolkit. EXACT is a software framework for describing, controlling, and analyzing computer experiments. It provides the experimentalist with convenient software tools to ease and organize the entire experimental process, including the description of factors and levels, the design of experiments, the control of experimental runs, the archiving of results, and analysis of results. As a case study for EXACT, we describe its interaction with FAST, the Sandia Framework for Agile Software Testing. EXACT and FAST now manage the nightly testing of several large software projects at Sandia. We also discuss EXACT's advanced features, which include a driver module that controls complex experiments such as comparisons of parallel algorithms. Copyright 2007 ACM.

Recently, extensive focus has been placed on determining the optimal locations of sensors within a distribution system to minimize the impact on public health from intentional intrusion events. Modified versions of these tools may have additional benefits for determining monitoring locations for other more common objectives associated with distribution systems. A modified Sensor Placement Optimization Tool (SPOT) is presented that can be used for satisfying more generic location problems such as determining monitoring locations for tracer tests or disinfectant byproduct sampling. The utility for the modified SPOT algorithm is discussed with respect to implementing a distribution system field-scale tracer study. © 2007 ASCE.

Since the first vector supercomputers in the mid-1970's, the largest scale applications have traditionally been floating point oriented numerical codes, which can be broadly characterized as the simulation of physics on a computer. Supercomputer architectures have evolved to meet the needs of those applications. Specifically, the computational work of the application tends to be floating point oriented, and the decomposition of the problem two or three dimensional. Today, an emerging class of critical applications may change those assumptions: they are combinatorial in nature, integer oriented, and irregular. The performance of both classes of applications is dominated by the performance of the memory system. This paper compares the memory performance sensitivity of both traditional and emerging HPC applications, and shows that the new codes are significantly more sensitive to memory latency and bandwidth than their traditional counterparts. Additionally, these codes exhibit lower base-line performance, which only exacerbates the problem. As a result, the construction of future supercomputer architectures to support these applications will most likely be different from those used to support traditional codes. Quantitatively understanding the difference between the two workloads will form the basis for future design choices. ©2007 IEEE.

An S-band image reject downconversion mixer with high intercept point and fully integrated single-ended ports, including a UHF output, is demonstrated using GaAs pHEMT technology. The image rejection is better than 20 dB across a wide IF bandwidth ranging from 400 MHz to beyond 1 GHz with a high input IP3 of 20 dBm. On-chip passive baluns are used to provide the single-ended to differential conversion in the RF and IF paths necessary for the resistive FET mixer core. A polyphase filter is used to generate the quadrature local oscillator (LO) components, while an integrated UHF lumped element quadrature hybrid combines the intermediate frequency (IF) components to achieve image rejection fully on-chip without the need for external components. © 2007 IEEE.

Abstract not provided.

Analysts, experimenters, and facilities have fallen into some poor practices in reporting many dosimetry metrics. While the experienced dosimetrist often knows the caveats that apply for a given dosimetry application, without proper reporting critical information is often lost before the data is received by the dosimetrist. In addition, the newcomers to the application of dosimetry are not being educated in the importance of a variation in the irradiation conditions. This paper captures some of the cases where care must be taken in expressing the proper context for a dosimetry metric. Examples focus on the interpretation of the response of a diamond photoconducting detector and a silicon transistor and highlight some common mistakes and some not-so-clear misinterpretations that even the experienced person often makes in this field. A careful study of the underlying physics reveals the non-intuitive trends in some metrics. Suggestions are made on how the community can minimize the chance of a dosimetry-related misinterpretation. © 2007 IEEE.

Si single junction photocells manufactured by Sandia National Laboratories and commercially available multiple junction GaAs photocells were tested in a pulsed high-dose mixed gamma-neutron environment. The Si photocells were also tested in steady state gamma environment at two different dose rates.

Lanthanide halide alloys have recently enabled scintillating gamma ray spectrometers comparable to room-temperature semiconductors (< 3% FWHM energy resolutions at 662keV). However brittle fracture of these materials hinders the growth of large volume crystals. Efforts to improve the strength through non-lanthanide alloy substitution, while preserving scintillation, are being pursued. Isovalent alloys nominal Ce0.9Al0.1Br 3, Ce0.9Ga0.1Br3, Ce 0.9Sc0.1Br3, Ce0.9In 0.1Br3 and Ce0.8Y0.2Br3, as well as aliovalent alloys nominal (CeBr3)0.99(CdCl 2)0.01, (CeBr3)0.99(CdBr 2)0.01, (CeBr3)0.99(ZnBr 2)0.01, (CeBr3)0.99(CaBr 2)0.01, (CeBr3)0.99(SrBr 2)0.01, (CeBr3)0.99(PbBr 2)0.01, (CeBr3)0.99(ZrBr 4)0.01, (CeBr3)0.99(HfBr 4)0.01 were prepared. All of these alloys exhibit bright fluorescence under UV excitation, with varying shifts in the spectral peaks and intensities relative to pure CeBr3. Further, these alloys scintillate when coupled to a photomultiplier tube (PMT) and exposed to 137Cs gamma rays. These data and the potential for improved crystal growth will be discussed.

The extended Finite Element Method (X-FEM) is a finite-element based discretization technique developed originally to model dynamic crack propagation [1]. Since that time the method has been used for modeling physics ranging from static meso-scale material failure to dendrite growth. Here we adapt the recent advances of Vitali and Benson [2] and Song et. al. [3] to model dynamic loading of a polycry stalline material. We use demonstration problems to examine the method's efficacy for modeling the dynamic response of polycrystalline materials at the meso-scale. Specifically, we use the X-FEM to model grain boundaries. This approach allows us to i) eliminate ad-hoc mixture rules for multi-material elements and ii) avoid explicitly meshing grain boundaries. © 2007 American Institute of Physics.

An analytical model is presented for predicting the critical air-flow velocity at the onset of water-droplet detachment from the GDL/channel interfaces in PEM fuel cells. Our model is based on the force balance between pressure drag that tends to detach the droplet and surface tension that tends to hold the droplet in place. In the present work, we consider the flow regime in which pressure drag, which arises from inertia effects, dominates over viscous shear - this is the flow regime of interest in real-world PEM fuel cell applications, both automotive and stationary. Our analytical model predicts that the critical air-flow velocity varies inversely (to the 2/3 power) with water-droplet size. It further predicts that making the GDL surface more hydrophobic, decreasing contact-angle hysteresis, and shrinking channel height reduce the critical air-flow velocity. Model predictions are compared with experimental data available from the literature and reasonably good agreement is obtained. © The Electrochemical Society.

Today's world demands new ways of thinking about security solutions. The problem space is complex and ambiguous. Solutions must be multidimensional, incorporating not only technology, but the social, economic, political, and religious dynamics of a security intervention. A facilitator-led experiential training program was designed for our technical staff that leads them out of the box. The course design is based upon the theories of cognitive flexibility and situated cognition, and uses a socio-constructivist approach. Participants are led by a senior systems engineer/facilitator through a series of exercises in which they observe contextually relevant right way/wrong way videos, engage in critical thinking assessments about what they observed, and solve logic puzzles. Group interaction and problem-solving is emphasized. As in the real world, there is no one "right" solution. Outcomes can include a broader understanding of the threat space, creative solutions that enable survival in spite of an evolving enemy, and a deeper sense of the complex dynamics involved in any security decision. Training impact is being evaluated using a mixed qualitative/quantitative approach. Survey data combined with ethnographic interviewing techniques will determine whether or not participants have transferred their new understandings to the work environment. ©2007 IEEE.

Soda-lime glass (SLG) is a highly available low cost glass formulation commonly used in window applications. Over the past decade, there have been a number of studies which have examined the Hugoniot elastic limit (HEL) of this material resulting in a wide range of values from 3.1 to 6.0 GPa. The determination of the HEL is complicated by many factors including ramp loading due to the convex downward curvature of the Hugoniot at low pressures. Results of transmitted wave profile experiments up to 20 GPa are presented and analyzed to determine the loading and release characteristics of SLG near the HEL. Results indicate a response that is more complex than the elastic - plastic response typical of many materials, possibly explaining the wide range in initially reported HEL values. © 2007 American Institute of Physics.

Attenuating wave profiles from shock experiments on tungsten carbide powder are compared to calculations from the continuum P-λ model and a 2-D mesoscale model to gain insight into the suitability of the two models. When calibrated, both models accurately capture the Hugoniot response of the powder and the arrival times of unattenuated steady waves. Their amplitudes are more accurately given by the mesoscale model since its reshock states are above the Hugoniot as seen experimentally; the P-λ model, in contrast, reshocks along the Hugoniot. When the attenuating wave is in the range of the Hugoniot data, the models predict attenuation correctly. However, when attenuation falls below the Hugoniot data both models are somewhat inaccurate, and the material response seems to lie between the two models. The final aspect considered is the wave rise time, which is qualitatively correct for the mesoscale model but completely inaccurate for the P-λ model. © 2007 American Institute of Physics.

Material heterogeneity appears to give rise to variability in the yield behavior of ceramics and metals under shock loading conditions. The line-imaging VISAR provides a way to measure this variability, which may then be quantified by Weibull statistics or other methods. Weibull methods assign a 2-parameter representation of failure phenomena and variability. We have conducted experiments with tantalum (25 and 40 μm grains) and silicon carbide (SiC-N with 5 μm grains). The tantalum HEL variability did not depend systematically on peak stress, grain size or sample thickness, although the previously observed precursor attenuation was present. SiC-N HEL variability within a single shot was approximately half that of single-point variability in a large family of shots; these results are more consistent with sample-to-sample variation than with variability due to changing shot parameters. © 2007 American Institute of Physics.

Fracture of materials has a huge consequence in our daily life ranging from structural damage to loss of life. Understanding the mechanism of crack initiation and propagation in materials is very important. Great effort, both theoretically and experimentally, has been made to understand the nature of crack propagation in crystalline materials. However, crack propagation in disordered systems such as highly cross-linked polymers (e.g. epoxies) is less understood. Many composites such as carbon fibers have an epoxy matrix, and thus it is important to understand the epoxy properties by themselves. We study fracture in highly cross-linked polymer networks bonded to a solid surface using large-scale molecular dynamics simulations. An initial crack is created by forbidding bonds to occur on a fraction of the solid surface up to a crack tip. The time and length scales involved in this process dictate the use of coarse grained bead-spring model of the epoxy network. In order to avoid unwanted boundary effects, large systems of up to 300 000 particles are used. Stress-strain curves are determined for each system from tensile pull molecular dynamics simulations. We found that crack propagation and also formation of voids ahead of the crack are directly related to the network structure.

A growing number of applications involve the transmission of high-intensity laser pulses through optical fibers. Previously, our particular interests led to a series of studies on single-fiber transmission of Q-switched, 1064 nm pulses from multimode Nd:YAG lasers through step-index, multimode, fused silica fibers. The maximum pulse energy that could be transmitted through a given fiber was limited by the onset of laser-induced breakdown or damage. Breakdown at the fiber entrance face was often the first limiting process encountered, but other mechanisms were observed that could result in catastrophic damage at either fiber face, within the initial "entry" segment of the fiber, and at other internal sites along the fiber path. These studies examined system elements that can govern the relative importance of different damage mechanisms, including laser characteristics, the design and alignment of laser-to-flber injection optics, fiber end-face preparation, and fiber routing. In particular, criteria were established for injection optics in order to maximize margins between transmission requirements and thresholds for laser-induced damage. Recent interests have led us to examine laser injection into multiple fibers. Effective methods for generating multiple beams are available, but the resulting beam geometry can lead to challenges in applying the criteria for optimum injection optics. To illustrate these issues, we have examined a three-fiber injection system consisting of a beam-shaping element, a primary injection lens, and a grating beamsplitter. Damage threshold characteristics were established by testing fibers using the injection geometry imposed by this system design.

Anhydrous cerium bromide (CeBr3) and cerium doped lanthanum bromide (Ce+3-LaBr3) were obtained by the dehydration of hydrates synthesized by a direct acidification process. The dehydration process involves heating in vacuum through three phase changes - hydrate, amorphous, and crystalline LaBr3. Incomplete removal of the bound water leads to the formation of oxybromides and the partial reduction of the lanthanum at high temperatures. It was found that upon the completion of dehydration (< 200°C) a complete solid solution can be formed between LaBr3 and CeBr3. These two compounds form a simple binary phase diagram. Challenges associated with the dehydration process are discussed.

We have synthesized and tested new highly fluorescent metal organic framework (MOF) materials based on stilbene dicarboxylic acid as a linker. The crystal structure and porosity of the product are dependent on synthetic conditions and choice of solvent and a low-density cubic form has been identified by x-ray diffraction. In this work we report experiments demonstrating scintillation properties of these crystals. Bright proton-induced luminescence with large shifts relative to the fluorescence excitation spectra were recorded, peaking near 475 nm. Tolerance to fast proton radiation was evaluated by monitoring this radio-luminescence to absorbed doses of several hundred MRAD.