The decomposition of unconfined rigid polyurethane foam has been modeled by a kinetic bond-breaking scheme describing degradation of a primary polymer and formation of a thermally stable secondary polymer. The bond-breaking scheme is resolved using percolation theory to describe evolving polymer fragments. The polymer fragments vaporize according to individual vapor pressures. Kinetic parameters for the model were obtained from Thermal Gravimetric Analysis (TGA). The chemical structure of the foam was determined from the preparation techniques and ingredients used to synthesize the foam. Scale-up effects were investigated by simulating the response of an incident heat flux of 25 W/cm{sup 2} on a partially confined 8.8-cm diameter by 15-cm long right circular cylinder of foam which contained an encapsulated component. Predictions of center, midradial, and component temperatures, as well as regression of the foam surface, were in agreement with measurements using thermocouples and X-ray imaging.
Experience has shown that the analyses of marine transport of spent fuel in the Environmental Impact Statement (EIS) were conservative. It is anticipated that for most shipments. The external dose rate for the loaded transportation cask will be more in line with recent shipments. At the radiation levels associated with these shipments, we would not expect any personnel to exceed radiation exposure limits for the public. Package dose rates usually well below the regulatory limits and personnel work practices following ALARA principles are keeping human exposures to minimal levels. However, the potential for Mure shipments with external dose rates closer to the exclusive-use regulatory limit suggests that DOE should continue to provide a means to assure that individual crew members do not receive doses in excess of the public dose limits. As a minimum, the program will monitor cask dose rates and continue to implement administrative procedures that will maintain records of the dose rates associated with each shipment, the vessel used, and the crew list for the vessel. DOE will continue to include a clause in the contract for shipment of the foreign research reactor spent nuclear fuel requiring that the Mitigation Action Plan be followed.
As model validation techniques gain more acceptance and increase in power, the demands on the modal parameter extractions increase. The estimation accuracy, the number of modes desired, and the data reduction efficiency are required features. An algorithm known as SMAC (Synthesize Modes And Correlate), based on principles of modal filtering, has been in development for a few years. SMAC has now been extended in two main areas. First, it has now been automated. Second, it has been extended to fit complex modes as well as real modes. These extensions have enhanced the power of modal extraction so that, typically, the analyst needs to manually fit only 10 percent of the modes in the desired bandwidth, whereas the automated routines will fit 90 percent of the modes. SMAC could be successfully automated because it generally does not produce computational roots.
Direct metal deposition technologies produce complex, near net shape components from Computer Aided Design (CAD) solid models. Most of these techniques fabricate a component by melting powder in a laser weld pool, rastering the weld bead to form a layer, and additively constructing subsequent layers. This report will describe anew direct metal deposition process, known as WireFeed, whereby a small diameter wire is used instead of powder as the feed material to fabricate components. Currently, parts are being fabricated from stainless steel alloys. Microscopy studies show the WireFeed parts to be filly dense with fine microstructural features. Mechanical tests show stainless steel parts to have high strength values with retained ductility. A model was developed to simulate the microstructural evolution and coarsening during the WireFeed process. Simulations demonstrate the importance of knowing the temperature distribution during fabrication of a WireFeed part. The temperature distribution influences microstructural evolution and, therefore, must be controlled to tailor the microstructure for optimal performance.
Development of bug-free, high-surety, complex software is quite difficult using current tools. The brittle nature of the programming constructs in popular languages such as C/C++ is one root cause. Brittle commands force the designer to rigidly specify the minutiae of tasks, e.g. using ''for(index=0;index>total;index++)'', rather than specifying the goals or intent of the tasks, e.g. ''ensure that all relevant data elements have been examined''. Specification of task minutiae makes code hard to comprehend, which in turn encourages design errors/limitations and makes future modifications quite difficult. This report describes an LDRD project to seed the development of a surety computer language, for stand-alone computing environments, to be implemented using the swarm intelligence of autonomous agents. The long term vision of this project was to develop a language with the following surety capabilities: (1) Reliability -- Autonomous agents can appropriate y decide when to act and when a task is complete, provide a natural means for avoiding brittle task specifications, and can overcome many hardware glitches. (2) Safety, security -- Watchdog safety and security agents can monitor other agents to prevent unauthorized or dangerous actions. (3) An immune system -- The small chunks of agent code can have an encryption scheme to enable detection and elimination of unauthorized and corrupted agents. This report describes the progress achieved during this small 9 month project and describes lessons learned.
Current supercomputers use large parallel arrays of tightly coupled processors to achieve levels of performance far surpassing conventional vector supercomputers. Shock-wave physics codes have been developed for these new supercomputers at Sandia National Laboratories and elsewhere. These parallel codes run fast enough on many simulations to consider using them to study the effects of varying design parameters on the performance of models of conventional munitions and other complex systems. Such studies maybe directed by optimization software to improve the performance of the modeled system. Using a shaped-charge jet design as an archetypal test case and the CTH parallel shock-wave physics code controlled by the Dakota optimization software, we explored the use of automatic optimization tools to optimize the design for conventional munitions. We used a scheme in which a lower resolution computational mesh was used to identify candidate optimal solutions and then these were verified using a higher resolution mesh. We identified three optimal solutions for the model and a region of the design domain where the jet tip speed is nearly optimal, indicating the possibility of a robust design. Based on this study we identified some of the difficulties in using high-fidelity models with optimization software to develop improved designs. These include developing robust algorithms for the objective function and constraints and mitigating the effects of numerical noise in them. We conclude that optimization software running high-fidelity models of physical systems using parallel shock wave physics codes to find improved designs can be a valuable tool for designers. While current state of algorithm and software development does not permit routine, ''black box'' optimization of designs, the effort involved in using the existing tools may well be worth the improvement achieved in designs.
This report documents the use of the FITS routine, which provides automated fits of various analytical, commonly used probability models from input data. It is intended to complement the previously distributed FITTING routine documented in RMS Report 14 (Winterstein et al., 1994), which implements relatively complex four-moment distribution models whose parameters are fit with numerical optimization routines. Although these four-moment fits can be quite useful and faithful to the observed data, their complexity can make them difficult to automate within standard fitting algorithms. In contrast, FITS provides more robust (lower moment) fits of simpler, more conventional distribution forms. For each database of interest, the routine estimates the distribution of annual maximum response based on the data values and the duration, T, over which they were recorded. To focus on the upper tails of interest, the user can also supply an arbitrary lower-bound threshold, {chi}{sub low}, above which a shifted distribution model--exponential or Weibull--is fit.
Path integral Monte Carlo simulations and calculations were performed on molecular hydrogen liquids. The equation-of-state, internal energies, and vapor liquid phase diagrams from simulation were found to be in quantitative agreement with experiments. Analytical calculations were performed on,H2 liquids using integral equation methods to study the degree of localization of the hydrogen molecules. Very little self-trapping or localization was found as a function of temperature and density. Good qualitative agreement was found between the integral equation calculations and the quantum Monte Carlo simulations for the radius of gyration of the hydrogen molecules. Path integral simulations were also performed on molecular hydrogen on graphite surfaces, slit pores, and in carbon nanotubes. Significant quantum effects on the adsorption of hydrogen were observed.
Thermionic energy conversion in a microminiature format shows potential as a viable, high efficiency, on-chip power source. Microminiature thermionic converters (MTC) with inter-electrode spacings on the order of microns are currently being prototyped and evaluated at Sandia. The remaining enabling technology is the development of low work function materials and processes than can be integrated into these converters. In this report, the authors demonstrate a method of incorporating thin film emitters into converters using rf sputtering. They find that the resultant films possess a minimum work function of 1.2 eV. Practical energy conversion is hindered by surface work function non-uniformity. They postulate the source of this heterogeneity to be a result of limited bulk and surface transport of barium. Several methods are proposed for maximizing transport, including increased film porosity and the use of metal terminating layers. They demonstrate a novel method for incorporating film porosity based on metal interlayer coalescence.
Monolithic, integrated acoustic wave chemical microsensors are being developed on gallium arsenide (GaAs) substrates. With this approach, arrays of microsensors and the high frequency electronic components needed to operate them reside on a single substrate, increasing the range of detectable analytes, reducing overall system size, minimizing systematic errors, and simplifying assembly and packaging. GaAs is employed because it is both piezoelectric, a property required to produce the acoustic wave devices, and a semiconductor with a mature microelectronics fabrication technology. Many aspects of integrated GaAs chemical sensors have been investigated, including: surface acoustic wave (SAW) sensors; monolithic SAW delay line oscillators; GaAs application specific integrated circuits (ASIC) for sensor operation; a hybrid sensor array utilizing these ASICS; and the fully monolithic, integrated SAW array. Details of the design, fabrication, and performance of these devices are discussed. In addition, the ability to produce heteroepitaxial layers of GaAs and aluminum gallium arsenide (AlGaAs) makes possible micromachined membrane sensors with improved sensitivity compared to conventional SAW sensors. Micromachining techniques for fabricating flexural plate wave (FPW) and thickness shear mode (TSM) microsensors on thin GaAs membranes are presented and GaAs FPW delay line and TSM resonator performance is described.
Inductive energy storage systems can have high energy density, lending to smaller, less expensive systems. The crucial element of an inductive energy storage system is the opening switch. This switch must conduct current while energy is stored in an inductor, then open quickly to transfer this energy to a load. Plasma can perform this function. The Plasma Opening Switch (POS) has been studied for more than two decades. Success with the conventional plasma opening switch has been limited. A system designed to significantly improve the performance of vacuum opening switches is described in this paper. The gap cleared of plasma is a rough figure-of-merit for vacuum opening switches. Typical opened gaps of 3 mm are reported for conventional switches. The goal for the system described in this paper is more than 3 cm. To achieve this, the command-triggered POS adds an active opening mechanism, which allows complete separation of conduction and opening. This separation is advantageous because of the widely different time scales of conduction and opening. The detrimental process of magnetic field penetration into the plasma during conduction is less important in this switch. The opening mechanism duration is much shorter than the conduction time, so penetration during opening is insignificant. Opening is accomplished with a fast magnetic field that pushes plasma out of the switch region. Plasma must be removed from the switch region to allow high voltage. This paper describes some processes important during conduction and opening, and show calculations on the trigger requirements. The design of the switch is shown. This system is designed to demonstrate both improved performance and nanosecond output jitter at levels greater than one terawatt. An amplification mechanism is described which reduces the trigger energy. Particle-in-cell simulations of the system are also shown.
Simple tangent, hard site chains near a hard wall are modeled with a Density Functional (DF) theory that uses the direct correlation function, c(r), as its ''input''. Two aspects of this DF theory are focused upon: (1) the consequences of variations in c(r)'s detailed form; and (2) the correct way to introduce c(r) into the DF formalism. The most important aspect of c(r) is found to be its integrated value, {cflx c}(0). Indeed, it appears that, for fixed {cflx c}(0), all reasonable guesses of the detailed shape of c(r) result in surprisingly similar density distributions, {rho}(r). Of course, the more accurate the c(r), the better the {rho}(r). As long as the length scale introduced by c(r) is roughly the hard site diameter and as long as the solution remains liquid-like, the {rho}(r) is found to be in good agreement with simulation results. The c(r) is used in DF theory to calculate the medium-induced-potential, U{sub M}(r) from the density distribution, {rho}(r). The form of U{sub M}(r) can be chosen to be one of a number of different forms. It is found that the forms for U{sub M}(r), which yield the most accurate results for the wall problem, are also those which were suggested as accurate in previous, related studies.
Solid freeform fabrication is the near-net-shape manufacturing of components by sequentially stacking thin layers of material until complicated three dimensional shapes are produced. The operation is computer controlled and requires no molds. This exciting new field of technology provides engineers with the ability to rapidly produce prototype parts directly from CAD drawings and oftentimes little or no machining is necessary after fabrication. Techniques for freeform fabrication with several types of plastics and metals are already quite advanced and maybe reviewed in references 1 and 2. Very complicated plastic models can be fabricated by stereolithography, selective laser sintering, fused deposition modeling, or three-dimensional ink jet printing. Metals may be freeformed by the LENS{trademark} technique and porous ceramic bodies by three dimensional printing into a porous powder bed. However, methods for freeform fabrication that utilize particulate slurries to build dense ceramics and composites are not as well developed. The techniques that are being developed for the freeform fabrication of dense structural ceramics primarily revolve around the sequential layering of ceramic loaded polymers or waxes. Laminated Object Manufacturing and CAM-LEM processing use controlled stacking and laser cutting of ceramic tapes [2,3]. Similar to fused deposition modeling, ceramic loaded polymer/wax filaments are being used for the fused deposition of ceramics [2,4]. Extrusion freeform fabrication uses high pressure extrusion to deposit layers of ceramic loaded polymer/wax systems[1]. Modified stereolithographic techniques are also being developed using ceramic loaded ultraviolet curable resins [2]. Pre-sintered parts made with any of these techniques typically have 40-55 vol.% polymeric binder. In this regard, these techniques are analogous to powder injection molding of ceramics. Very long and complicated burnout heat treatments are necessary to produce a dense ceramic, free of organics. Heating rates of 0.2 degrees Celsius per minute are common. [5] Thus, while a part maybe rapidly prototype within a few hours, it takes several days to densify. In contrast, robocasting is a freeform fabrication technique developed at Sandia National Labs that utilizes particulate slurries but does not require organic binders. Since binder burnout is not an issue, a dense ceramic part maybe freeformed, dried, and sintered in less than 24 hours. In some regards, robocasting is analogous to the ceramic near-net-shape processing techniques, slip casting and gel casting [6]; however, robocasting is moldless and fabrication times can be quicker.
Device penetration into media such as metal and soil is an application of some engineering interest. Often, these devices contain internal components and it is of paramount importance that all significant components survive the severe environment that accompanies the penetration event. In addition, the system must be robust to perturbations in its operating environment, some of which exhibit behavior which can only be quantified to within some level of uncertainty. In the analysis discussed herein, methods to address the reliability of internal components for a specific application system are discussed. The shock response spectrum (SRS) is utilized in conjunction with the Advanced Mean Value (AMV) and Response Surface methods to make probabilistic statements regarding the predicted reliability of internal components. Monte Carlo simulation methods are also explored.
Previous applications of DF theory required a single chain Monte Carlo simulation to be performed within a self-consistent loop. In the current work, a methodology is developed which permits the simulation to be taken out of the iterative loop. Consequently, the calculation of the self-consistent, medium-induced-potential, or field, is decoupled from the simulation. This approach permits different densities, different forms of U{sub M}(r), and different wall-polymer interactions to be investigated from a single Monte Carlo simulation. The increase in computational efficiency is immense.
In the world of computers a trusted object is a collection of possibly-sensitive data and programs that can be allowed to reside and execute on a computer, even on an adversary's machine. Beyond the scope of one computer we believe that network-based agents in high-consequence and highly reliable applications will depend on this approach, and that the basis for such objects is what we call ''faithful execution.''
Minimum-time trajectory tracking of an under-actuated mechanical system called the Acrobot is presented. The success of the controller is demonstrated by the fact that the tracking error is reduced by more than an order of magnitude when compared to the open-loop system response. The control law is obtained by linearizing the system about the nominal trajectory and applying differential dynamic programming to the resulting linear time-varying system, while using a weighted sum of the state-deviation and input-deviation as the cost function.
A study to characterize the low-temperature reactive processes for o-AP and an AP/HTPB-based propellant (class 1.3) is being conducted in the laboratory using the techniques of simultaneous thermogravimetric modulated beam mass spectrometry (STMBMS) and scanning electron microscopy (SEM). The results presented in this paper are a follow up of the previous work that showed the overall decomposition to be complex and controlled by both physical and chemical processes. The decomposition is characterized by the occurrence of one major event that consumes up to {approx}35% of the AP, depending upon particle size, and leaves behind a porous agglomerate of AP. The major gaseous products released during this event include H{sub 2}O, O{sub 2}, Cl{sub 2}, N{sub 2}O and HCl. The recent efforts provide further insight into the decomposition processes for o-AP. The temporal behaviors of the gas formation rates (GFRs) for the products indicate that the major decomposition event consists of three chemical channels. The first and third channels are affected by the pressure in the reaction cell and occur at the surface or in the gas phase above the surface of the AP particles. The second channel is not affected by pressure and accounts for the solid-phase reactions characteristic of o-AP. The third channel involves the interactions of the decomposition products with the surface of the AP. SEM images of partially decomposed o-AP provide insight to how the morphology changes as the decomposition progresses. A conceptual model has been developed, based upon the STMBMS and SEM results, that provides a basic description of the processes. The thermal decomposition characteristics of the propellant are evaluated from the identities of the products and the temporal behaviors of their GFRs. First, the volatile components in the propellant evolve from the propellant as it is heated. Second, the hot AP (and HClO{sub 4}) at the AP-binder interface oxidize the binder through reactions that preferentially strip the hydrogen from the binder and yield HCl and H{sub 2}O with some oxidation of the carbon from the binder. Third, the o-AP in the propellant decomposes in the same manner as in neat o-AP. Finally, AP-derived gaseous products interact with other ingredients in the propellant.
A set of linear and nonlinear stability analysis tools have been developed to analyze steady state incompressible flows in 3D geometries. The algorithms have been implemented to be scalable to hundreds of parallel processors. The linear stability of steady state flows are determined by calculating the rightmost eigenvalues of the associated generalize eigenvalue problem. Nonlinear stability is studied by bifurcation analysis techniques. The boundaries between desirable and undesirable operating conditions are determined for buoyant flow in the rotating disk CVD reactor.
Mesa and planar GaN Schottky diode rectifiers with reverse breakdown voltages (V{sub RB}) up to 550V and >2000V, respectively, have been fabricated. The on-state resistance, R{sub ON}, was 6m{Omega}{center_dot} cm{sup 2} and 0.8{Omega}cm{sup 2}, respectively, producing figure-of-merit values for (V{sub RB}){sup 2}/R{sub ON} in the range 5-48 MW{center_dot}cm{sup -2}. At low biases the reverse leakage current was proportional to the size of the rectifying contact perimeter, while at high biases the current was proportional to the area of this contact. These results suggest that at low reverse biases, the leakage is dominated by the surface component, while at higher biases the bulk component dominates. On-state voltages were 3.5V for the 550V diodes and {ge}15 for the 2kV diodes. Reverse recovery times were <0.2{micro}sec for devices switched from a forward current density of {approx}500A{center_dot}cm{sup -2} to a reverse bias of 100V.
The topic of Privacy is complex, multi-faceted, and often emotionally laden. This paper will cover the following topics, in an effort to further understanding of federal regulations and activities, the balancing act that necessarily occurs in business, and what role a records manager can play. The topics are: Definitions; The Privacy Act; ''Private'' companies; Potential areas of concern; Expectations; Corporate responsibilities; Case studies; and Records Manager's role.
The non-linear stress-strain relation for crosslinked polymer networks is studied using molecular dynamics simulations. Previously we demonstrated the importance of trapped entanglements in determining the elastic and relaxational properties of networks. Here we present new results for the stress versus strain for both dry and swollen networks. Models which limit the fluctuations of the network strands like the tube model are shown to describe the stress for both elongation and compression. For swollen networks, the total modulus is found to decrease like (V{sub o}/V){sup 2/3} and goes to the phantom model result only for short strand networks.
Predicting the properties of nonequilibrium systems from molecular simulations is a growing area of interest. One important class of problems involves steady state diffusion. To study these cases, a grand canonical molecular dynamics approach has been developed by Heffelfinger and van Swol [J. Chem. Phys., 101, 5274 (1994)]. With this method, the flux of particles, the chemical potential gradients, and density gradients can all be measured in the simulation. In this paper, we present a complementary approach that couples a nonlocal density functional theory (DFT) with a transport equation describing steady-state flux of the particles. We compare transport-DFT predictions to GCMD results for a variety of ideal (color diffusion), and nonideal (uphill diffusion and convective transport) systems. In all cases excellent agreement between transport-DFT and GCMD calculations is obtained with diffusion coefficients that are invariant with respect to density and external fields.
Magnesium vanadates are potentially important catalytic materials for the conversion of alkanes to alkenes via oxidative dehydrogenation. However, little is known about the active sites at which the catalytic reactions take place. It may be possible to obtain a significant increase in the catalytic efficiency if the effects of certain material properties on the surface reactions could be quantified and optimized through the use of appropriate preparation techniques. Given that surface reactivity is often dependent upon surface structure and that the atomic level structure of the active sites in these catalysts is virtually unknown, we desire thin film samples consisting of a single magnesium vanadate phase and a well defined crystallographic orientation in order to reduce complexity and simplify the study of active sites. We report on the use of reactive RF sputter deposition to fabricate very highly oriented, stoichiometric Mg{sub 3}(VO{sub 4}){sub 2} thin films for use in these surface analysis studies. Deposition of samples onto amorphous substrates resulted in very poor crystallinity. However, deposition of Mg{sub 3}(VO{sub 4}){sub 2} onto well-oriented, lattice-matched thin film ''seed'' layers such as Ti(0001), Au(111), or Pt(111) resulted in very strong preferential (042) crystallographic orientation (pseudo-hexagonal oxygen planes parallel to the substrate). This strong preferential growth of the Mg{sub 3}VO{sub 4}{sub 2} suggests epitaxial (single-crystal) growth of this mixed metal oxide on the underlying metal seed layer. The effects of the seed layer material, deposition temperature, and post-deposition reactive treatments on thin film properties such as stoichiometry, crystallographic orientation, and chemical interactions will be discussed.