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.
The progress in the fabrication of high voltage GaN and AlGaN rectifiers, GaN/AlGaN HBT and GaN MOSFET is reviewed. Improvements in epitaxial layer quality are studied. The advances in fabrication techniques that led to the improvement of device performance are discussed.
This paper documents work performed to convert scanned range data to CAD solid model representation. The work successfully developed surface fitting algorithms for quadric surfaces (e.g. plane, cone, cylinder, and sphere), and a segmentation algorithm based entirely on surface type, rather than on a differential metric like Gaussian curvature. Extraction of all CAD-required parameters for quadric surface representation was completed. Approximate face boundaries derived from the original point cloud were constructed. Work to extrapolate surfaces, compute exact edges and solid connectivity was begun, but left incomplete due to funding reductions. The surface fitting algorithms are robust in the face of noise and degenerate surface forms.
As part of a project with SEMATECH, detailed chemical reaction mechanisms have been developed that describe the gas-phase and surface chemistry occurring during the fluorocarbon plasma etching of silicon dioxide and related materials. The fluorocarbons examined are C{sub 2}F{sub 6}, CHF{sub 3} and C{sub 4}F{sub 8}, while the materials studied are silicon dioxide, silicon, photoresist, and silica-based low-k dielectrics. These systems were examined at different levels, ranging from in-depth treatment of C{sub 2}F{sub 6} plasma etch of oxide, to a fairly cursory examination of C{sub 4}F{sub 8} etch of the low-k dielectric. Simulations using these reaction mechanisms and AURORA, a zero-dimensional model, compare favorably with etch rates measured in three different experimental reactors, plus extensive diagnostic absolute density measurements of electron and negative ions, relative density measurements of CF, CF{sub 2}, SiF and SiF{sub 2} radicals, ion current densities, and mass spectrometric measurements of relative ion densities.
As computers become faster, have more memory, and use multiple parallel processors, large, complex codes that more accurately simulate physical phenomena have emerged to utilize this capability. Most problems can benefit from this approach and many require it. But not all! There are problems for which simpler methods on more modest computers still work. The trick is to identify those problems, write the codes, and make their implementation sufficiently simple that they can be used conveniently by those who could profit from them. A Simple Plasma Code has been written with this philosophy in mind. It retains just enough physics to allow realistic simulations to be formulated and run quickly, even on a personal computer. This paper describes the physical model, its numerical implementation, and presents a sample simulation.
Window and free surface interfaces perturb the flow in compression wave experiments. The velocity of these interfaces is routinely measured in shock-compression experiments using interferometry (i.e., VISAR). Interface perturbations often must be accounted for before meaningful material property results can be obtained. For shockless experiments when stress is a single valued function of strain, the governing equations of motion are hyperbolic and can be numerically integrated forward or backward in either time or space with assured stability. Using the VISAR results as ''initial conditions'' the flow fields are integrated backward in space to the interior of the specimen where the VISAR interface has not perturbed the flow at earlier times and results can be interpreted as if the interface had not been present. This provides a rather exact correction for free surface perturbations. The method can also be applied to window interfaces by selecting the appropriate initial conditions. Applications include interpreting Z-accelerator ramp wave experiments. The method applies to multiple layers and multiple reverberations. For an elastic-plastic material model the flow is dissipative and the governing equations are parabolic. When the parabolic terms are small, the equations also can be successfully integrated backward in space. This is verified by using a traditional elastic-plastic wave propagation code with a backward-derived stress history as the boundary condition for a forward calculation. Calculated free surface histories match the starting VISAR record verifying that the backward method produced an accurate solution to the governing equations. With our cooperation, workers at Los Alamos have successfully applied the Sandia-developed backward technique for the time-dependent quasielastic material model and are analyzing stress histories at a spall plane using the VISAR free surface velocity measurement from a ''pullback'' experiment.
Superresolution concepts offer the potential of resolution beyond the classical limit. This great promise has not generally been realized. In this study we investigate the potential application of superresolution concepts to synthetic aperture radar. The analytical basis for superresolution theory is discussed. The application of the concept to synthetic aperture radar is investigated as an operator inversion problem. Generally, the operator inversion problem is ill posed. A criterion for judging superresolution processing of an image is presented.
This report is a supplement to ''The Unique Signal Concept for Detonation Safety in Nuclear Weapons,'' SAND91-1269, which provides a prerequisite fundamental background on the unique signal (UQS) concept. The UQS is one of the key constituents of Enhanced Nuclear Detonation Safety (ENDS), as outlined in Section 1 of that report. There have been many documents written over the past quarter of a century describing various aspects of the UQS, but none of these emphasized the mathematical approaches that help explain why the UQS is effective in resisting inadvertent pre-arming, even in abnormal environments and how UQS implementations can be quantitatively assessed. The intent of this report is to describe various pertinent mathematical methodologies (many of which have not been previously reported) without duplicating, any more than necessary, background information available in other reports. Mathematical UQS analysis is needed because of quantitative requirements associated with ENDS, and because limited comparisons of various implementation approaches can be quantified under mathematical modeling assumptions. Some of the mathematics-based results shown in this report are presented to explain: (1) The reasons that the UQS methodology can provide greater protection against accident environments than could combinational techniques (Sections 2.1 through 2.4); (2) The reason that the probability of inadvertently duplicating a UQS comprising n bivalued events cannot be estimated as low as (1/2) inches (Section 2.4); (3) The value of, and the Sandia National Laboratories policy on independent sequential communication of UQS events (Section 3.4); and (4) The care that must be exercised if any signal processing is necessary (Section 4). There are also numerous examples (e.g., in Appendices A and B) of ill-advised deviations from UQS methodology that can seriously degrade safety. These examples help demonstrate that the UQS methodology should not be compromised.
Dual control volume molecular dynamics was employed to study the flux of methane through channels of thin silicalite membranes. The DCANIS force field was analyzed to describe the adsorption isotherms of methane and ethane in silicalite. The alkane parameters and silicalite parameters were determined by fiiting the DCANIS force field to single-component vapor-liquid coexistence curves (VLCC) and adsorption isotherms respectively. The adsorption layers on the surfaces of thin silicalite membranes showed a sifnificant resistance to the flux of methane. The results depicted the insensitivity of permeance to both the average pressure and pressure drop.
For highly cross-linked polymer networks bonded to a solid surface, the effect of interfacial bond density and system size on interfacial fracture is studied using molecular dynamics simulations. Results for tensile and shear mode simulations are given. The correspondence between the stress-strain curve and the sequence of molecular deformations is obtained. The failure strain for a fully bonded surface is equal to the strain necessary to make taut the average of the minimal paths through the network from a bonded site on the bottom solid surface to a bonded site on the top surface. At fractional interfacial bond densities, cavities form above the nonbonded surface, yielding an inhomogeneous strain profile and a smaller failure strain. The failure strain and stress are linearly proportional to the number of bonds at the interface except in the tensile mode when number of bonds is so few that van der Waals interactions dominate. The failure mode is successfully constructed to be interracial by limiting the interfacial bond density to be less than the bulk bond density.
Advanced thin-film processing and packaging technologies are employed in the fabrication of new planar thin-film multijunction thermal converters (MJTCs). The processing, packaging, and design features build on experience gained from prior NIST demonstrations of thin-film converters and are optimized for improved sensitivity, bandwidth, manufacturability, and reliability.