In this paper, we discuss several specific threats directed at the routing data of an ad hoc network. We address security issues that arise from wrapping authentication mechanisms around ad hoc routing data. We show that this bolt-on approach to security may make certain attacks more difficult, but still leaves the network routing data vulnerable. We also show that under a certain adversarial model, most existing routing protocols cannot be secured with the aid of digital signatures.
The recent unprecedented growth of global network (Internet) usage has created an ever-increasing amount of congestion. Telecommunication companies (Telco) and Internet Service Providers (ISP's), which provide access and distribution through the network, are increasingly more aware of the need to manage this growth. Congestion, if left unmanaged, will result in a degradation of the over-all network. These access and distribution networks currently lack formal mechanisms to select Quality of Service (QoS) attributes for data transport. Network services with a requirement for expediency or consistent amounts of bandwidth cannot function properly in a communication environment without the implementation of a QoS structure. This report describes and implements such a structure that results in the ability to identify, prioritize, and police critical application flows.
We investigate the rate of convergence of stochastic basis elements to the solution of a stochastic operator equation. As in deterministic finite elements, the solution may be approximately represented as the linear combination of basis elements. In the stochastic case, however, the solution belongs to a Hilbert space of functions defined on a cross product domain endowed with the product of a deterministic and probabilistic measure. We show that if the dimension of the stochastic space is n, and the desired accuracy is of order {var_epsilon}, the number of stochastic elements required to achieve this level of precision, in the Galerkin method, is on the order of | ln {var_epsilon} |{sup n}.
The GEO-SEQ Project is investigating methods for geological sequestration of CO{sub 2}. This project, which is directed by LBNL and includes a number of other industrial, university, and national laboratory partners, is evaluating computer simulation methods including TOUGH2 for this problem. The TOUGH2 code, which is a widely used code for flow and transport in porous and fractured media, includes simplified methods for gas diffusion based on a direct application of Fick's law. As shown by Webb (1998) and others, the Dusty Gas Model (DGM) is better than Fick's Law for modeling gas-phase diffusion in porous media. In order to improve gas-phase diffusion modeling for the GEO-SEQ Project, the EOS7R module in the TOUGH2 code has been modified to include the Dusty Gas Model as documented in this report. In addition, the liquid diffusion model has been changed from a mass-based formulation to a mole-based model. Modifications for separate and coupled diffusion in the gas and liquid phases have also been completed. The results from the DGM are compared to the Fick's law behavior for TCE and PCE diffusion across a capillary fringe. The differences are small due to the relatively high permeability (k = 10{sup -11} m{sup 2}) of the problem and the small mole fraction of the gases. Additional comparisons for lower permeabilities and higher mole fractions may be useful.
Measurement and signal intelligence demands has created new requirements for information management and interoperability as they affect surveillance and situational awareness. Integration of on-board autonomous learning and adaptive control structures within a remote sensing platform architecture would substantially improve the utility of intelligence collection by facilitating real-time optimization of measurement parameters for variable field conditions. A problem faced by conventional digital implementations of intelligent systems is the conflict between a distributed parallel structure on a sequential serial interface functionally degrading bandwidth and response time. In contrast, optically designed networks exhibit the massive parallelism and interconnect density needed to perform complex cognitive functions within a dynamic asynchronous environment. Recently, all-optical self-organizing neural networks exhibiting emergent collective behavior which mimic perception, recognition, association, and contemplative learning have been realized using photorefractive holography in combination with sensory systems for feature maps, threshold decomposition, image enhancement, and nonlinear matched filters. Such hybrid information processors depart from the classical computational paradigm based on analytic rules-based algorithms and instead utilize unsupervised generalization and perceptron-like exploratory or improvisational behaviors to evolve toward optimized solutions. These systems are robust to instrumental systematics or corrupting noise and can enrich knowledge structures by allowing competition between multiple hypotheses. This property enables them to rapidly adapt or self-compensate for dynamic or imprecise conditions which would be unstable using conventional linear control models. By incorporating an intelligent optical neuroprocessor in the back plane of an imaging sensor, a broad class of high-level cognitive image analysis problems including geometric change detection, pattern recognition, and correlated feature extraction can be realized in an inherently parallel fashion without information bottlenecking or external supervision, Using this approach, we believe that autonomous control systems embodied with basic adaptive decision-theoretic capabilities can be developed for imaging and surveillance sensors to improve discrimination in stressing operational environments.
Conventional systems surety analysis is basically restricted to measurable or physical-model-derived data. However, most analyses, including high-consequence system surety analysis, must also utilize subjective information. In order to address this need, there has been considerable effort on analytically incorporating engineering judgment. For example, Dempster-Shafer theory establishes a framework in which frequentist probability and Bayesian incorporation of new data are subsets. Although Bayesian and Dempster-Shafer methodology both allow judgment, neither derives results that can indicate the relative amounts of subjective judgment and measurable data in the results. The methodology described in this report addresses these problems through a hybrid-mathematics-based process that allows tracking of the degree of subjective information in the output, thereby providing more informative (as well as more appropriate) results. In addition, most high consequence systems offer difficult-to-analyze situations. For example, in the Sandia National Laboratories nuclear weapons program, the probability that a weapon responds safely when exposed to an abnormal environment (e.g., lightning, crush, metal-melting temperatures) must be assured to meet a specific requirement. There are also non-probabilistic DOE and DoD requirements (e.g., for determining the adequacy of positive measures). The type of processing required for these and similar situations transcends conventional probabilistic and human factors methodology. The results described herein address these situations by efficiently utilizing subjective and objective information in a hybrid mathematical structure in order to directly apply to the surety assessment of high consequence systems. The results can also improve the quality of the information currently provided to decision-makers. To this end, objective inputs are processed in a conventional manner; while subjective inputs are derived from the combined engineering judgment of experts in the appropriate disciplines. In addition to providing output constituents (including portrayal of uncertainty) corresponding to combination of these input types, their individual contributions to the resultant uncertainty are determined and provided as part of the output information. Finally, the safety assessment is complemented by a latent effects analysis, facilitated by soft-aggregation accumulation of observed operational constituents.
An effort is underway at Sandia National Laboratories to develop a library of algorithms to search for potential interactions between surfaces represented by analytic and discretized topological entities. This effort is also developing algorithms to determine forces due to these interactions for transient dynamics applications. This document describes the Application Programming Interface (API) for the ACME (Algorithms for Contact in a Multiphysics Environment) library.
We study the ballistic and diffusive magnetoquantum transport using a typical quantum point contact geometry for single and tunnel-coupled double wires that are wide (less than or similar to1 mum) in one perpendicular direction with densely populated sublevels and extremely confined in the other perpendicular (i.e., growth) direction. A general analytic solution to the Boltzmann equation is presented for multisublevel elastic scattering at low temperatures. The solution is employed to study interesting magnetic-field dependent behavior of the conductance such as a large enhancement and quantum oscillations of the conductance for various structures and field orientations. These phenomena originate from the following field-induced properties: magnetic confinement, displacement of the initial- and final-state wave functions for scattering, variation of the Fermi velocities, mass enhancement, depopulation of the sublevels and anticrossing (in double quantum wires). The magnetoconductance is strikingly different in long diffusive (or rough. dirty) wires from the quantized conductance in short ballistic (or clean) wires. Numerical results obtained for the rectangular confinement potentials in the growth direction are satisfactorily interpreted in terms of the analytic solutions based on harmonic confinement potentials. Some of the predicted features of the field-dependent diffusive and quantized conductances are consistent with recent data from GaAs/AlxGa1-xAs double quantum wires.
Reduced prestressing and degradation of prestressing tendons in concrete containment vessels were investigated using finite element analysis of a typical prestressed containment vessel. The containment was analyzed during a loss of coolant accident (LOCA) with varying levels of prestress loss and with reduced tendon area. It was found that when selected hoop prestressing tendons were completely removed (as if broken) or when the area of selected hoop tendons was reduced, there was a significant impact on the ultimate capacity of the containment vessel. However, when selected hoop prestressing tendons remained, but with complete loss of prestressing, the predicted ultimate capacity was not significantly affected for this specific loss of coolant accident. Concrete cracking occurred at much lower levels for all cases. For cases where selected vertical tendons were analyzed with reduced prestressing or degradation of the tendons, there also was not a significant impact on the ultimate load carrying capacity for the specific accident analyzed. For other loading scenarios (such as seismic loading) the loss of hoop prestressing with the tendons remaining could be more significant on the ultimate capacity of the containment vessel than found for the accident analyzed. A combination of loss of prestressing and degradation of the vertical tendons could also be more critical during other loading scenarios.
The Sandia coilgun [1,2,3,4,5] is an inductive electromagnetic launcher. It consists of a sequence of powered, multi-turn coils surrounding a flyway of circular cross-section through which a conducting armature passes. When the armature is properly positioned with respect to a coil, a charged capacitor is switched into the coil circuit. The rising coil currents induce a current in the armature, producing a repulsive accelerating force. The basic numerical tool for modeling the coilgun is the SLINGSHOT code, an expanded, user-friendly successor to WARP-10 [6]. SLINGSHOT computes the currents in the coils and armature, finds the forces produced by those currents, and moves the armature through the array of coils. In this approach, the cylindrically symmetric coils and armature are subdivided into concentric hoops with rectangular cross-section, in each of which the current is assumed to be uniform. The ensemble of hoops are treated as coupled circuits. The specific heats and resistivities of the hoops are found as functions of temperature and used to determine the resistive heating. The code calculates the resistances and inductances for all hoops, and the mutual inductances for all hoop pairs. Using these, it computes the hoop currents from their circuit equations, finds the forces from the products of these currents and the mutual inductance gradient, and moves the armature. Treating the problem as a set of coupled circuits is a fast and accurate approach compared to solving the field equations. Its use, however, is restricted to problems in which the symmetry dictates the current paths. This paper is divided into three parts. The first presents a demonstration of the code. The second describes the input and output. The third part describes the physical models and numerical methods used in the code. It is assumed that the reader is familiar with coilguns.
The methodology in this report addresses the safety effects of organizational and operational factors that can be measured through ''inspection.'' The investigation grew out of a preponderance of evidence that the safety ''culture'' (attitude of employees and management toward safety) was frequently one of the major root causes behind accidents or safety-relevant failures. The approach is called ''Markov latent effects'' analysis. Since safety also depends on a multitude of factors that are best measured through well known risk analysis methods (e.g., fault trees, event trees, FMECA, physical response modeling, etc.), the Markov latent effects approach supplements conventional safety assessment and decision analysis methods. A top-down mathematical approach is developed for decomposing systems, for determining the most appropriate items to be measured, and for expressing the measurements as imprecise subjective metrics through possibilistic or fuzzy numbers. A mathematical model is developed that facilitates combining (aggregating) inputs into overall metrics and decision aids, also portraying the inherent uncertainty. A major goal of the modeling is to help convey the top-down system perspective. Metrics are weighted according to significance of the attribute with respect to subsystems and are aggregated nonlinearly. Since the accumulating effect responds less and less to additional contribution, it is termed ''soft'' mathematical aggregation, which is analogous to how humans frequently make decisions. Dependence among the contributing factors is accounted for by incorporating subjective metrics on commonality and by reducing the overall contribution of these combinations to the overall aggregation. Decisions derived from the results are facilitated in several ways. First, information is provided on input ''Importance'' and ''Sensitivity'' (both Primary and Secondary) in order to know where to place emphasis on investigation of root causes and in considering new controls that may be necessary. Second, trends in inputs and outputs are tracked in order to obtain significant information, including cyclic information, for the decision process. Third, Early Alerts are provided in order to facilitate pre-emptive action. Fourth, the outputs are compared to soft thresholds provided by sigmoid functions. The methodology has been implemented in a software tool.
The unique properties of carbon have made it both a fascinating and an important subject of experimental and theoretical studies for many years [1]-[4]. The contrast between its best-known elemental forms, graphite and diamond, is particularly striking. Graphite is black, has a rather low density and high compressibility (close to that of magnesium), and is greasy enough to be useful as a lubricant and in pencil leads. Diamond is brilliantly translucent, 60% more dense than graphite, less compressible than either tungsten or corundum, and its hardness makes it useful for polishing and cutting. This variability in properties, as well as that observed among the many classes of carbon compounds, arises because of profound differences in electronic structure of the carbon bonds [5]. A number of other solid forms of carbon are known. Pyrolytic graphite [6] is a polycrystalline material in which the individual crystallites have a structure quite similar to that of natural graphite. Fullerite (solid C 60), discovered only ten years ago [7], consists of giant molecules in which the atoms are arranged into pentagons and hexagons on the surface of a spherical cage. Amorphous carbon [8][9], including carbon black and ordinary soot, is a disordered form of graphite in which the hexagonally bonded layers are randomly oriented. Glassy carbons [9][10], on the other hand, have more random structures. Many other structures have been discussed [1][9].
This report describes least squares solution methods and linearized estimates of solution errors caused by data errors. These methods are applied to event locating systems which use time-of-arrival (TOA) data. Analyses are presented for algorithms that use the TOA data in a ''direct'' manner and for algorithms utilizing Time-of-arrival Squared (TSQ) methods. Location and error estimation results were applied to a ''typical'' satellite TOA detecting system. Using Monte Carlo methods, it was found that the linearized location error estimates were valid for random data errors with relatively large variances and relatively poor event/sensor geometries. In addition to least squares methods, which use an L{sub 2} norm, methods were described for L{sub 1} and L{sub {infinity}} norms. In general, these latter norms offered little improvement over least squares methods. Reduction of the location error variances can be effected by using information in addition to the TOA data themselves by adding judiciously chosen ''conditioning'' equation(s) to the least squares system. However, the added information can adversely affect the mean errors. Also, conditioned systems may offer location solutions where nonconditioned scenarios may not be solvable. Solution methods and linearized error estimates are given for ''conditioned'' systems. It was found that for significant data errors, the linearized estimates were also close to the Monte Carlo results.
Aluminum oxide (ALOX) filled epoxy is the dielectric encapsulant in shock driven high-voltage power supplies. ALOX encapsulants display a high dielectric strength under purely electrical stress, but minimal information is available on the combined effects of high voltage and mechanical shock. We report breakdown results from applying electrical stress in the form of a unipolar high-voltage pulse of the order of 10-{micro}s duration, and our findings may establish a basis for understanding the results from proposed combined-stress experiments. A test specimen geometry giving approximately uniform fields is used to compare three ALOX encapsulant formulations, which include the new-baseline 459 epoxy resin encapsulant and a variant in which the Alcoa T-64 alumina filler is replaced with Sumitomo AA-10 alumina. None of these encapsulants show a sensitivity to ionizing radiation. We also report results from specimens with sharp-edged electrodes that cause strong, localized field enhancement as might be present near electrically-discharged mechanical fractures in an encapsulant. Under these conditions the 459-epoxy ALOX encapsulant displays approximately 40% lower dielectric strength than the older Z-cured Epon 828 formulation. An investigation of several processing variables did not reveal an explanation for this reduced performance. The 459-epoxy encapsulant appears to suffer electrical breakdown if the peak field anywhere reaches a critical level. The stress-strain characteristics of Z-cured ALOX encapsulant are measured under high triaxial pressure and we find that this stress causes permanent deformation and a network of microscopic fractures. Recommendations are made for future experimental work.
Hafnium diboride-silicon carbide (HS) and zirconium diboride-silicon carbide (ZS) composites are potential materials for high temperature, thermal shock applications such as for components on re-entry vehicles. In order to establish material constants necessary for evaluation of in situ fracture, bars fractured in four-point flexure were examined using fractographic principles. The fracture toughness was determined from measurements of the critical crack sizes and the strength values and the crack branching constants were established to use in forensic fractography for future in-flight tests. The fracture toughnesses range from about 13 MPam{sup 1/2} at room temperature to about 6 MPam{sup 1/2} at 1400 C for ZrB{sub 2}-Sic composites and from about 13 MPam{sup 1/2} at room temperature to about 4 MPam{sup 1/2} at 1400 C for HfB{sub 2}-SiC composites. Thus, the toughnesses of either the HS or ZS composites have the potential for use in thermal shock applications. Processing and manufacturing defects limited the strength of the test bars. However, examination of the microstructure on the fracture surfaces shows that the processing of these composites can be improved. There is potential for high toughness composites with high strength to be used in thermal shock conditions if the processing and handling are controlled.
One of the tasks performed routinely by the Electromagnetics and Plasma Physics Analysis Department at Sandia National Laboratories is analyzing the effects of direct-strike lightning on Faraday cages that protect sensitive items. The Faraday cages analyzed thus far have many features in common. This report is an attempt to collect equations and other information that have been routinely used in the past in order to facilitate future analysis.
The design of experiments (DOEx) approach was used to characterize the Precision Laser Beam Welding Process with respect to four processing factors: Angle of Attack, Volts, Pulse Length, and Focus. The experiment was performed with Lap Joints, Nickel-Wire Joints, and Kovar-Wire Joints. The laser welding process and these types of welds are used in the manufacture of MC4368A Neutron Generators. For each weld type an individual optimal condition and operating window was identified. The widths of the operating windows that were identified by experimentation indicate that the laser weld process is very robust. This is highly desirable because it means that the quality of the resulting welds is not sensitive to the exact values of the processing factors within the operating windows. Statistical process control techniques can be used to ensure that the processing factors stay well within the operating window.
An optical sensor system has been developed for the autonomous monitoring of NO{sub 2} evolution in energetic material aging studies. The system is minimally invasive, requiring only the presence of a small sensor film within the aging chamber. The sensor material is a perylene/PMMA film that is excited by a blue LED light source and the fluorescence detected with a CCD spectrometer. Detection of NO{sub 2} gas is done remotely through the glass window of the aging chamber. Irreversible reaction of NO{sub 2} with perylene, producing the non-fluorescent nitroperylene, provides the optical sensing scheme. The rate of fluorescence intensity loss over time can be modeled using a numerical solution to the coupled diffusion and a nonlinear chemical reaction problem to evaluate NO{sub 2} concentration levels. The light source, spectrometer, spectral acquisition, and data processing were controlled through a Labivew program run by a laptop PC. Due to the long times involved with materials aging studies the system was designed to turn on, warm up, acquire data, power itself off, then recycle at a specific time interval. This allowed the monitoring of aging HE material over the period of several weeks with minimal power consumption and stable LED light output. Despite inherent problems with gas leakage of the aging chamber they were able to test the sensor system in the field under an accelerated aging study of rocket propellant. They found that the propellant evolved NO{sub 2} at a rate that yielded a concentration of between 10 and 100 ppm. The sensor system further revealed that the propellant, over an aging period of 25 days, evolves NO{sub 2} with cyclic behavior between active and dormant periods.
Dynamic thermography is a promising technology for inspecting metallic and composite structures used in high-consequence industries. However, the reliability and inspection sensitivity of this technology has historically been limited by the need for extensive operator experience and the use of human judgment and visual acuity to detect flaws in the large volume of infrared image data collected. To overcome these limitations new automated data analysis algorithms and software is needed. The primary objectives of this research effort were to develop a data processing methodology that is tied to the underlying physics, which reduces or removes the data interpretation requirements, and which eliminates the need to look at significant numbers of data frames to determine if a flaw is present. Considering the strengths and weakness of previous research efforts, this research elected to couple both the temporal and spatial attributes of the surface temperature. Of the possible algorithms investigated, the best performing was a radiance weighted root mean square Laplacian metric that included a multiplicative surface effect correction factor and a novel spatio-temporal parametric model for data smoothing. This metric demonstrated the potential for detecting flaws smaller than 0.075 inch in inspection areas on the order of one square foot. Included in this report is the development of a thermal imaging model, a weighted least squares thermal data smoothing algorithm, simulation and experimental flaw detection results, and an overview of the ATAC (Automated Thermal Analysis Code) software that was developed to analyze thermal inspection data.
This report documents how active structural control was used to significantly enhance the metal removal rate of a milling machine. An active structural control system integrates actuators, sensors, a control law and a processor into a structure for the purpose of improving the dynamic characteristics of the structure. Sensors measure motion, and the control law, implemented in the processor, relates this motion to actuator forces. Closed-loop dynamics can be enhanced by proper control law design. Actuators and sensors were imbedded within a milling machine for the purpose of modifying dynamics in such a way that mechanical energy, produced during cutting, was absorbed. This limited the on-set of instabilities and allowed for greater depths of cut. Up to an order of magnitude improvement in metal removal rate was achieved using this system. Although demonstrations were very successful, the development of an industrial prototype awaits improvements in the technology. In particular, simpler system designs that assure controllability and observability and control algorithms that allow for adaptability need to be developed.
Many problems in aeronautics can be described in terms of nonlinear systems of equations. Carleman developed a technique to linearize such equations that could lead to analytical solutions of nonlinear problems. Nonlinear problems are difficult to solve in closed form and therefore the construction of such solutions is usually nontrivial. This research will apply the Carleman linearization technique to three model problems: a two-degree-of-freedom (2DOF) ballistic trajectory, Blasius' boundary layer, and Van der Pol's equation and evaluate how well the technique can adequately approximate the solutions of these ordinary differential equations.
This research explores four experiments of adaptive host-based intrusion detection (ID) techniques in an attempt to develop systems that can detect novel exploits. The technique considered to have the most potential is adaptive critic designs (ACDs) because of their utilization of reinforcement learning, which allows learning exploits that are difficult to pinpoint in sensor data. Preliminary results of ID using an ACD, an Elman recurrent neural network, and a statistical anomaly detection technique demonstrate an ability to learn to distinguish between clean and exploit data. We used the Solaris Basic Security Module (BSM) as a data source and performed considerable preprocessing on the raw data. A detection approach called generalized signature-based ID is recommended as a middle ground between signature-based ID, which has an inability to detect novel exploits, and anomaly detection, which detects too many events including events that are not exploits. The primary results of the ID experiments demonstrate the use of custom data for generalized signature-based intrusion detection and the ability of neural network-based systems to learn in this application environment.
A new scheme to simulate elastic collisions in particle simulation codes is presented. The new scheme aims at simulating the collisions in the highly collisional regime, in which particle simulation techniques typically become computationally expensive. The new scheme is based on the concept of a grid-based collision field. According to this scheme, the particles perform a single collision with the background grid during a time step. The properties of the background field are calculated from the moments of the distribution function accumulated on the grid. The collision operator is based on the Langevin equation. Based on comparisons with other methods, it is found that the Langevin method overestimates the collision frequency for dilute gases.
A pointing control system is developed and tested for a flying gimbaled telescope. The two-axis pointing system is capable of sub-microradian pointing stability and high accuracy in the presence of large host vehicle jitter. The telescope also has high agility - it is capable of a 50° retarget (in both axes simultaneously) in less than 2 s. To achieve the design specifications, high-accuracy, high-resolution, two-speed resolvers were used, resulting in gimbal-angle measurements stable to 1.5 μrad. In addition, on-axis inertial angle displacement sensors were mounted on the telescope to provide host-vehicle jitter cancellation. The inertial angle sensors are accurate to about 100 nrad, but do not measure low-frequency displacements below 2 Hz. The gimbal command signal includes host-vehicle attitude information, which is band-limited. This provides jitter data below 20 Hz, but includes a variable latency between 15 and 25 ms. One of the most challenging aspects of this design was to combine the inertial-angle-sensor data with the less perfect information in the command signal to achieve maximum jitter reduction. The optimum blending of these two signals, along with the feedback compensation were designed using Quantitative Feedback Theory.
This report describes the underlying principles and goals of the Sandia ASCI Verification and Validation Program Validation Metrics Project. It also gives a technical description of two case studies, one in structural dynamics and the other in thermomechanics, that serve to focus the technical work of the project in Fiscal Year 2001.
This report documents an exploratory FY 00 LDRD project that sought to demonstrate the first steps toward a realistic computational representation of the variability encountered in individual human behavior. Realism, as conceptualized in this project, required that the human representation address the underlying psychological, cultural, physiological, and environmental stressors. The present report outlines the researchers' approach to representing cognitive, cultural, and physiological variability of an individual in an ambiguous situation while faced with a high-consequence decision that would greatly impact subsequent events. The present project was framed around a sensor-shooter scenario as a soldier interacts with an unexpected target (two young Iraqi girls). A software model of the ''Sensor Shooter'' scenario from Desert Storm was developed in which the framework consisted of a computational instantiation of Recognition Primed Decision Making in the context of a Naturalistic Decision Making model [1]. Recognition Primed Decision Making was augmented with an underlying foundation based on our current understanding of human neurophysiology and its relationship to human cognitive processes. While the Gulf War scenario that constitutes the framework for the Sensor Shooter prototype is highly specific, the human decision architecture and the subsequent simulation are applicable to other problems similar in concept, intensity, and degree of uncertainty. The goal was to provide initial steps toward a computational representation of human variability in cultural, cognitive, and physiological state in order to attain a better understanding of the full depth of human decision-making processes in the context of ambiguity, novelty, and heightened arousal.
Virtual Private Networking is a new communications technology that promises lower cost, more secure wide area communications by leveraging public networks such as the Internet. Sandia National Laboratories has embraced the technology for interconnecting remote sites to Sandia's corporate network, and for enabling remote access users for both dial-up and broadband access.
Public concern regarding the effects of noise generated by the detonation of excess and obsolete explosive munitions at U.S. Army demolition ranges is a continuing issue for the Army's demilitarization and disposal groups. Recent concerns of citizens living near the McAlester Army Ammunition Plant (MCAAP) in Oklahoma have lead the U.S. Army Defense Ammunition Center (DAC) to conduct a demonstration and evaluation of noise abatement techniques that could be applied to the MCAAP demolition range. With the support of the DAC, MCAAP, and Sandia National Laboratories (SNL), three types of noise abatement techniques were applied: aqueous foams, overburden (using combinations of sand beds and dirt coverings), and rubber or steel blast mats. Eight test configurations were studied and twenty-four experiments were conducted on the MCAAP demolition range in July of 2000. Instrumentation and data acquisition systems were fielded for the collection of near-field blast pressures, far-field acoustic pressures, plant boundary seismic signals, and demolition range meteorological conditions. The resulting data has been analyzed and reported, and a ranking of each technique's effects has been provided to the DAC.
A simplified and bounding methodology for analyzing the pressure buildup and hydrogen concentration within an unvented 2R container was developed (the 2R is a sealed container within a 6M package). The specific case studied was the gas buildup due to alpha radiolysis of water moisture sorbed on small quantities (less than 20 Ci per package) of plutonium oxide. Analytical solutions for gas pressure buildup and hydrogen concentration within the unvented 2R container were developed. Key results indicated that internal pressure buildup would not be significant for a wide range of conditions. Hydrogen concentrations should also be minimal but are difficult to quantify due to a large variation/uncertainty in model parameters. Additional assurance of non-flammability can be obtained by the use of an inert backfill gas in the 2R container.
The efficiency of the design-to-analysis process for translating solid-model-based design data to computational analysis model data plays a central role in the application of computational analysis to engineering design and certification. A review of the literature from within Sandia as well as from industry shows that the design-to-analysis process involves a number of complex organizational and technological issues. This study focuses on the design-to-analysis process from a business process standpoint and is intended to generate discussion regarding this important issue. Observations obtained from Sandia staff member and management interviews suggest that the current Sandia design-to-analysis process is not mature and that this cross-organizational issue requires committed high-level ownership. A key recommendation of the study is that additional resources should be provided to the computer aided design organizations to support design-to-analysis. A robust community of practice is also needed to continuously improve the design-to-analysis process and to provide a corporate perspective.
A Synthetic Aperture Radar (SAR) image is a two-dimensional projection of the radar reflectivity from a 3-dimensional object or scene. Stereoscopic SAR employs two SAR images from distinct flight paths that can be processed together to extract information of the third collapsed dimension (typically height) with some degree of accuracy. However, more than two SAR images of the same scene can similarly be processed to further improve height accuracy, and hence 3-dimensional position accuracy. This report shows how.
This work points out that the costates are actually discontinuous functions of time for optimal control problems with Coulomb friction. In particular these discontinuities occur at the time points where the velocity of the system changes sign. To our knowledge, this has not been noted before. This phenomenon is demonstrated on a minimum-time problem with Coulomb friction and the consistency of discontinuous costates and switching functions with respect to the input switches is shown.
In order to achieve higher rendering performance, the use of parallel sort-last architecture on a PC cluster is presented. The sort-last library (libpglc) can be linked to an existing parallel application to achieve high rendering rates. The efficient use of 64 commodity graphics cards enables to establish pace-setting rendering performance of 300 million triangles per second on extremely large data.
This report is a summary of the work completed in FY00 for science-based characterization of the processes used to fabricate cermet vias in source feedthrus. In particular, studies were completed to characterize the CND50 cermet slurry, characterize solvent imbibition, and identify critical via filling variables. These three areas of interest are important to several processes pertaining to the production of neutron generator tubes. Rheological characterization of CND50 slurry prepared with 94ND2 and Sandi94 primary powders were also compared. The 94ND2 powder was formerly produced at the GE Pinellas Plant and the Sandi94 is the new replacement powder produced at CeramTec. Processing variables that may effect the via-filling process were also studied and include: the effect of solids loading in the CND50 slurry; the effect of milling time; and the effect of Nuosperse (a slurry ''conditioner''). Imbibition characterization included a combination of experimental, theoretical, and computational strategies to determine solvent migration though complex shapes, specifically vias in the source feedthru component. Critical factors were determined using a controlled set of experiments designed to identify those variables that influence the occurrence of defects within the cermet filled via. These efforts were pursued to increase part production reliability, understand selected fundamental issues that impact the production of slurry-filled parts, and validate the ability of the computational fluid dynamics code, GOMA, to simulate these processes. Suggestions are made for improving the slurry filling of source feedthru vias.
As part of the full scale fuel fire experimental program, a series of JP-8 pool fire experiments with a large cylindrical calorimeter (3.66 m diameter), representing a C-141 aircraft fuselage, at the lee end of the fuel pool were performed at Naval Air Warfare Center, Weapons Division (NAWCWPNS). The series was designed to support Weapon System Safety Assessment (WSSA) needs by addressing the case of a transport aircraft subjected to a large fuel fire. The data collected from this mock series will allow for characterization of the fire environment via a survivable test fixture. This characterization will provide important background information for a future test series utilizing the same fuel pool with an actual C-141 aircraft in place of the cylindrical calorimeter.
This report provides an independent assessment of information on mixed waste streams, chemical compatibility information on polymers, and standard test methods for polymer properties. It includes a technology review of mixed low-level waste (LLW) streams and material compatibilities, validation for the plan to test the compatibility of simulated mixed wastes with potential seal and liner materials, and the test plan itself. Potential packaging materials were reviewed and evaluated for compatibility with expected hazardous wastes. The chemical and physical property measurements required for testing container materials were determined. Test methodologies for evaluating compatibility were collected and reviewed for applicability. A test plan to meet US Department of Energy and Environmental Protection Agency requirements was developed. The expected wastes were compared with the chemical resistances of polymers, the top-ranking polymers were selected for testing, and the most applicable test methods for candidate seal and liner materials were determined. Five recommended solutions to simulate mixed LLW streams are described. The test plan includes descriptions of test materials, test procedures, data collection protocols, safety and environmental considerations, and quality assurance procedures. The recommended order of testing to be conducted is specified.
In this study, the erosion properties of four sediments related to the Canaveral Ocean Dredged Material Disposal Site have been determined as a function of density, consolidation, and shear stress by means of a high shear stress sediment erosion flume at Sandia National Laboratories. Additional analysis was completed for each sediment to determine mineralogy, particle size, and organic content. This was done to support numerical modeling efforts, aid in effective management, and minimize environmental impact. The motivation for this work is based on concerns of dredged material transporting beyond the designated site and estimates of site capacity.
This paper provides an overview John Holland's Echo model, describes an implementation of the model, documents results from preliminary experiments using the model, and proposes further research in using Echo to study complex adaptive systems. Echo simulates the behavior of complex adaptive systems and can provide an experimental testbed for exploring theories of, and developing tools useful for analyzing these systems. Preliminary results indicate that the dynamic behavior of Echo can be used to generate interesting, time-series data that will be useful for evaluating the applicability of and developing tools, techniques, and possibly general theories, for the analysis of specific complex adaptive systems.
This report contains the design basis for a generic molten-salt solar power tower. A solar power tower uses a field of tracking mirrors (heliostats) that redirect sunlight on to a centrally located receiver mounted on top a tower, which absorbs the concentrated sunlight. Molten nitrate salt, pumped from a tank at ground level, absorbs the sunlight, heating it up to 565 C. The heated salt flows back to ground level into another tank where it is stored, then pumped through a steam generator to produce steam and make electricity. This report establishes a set of criteria upon which the next generation of solar power towers will be designed. The report contains detailed criteria for each of the major systems: Collector System, Receiver System, Thermal Storage System, Steam Generator System, Master Control System, and Electric Heat Tracing System. The Electric Power Generation System and Balance of Plant discussions are limited to interface requirements. This design basis builds on the extensive experience gained from the Solar Two project and includes potential design innovations that will improve reliability and lower technical risk. This design basis document is a living document and contains several areas that require trade-studies and design analysis to fully complete the design basis. Project- and site-specific conditions and requirements will also resolve open To Be Determined issues.
This report documents an investigation of irreversible electrical breakdown in ZnO varistors due to short pulses of high electric field and current density. For those varistors that suffer breakdown, there is a monotonic, pulse-by-pulse degradation in the switching electric field. The electrical and structural characteristics of varistors during and after breakdown are described qualitatively and quantitatively. Once breakdown is nucleated, the degradation typically follows a well-defined relationship between the number of post-initiation pulses and the degraded switching voltage. In some cases the degraded varistor has a remnant 20 {micro}m diameter hollow track showing strong evidence of once-molten ZnO. A model is developed for both electrical and thermal effects during high energy pulsing. The breakdown is assumed to start at one electrode and advance towards the other electrode as a thin filament of conductive material that grows incrementally with each successive pulse. The model is partially validated by experiments in which the varistor rod is cut at several different lengths from the electrode. Invariably one section of the cut varistor has a switching field that is not degraded while the other section(s) are heavily degraded. Based on the experiments and models of behavior during breakdown, some speculations about the nature of the nucleating mechanism are offered in the last section.
This report describes the 19-foot diameter blast tunnel at Sandia National Laboratories. The blast tunnel configuration consists of a 6 foot diameter by 200 foot long shock tube, a 6 foot diameter to 19 foot diameter conical expansion section that is 40 feet long, and a 19 foot diameter test section that is 65 feet long. Therefore, the total blast tunnel length is 305 feet. The development of this 19-foot diameter blast tunnel is presented. The small scale research test results using 4 inch by 8 inch diameter and 2 foot by 6 foot diameter shock tube facilities are included. Analytically predicted parameters are compared to experimentally measured blast tunnel parameters in this report. The blast tunnel parameters include distance, time, static, overpressure, stagnation pressure, dynamic pressure, reflected pressure, shock Mach number, flow Mach number, shock velocity, flow velocity, impulse, flow duration, etc. Shadowgraphs of the shock wave are included for the three different size blast tunnels.
Strain-induced self-assembly during semiconductor heteroepitaxy offers a promising approach to produce quantum nanostructures for nanologic and optoelectronics applications. Our current research direction aims to move beyond self-assembly of the basic quantum dot towards the fabrication of more complex, potentially functional structures such as quantum dot molecules and quantum wires. This report summarizes the steps taken to improve the growth quality of our GeSi molecular beam epitaxy process, and then highlights the outcomes of this effort.
Methods to determine unsaturated hydraulic properties can exhibit random and nonunique behavior. We assess the causes for these behaviors by visualizing microscale phase displacement processes that occur during equilibrium retention and transient outflow experiments. For both types of experiments we observe the drainage process to be composed of a mixture of fast air fingering and slower air back-filling. The influence of each of these microscale processes is controlled by a combination of the size and the speed of the applied boundary step, the initial saturation and its structure, and small-scale heterogeneities. Because the mixture of these microscale processes yields macroscale effective behavior, measured unsaturated flow properties are also a function of these controls. Such results suggest limitations on the current definitions and uniqueness of unsaturated hydraulic properties.
As a participating national lab in the inter-institutional effort to resolve performance issues of the non-elutable ion exchange technology for Cs extraction, they have carried out a series of characterization studies of UOP IONSIV{reg_sign} IE-911 and its component parts. IE-911 is a bound form (zirconium hydroxide-binder) of crystalline silicotitanate (CST) ion exchanger. The crystalline silicotitanate removes Cs from solutions by selective ion exchange. The performance issues of primary concern are: (1) excessive Nb leaching and subsequent precipitation of column-plugging Nb-oxide material, and (2) precipitation of aluminosilicate on IE-911 pellet surfaces, which may be initiated by dissolution of Si from the IE-911, thus creating a supersaturated solution with respect to silica. In this work, they have identified and characterized Si- and Nb-oxide based impurity phases in IE-911, which are the most likely sources of leachable Si and Nb, respectively. Furthermore, they have determined the criteria and mechanism for removal from IE-911 of the Nb-based impurity phase that is responsible for the Nb-oxide column plugging incidents.
This report summarizes progress from the Laboratory Directed Research and Development (LDRD) program during fiscal year 2000. In addition to a programmatic and financial overview, the report includes progress reports from 244 individual R and D projects in 13 categories.
Umbra is a new Sandia-developed modeling and simulation framework. The Umbra framework allows users to quickly build models and simulations for intelligent system development, analysis, experimentation, and control and supports tradeoff analyses of complex robotic systems, device, and component concepts. Umbra links together heterogeneous collections of modeling tools. The models in Umbra include 3D geometry and physics models of robots, devices and their environments. Model components can be built with varying levels of fidelity and readily switched to allow models built with low fidelity for conceptual analysis to be gradually converted to high fidelity models for later phase detailed analysis. Within control environments, the models can be readily replaced with actual control elements. This paper describes Umbra at a functional level and describes issues that Sandia uses Umbra to address.
The construction of inverse states in a finite field F{sub P{sub P{alpha}}} enables the organization of the mass scale by associating particle states with residue class designations. With the assumption of perfect flatness ({Omega}total = 1.0), this approach leads to the derivation of a cosmic seesaw congruence which unifies the concepts of space and mass. The law of quadratic reciprocity profoundly constrains the subgroup structure of the multiplicative group of units F{sub P{sub {alpha}}}* defined by the field. Four specific outcomes of this organization are (1) a reduction in the computational complexity of the mass state distribution by a factor of {approximately}10{sup 30}, (2) the extension of the genetic divisor concept to the classification of subgroup orders, (3) the derivation of a simple numerical test for any prospective mass number based on the order of the integer, and (4) the identification of direct biological analogies to taxonomy and regulatory networks characteristic of cellular metabolism, tumor suppression, immunology, and evolution. It is generally concluded that the organizing principle legislated by the alliance of quadratic reciprocity with the cosmic seesaw creates a universal optimized structure that functions in the regulation of a broad range of complex phenomena.
This report describes the use of PorSalsa, a parallel-processing, finite-element-based, unstructured-grid code for the simulation of subsurface nonisothermal two-phase, two component flow through heterogeneous porous materials. PorSalsa can also model the advective-dispersive transport of any number of species. General source term and transport coefficient implementation greatly expands possible applications. Spatially heterogeneous flow and transport data are accommodated via a flexible interface. Discretization methods include both Galerkin and control volume finite element methods, with various options for weighting of nonlinear coefficients. Time integration includes both first and second-order predictor/corrector methods with automatic time step selection. Parallel processing is accomplished by domain decomposition and message passing, using MPI, enabling seamless execution on single computers, networked clusters, and massively parallel computers.
Arithmetic conditions relating particle masses can be defined on the basis of (1) the supersymmetric conservation of congruence and (2) the observed characteristics of particle reactions and stabilities. Stated in the form of common divisors, these relations can be interpreted as expressions of genetic elements that represent specific particle characteristics. In order to illustrate this concept, it is shown that the pion triplet ({pi}{sup {+-}}, {pi}{sup 0}) can be associated with the existence of a greatest common divisor d{sub 0{+-}} in a way that can account for both the highly similar physical properties of these particles and the observed {pi}{sup {+-}}/{pi}{sup 0} mass splitting. These results support the conclusion that a corresponding statement holds generally for all particle multiplets. Classification of the respective physical states is achieved by assignment of the common divisors to residue classes in a finite field F{sub P{sub {alpha}}} and the existence of the multiplicative group of units F{sub P{sub {alpha}}} enables the corresponding mass parameters to be associated with a rich subgroup structure. The existence of inverse states in F{sub P{sub {alpha}}} allows relationships connecting particle mass values to be conveniently expressed in a form in which the genetic divisor structure is prominent. An example is given in which the masses of two neutral mesons (K{degree} {r_arrow} {pi}{degree}) are related to the properties of the electron (e), a charged lepton. Physically, since this relationship reflects the cascade decay K{degree} {r_arrow} {pi}{degree} + {pi}{degree}/{pi}{degree} {r_arrow} e{sup +} + e{sup {minus}}, in which a neutral kaon is converted into four charged leptons, it enables the genetic divisor concept, through the intrinsic algebraic structure of the field, to provide a theoretical basis for the conservation of both electric charge and lepton number. It is further shown that the fundamental source of supersymmetry can be expressed in terms of hierarchical relationships between odd and even order subgroups of F{sub P{sub {alpha}}}, an outcome that automatically reflects itself in the phenomenon of fermion/boson pairing of individual particle systems. Accordingly, supersymmetry is best represented as a group rather than a particle property. The status of the Higgs subgroup of order 4 is singular; it is isolated from the hierarchical pattern and communicates globally to the mass scale through the seesaw congruence by (1) fusing the concepts of mass and space and (2) specifying the generators of the physical masses.
The stress of scandium dideuteride, ScD{sub 2}, thin films is investigated during each stage of vacuum processing including metal deposition via evaporation, reaction and cooldown. ScD{sub 2} films with thin Cr underlayers are fabricated on three different substrate materials: molybdenum-alumina cermet, single crystal sapphire and quartz. In all experiments, the evaporated Cr and Sc metal is relatively stress-free. However, reaction of scandium metal with deuterium at elevated temperature to form a stoichiometric dideuteride phase leads to a large compressive in-plane film stress. Compression during hydriding results from an increased atomic density compared with the as-deposited metal film. After reaction with deuterium, samples are cooled to ambient temperature, and a tensile stress develops due to mismatched coefficients of thermal expansion (CTE) of the substrate-film couple. The residual film stress and the propensity for films to crack during cooldown depends principally on the substrate material when using identical process parameters. Films deposited onto quartz substrates show evidence of stress relief during cooldown due to a large CTE misfit; this is correlated with crack nucleation and propagation within films. All ScD{sub 2} layers remain in a state of tension when cooled to 30 C. An in-situ, laser-based, wafer curvature sensor is designed and implemented for studies of ScD{sub 2} film stress during processing. This instrument uses a two-dimensional array of laser beams to noninvasively monitor stress during sample rotation and with samples stationary. Film stress is monitored by scattering light off the backside of substrates, i.e., side opposite of the deposition flux.
An experiment to measure surface pressure data on a series of three stainless steel simulated parachute ribbons was conducted. During the first phase of the test, unsteady pressure measurements were made on the windward and leeward sides of the ribbons to determine the statistical properties of the surface pressures. Particle Image Velocimetry (PIV) measurements were simultaneously made to establish the velocity field in the wake of the ribbons and its correlation with the pressure measurements. In the second phase of the test, steady-state pressure measurements were made to establish the pressure distributions. In the third phase, the stainless steel ribbons were replaced with nylon ribbons and PIV measurements were made in the wake. A detailed error analysis indicates that the accuracy of the pressure measurements was very good. However, an anomaly in the flow field caused the wake behind the stainless steel ribbons to establish itself in a stable manner on one side of the model. This same stability was not present for the nylon ribbon model although an average of the wake velocity data indicated an apparent 2{degree} upwash in the wind tunnel flow field. Since flow angularity upstream of the model was not measured, the use of the data for code validation is not recommended without a second experiment to establish that upstream boundary condition.
In support of two major SNL programs, the Long-term Inflow and Structural Test (LIST) program and the Blade Manufacturing Initiative (BMI), three Micon 65/13M wind turbines have been erected at the USDA Agriculture Research Service (ARS) center in Bushland, Texas. The inflow and structural response of these turbines are being monitored with an array of 60 instruments: 34 to characterize the inflow, 19 to characterize structural response and 7 to characterize the time-varying state of the turbine. The primary characterization of the inflow into the LIST turbine relies upon an array of five sonic anemometers. Primary characterization of the structural response of the turbine uses several sets of strain gauges to measure bending loads on the blades and the tower and two accelerometers to measure the motion of the nacelle. Data are sampled at a rate of 30 Hz using a newly developed data acquisition system. The system features a time-synchronized continuous data stream and telemetered data from the turbine rotor. This paper documents the instruments and infrastructure that have been developed to monitor these turbines and their inflow.
Numerical methods may require derivatives of functions whose values are known only on irregularly spaced calculation points. This document presents and quantifies the performance of Moving Least-Squares (MLS), a method of derivative evaluation on irregularly spaced points that has a number of inherent advantages. The user selects both the spatial dimension of the problem and order of the highest conserved moment. The accuracy of calculations is maintained on highly irregularly spaced points. Not required are creation of additional calculation points or interpolation of the calculation points onto a regular grid. Implementation of the method requires the use of only a relatively small number of calculation points. The method is fast, robust and provides smooth results even as the order of the derivative increases.
Meso-scale manufacturing processes are bridging the gap between silicon-based MEMS processes and conventional miniature machining. These processes can fabricate two and three-dimensional parts having micron size features in traditional materials such as stainless steels, rare earth magnets, ceramics, and glass. Meso-scale processes that are currently available include, focused ion beam sputtering, micro-milling, micro-turning, excimer laser ablation, femtosecond laser ablation, and micro electro discharge machining. These meso-scale processes employ subtractive machining technologies (i.e., material removal), unlike LIGA, which is an additive meso-scale process. Meso-scale processes have different material capabilities and machining performance specifications. Machining performance specifications of interest include minimum feature size, feature tolerance, feature location accuracy, surface finish, and material removal rate. Sandia National Laboratories is developing meso-scale mechanical components and actuators which require meso-scale parts fabricated in a variety of materials. Subtractive meso-scale manufacturing processes expand the functionality of meso-scale components and complement silicon based MEMS and LIGA technologies.
This technical report presents the initial proposal and renewable proposals for an LDRD project whose intended goal was to enable applications to take full advantage of the hardware available on Sandia's current and future massively parallel supercomputers by analyzing various ways of combining distributed-memory and shared-memory programming models. Despite Sandia's enormous success with distributed-memory parallel machines and the message-passing programming model, clusters of shared-memory processors appeared to be the massively parallel architecture of the future at the time this project was proposed. They had hoped to analyze various hybrid programming models for their effectiveness and characterize the types of application to which each model was well-suited. The report presents the initial research proposal and subsequent continuation proposals that highlight the proposed work and summarize the accomplishments.
This paper studies the implementation of polar format, synthetic aperture radar image formation in modern Field Programmable Gate Arrays (FPGA's). The polar format algorithm is described in rough terms and each of the processing steps is mapped to FPGA logic. This FPGA logic is analyzed with respect to throughput and circuit size for compatibility with airborne image formation.