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.
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.
