A novel dual stage chemiluminescence detection system incorporating individually controlled hot stages has been developed and applied to probe for material interaction effects during polymer degradation. Utilization of this system has resulted in experimental confirmation for the first time that in an oxidizing environment a degrading polymer A (in this case polypropylene, PP) is capable of infecting a different polymer B (in this case polybutadiene, HTPB) over a relatively large distance. In the presence of the infectious degrading polymer A, the thermal degradation of polymer B is observed over a significantly shorter time period. Consistent with infectious volatiles from material A initiating the degradation process in material B it was demonstrated that traces (micrograms) of a thermally sensitive peroxide in the vicinity of PP could induce degradation remotely. This observation documents cross-infectious phenomena between different polymers and has major consequences for polymer interactions, understanding fundamental degradation processes and long-term aging effects under combined material exposures.
Photocatalytic porphyrins are used to reduce metal complexes from aqueous solution and, further, to control the deposition of metals onto porphyrin nanotubes and surfactant assembly templates to produce metal composite nanostructures and nanodevices. For example, surfactant templates lead to spherical platinum dendrites and foam-like nanomaterials composed of dendritic platinum nanosheets. Porphyrin nanotubes are reported for the first time, and photocatalytic porphyrin nanotubes are shown to reduce metal complexes and deposit the metal selectively onto the inner or outer surface of the tubes, leading to nanotube-metal composite structures that are capable of hydrogen evolution and other nanodevices.
This report summarizes the results of a five-year effort to understand the mechanisms and develop models that predict the corrosion of refractories in oxygen-fuel glass-melting furnaces. Thermodynamic data for the Si-O-(Na or K) and Al-O-(Na or K) systems are reported, allowing equilibrium calculations to be performed to evaluate corrosion of silica- and alumina-based refractories under typical furnace operating conditions. A detailed analysis of processes contributing to corrosion is also presented. Using this analysis, a model of the corrosion process was developed and used to predict corrosion rates in an actual industrial glass furnace. The rate-limiting process is most likely the transport of NaOH(gas) through the mass-transport boundary layer from the furnace atmosphere to the crown surface. Corrosion rates predicted on this basis are in better agreement with observation than those produced by any other mechanism, although the absolute values are highly sensitive to the crown temperature and the NaOH(gas) concentration at equilibrium and at the edge of the boundary layer. Finally, the project explored the development of excimer laser induced fragmentation (ELIF) fluorescence spectroscopy for the detection of gas-phase alkali hydroxides (e.g., NaOH) that are predicted to be the key species causing accelerated corrosion in these furnaces. The development of ELIF and the construction of field-portable instrumentation for glass furnace applications are reported and the method is shown to be effective in industrial settings.
Proposed for publication in Journal of Applied Physics.
Computation of space-charge current-limiting effects across a vacuum cavity between parallel electrodes has previously been carried out only for thermionic emission spectra. In some applications, where the current arises from an injected electron beam or photo-Compton emission from electrode walls, the electron energy spectra may deviate significantly from Maxwellian. Considering the space charge as a collisionless plasma, we derive an implicit equation for the peak cavity potential assuming steady-state currents. For the examples of graphite, nickel, and gold electrodes exposed to x rays, we find that cavity photoemission currents are typically more severely space-charge limited than they would be with the assumption of a purely Maxwellian energy distribution.
Network administrators and security analysts often do not know what network services are being run in every corner of their networks. If they do have a vague grasp of the services running on their networks, they often do not know what specific versions of those services are running. Actively scanning for services and versions does not always yield complete results, and patch and service management, therefore, suffer. We present Net-State, a system for monitoring, storing, and reporting application and operating system version information for a network. NetState gives security and network administrators the ability to know what is running on their networks while allowing for user-managed machines and complex host configurations. Our architecture uses distributed modules to collect network information and a centralized server that stores and issues reports on that collected version information. We discuss some of the challenges to building and operating NetState as well as the legal issues surrounding the promiscuous capture of network data. We conclude that this tool can solve some key problems in network management and has a wide range of possibilities for future uses.
The large number of government and industry activities supporting the Unit of Action (UA), with attendant documents, reports and briefings, can overwhelm decision-makers with an overabundance of information that hampers the ability to make quick decisions often resulting in a form of gridlock. In particular, the large and rapidly increasing amounts of data and data formats stored on UA Advanced Collaborative Environment (ACE) servers has led to the realization that it has become impractical and even impossible to perform manual analysis leading to timely decisions. UA Program Management (PM UA) has recognized the need to implement a Decision Support System (DSS) on UA ACE. The objective of this document is to research the commercial Knowledge Discovery and Data Mining (KDDM) market and publish the results in a survey. Furthermore, a ranking mechanism based on UA ACE-specific criteria has been developed and applied to a representative set of commercially available KDDM solutions. In addition, an overview of four R&D areas identified as critical to the implementation of DSS on ACE is provided. Finally, a comprehensive database containing detailed information on surveyed KDDM tools has been developed and is available upon customer request.
The stainless steel alloy 17-4PH contains a martensitic microstructure and second phase delta ({delta}) ferrite. Strengthening of 17-4PH is attributed to Cu-rich precipitates produced during age hardening treatments at 900-1150 F (H900-H1150). For wrought 17-4PH, the effects of heat treatment and microstructure on mechanical properties are well-documented [for example, Ref. 1]. Fewer studies are available on cast 17-4PH, although it has been a popular casting alloy for high strength applications where moderate corrosion resistance is needed. Microstructural features and defects particular to castings may have adverse effects on properties, especially when the alloy is heat treated to high strength. The objective of this work was to outline the effects of microstructural features specific to castings, such as shrinkage/solidification porosity, on the mechanical behavior of investment cast 17-4PH. Besides heat treatment effects, the results of metallography and SEM studies showed that the largest effect on mechanical properties is from shrinkage/solidification porosity. Figure 1a shows stress-strain curves obtained from samples machined from castings in the H925 condition. The strength levels were fairly similar but the ductility varied significantly. Figure 1b shows an example of porosity on a fracture surface from a room-temperature, quasi-static tensile test. The rounded features represent the surfaces of dendrites which did not fuse or only partially fused together during solidification. Some evidence of local areas of fracture is found on some dendrite surfaces. The shrinkage pores are due to inadequate backfilling of liquid metal and simultaneous solidification shrinkage during casting. A summary of percent elongation results is displayed in Figure 2a. It was found that higher amounts of porosity generally result in lower ductility. Note that the porosity content was measured on the fracture surfaces. The results are qualitatively similar to those found by Gokhale et al. and Surappa et al. in cast A356 Al and by Gokhale et al. for a cast Mg alloys. The quantitative fractography and metallography work by Gokhale et al. illustrated the strong preference for fracture in regions of porosity in cast material. That is, the fracture process is not correlated to the average microstructure in the material but is related to the extremes in microstructure (local regions of high void content). In the present study, image analysis on random cross-sections of several heats indicated an overall porosity content of 0.03%. In contrast, the area % porosity was as high as 16% when measured on fracture surfaces of tensile specimens using stereology techniques. The results confirm that the fracture properties of cast 17-4PH cannot be predicted based on the overall 'average' porosity content in the castings.
This paper describes the development of a set of software tools useful for analyzing ultra-wideband (UWB) antennas and structures. These tools are used to perform finite difference time domain (FDTD) simulation of a conical antenna with continuous wave (CW) and UWB pulsed excitations. The antenna is analyzed using spherical coordinate-based FDTD equations that are derived from first principles. The simulation results for CW excitation are compared to simulation and measured results from published sources; the results for UWB excitation are new.
This document contains a description of the verification and validation process used for the RADTRAN 5.5 code. The verification and validation process ensured the proper calculational models and mathematical and numerical methods were used in the RADTRAN 5.5 code for the determination of risk and consequence assessments. The differences between RADTRAN 5 and RADTRAN 5.5 are the addition of tables, an expanded isotope library, and the additional User-Defined meteorological option for accident dispersion. 3
A continuum-scale, evolutionary model of bubble nucleation, growth and He release for aging metal tritides is described which accounts for major features of the tritide database. Bubble nucleation, modeled as self-trapping of interstitially diffusing He atoms, occurs during the first few days following tritium introduction into the metal. Bubble growth by dislocation loop punching yields good agreement between He atomic volumes and bubble pressures determined from bulk swelling and 3He NMR data. The bubble spacing distribution determined from NMR is shown to remain fixed with age, justifying the separation of nucleation and growth phases and providing a sensitive test of the growth formulation. Late in life, bubble interactions are proposed to produce cooperative stress effects, which lower the bubble pressure. Helium generated near surfaces and surface-connected porosity accounts for the low-level early helium release. Use of an average ligament stress criterion predicts an onset of inter-bubble fracture in good agreement with the He/Metal ratio observed for rapid He release. From the model, it is concluded that He retention can be controlled through control of bubble nucleation.
As interest in 3D face recognition increases the importance of the initial alignment problem does as well. In this paper we present a method utilizing the registered 2D color and range image of a face to automatically identify the eyes, nose, and mouth. These features are important to initially align faces in both standard 2D and 3D face recognition algorithms. For our algorithm to run as fast as possible, we focus on the 2D color information. This allows the algorithm to run in approximately 4 seconds on a 640X480 image with registered range data. On a database of 1,500 images the algorithm achieved a facial feature detection rate of 99.6% with 0.4% of the images skipped due to hair obstruction of the face.
A continuum-scale, evolutionary model of bubble nucleation, growth and He release for aging metal tritides is described which accounts for major features of the tritide database. Bubble nucleation, modeled as self-trapping of interstitially diffusing He atoms, occurs during the first few days following tritium introduction into the metal. Bubble growth by dislocation loop punching yields good agreement between He atomic volumes and bubble pressures determined from bulk swelling and 3He NMR data. The bubble spacing distribution determined from NMR is shown to remain fixed with age, justifying the separation of nucleation and growth phases and providing a sensitive test of the growth formulation. Late in life, bubble interactions are proposed to produce cooperative stress effects, which lower the bubble pressure. Helium generated near surfaces and surface-connected porosity accounts for the low-level early helium release. Use of an average ligament stress criterion predicts an onset of inter-bubble fracture in good agreement with the He/Metal ratio observed for rapid He release. From the model, it is concluded that He retention can be controlled through control of bubble nucleation.
The Z-Pinch Power Plant uses the results from Sandia National Laboratories' Z accelerator in a power plant application to generate energy pulses using inertial confinement fusion. A collaborative project has been initiated by Sandia to investigate the scientific principles of a power generation system using this technology. Research is under way to develop an integrated concept that describes the operational issues of a 1000 MW electrical power plant. Issues under consideration include: 1-20 gigajoule fusion pulse containment, repetitive mechanical connection of heavy hardware, generation of terawatt pulses every 10 seconds, recycling often thousand tons of steel, and manufacturing of millions of hohlraums and capsules per year. Additionally, waste generation and disposal issues are being examined. This paper describes the current concept for the plant and also the objectives for future research.
The diminished response of thermoluminescent phosphors over time is a well-documented challenge to thermoluminescent dosimetry. Wide ranges in fading rates for various phosphor types have been reported, making it necessary for many external dosimetry programs to perform individual studies on thermoluminescent fade. Sandia National Laboratories currently uses the Thermo 8802 LiF:Mg,Ti thermoluminescent dosimeter (TLD) in its personnel external dosimetry program. Doses received in the field are calculated by applying a fade algorithm published by the manufacturer to TLD readings. Since the algorithm was established by characterizing the diminished response of a TLD similar to the 8802, Sandia chose to model its fade study after the analysis done by Thermo. As a result, the parameters of each experiment were comparable, and data from the two studies were compared to determine whether or not the current algorithm should be modified specifically for use at Sandia. Cards were irradiated using an internal 90Sr/90Y source, and pre- and post-irradiation fading rates were monitored over a period of 18 wk. While significant fading was demonstrated, results closely matched those found in the original Thermo study.
The molten salt Flibe, a combination of lithium and beryllium fluorides studied for molten salt fission reactors, has been proposed as a breeder and coolant for fusion applications. The melting points of 2LiF-BeF2 and LiF-BeF2 are 460°C and 363°C, but LiF-BeF2 is rather viscous and has less lithium for breeding. In the Advanced Power Extraction (APEX) Program, concepts with a free flowing liquid for the first wall and blanket were investigated. Flinabe (a mixture of LiF, BeF2 and NaF) was selected for a molten salt design because a melting temperature below 350°C appeared possible and this provided an attractive operating temperature window for a reactor. To confirm that a ternary salt with a low melting temperature existed, several combinations of the fluoride salts, LiF, NaF and BeF2, were melted in a stainless steel crucible under vacuum. One had an apparent melting temperature of 305°C. The test system, preparation of the mixtures, melting procedures and temperature curves for the melting and cooling are presented along with the apparent melting points. Thermal modeling of the salt pool and crucible is reported in an accompanying paper.
Ion mobility spectrometry (IMS) is recognized as one of the most sensitive and versatile techniques for the detection of trace levels of organic vapors. IMS is widely used for detecting contraband narcotics, explosives, toxic industrial compounds and chemical warfare agents. Increasing threat of terrorist attacks, the proliferation of narcotics, Chemical Weapons Convention treaty verification as well as humanitarian de-mining efforts has mandated that equal importance be placed on the analysis time as well as the quality of the analytical data. (1) IMS is unrivaled when both speed of response and sensitivity has to be considered. (2) With conventional (signal averaging) IMS systems the number of available ions contributing to the measured signal to less than 1%. Furthermore, the signal averaging process incorporates scan-to-scan variations decreasing resolution. With external second gate Fourier Transform ion mobility spectrometry (FT-IMS), the entrance gate frequency is variable and can be altered in conjunction with other data acquisition parameters to increase the spectral resolution. The FT-IMS entrance gate operates with a 50% duty cycle and so affords a 7 to 10-fold increase in sensitivity. Recent data on high explosives are presented to demonstrate the parametric optimization in sensitivity and resolution of our system.
This report summarizes the work performed as part of a one-year LDRD project, 'Evolutionary Complexity for Protection of Critical Assets.' A brief introduction is given to the topics of genetic algorithms and genetic programming, followed by a discussion of relevant results obtained during the project's research, and finally the conclusions drawn from those results. The focus is on using genetic programming to evolve solutions for relatively simple algebraic equations as a prototype application for evolving complexity in computer codes. The results were obtained using the lil-gp genetic program, a C code for evolving solutions to user-defined problems and functions. These results suggest that genetic programs are not well-suited to evolving complexity for critical asset protection because they cannot efficiently evolve solutions to complex problems, and introduce unacceptable performance penalties into solutions for simple ones.
The Sandia National Laboratories Corporate Mentor Program provides a mechanism for the development and retention of Sandia's people and knowledge. The relationships formed among staff members at different stages in their careers offer benefits to all. These relationships can provide experienced employees with new ideas and insight and give less experienced employees knowledge of Sandia's culture, strategies, and programmatic direction. The program volunteer coordinators are dedicated to the satisfaction of the participants, who come from every area of Sandia. Since its inception in 1995, the program has sustained steady growth and excellent customer satisfaction. This report summarizes the accomplishments, activities, enhancements, and evaluation data for the Corporate Mentor Program for the 2003/2004 program year ending May 1, 2004.
This SAND report provides the technical progress through June 2004 of the Sandia-led project, ''Carbon Sequestration in Synechococcus Sp.: From Molecular Machines to Hierarchical Modeling'', funded by the DOE Office of Science Genomes to Life Program. Understanding, predicting, and perhaps manipulating carbon fixation in the oceans has long been a major focus of biological oceanography and has more recently been of interest to a broader audience of scientists and policy makers. It is clear that the oceanic sinks and sources of CO{sub 2} are important terms in the global environmental response to anthropogenic atmospheric inputs of CO{sub 2} and that oceanic microorganisms play a key role in this response. However, the relationship between this global phenomenon and the biochemical mechanisms of carbon fixation in these microorganisms is poorly understood. In this project, we will investigate the carbon sequestration behavior of Synechococcus Sp., an abundant marine cyanobacteria known to be important to environmental responses to carbon dioxide levels, through experimental and computational methods. This project is a combined experimental and computational effort with emphasis on developing and applying new computational tools and methods. Our experimental effort will provide the biology and data to drive the computational efforts and include significant investment in developing new experimental methods for uncovering protein partners, characterizing protein complexes, identifying new binding domains. We will also develop and apply new data measurement and statistical methods for analyzing microarray experiments. Computational tools will be essential to our efforts to discover and characterize the function of the molecular machines of Synechococcus. To this end, molecular simulation methods will be coupled with knowledge discovery from diverse biological data sets for high-throughput discovery and characterization of protein-protein complexes. In addition, we will develop a set of novel capabilities for inference of regulatory pathways in microbial genomes across multiple sources of information through the integration of computational and experimental technologies. These capabilities will be applied to Synechococcus regulatory pathways to characterize their interaction map and identify component proteins in these pathways. We will also investigate methods for combining experimental and computational results with visualization and natural language tools to accelerate discovery of regulatory pathways. The ultimate goal of this effort is develop and apply new experimental and computational methods needed to generate a new level of understanding of how the Synechococcus genome affects carbon fixation at the global scale. Anticipated experimental and computational methods will provide ever-increasing insight about the individual elements and steps in the carbon fixation process, however relating an organism's genome to its cellular response in the presence of varying environments will require systems biology approaches. Thus a primary goal for this effort is to integrate the genomic data generated from experiments and lower level simulations with data from the existing body of literature into a whole cell model. We plan to accomplish this by developing and applying a set of tools for capturing the carbon fixation behavior of complex of Synechococcus at different levels of resolution. Finally, the explosion of data being produced by high-throughput experiments requires data analysis and models which are more computationally complex, more heterogeneous, and require coupling to ever increasing amounts of experimentally obtained data in varying formats. These challenges are unprecedented in high performance scientific computing and necessitate the development of a companion computational infrastructure to support this effort.
This manual describes the input syntax to the ALEGRA radiation transport package. All input and output variables are defined, as well as all algorithmic controls. This manual describes the radiation input syntax for ALEGRA-HEDP. The ALEGRA manual[2] describes how to run the code and general input syntax. The ALEGRA-HEDP manual[13] describes the input for other physics used in high energy density physics simulations, as well as the opacity models used by this radiation package. An emission model, which is the lowest order radiation transport approximation, is also described in the ALEGRA-HEDP manual. This document is meant to be used with these other manuals.
PDC drill bit performance has been greatly improved over the past three decades by innovations in bit design and how these designs are applied. The next leap forward is most likely to come from using high-speed, real-time downhole data to optimize the performance of these sophisticated bits on an application-by-application basis. By effectively managing conditions of lateral, axial and torsional acceleration, damage to cutting structures can be minimized for improved penetration rates. Avoiding these damaging vibrations is essential to increasing bit durability and improving overall drilling economics. This paper describes one set of independent drilling optimization results obtained as part of a series of controlled demonstrations of PDC bits. Sandia National Laboratories on behalf of the U. S. Department of Energy (DOE) managed this work. The effort was organized as a Cooperative Research and Development Agreement (CRADA) established between Sandia and four bit manufacturers with DOE funding and in-kind contributions by the industry partners. The goal of this CRADA was to demonstrate drag bit performance in formations with degrees of hardness typical of those encountered while drilling geothermal wells. The test results indicate that the surface weight-on-bit (WOB), revolutions per minute (RPM) and torque readings traditionally used to guide adjustments in the drilling parameters do not always provide the true picture of what is really taking place at the bit. Instead, a holistic approach combining traditional methods of optimization together with high-speed, real-time data enables far better decisions for improving bit performance and avoiding damaging situations that may have otherwise gone unnoticed.
ALEGRA is an arbitrary Lagrangian-Eulerian finite element code that emphasizes large distortion and shock propagation in inviscid fluids and solids. This document describes user options for modeling resistive magnetohydrodynamic, thermal conduction, and radiation emission effects.
An experimental investigation is made into the fluid mechanics and heat transfer of a circular cylinder immersed in a wall-bounded turbulent mixed-convection flow of water. The cylinder is oriented spanwise to the forced channel flow and within the thermal boundary layer of the heated lower wall. The flow channel is capped with a cold, near-adiabatic upper wall producing a fully turbulent gap Rayleigh number of 108. A low-speed crossflow is applied to advect the turbulent thermal plumes over the cylinder surface. We present spatially resolved cylinder-surface heat-flux data alongside 2-D PIV imaging of the streamwise and wall-normal velocity components for two flow conditions in the mixed-convection heat-transfer regime. The measured cylinder-wake flowfield reflects the complex coupling between the separated wake flow, the highly turbulent freestream and the buoyant wall and cylinder boundary layers. A method for measurement of spatially resolved surface heat fluxes based on the measured cylinder-surface temperature distribution and a well-posed two-dimensional solution to the conduction problem in the cylinder wall is presented. The resulting spatially resolved flux measurements show enhanced surface heat transfer, which results from the intense buoyancy generated free-stream turbulence and mixing in the cylinder wake. This work extends the literature on thermal convection with crossflow well into the turbulent regime and is, to our knowledge, the first investigation of surface heat-transfer to an object of engineering importance placed in this type of turbulent mixed-convection flowfield. The data are currently being utilized for validation of mixed-convection turbulence models at Sandia and comparisons between the computational and experimental results are presented.
Mobile manipulator systems used by emergency response operators consist of an articulated robot arm, a remotely driven base, a collection of cameras, and a remote communications link. Typically the system is completely teleoperated, with the operator using live video feedback to monitor and assess the environment, plan task activities, and to conduct the operations via remote control input devices. The capabilities of these systems are limited, and operators rarely attempt sophisticated operations such as retrieving and utilizing tools, deploying sensors, or building up world models. This project has focused on methods to utilize this video information to enable monitored autonomous behaviors for the mobile manipulator system, with the goal of improving the overall effectiveness of the human/robot system. Work includes visual servoing, visual targeting, utilization of embedded video in 3-D models, and improved methods of camera utilization and calibration.
Treatment systems that can neutralize biological agents are needed to mitigate risks from novel and legacy biohazards. Tests with Bacillus thuringiensis and Bacillus steurothemophilus spores were performed in a 190-liter, 1-112 lb TNT equivalent rated Explosive Destruction System (EDS) system to evaluate its capability to treat and destroy biological agents. Five tests were conducted using three different agents to kill the spores. The EDS was operated in steam autoclave, gas fumigation and liquid decontamination modes. The first three tests used EDS as an autoclave, which uses pressurized steam to kill the spores. Autoclaving was performed at 130-140 deg C for up to 2-hours. Tests with chlorine dioxide at 750 ppm concentration for 1 hour and 10% (vol) aqueous chlorine bleach solution for 1 hour were also performed. All tests resulted in complete neutralization of the bacterial spores based on no bacterial growth in post-treatment incubations. Explosively opening a glass container to expose the bacterial spores for treatment with steam was demonstrated and could easily be done for chlorine dioxide gas or liquid bleach.
We have developed a novel approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only a sparse set of distance constraints, such as those derived from MS3-D, dipolar-EPR and FRET experiments. Algorithms have been written for searching the conformational space of membrane protein folds matching the set of distance constraints, which provides initial structures for local conformational searches. Local conformation search is achieved by optimizing these candidates against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. This results in refined helical bundles to which the interhelical loops and amino acid side-chains are added. Using a set of only 27 distance constraints extracted from the literature, our methods successfully recover the structure of dark-adapted rhodopsin to within 3.2 {angstrom} of the crystal structure.
We report a new nanolaser technique for measuring characteristics of human mitochondria. Because mitochondria are so small, it has been difficult to study large populations using standard light microscope or flow cytometry techniques. We recently discovered a nano-optical transduction method for high-speed analysis of submicron organelles that is well suited to mitochondrial studies. This ultrasensitive detection technique uses nano-squeezing of light into photon modes imposed by the ultrasmall organelle dimensions in a semiconductor biocavity laser. In this paper, we use the method to study the lasing spectra of normal and diseased mitochondria. We find that the diseased mitochondria exhibit larger physical diameter and standard deviation. This morphological differences are also revealed in the lasing spectra. The diseased specimens have a larger spectral linewidth than the normal, and have more variability in their statistical distributions.
The brain is often identified with decision-making processes in the biological world. In fact, single cells, single macromolecules (proteins) and populations of molecules also make simple decisions. These decision processes are essential to survival and to the biological self-assembly and self-repair processes that we seek to emulate. How do these tiny systems make effective decisions? How do they make decisions in concert with a cooperative network of other molecules or cells? How can we emulate the decision-making behaviors of small-scale biological systems to program and self-assemble microsystems? This LDRD supported research to answer these questions. Our work included modeling and simulation of protein populations to help us understand, mimic, and categorize molecular decision-making mechanisms that nonequilibrium systems can exhibit. This work is an early step towards mimicking such nanoscale and microscale biomolecular decision-making processes in inorganic systems.
Chemically prepared zinc oxide powders are fabricated for the production of high aspect ratio varistor components. Colloidal processing in water was performed to reduce agglomerates to primary particles, form a high solids loading slurry, and prevent dopant migration. The milled and dispersed powder exhibited a viscoelastic to elastic behavioral transition at a volume loading of 43-46%. The origin of this transition was studied using acoustic spectroscopy, zeta potential measurements and oscillatory rheology. The phenomenon occurs due to a volume fraction solids dependent reduction in the zeta potential of the solid phase. It is postulated to result from divalent ion binding within the polyelectrolyte dispersant chain, and was mitigated using a polyethylene glycol plasticizing additive. Chemically prepared zinc oxide powders were processed for the production of high aspect ratio varistor components. Near net shape casting methods including slip casting and agarose gelcasting were evaluated for effectiveness in achieving a uniform green microstructure achieving density values near the theoretical maximum during sintering. The structure of the green parts was examined by mercury porisimetry. Agarose gelcasting produced green parts with low solids loading values and did not achieve high fired density. Isopressing the agarose cast parts after drying raised the fired density to greater than 95%, but the parts exhibited catastrophic shorting during electrical testing. Slip casting produced high green density parts, which exhibited high fired density values. The electrical characteristics of slip cast parts are comparable with dry pressed powder compacts. Alternative methods for near net shape forming of ceramic dispersions were investigated for use with the chemically prepared ZnO material. Recommendations for further investigation to achieve a viable production process are presented.
ALEGRA is an arbitrary Lagrangian-Eulerian multi-material finite element code used for modeling solid dynamics problems involving large distortion and shock propagation. This document describes the basic user input language and instructions for using the software.
An exploratory effort in the application of carbon epoxy composite structural materials to a multi-axis gimbal arm design is described. An existing design in aluminum was used as a baseline for a functionally equivalent redesigned outer gimbal arm using a carbon epoxy composite material. The existing arm was analyzed using finite element techniques to characterize performance in terms of strength, stiffness, and weight. A new design was virtually prototyped. using the same tools to produce a design with similar stiffness and strength, but reduced overall weight, than the original arm. The new design was prototyped using Rapid Prototyping technology, which was subsequently used to produce molds for fabricating the carbon epoxy composite parts. The design tools, process, and results are discussed.
Sensitivity analysis is critically important to numerous analysis algorithms, including large scale optimization, uncertainty quantification,reduced order modeling, and error estimation. Our research focused on developing tools, algorithms and standard interfaces to facilitate the implementation of sensitivity type analysis into existing code and equally important, the work was focused on ways to increase the visibility of sensitivity analysis. We attempt to accomplish the first objective through the development of hybrid automatic differentiation tools, standard linear algebra interfaces for numerical algorithms, time domain decomposition algorithms and two level Newton methods. We attempt to accomplish the second goal by presenting the results of several case studies in which direct sensitivities and adjoint methods have been effectively applied, in addition to an investigation of h-p adaptivity using adjoint based a posteriori error estimation. A mathematical overview is provided of direct sensitivities and adjoint methods for both steady state and transient simulations. Two case studies are presented to demonstrate the utility of these methods. A direct sensitivity method is implemented to solve a source inversion problem for steady state internal flows subject to convection diffusion. Real time performance is achieved using novel decomposition into offline and online calculations. Adjoint methods are used to reconstruct initial conditions of a contamination event in an external flow. We demonstrate an adjoint based transient solution. In addition, we investigated time domain decomposition algorithms in an attempt to improve the efficiency of transient simulations. Because derivative calculations are at the root of sensitivity calculations, we have developed hybrid automatic differentiation methods and implemented this approach for shape optimization for gas dynamics using the Euler equations. The hybrid automatic differentiation method was applied to a first order approximation of the Euler equations and used as a preconditioner. In comparison to other methods, the AD preconditioner showed better convergence behavior. Our ultimate target is to perform shape optimization and hp adaptivity using adjoint formulations in the Premo compressible fluid flow simulator. A mathematical formulation for mixed-level simulation algorithms has been developed where different physics interact at potentially different spatial resolutions in a single domain. To minimize the implementation effort, explicit solution methods can be considered, however, implicit methods are preferred if computational efficiency is of high priority. We present the use of a partial elimination nonlinear solver technique to solve these mixed level problems and show how these formulation are closely coupled to intrusive optimization approaches and sensitivity analyses. Production codes are typically not designed for sensitivity analysis or large scale optimization. The implementation of our optimization libraries into multiple production simulation codes in which each code has their own linear algebra interface becomes an intractable problem. In an attempt to streamline this task, we have developed a standard interface between the numerical algorithm (such as optimization) and the underlying linear algebra. These interfaces (TSFCore and TSFCoreNonlin) have been adopted by the Trilinos framework and the goal is to promote the use of these interfaces especially with new developments. Finally, an adjoint based a posteriori error estimator has been developed for discontinuous Galerkin discretization of Poisson's equation. The goal is to investigate other ways to leverage the adjoint calculations and we show how the convergence of the forward problem can be improved by adapting the grid using adjoint-based error estimates. Error estimation is usually conducted with continuous adjoints but if discrete adjoints are available it may be possible to reuse the discrete version for error estimation. We investigate the advantages and disadvantages of continuous and discrete adjoints through a simple example.
The purpose of the Sandia National Laboratories Advanced Simulation and Computing (ASC) Software Quality Plan is to clearly identify the practices that are the basis for continually improving the quality of ASC software products. The plan defines the ASC program software quality practices and provides mappings of these practices to Sandia Corporate Requirements CPR 1.3.2 and 1.3.6 and to a Department of Energy document, 'ASCI Software Quality Engineering: Goals, Principles, and Guidelines'. This document also identifies ASC management and software project teams responsibilities in implementing the software quality practices and in assessing progress towards achieving their software quality goals.
The purpose of the Sandia National Laboratories (SNL) Advanced Simulation and Computing (ASC) Software Quality Plan is to clearly identify the practices that are the basis for continually improving the quality of ASC software products. Quality is defined in DOE/AL Quality Criteria (QC-1) as conformance to customer requirements and expectations. This quality plan defines the ASC program software quality practices and provides mappings of these practices to the SNL Corporate Process Requirements (CPR 1.3.2 and CPR 1.3.6) and the Department of Energy (DOE) document, ASCI Software Quality Engineering: Goals, Principles, and Guidelines (GP&G). This quality plan identifies ASC management and software project teams' responsibilities for cost-effective software engineering quality practices. The SNL ASC Software Quality Plan establishes the signatories commitment to improving software products by applying cost-effective software engineering quality practices. This document explains the project teams opportunities for tailoring and implementing the practices; enumerates the practices that compose the development of SNL ASC's software products; and includes a sample assessment checklist that was developed based upon the practices in this document.
Macroscopic quantum states such as superconductors, Bose-Einstein condensates and superfluids are some of the most unusual states in nature. In this project, we proposed to design a semiconductor system with a 2D layer of electrons separated from a 2D layer of holes by a narrow (but high) barrier. Under certain conditions, the electrons would pair with the nearby holes and form excitons. At low temperature, these excitons could condense to a macroscopic quantum state either through a Bose-Einstein condensation (for weak exciton interactions) or a BCS transition to a superconductor (for strong exciton interactions). While the theoretical predictions have been around since the 1960's, experimental realization of electron-hole bilayer systems has been extremely difficult due to technical challenges. We identified four characteristics that if successfully incorporated into a device would give the best chances for excitonic condensation to be observed. These characteristics are closely spaced layers, low disorder, low density, and independent contacts to allow transport measurements. We demonstrated each of these characteristics separately, and then incorporated all of them into a single electron-hole bilayer device. The key to the sample design is using undoped GaAs/AlGaAs heterostructures processed in a field-effect transistor geometry. In such samples, the density of single 2D layers of electrons could be varied from an extremely low value of 2 x 10{sup 9} cm{sup -2} to high values of 3 x 10{sup 11} cm{sup -2}. The extreme low values of density that we achieved in single layer 2D electrons allowed us to make important contributions to the problem of the metal insulator transition in two dimensions, while at the same time provided a critical base for understanding low density 2D systems to be used in the electron-hole bilayer experiments. In this report, we describe the processing advances to fabricate single and double layer undoped samples, the low density results on single layers, and evidence for gateable undoped bilayers.
The objective of this LDRD project was to develop a programmable diffraction grating fabricated in SUMMiT V{trademark}. Two types of grating elements (vertical and rotational) were designed and demonstrated. The vertical grating element utilized compound leveraged bending and the rotational grating element used vertical comb drive actuation. This work resulted in two technical advances and one patent application. Also a new optical configuration of the Polychromator was demonstrated. The new optical configuration improved the optical efficiency of the system without degrading any other aspect of the system. The new configuration also relaxes some constraint on the programmable diffraction grating.
We summarize the results of a project to develop evolutionary computing methods for the design of behaviors of embodied agents in the form of autonomous vehicles. We conceived and implemented a strategy called graduated embodiment. This method allows high-level behavior algorithms to be developed using genetic programming methods in a low-fidelity, disembodied modeling environment for migration to high-fidelity, complex embodied applications. This project applies our methods to the problem domain of robot navigation using adaptive waypoints, which allow navigation behaviors to be ported among autonomous mobile robots with different degrees of embodiment, using incremental adaptation and staged optimization. Our approach to biomimetic behavior engineering is a hybrid of human design and artificial evolution, with the application of evolutionary computing in stages to preserve building blocks and limit search space. The methods and tools developed for this project are directly applicable to other agent-based modeling needs, including climate-related conflict analysis, multiplayer training methods, and market-based hypothesis evaluation.
Forward-to-reverse bias step-recovery measurements were performed on In.07Ga.93N/GaN and Al.36Ga.64N/Al.46Ga.54N quantum-well (QW) light-emitting diodes grown on sapphire. With the QW sampling the minority-carrier hole density at a single position, distinctive two-phase optical decay curves were observed. Using diffusion equation solutions to self-consistently model both the electrical and optical responses, hole transport parameters tp = 758 {+-} 44 ns, Lp = 588 {+-} 45 nm, and up = 0.18 {+-} 0.02 cm2/Vs were obtained for GaN. The mobility was thermally activated with an activation energy of 52 meV, suggesting trap-modulated transport. Optical measurements of sub-bandgap peaks exhibited slow responses approaching the bulk lifetime. For Al.46Ga.54N, a longer lifetime of tp = 3.0 us was observed, and the diffusion length was shorter, Lp = 280 nm. Mobility was an order of magnitude smaller than in GaN, up = 10-2 cm2/Vs, and was insensitive to temperature, suggesting hole transport through a network of defects.
This report documents the results of an LDRD program entitled 'System of Systems Modeling and Analysis' that was conducted during FY 2003 and FY 2004. Systems that themselves consist of multiple systems (referred to here as System of Systems or SoS) introduce a level of complexity to systems performance analysis and optimization that is not readily addressable by existing capabilities. The objective of the 'System of Systems Modeling and Analysis' project was to develop an integrated modeling and simulation environment that addresses the complex SoS modeling and analysis needs. The approach to meeting this objective involved two key efforts. First, a static analysis approach, called state modeling, has been developed that is useful for analyzing the average performance of systems over defined use conditions. The state modeling capability supports analysis and optimization of multiple systems and multiple performance measures or measures of effectiveness. The second effort involves time simulation which represents every system in the simulation using an encapsulated state model (State Model Object or SMO). The time simulation can analyze any number of systems including cross-platform dependencies and a detailed treatment of the logistics required to support the systems in a defined mission.
We report micro-Raman studies of self-heating in an AlGaN/GaN heterostructure field-effect transistor using below (visible 488.0 nm) and near (UV 363.8 nm) GaN band-gap excitation. The shallow penetration depth of the UV light allows us to measure temperature rise ({Delta}T) in the two-dimensional electron gas (2DEG) region of the device between drain and source. Visible light gives the average {Delta}T in the GaN layer, and that of the SiC substrate, at the same lateral position. Combined, we depth profile the self-heating. Measured {Delta}T in the 2DEG is consistently over twice the average GaN-layer value. Electrical and thermal transport properties are simulated. We identify a hotspot, located at the gate edge in the 2DEG, as the prevailing factor in the self-heating.
Single molecule fluorophores were studied for the first time with a new confocal fluorescence microscope that allows the wavelength and emission time to be simultaneously measured with single molecule sensitivity. In this apparatus, the photons collected from the sample are imaged through a dispersive optical system onto a time and position sensitive detector. This detector records the wavelength and emission time of each detected photon relative to an excitation laser pulse. A histogram of many events for any selected spatial region or time interval can generate a full fluorescence spectrum and correlated decay plot for the given selection. At the single molecule level, this approach makes entirely new types of temporal and spectral correlation spectroscopy of possible. This report presents the results of simultaneous time- and frequency-resolved fluorescence measurements of single rhodamine 6G (R6G), tetramethylrhodamine (TMR), and Cy3 embedded in thin films of polymethylmethacrylate (PMMA).
Order-of-accuracy verification is necessary to ensure that software correctly solves a given set of equations. One method to verify the order of accuracy of a code is the method of manufactured solutions. In this study, a manufactured solution has been derived and implemented that allows verification of not only the Euler, Navier-Stokes, and Reynolds-Averaged Navier-Stokes (RANS) equation sets, but also some of their associated boundary conditions (BC's): slip, no-slip (adiabatic and isothermal), and outflow (subsonic, supersonic, and mixed). Order-of-accuracy verification has been performed for the Euler and Navier-Stokes equations and these BC's in a compressible computational fluid dynamics code. All of the results shown are on skewed, non-uniform meshes. RANS results will be presented in a future paper. The observed order of accuracy was lower than the expected order of accuracy in two cases. One of these cases resulted in the identification and correction of a coding mistake in the CHAD gradient correction that was reducing the observed order of accuracy. This mistake would have been undetectable on a Cartesian mesh. During the search for the CHAD gradient correction problem, an unrelated coding mistake was found and corrected. The other case in which the observed order of accuracy was less than expected was a test of the slip BC; although no specific coding or formulation mistakes have yet been identified. After the correction of the identified coding mistakes, all of the aforementioned equation sets and BC's demonstrated the expected (or at least acceptable) order of accuracy except the slip condition.
Rate constants for the thermal dissociation of Si{sub 2}H{sub 6} are predicted with a novel transition state model. The saddle points for dissociation on the Si{sub 2}H{sub 6} potential energy surface are lower in energy than the corresponding separated products, as confirmed by high level ab initio quantum mechanical calculations. Thus, the dissociations of Si{sub 2}H{sub 6} to produce SiH{sub 2} + SiH{sub 4} (R1) and H{sub 3}SiSiH + H{sub 2} (R2) both proceed through tight inner transition states followed by loose outer transition states. The present 'dual' transition state model couples variational phase space theory treatments of the outer transition states with ab initio based fixed harmonic vibrator treatments of the inner transition states to obtain effective numbers of states for the two transition states acting in series. It is found that, at least near room temperature, such a dual transition state model is generally required for the proper description of each of the dissociations. Only at quite high temperatures, i.e., above 2000 K for (R1) and 600 K for (R2), does a single fixed inner transition state provide an adequate description. Similarly, only at quite low temperatures (below 100 and 10 K for (R1) and (R2), respectively) does a single outer transition state provide an adequate description. Pressure dependent rate constants are obtained from solutions to the multichannel master equation. These calculations confirm that dissociation channel (R2) is negligible under conditions relevant to the thermal chemical vapor deposition (CVD) processes. Rate constants for the chemical activation reactions, SiH{sub 2} + SiH{sub 4} {yields} Si{sub 2}H{sub 6} (R-1) and SiH{sub 2} + SiH{sub 4} {yields} H{sub 3}SiSiH + H{sub 2} (R3), are also evaluated within the dual transition state model. It is found that reaction R3 is the dominant channel for low pressures and high temperatures, i.e., below 100 Torr for temperatures above 1100 K.
Radio frequency microelectromechanical systems (RF MEMS) are an enabling technology for next-generation communications and radar systems in both military and commercial sectors. RF MEMS-based reconfigurable circuits outperform solid-state circuits in terms of insertion loss, linearity, and static power consumption and are advantageous in applications where high signal power and nanosecond switching speeds are not required. We have demonstrated a number of RF MEMS switches on high-resistivity silicon (high-R Si) that were fabricated by leveraging the volume manufacturing processes available in the Microelectronics Development Laboratory (MDL), a Class-1, radiation-hardened CMOS manufacturing facility. We describe novel tungsten and aluminum-based processes, and present results of switches developed in each of these processes. Series and shunt ohmic switches and shunt capacitive switches were successfully demonstrated. The implications of fabricating on high-R Si and suggested future directions for developing low-loss RF MEMS-based circuits are also discussed.
The focus of this paper is a penalty-based strategy for preconditioning elliptic saddle point systems. As the starting point, we consider the regularization approach of Axelsson in which a related linear system, differing only in the (2,2) block of the coefficient matrix, is introduced. By choosing this block to be negative definite, the dual unknowns of the related system can be eliminated resulting in a positive definite primal Schur complement. Rather than solving the Schur complement system exactly, an approximate solution is obtained using a substructuring preconditioner. The approximate primal solution together with the recovered dual solution then define the preconditioned residual for the original system. The approach can be applied to a variety of different saddle point problems. Although the preconditioner itself is symmetric and indefinite, all the eigenvalues of the preconditioned system are real and positive if certain conditions hold. Stronger conditions also ensure that the eigenvalues are bounded independently of mesh parameters. An interesting feature of the approach is that conjugate gradients can be used as the iterative solution method rather than GMRES. The effectiveness of the overall strategy hinges on the preconditioner for the primal Schur complement. Interestingly, the primary condition ensuring real and positive eigenvalues is satisfied automatically in certain instances if a Balancing Domain Decomposition by Constraints (BDDC) preconditioner is used. Following an overview of BDDC, we show how its constraints can be chosen to ensure insensitivity to parameter choices in the (2,2) block for problems with a divergence constraint. Examples for different saddle point problems are presented and comparisons made with other approaches.
Confinement within the nanoscale pores of a zeolite strongly modifies the behavior of small molecules. Typical of many such interesting and important problems, realistic modeling of this phenomena requires simultaneously capturing the detailed behavior of chemical bonds and the possibility of collective dynamics occurring in a complex unit cell (672 atoms in the case of Zeolite-4A). Classical simulations alone cannot reliably model the breaking and formation of chemical bonds, while quantum methods alone are incapable of treating the extended length and time scales characteristic of complex dynamics. We have developed a robust and efficient model in which a small region treated with the Kohn-Sham density functional theory is embedded within a larger system represented with classical potentials. This model has been applied in concert with first-principles electronic structure calculations and classical molecular dynamics and Monte Carlo simulations to study the behavior of water, ammonia, the hydroxide ion, and the ammonium ion in Zeolite-4a. Understanding this behavior is important to the predictive modeling of the aging of Zeolite-based desiccants. In particular, we have studied the absorption of these molecules, interactions between water and the ammonium ion, and reactions between the hydroxide ion and the zeolite cage. We have shown that interactions with the extended Zeolite cage strongly modifies these local chemical phenomena, and thereby we have proven out hypothesis that capturing both local chemistry and collective phenomena is essential to realistic modeling of this system. Based on our results, we have been able to identify two possible mechanisms for the aging of Zeolite-based desiccants.
Laser diode ignition experiments were conducted in an effort to characterize the effects of scale and heating rate on micro-scale explosive ignition criteria. Over forty experiments were conducted with various laser power densities and laser spot sizes. In addition, relatively simple analytical and numerical calculations were performed to assist with interpretation of the experimental data and characterization of the explosive ignition criteria.
In this article, we discuss in detail the addition of hydrogen atoms to diacetylene and the reverse dissociation reactions, H + C{sub 4}H{sub 2} {leftrightarrow} i-C{sub 4}H{sub 3} (R1) and H + C{sub 4}H{sub 2} n-C{sub 4}H{sub 3} (R2). The theory utilizes high-level electronic structure methodology to characterize the potential energy surface, Rice-Ramsperger-Kassel-Marcus (RRKM) theory to calculate microcanonical/J-resolved rate coefficients, and a two-dimensional master-equation approach to extract phenomenological (thermal) rate coefficients. Comparison is made with experimental results where they are available. The rate coefficients k{sub 1}(T, p) and k{sub 2}(T, p) are cast in forms that can be used in chemical kinetic modeling. In addition, we predict values of the heats of formation of i-C{sub 4}H{sub 3} and n-C{sub 4}H{sub 3} and discuss their importance in flame chemistry. Our basis-set extrapolated, quadratic-configuration-interaction with single and double excitations (and triple excitations added perturbatively), QCISD(T), predictions of these heats of formation at 298 K are 130.8 kcal/mol for n-C{sub 4}H{sub 3} and 119.3 kcal/mol for the i-isomer; multireference CI calculations with a nine-electron, nine-orbital, complete-active-space (CAS) reference wavefunction give just slightly larger values for these parameters. Our results are in good agreement with the recent focal-point analysis of Wheeler et al. (J. Chem. Phys. 2004, 121, 8800-8813), but they differ substantially for {Delta} H{sub f 298}{sup 0}(n-C{sub 4}H{sub 3}) with the earlier diffusion Monte Carlo predictions of Krokidis et al.
A series of numerical simulations have been performed to determine scaling laws for fast ignition break even of a hot spot formed by energetic particles created by a short pulse laser. Hot spot break even is defined to be when the fusion yield is equal to the total energy deposited in the hot spot through both the initial compression and the subsequent heating. In these simulations, only a small portion of a previously compressed mass of deuterium-tritium fuel is heated on a short time scale, i.e., the hot spot is tamped by the cold dense fuel which surrounds it. The hot spot tamping reduces the minimum energy required to obtain break even as compared to the situation where the entire fuel mass is heated, as was assumed in a previous study [S. A. Slutz, R. A. Vesey, I. Shoemaker, T. A. Mehlhorn, and K. Cochrane, Phys. Plasmas 7, 3483 (2004)]. The minimum energy required to obtain hot spot break even is given approximately by the scaling law E{sub T} = 7.5({rho}/100){sup -1.87} kJ for tamped hot spots, as compared to the previously reported scaling of E{sub UT} = 15.3({rho}/100){sup -1.5} kJ for untamped hotspots. The size of the compressed fuel mass and the focusability of the particles generated by the short pulse laser determines which scaling law to use for an experiment designed to achieve hot spot break even.
Many accelerators at Sandia National Laboratories utilize the Rimfire gas switch for high-voltage, high-power switching. Future accelerators will have increased performance requirements for switching elements. When designing improved versions of the Rimfire switch, there is a need for quick and accurate simulation of the electrical effects of geometry changes. This paper presents an advanced circuit model of the Rimfire switch that can be used for these simulations. The development of the model is shown along with comparisons to past models and experimental results.
The Sandia Lightning Simulator at Sandia National Laboratories can provide up to 200 kA for a simulated single lightning stroke, 100 kA for a subsequent stroke, and hundreds of Amperes of continuing current. It has recently been recommissioned after a decade of inactivity and the single-stroke capability demonstrated. The simulator capabilities, basic design components, upgrades, and diagnostic capabilities are discussed in this paper.
An experimental investigation is made into the fluid mechanics and heat transfer of a circular cylinder immersed in a wall-bounded turbulent mixed-convection flow of water. The cylinder is oriented spanwise to the forced channel flow and within the thermal boundary layer of the heated lower wall. The flow channel is capped with a cold, near-adiabatic upper wall producing a fully turbulent gap Rayleigh number of 10{sup 8}. A low-speed crossflow is applied to advect the turbulent thermal plumes over the cylinder surface. We present spatially resolved cylinder-surface heat-flux data alongside 2-D PIV imaging of the streamwise and wall-normal velocity components for two flow conditions in the mixed-convection heat-transfer regime. The measured cylinder-wake flowfield reflects the complex coupling between the separated wake flow, the highly turbulent freestream and the buoyant wall and cylinder boundary layers. A method for measurement of spatially resolved surface heat fluxes based on the measured cylinder-surface temperature distribution and a well-posed two-dimensional solution to the conduction problem in the cylinder wall is presented. The resulting spatially resolved flux measurements show enhanced surface heat transfer, which results from the intense buoyancy generated free-stream turbulence and mixing in the cylinder wake. This work extends the literature on thermal convection with crossflow well into the turbulent regime and is, to our knowledge, the first investigation of surface heat-transfer to an object of engineering importance placed in this type of turbulent mixed-convection flowfield. The data are currently being utilized for validation of mixed convection turbulence models at Sandia and comparisons between the computational and experimental results are presented.
This article provides a brief review of the field of electroporation and introduces a new microdevice that facilitates studies to test theories, gain understanding, and control this important biomedical technology. Electroporation, a bio-electrochemical process whose fundamentals are not yet understood, is a means of permeating the cell membrane by applying a voltage across the cell and forming nano-scale pores in the membrane. It has become an important field in biotechnology and medicine for the controlled introduction of macromolecules, such as gene constructs and drugs, into various cells. It is viewed as an engineering alternative to biological techniques for the genetic engineering of cells. To study and control electroporation, we have created a low-cost microelectroporation chip that incorporates a live biological cell with an electric circuit. The device revealed an important behavior of cells in electrical fields. They produce measurable electrical information about the electroporation state of the cell that may enable precise control of the process. The device can be used to facilitate fundamental studies of electroporation and can become useful in providing precise control over biotechnological processes.
The WLF equation is typically used to describe the dependence of polymer mobility on temperature at atmospheric pressure. Tests at different pressures would at least require different WLF parameterization. Completely different tests, for example, probing the temperature dependence of mobility at constant density, would require even greater modifications. By performing molecular dynamics simulations on simple chain molecules equilibrated at different thermodynamic states, we have shown that the mobility depends in a more general sense on the potential energy density of the system. That is, mobilities for any equilibrated state collapse onto one master curve when plotted against the potential energy density. Moreover, this relationship can be fit by either a 'generalized' WLF equation or by a power-law relationship observed in critical phenomena. When this mobility relationship is used within a rheologically simple, thermodynamically consistent, viscoelastic framework, quantitative agreement is seen between experimental data and theoretical predictions on a range of tests covering enthalpy relaxation to mechanical yield to physical aging.
This document provides a detailed discussion and a guide for the use of the RadCat 2.0 Graphical User Interface input file generator for the RADTRAN 5.5 code. The differences between RadCat 2.0 and RadCat 1.0 can be attributed to the differences between RADTRAN 5 and RADTRAN 5.5 as well as clarification for some of the input parameters. 3
The synthesis of a photoswitchable polymer by grafting an azobenzene dye to methacrylate followed by polymerization is presented. The azobenzene dye undergoes a trans-cis photoisomerization that causes a persistent change in the refractive index of cast polymer films. This novel polymer was incorporated into superlattices prepared by spin casting and the optical activity of the polymer was maintained. A modified coextruder that allows the rapid production of soft matter superlattices was designed and fabricated.
Two methods for creating a hybrid level-set (LS)/particle method for modeling surface evolution during feature-scale etching and deposition processes are developed and tested. The first method supplements the LS method by introducing Lagrangian marker points in regions of high curvature. Once both the particle set and the LS function are advanced in time, minimization of certain objective functions adjusts the LS function so that its zero contour is in closer alignment with the particle locations. It was found that the objective-minimization problem was unexpectedly difficult to solve, and even when a solution could be found, the acquisition of it proved more costly than simply expanding the basis set of the LS function. The second method explored is a novel explicit marker-particle method that we have named the grid point particle (GPP) approach. Although not a LS method, the GPP approach has strong procedural similarities to certain aspects of the LS approach. A key aspect of the method is a surface rediscretization procedure--applied at each time step and based on a global background mesh--that maintains a representation of the surface while naturally adding and subtracting surface discretization points as the surface evolves in time. This method was coded in 2-D, and tested on a variety of surface evolution problems by using it in the ChISELS computer code. Results shown for 2-D problems illustrate the effectiveness of the method and highlight some notable advantages in accuracy over the LS method. Generalizing the method to 3D is discussed but not implemented.
In this study, we describe the extension of the 2-d preliminary design bluff body drag estimation tool developed by De Chant1 to apply for 3-d flows. As with the 2-d method, the 3-d extension uses a combined approximate Green's function/Gram-Charlier series approach to retain the body geometry information. Whereas, the 2-d methodology relied solely upon the use of small disturbance theory for the inviscid flow field associated with the body of interest to estimate the near-field initial conditions, e.g. velocity defect, the 3-d methodology uses both analytical (where available) and numerical inviscid solutions. The defect solution is then used as an initial condition in an approximate 3-d Green's function solution. Finally, the Green's function solution is matched to the 3-d analog of the classical 2-d Gram-Charlier series and then integrated to yield the net form drag on the bluff body. Preliminary results indicate that drag estimates computed are of accuracy equivalent to the 2-d method for flows with large separation, i.e. less than 20% relative error. As was the lower dimensional method, the 3-d concept is intended to be a supplement to turbulent Navier-Stokes and experimental solution for estimating drag coefficients over blunt bodies.
The Method of Nearby Problems is employed to generate exact solutions to equations 'nearby' the steady and unsteady Burgers equation. Burgers equation is chosen because of the existence of exact solutions, and these exact solutions are discussed. Legendre polynomials are used to derive the exact solutions to the nearby problems, and the application of Legendre polynomials for both 1D and 2D problems is also discussed. Results are presented for the steady-state Burgers equation corresponding to a viscous shock wave for Reynolds numbers of 8, 16, and 512. The low Reynolds number cases are well approximated by 10th order Legendre polynomial fits, while the high Reynolds number case is not. The unsteady Burgers equation corresponding to coalescence of two viscous shock waves at a Reynolds number of 8 is also examined. Preliminary results indicate that further investigation is required to accurately capture this 2D solution.
We performed calculations to investigate the classical theories of chain branching and thermal--run--away that lead to the rapid oxidation of fuels. Mathematically, both theories infer the existence of eigenvalues with positive real parts i.e., explosive modes. We found in studies of homogeneous hydrogen--air and the methane--air mixtures that when ignition is initiated by a sufficiently high initial temperature, the transient response of the system exhibits two stages. The first stage is characterized by the existence of explosive modes. The ensuing second stage consists of fast exponential decay modes that bring the system to its equilibrium point. We demonstrated with two examples that the existence of explosive modes is not a necessary condition for the existence of a premixed flame. Homogeneous ignition calculations for mixtures with an initial concentration of radical species suggest that the diffusive transport of radical species is probably responsible for the lack of explosive modes in premixed flames.
Modal analysis of three-dimensional structures frequently involves finite element discretizations with millions of unknowns and requires computing hundreds or thousands of eigenpairs. In this presentation we review methods based on domain decomposition for such eigenspace computations in structural dynamics. We distinguish approaches that solve the eigenproblem algebraically (with minimal connections to the underlying partial differential equation) from approaches that tightly couple the eigensolver with the partial differential equation.
Low-temperature combustion concepts that utilize cooled EGR, early/retarded injection, high swirl ratios, and modest compression ratios have recently received considerable attention. To understand the combustion and, in particular, the soot formation process under these operating conditions, a modeling study was carried out using the KIVA-3V code with an improved phenomenological soot model. This multi-step soot model includes particle inception, surface growth, surface oxidation, and particle coagulation. Additional models include a piston-ring crevice model, the KH/RT spray breakup model, a droplet wall impingement model, a wall heat transfer model, and the RNG k-{var_epsilon} turbulence model. The Shell model was used to simulate the ignition process, and a laminar-and-turbulent characteristic time combustion model was used for the post-ignition combustion process. A low-load (IMEP=3 bar) operating condition was considered and the predicted in-cylinder pressures and heat release rates were compared with measurements. Predicted soot mass, soot particle size, soot number density distributions and other relevant quantities are presented and discussed. The effects of variable EGR rate (0-68%), injection pressure (600-1200 bar), and injection timing were studied. The predictions demonstrate that both EGR and retarded injection are beneficial for reducing NO{sub x} emissions, although the former has a more pronounced effect. Additionally, higher soot emissions are typically predicted for the higher EGR rates. However, when the EGR rate exceeds a critical value (over 65% in this study), the soot emissions decrease. Reduced soot emissions are also predicted when higher injection pressures or retarded injection timings are employed. The reduction in soot with retarded injection is less than what is observed experimentally, however.
Understanding and characterizing the electrical properties of multi-conductor shielded and unshielded cables is an important endeavor for many diverse applications, including airlines, land based communications, nuclear weapons, and any piece of hardware containing multi-conductor cabling. Determining the per unit length capacitance and inductance based on the geometry of the conductors, number of conductors, and characteristics of the shield can prove quite difficult. Relating the inductance and capacitance to shielding effectiveness can be even more difficult. An exceedingly large number of measurements were taken to characterize eight multi-conductor cables, of which four were 3-conductor cables and four were 18-conductor cables. Each set of four cables contained a shielded cable and an unshielded cable with the inner conductors twisted together and a shielded cable and an unshielded cable with the inner conductors not twisted together (or straight). Male LJT connectors were attached on either end of the cable and each cable had a finished length of 22.5 inches. The measurements performed were self and mutual inductance, self and mutual capacitance, and effective height. For the 18 conductor case there ended up being an 18 by 18 element matrix for inductance (with the self inductance terms lying on the diagonal) and an 18 by 18 matrix for capacitance. The effective height of each cable was measured over a frequency range from 220 MHz to 18 GHz in a Mode-Stirred Chamber. The effective height of each conductor of each cable was measured individually and all shorted together, producing 19 frequency responses for each 18 conductor cable. Shielding effectiveness was calculated using the effective heights from the shielded and unshielded cables. The results of these measurements and the statistical analysis of the data will be presented. There will also be a brief presentation of comparison with numerical models.
This report describes an approach for extending the one-dimensional turbulence (ODT) model of Kerstein [6] to treat turbulent flow in three-dimensional (3D) domains. This model, here called ODTLES, can also be viewed as a new LES model. In ODTLES, 3D aspects of the flow are captured by embedding three, mutually orthogonal, one-dimensional ODT domain arrays within a coarser 3D mesh. The ODTLES model is obtained by developing a consistent approach for dynamically coupling the different ODT line sets to each other and to the large scale processes that are resolved on the 3D mesh. The model is implemented computationally and its performance is tested and evaluated by performing simulations of decaying isotropic turbulence, a standard turbulent flow benchmarking problem.
Two-color resonant four-wave-mixing spectroscopy (TC-RFWM) is used to investigate ground-state energy transfer of hydroxyl radical in atmospheric-pressure flames. Two amplified distributed-feedback dye lasers produce 50-ps, nearly transform-limited, infrared (IR) and ultraviolet pulses. The infrared pump laser is tuned to individual rovibrational transitions of OH X {sup 2}{pi}{sub 3/2} (v{prime}=1, N{prime}) {l_arrow} X {sup 2}{pi}{sub 3/2} (v{double_prime}=0, N{double_prime}), and the ultraviolet pulse probes either the directly pumped or collisionally populated intermediate levels via A{sup 2}{Sigma}{sup +} (v*=1, N*) {l_arrow} X{sup 2}{pi}{sub 3/2}(v{prime}=1, N{prime}). By time-delaying the probe pulse with respect to the pump pulse, and appropriately constraining the polarizations of each of the four fields taking part in the wave-mixing process, we are able to independently and unambiguously measure the moments of the rotational angular momentum distribution in single rotational levels of the ground state. We present measurements of population, alignment, and orientation decay in X {sup 2}{pi}{sub 3/2} for several flame conditions. These experiments provide data necessary for the development of accurate models for diagnostic techniques using saturating laser pulses.
An alternative theory of solid mechanics, known as the peridynamic theory, formulates problems in terms of integral equations rather than partial differential equations. This theory assumes that particles in a continuum interact with each other across a finite distance, as in molecular dynamics. Damage is incorporated in the theory at the level of these two-particle interactions, so localization and fracture occur as a natural outgrowth of the equation of motion and constitutive models. A numerical method for solving dynamic problems within the peridynamic theory is described. Accuracy and numerical stability are discussed. Examples illustrate the properties of the method for modeling brittle dynamic crack growth.
Three methods that were used to measure the chemical changes associated with oxidative degradation of polymeric materials are presented. The first method is based on the nuclear activation of {sup 18}O in an elastomer that was thermally aged in an {sup 18}O{sub 2} atmosphere. Second, the alcohol groups in a thermally aged elastomer were derivatized with trifluoroacetic anhydride and their concentration measured via {sup 19}F NMR spectroscopy. Finally, a respirometer was used to directly measure the oxidative rates of a polyurethane foam as a function of aging temperature. The measurement of the oxidation rates enabled acceleration factors for oxidative degradation of these materials to be calculated.
An environmentally friendly method and materials study for desalinating inland brackish waters (i.e., coal bed methane produced waters) using a set of ion-exchange materials is presented. This desalination process effectively removes anions and cations in separate steps with minimal caustic waste generation. The anion-exchange material, hydrotalcite (HTC), exhibits an ion-exchange capacity (IEC) of {approx} 3 mequiv g{sup -1}. The cation-exchange material, an amorphous aluminosilicate permutite-like material, (Na{sub x+2y}Al{sub x}Si{sub 1-x}O{sub 2+y}), has an IEC of {approx}2.5 mequiv g{sup -1}. These ion-exchange materials were studied and optimized because of their specific ion-exchange capacity for the ions of interest and their ability to function in the temperature and pH regions necessary for cost and energy effectiveness. Room temperature, minimum pressure column studies (once-pass through) on simulant brackish water (total dissolved solids (TDS) = 2222 ppm) resulted in water containing TDS = 25 ppm. A second once-pass through column study on actual produced water (TDS = {approx}11,000) with a high carbonate concentration used an additional lime softening step and resulted in a decreased TDS of 600 ppm.
The Arsenic Water Technology Partnership program is a multi-year program funded by a congressional appropriation through the Department of Energy. The program is designed to move technologies from benchscale tests to field demonstrations. It will enable water utilities, particularly those serving small, rural communities and Indian tribes, to implement the most cost-effective solutions to their arsenic treatment needs. As part of the Arsenic Water Technology Partnership program, Sandia National Laboratories is carrying out field demonstration testing of innovative technologies that have the potential to substantially reduce the costs associated with arsenic removal from drinking water. The scope for this work includes: (1) Selection of sites and identification of technologies for pilot demonstrations; (2) Laboratory studies to develop rapid small-scale test methods; and (3) Pilot-scale studies at community sites involving side-by-side tests of innovative technologies. The goal of site selection is to identify sites that allow examination of treatment processes and systems under conditions that are relevant to different geochemical settings throughout the country. A number of candidate sites have been identified through reviews of groundwater quality databases, conference proceedings and discussions with state and local officials. These include sites in New Mexico, Arizona, Colorado, Oklahoma, Michigan, and California. Candidate technologies for the pilot tests are being reviewed through vendor forums, proof-of-principle benchscale studies managed by the American Water Works Association Research Foundation (AwwaRF) and the WERC design contest. The review considers as many potential technologies as possible and screens out unsuitable ones by considering data from past performance testing, expected costs, complexity of operation and maturity of the technology. The pilot test configurations will depend on the site-specific conditions such as access, power availability, waste disposal options and availability of permanent structures to house the test. Conducting pilot tests for media comparison at all sites in need of arsenic treatment would be extremely time consuming and costly. Laboratory studies are being conducted using rapid small-scale column tests (RSSCTs) to predict the performance of pilot-scale adsorption columns. RSSCTs are a rapid and inexpensive method of investigating innovative technologies while varying water quality and/or system design. RSSCTs are scaled down columns packed with smaller diameter adsorption media that receive higher hydraulic loading rates to significantly reduce the duration of experiments. Results for RSSCTs can be obtained in a matter of days to a few weeks, whereas pilot tests can take a number of months to over a year. In the pilot tests, the innovative technologies will be evaluated in terms of adsorptive capacity for arsenic; robustness of performance with respect to water quality parameters including pH, TDS, foulants such as Fe, Mn, silica, and organics, and other metals and radionuclides; and potentially deleterious effects on the water system such as pipe corrosion from low pH levels, fluoride removal, and generation of disinfection by-products. The new arsenic MCL will result in modification of many rural water systems that otherwise would not require treatment. Simultaneous improvement of water quality in systems that will require treatment for other contaminants such as uranium, radon and radium would be an added benefit of this program.
Military test and training ranges operate with live fire engagements to provide realism important to the maintenance of key tactical skills. Ordnance detonations during these operations typically produce minute residues of parent explosive chemical compounds. Occasional low order detonations also disperse solid phase energetic material onto the surface soil. These detonation remnants are implicated in chemical contamination impacts to groundwater on a limited set of ranges where environmental characterization projects have occurred. Key questions arise regarding how these residues and the environmental conditions (e.g., weather and geostratigraphy) contribute to groundwater pollution impacts. This report documents interim results of a mass transfer model evaluating mass transfer processes from solid phase energetics to soil pore water based on experimental work obtained earlier in this project. This mass transfer numerical model has been incorporated into the porous media simulation code T2TNT. Next year, the energetic material mass transfer model will be developed further using additional experimental data.
The use of rotational echo adiabatic passage double resonance (REAPDOR) solid-state nuclear magnetic resonance (NMR) to determine the site location of an adsorbed polar molecule in a zeolite cage is presented. Nitrogen-15 labeled ammonia is used as a probe molecule to characterize the initial adsorption site in 3A zeolite molecular sieves. The relative position of the ammonia adsorption site in the cage is determined by measuring the internuclear distance between the N on ammonia and both a Na cation site and an Al framework environment using {sup 15}N/{sup 23}Na and {sup 15}N/{sup 27}Al REAPDOR NMR experiments, respectively. The measured internuclear distances are similar to a specific ammonia adsorption site for the zeolite 4A ammonia sorption complex located using X-ray diffraction. Additional details regarding the ammonia hydrogen-bonding environment can be extracted from {sup 1}H/{sup 23}Na and {sup 1}H/{sup 27}Al REAPDOR NMR measurements.
The presentation outline of this paper is: (1) How identification of chemical hazards fits into a security risk analysis approach; (2) Techniques for target identification; and (3) Identification of chemical hazards by different organizations. The summary is: (1) There are a number of different methodologies used within the chemical industry which identify chemical hazards: (a) Some develop a manual listing of potential targets based on published lists of hazardous chemicals or chemicals of concern, 'expert opinion' or known hazards. (b) Others develop a prioritized list based on chemicals found at a facility and consequence analysis (offsite release affecting population, theft of material, product tampering). (2) Identification of chemical hazards should include not only intrinsic properties of the chemicals but also potential reactive chemical hazards and potential use for activities off-site.
Sandia National Laboratories has previously tested a capability to impose a 7.5 g-rms (30 g peak) radial vibration load up to 2 kHz on a 25 lb object with superimposed 50 g acceleration at its centrifuge facility. This was accomplished by attaching a 3,000 lb Unholtz-Dickie mechanical shaker at the end of the centrifuge arm to create a 'Vibrafuge'. However, the combination of non-radial vibration directions, and linear accelerations higher than 50g's are currently not possible because of the load capabilities of the shaker and the stresses on the internal shaker components due to the combined centrifuge acceleration. Therefore, a new technique using amplified piezo-electric actuators has been developed to surpass the limitations of the mechanical shaker system. They are lightweight, modular and would overcome several limitations presented by the current shaker. They are 'scalable', that is, adding more piezo-electric units in parallel or in series can support larger-weight test articles or displacement/frequency regimes. In addition, the units could be mounted on the centrifuge arm in various configurations to provide a variety of input directions. The design along with test results will be presented to demonstrate the capabilities and limitations of the new piezo-electric Vibrafuge.
The influence of thermal stratification on auto-ignition at constant volume and high pressure is studied by Direct Numerical Simulation (DNS) with complex H{sub 2}/air chemistry with a view to providing better understanding of combustion processes in homogeneous charge compression ignition engines. In particular the dependence of overall ignition progress on initial mixture conditions is determined. The propagation speed of ignition fronts that emanate from 'hot spots' given by a temperature spectrum is monitored by using the displacement velocity of a scalar that tracks the location of maximum heat release. The evolution of the front velocity is compared for different initial temperature distributions and the role of scalar dissipation of heat and mass is identified. It is observed that both deagrative as well as spontaneous ignition front propagation occur depending upon the local temperature gradient. It is found that the ratio of the instantaneous front speed to the deflagrative speed is a good measure of the local mode of propagation. This is verified by examining the energy and species balances. A parametric study in the amplitudes of the initial temperature fluctuation is performed and shows that this parameter has a significant influence on the observed combustion mode. Higher levels of stratification lead to more front-like structures. Predictions of the multi-zone model are presented and explained using the diagnostics developed.
The impacts of small niobium additions to processing, microstructure, and electrical properties in the Zr-rich lead zirconate titanate ceramics (PZT 95/5) were investigated. The influence of niobium content on dielectric responses and the characteristics of ferroelectric behaviors, as well as the relative phase stability and the hydrostatic pressure induced ferroelectric-to- antiferroelectric phase transformation are reported. Results indicate that increasing the niobium concentration in the solid solutions enhances densification, refines the microstructure, decreases dielectric constant and spontaneous polarization, and stabilizes the ferroelectric phase. The stabilization of ferroelectric phase with respect to the antiferroelectric phase near PZT 95/5 composition dramatically increases the pressure required for the ferroelectric-to-antiferroelectric phase transformation. These observations were correlated to the creation of A-site vacancies and a slight modification of the crystal structure. The importance of these composition-property relationships on device application will be presented.
Optical actuation of microelectromechanical systems (MEMS) is advantageous for applications for which electrical isolation is desired. Thirty-two polycrystalline silicon opto-thermal actuators, optically-powered MEMS thermal actuators, were designed, fabricated, and tested. The design of the opto-thermal actuators consists of a target for laser illumination suspended between angled legs that expand when heated, providing the displacement and force output. While the amount of displacement observed for the opto-thermal actuators was fairly uniform for the actuators, the amount of damage resulting from the laser heating ranged from essentially no damage to significant amounts of damage on the target. The likelihood of damage depended on the target design with two of the four target designs being more susceptible to damage. Failure analysis of damaged targets revealed the extent and depth of the damage.
A three-dimensional tungsten photonic crystal is thermally excited and shown to emit light at a narrow band, {lambda} = 3.3-4.25 {micro}m. The emission is experimentally observed to exceed that of the free-space Planck radiation over a wide temperature range, T = 475-850 K. it is proposed that an enhanced density of state associated with the propagating electromagnetic Bloch waves in the photonic crystal is responsible for this experimental finding.
This paper provides a brief overview of the fields of biological micro-electromechanical systems (bioMEMs) and associated nanobiotechnologies, collectively denoted as BioMicroNano. Although they are developing at a very rapid pace and still redefining themselves, several stabilized areas of research and development can be identified. Six major areas are delineated, and specific examples are discussed and illustrated. Various applications of the technologies are noted, and potential market sizes are compared.
Simplified models that are based on macroscopic force balances and droplet-geometry approximations are presented for predicting the onset of instability leading to removal of water droplets at the gas diffusion layer (GDL)/gas flow channel (GFC) interface. Visualization experiments are carried out to observe the formation, growth, and removal or instability of the water droplets at the GDL/GFC interface of a simulated polymer electrolyte fuel cell cathode. Droplet-instability diagrams or windows computed by the simplified models are compared with those measured experimentally, and good agreement is obtained. Two-dimensional flow simulations employing the finite element method coupled with an arbitrary Lagrangian-Eulerian formulation for determining the liquid/gas interface position are also performed to assess the simplified cylindrical-droplet model. Necessary conditions for preventing fully grown droplets from lodging in the flow channel are derived using the simplified models. It is found that droplet removal can be enhanced by increasing flow channel length or mean gas flow velocity, decreasing channel height or contact angle hysteresis, or making the GDL/GFC interface more hydrophobic.
Concerns over the illicit trafficking of radiological and nuclear materials were focused originally on the lack of security and accountability of such material throughout the former Soviet states. This is primarily attributed to the frequency of events that have occurred involving the theft and trafficking of critical material components that could be used to construct a Radiological Dispersal Device (RDD) or even a rudimentary nuclear device. However, with the continued expansion of nuclear technology and the deployment of a global nuclear fuel cycle these materials have become increasingly prevalent, affording a more diverse inventory of dangerous materials and dual-use items. To further complicate the matter, the list of nuclear consumers has grown to include: (1) Nation-states that have gone beyond the IAEA agreed framework and additional protocols concerning multiple nuclear fuel cycles and processes that reuse the fuel through reprocessing to exploit technologies previously confined to the more industrialized world; (2) Terrorist organizations seeking to acquire nuclear and radiological material due to the potential devastation and psychological effect of their use; (3) Organized crime, which has discovered a lucrative market in trafficking of illicit material to international actors and/or countries; and (4) Amateur smugglers trying to feed their families in a post-Soviet era. An initial look at trafficking trends of this type seems scattered and erratic, localized primarily to a select group of countries. This is not necessarily the case. The success with which other contraband has been smuggled throughout the world suggests that nuclear trafficking may be carried out with relative ease along the same routes by the same criminals or criminal organizations. Because of the inordinately high threat posed by terrorist or extremist groups acquiring the ingredients for unconventional weapons, it is necessary that illicit trafficking of these materials be better understood as to prepare for the sustained global development of the nuclear fuel cycle. Conversely, modeling and analyses of this activity must not be limited in their scope to loosely organized criminal smuggling, but address the problem as a commercial, industrial project for the covert development of nuclear technologies and unconventional weapon development.
The nation's electric power sector is highly interdependent with the economic sectors it serves; electric power needs are driven by economic activity while the economy itself depends on reliable and sustainable electric power. To advance higher level understandings of the vulnerabilities that result from these interdependencies and to identify the loss prevention and loss mitigation policies that best serve the nation, the National Infrastructure Simulation and Analysis Center is developing and using N-ABLE{trademark}, an agent-based microeconomic framework and simulation tool that models these interdependencies at the level of collections of individual economic firms. Current projects that capture components of these electric power and economic sector interdependencies illustrate some of the public policy issues that should be addressed for combined power sector reliability and national economic security.
Technologies that could quickly detect and identify virus particles would play a critical role in fighting bioterrorism and help to contain the rapid spread of disease. Of special interest is the ability to detect the presence and movement of virions without chemically modifying them by attaching molecular probes. This would be useful for rapid detection of pathogens in food or water supplies without the use of expensive chemical reagents. Such detection requires new devices to quickly screen for the presence of tiny pathogens. To develop such a device, we fabricated nanochannels to transport virus particles through ultrashort laser cavities and measured the lasing output as a sensor for virions. To understand this transduction mechanism, we also investigated light scattering from virions, both to determine the magnitude of the scattered signal and to use it to investigate the motion of virions.
The thermal interdiffusion of AlSb/GaSb multiquantum wells was measured and the intrinsic diffusivities of Al and Ga determined over a temperature range of 823-948 K for 30-9000 s. The 77-K photoluminescence (PL) was used to monitor the extent of interdiffusion through the shifts in the superlattice luminescence peaks. The chemical diffusion coefficient was quantitatively determined by fitting the observed PL peak shifts to the solution of the Schroedinger equation, using a potential derived from the solution of the diffusion equation. The value of the interdiffusion coefficient ranged from 5.2 x 10{sup -4} to 0.06 nm{sup 2}/s over the conditions studied and was characterized by an activation energy of 3.0 {+-} 0.1 eV. The intrinsic diffusion coefficients for Al and Ga were also determined with higher values for Al than for Ga, described by activation energies of 2.8 {+-} 0.4 and 1.1 {+-} 0.1 eV, respectively.
A number of industrial combustion systems are adopting oxygen-enhanced firing to improve heat transfer characteristics and reduce emissions. The exhaust gas from these systems is dominated by H2O and CO2 and therefore has substantially different gas properties from traditional combustion exhaust. In the past, laser-induced breakdown spectroscopy (LIBS) has been successfully used for the evaluation of alkali aerosol concentrations in air-based combustion systems. This paper presents results of LIBS measurements of alkali concentrations in a laboratory calibration setup and in an oxygen/natural gas container glass furnace. It shows how both gas conditions (composition and temperature) and the molecular form of the alkali species affect the LIBS signals. The paper proposes strategies for mitigating these effects in future applications of LIBS in oxygen-enhanced combustion systems.
The dimensionless extinction coefficient (K{sub e}) of soot must be known to quantify laser extinction measurements of soot concentration and to predict optical attenuation through smoke clouds. Previous investigations have measured K{sub e} for post-flame soot emitted from laminar and turbulent diffusion flames and smoking laminar premixed flames. This paper presents the first measurements of soot K{sub e} from within laminar diffusion flames, using a small extractive probe to withdraw the soot from the flame. To measure K{sub e}, two laser sources (635 nm and 1310 nm) were coupled to a transmission cell, followed by gravimetric sampling. Coannular diffusion flames of methane, ethylene and nitrogen-diluted kerosene burning in air were studied, together with slot flames of methane and ethylene. K{sub e} was measured at the radial location of maximum soot volume fraction at several heights for each flame. Results for K{sub e} at both 635 nm and 1310 nm for ethylene and kerosene coannular flames were in the range of 9-10, consistent with the results from previous studies of post-flame soot. The ethylene slot flame and the methane flames have lower K{sub e} values, in some cases as low as 2.0. These lower values of K{sub e} are found to result from the contributions of (a) the condensation of PAH species during the sampling of soot, (b) the wavelength-dependent absorptivity of soot precursor particles, and, in the case of methane, (c) the negligible contribution of soot scattering to the extinction coefficient. RDG calculations of soot scattering, in combination with the measured K{sub e} values, imply that the soot refractive index is in the vicinity of 1.75-1.03i at 635 nm.
A severe fire and explosion occurred at a propane storage yard in Truth or Consequences, N.M., when a truck ran into the pumping and plumbing system beneath a large propane tank. The storage tank emptied when the liquid-phase excess flow valve tore out of the tank. The ensuing fire engulfed several propane delivery trucks, causing one of them to explode. A series of elevated-temperature stress-rupture tears developed along the top of a 9800 L (2600 gal) truck-mounted tank as it was heated by the fire. Unstable fracture then occurred suddenly along the length of the tank and around both end caps, along the girth welds connecting the end caps to the center portion of the tank. The remaining contents of the tank were suddenly released, aerosolized, and combusted, creating a powerful boiling liquid expanding vapor explosion (BLEVE). Based on metallography of the tank pieces, the approximate tank temperature at the onset of the BLEVE was determined. Metallurgical analysis of the ruptured tank also permitted several hypotheses regarding BLEVE mechanisms to be evaluated. Suggestions are made for additional work that could provide improved predictive capabilities regarding BLEVEs and for methods to decrease the susceptibility of propane tanks to BLEVEs.
This overview is intended to provide the reader with insight into basic reliability issues often confronted when designing long-term geothermal well monitoring equipment. No single system is looked at. General examples of the long-term reliability of other industries are presented. Examples of reliability issues involving electronic components and sensors along with fiber optic sensors and cables are given. This paper will aid in building systems where a long operating life is required. However, as no introductory paper can cover all reliability issues, basic assembly practices and testing concepts are presented.