A probabilistic, transient, three-phase model of chemical transport through human skin has been developed to assess the relative importance of uncertain parameters and processes during chemical exposure assessments and transdermal drug delivery. Penetration routes through the skin that were modeled include the following: (1) intercellular diffusion through the multiphase stratum corneum; (2) aqueous-phase diffusion through sweat ducts; and (3) oil-phase diffusion through hair follicles. Uncertainty distributions were developed for the model parameters, and a Monte Carlo analysis was performed to simulate probability distributions of mass fluxes through each of the routes for a hypothetical scenario of chemical transport through the skin. At early times (60 seconds), the sweat ducts provided a significant amount of simulated mass flux into the bloodstream. At longer times (1 hour), diffusion through the stratum corneum became important because of its relatively large surface area. Sensitivity analyses using stepwise linear regression were also performed to identify model parameters that were most important to the simulated mass fluxes at different times.
We have developed a new inchworm actuator, consisting of a plate and two frictional clamps, which utilizes leveraged bending for improved amplitude control. Here we investigate its friction characteristics. We measure its average slip per cycle as a function of friction load and independently measure the clamp friction coefficients. A model is developed that takes into account the electromechanics of the actuation plate, boundary conditions and clamp friction. We find that the model does not satisfactorily describe the operation of the actuator. We attribute this to pre-sliding tangential deflections and devise a test whose results are consistent with this phenomenon. This suggests that stable pre-sliding deflections control the behavior of nanometer-scale slip events in MEMS.
A series of experiments has been performed in the Sandia National Laboratories FLAME facility with a 2-meter diameter JP-8 fuel pool fire. Sandia heat flux gages were employed to measure the incident flux at 8 locations outside the flame. Experiments were repeated to generate sufficient data for accurate confidence interval analysis. Additional sources of error are quantified and presented together with the data. The goal of this paper is to present these results in a way that is useful for validation of computer models that are capable of predicting heat flux from large fires. We anticipate using these data for comparison to validate models within the Advanced Simulation and Computing (ASC, formerly ASCI) codes FUEGO and SYRINX that predict fire dynamics and radiative transport through participating media. We present preliminary comparisons between existing models and experimental results.
An improved model for the gas damping of out-of-plane motion of a microbeam is developed based on the Reynolds equation (RE). A boundary condition for the RE is developed that relates the pressure at the beam perimeter to the beam motion. The two coefficients in this boundary condition are determined from Navier-Stokes (NS) simulations with the slip boundary condition for small slip lengths (relative to the gap height) and from Direct Simulation Monte Carlo (DSMC) molecular gas dynamics simulations for larger slip lengths. This boundary condition significantly improves the accuracy of the RE for cases where the beam width is only slightly greater than the gap height.
RMPP (reliable message passing protocol) is a lightweight transport protocol designed for clusters that provides end-to-end flow control and fault tolerance. In this article, presentations were made that compares RMPP to TCP, UDP, and "Utopia". The article compared the protocols on four benchmarks: bandwidth, latency, all-to-all, and communication-computation overlap. The results have shown that message-based protocols like RMPP have several advantages over TCP including ease of implementation, support for computation/communication overlap, and low CPU overhead.
The general problem considered is an optimization problem involving product design where some initial data are available and computer simulation is to be used to obtain more information. Resources and system complexity together restrict the number of simulations that can be performed in search of optimal settings for the product parameters. Consequently levels of these parameters, used in the simulations, (the experimental design) must be selected in an efficient way. We describe an algorithmic 'response-modeling' approach for performing this selection. The algorithm is illustrated using a rolamite design application. We provide (as examples) optimal one, two and three-point experimental designs for the rolamite computational analyses.
In many micro-scale fluid dynamics problems, molecular-level processes can control the interfacial energy and viscoelastic properties at a liquid-solid interface. This leads to a flow behavior that is very different from those similar fluid dynamics problems at the macro-scale. Presently, continuum modeling fails to capture this flow behavior. Molecular dynamics simulations have been applied to investigate these complex fluid-wall interactions at the nano-scale. Results show that the influence of the wall crystal lattice orientation on the fluid-wall interactions can be very important. To address those problems involving interactions of multiple length scales, a coupled atomistic-continuum model has been developed and applied to analyze flow in channels with atomically smooth walls. The present coupling strategy uses the molecular dynamics technique to probe the non-equilibrium flow near the channel walls and applies constraints to the fluid particle motion, which is coupled to the continuum flow modeling in the interior region. We have applied this new methodology to investigate Couette flow in micro-channels.
Existing paper-based site characterization models of salt domes at the four active U.S. Strategic Petroleum Reserve sites have been converted to digital format and visualized using modern computer software. The four sites are the Bayou Choctaw dome in Iberville Parish, Louisiana; the Big Hill dome in Jefferson County, Texas; the Bryan Mound dome in Brazoria County, Texas; and the West Hackberry dome in Cameron Parish, Louisiana. A new modeling algorithm has been developed to overcome limitations of many standard geological modeling software packages in order to deal with structurally overhanging salt margins that are typical of many salt domes. This algorithm, and the implementing computer program, make use of the existing interpretive modeling conducted manually using professional geological judgement and presented in two dimensions in the original site characterization reports as structure contour maps on the top of salt. The algorithm makes use of concepts of finite-element meshes of general engineering usage. Although the specific implementation of the algorithm described in this report and the resulting output files are tailored to the modeling and visualization software used to construct the figures contained herein, the algorithm itself is generic and other implementations and output formats are possible. The graphical visualizations of the salt domes at the four Strategic Petroleum Reserve sites are believed to be major improvements over the previously available two-dimensional representations of the domes via conventional geologic drawings (cross sections and contour maps). Additionally, the numerical mesh files produced by this modeling activity are available for import into and display by other software routines. The mesh data are not explicitly tabulated in this report; however an electronic version in simple ASCII format is included on a PC-based compact disk.
A major cause of failures in heat exchangers and steam generators in nuclear power plants is degradation of the tubes within them. The tube failure is often caused by the development of cracks that begin on the outer surface of the tube and propagate both inwards and laterally. A new technique was researched for detection of defects using a continuous-wave radar method within metal tubing. The experimental program resulted in a completed product development schedule and the design of an experimental apparatus for studying handling of the probe and data acquisition. These tests were completed as far as the prototypical probe performance allowed. The prototype probe design did not have sufficient sensitivity to detect a defect signal using the defined radar technique and did not allow successful completion of all of the project milestones. The best results from the prototype probe could not detect a tube defect using the radar principle. Though a more precision probe may be possible, the cost of design and construction was beyond the scope of the project. This report describes the probe development and the status of the design at the termination of the project.
This paper presents an analysis of utilizing unused cycles on supercomputers through the use of many small jobs. What we call "interstitial computing," is important to supercomputer centers for both productivity and political reasons. Interstitial computing makes use of the fact that small jobs are more or less fungible consumers of compute cycles that are more efficient for bin packing than the typical jobs on a supercomputer. An important feature of interstitial computing is that it not have a significant impact on the makespan of native jobs on the machine. Also, a facility can obtain higher utilizations that may only be otherwise possible with more complicated schemes or with very long wait times. The key contribution of this paper is that it provides theoretical and empirical guidelines for users and administrators for how currently unused supercomputer cycles may be exploited. We find that that interstitial computing is a more effective means for increasing machine utilization than increasing native job run times or size.
Improvements have been made at TRIUMF to permit higher proton intensities of up to 1010 cm-2s-1 over the energy range 20-500 MeV. This improved capability enables the study of displacement damage effects that require higher fluence irradiations. In addition, a high energy neutron irradiation capability has been developed for terrestrial cosmic ray soft error rate (SER) characterization of integrated circuits. The neutron beam characteristics of this facility are similar to those currently available at the Los Alamos National Laboratory WNR test facility. SER data measured on several SRAMs using the TRIUMF neutron beam are in good agreement with the results obtained on the same devices using the WNR facility. The TRIUMF neutron beam also contains thermal neutrons that can be easily removed by a sheet of cadmium. The ability to choose whether thermal neurons are present is a useful attribute not possible at the WNR.
A sample of polymeric propellant binder was aged from 0 to 60 days at 95 C and analyzed using FT-IR step scan photoacoustic spectroscopy. This technique has the ability of to obtain spectra of the polymer as a function of depth into the polymer material. Multivariate curve resolution was applied to the spectra data obtained to extract the contributions of the aged and un-aged spectral components from the spectra. It was found that multivariate curve resolution could efficiently separate highly overlapped spectra and yielded insights into the aging process.
We use a Monte Carlo approach to explore the potential impact of observation and inversion model errors on the spatial statistics of field-estimated unsaturated hydraulic properties. For this analysis we simulate tension infiltrometer measurements in a series of idealized realities, each consisting of spatially correlated random property fields. We consider only simple measurement errors that can be easily modeled. We show that estimated hydraulic properties are strongly biased by small, simple observation and inversion model errors. This bias can lead to order-of-magnitude errors in spatial statistics and artificial cross correlation between measured properties. The magnitude of bias varies with the true mean of the property field, the type of error considered, and the type of spatial statistic. We find no unique indicators of bias as property values may appear reasonable and spatial statistics may look realistic. Our results suggest new concerns for geostatisticians, stochastic modelers, and unsaturated zone practitioners who are unaware of the potential impact of spatial bias in field-estimated properties.
Shear-induced migration of particles is studied during the slow flow of suspensions of neutrally buoyant spheres, at 50% particle volume fraction, in an inelastic but shear-thinning, suspending fluid. The suspension is flowing in between a rotating inner cylinder and a stationary outer cylinder. The conditions are such that nonhydrodynamic effects are negligible. Nuclear magnetic resonance (NMR) imaging demonstrates that the movement of particles away from the high shear rate region is more pronounced than for a Newtonian suspending liquid. We test a continuum constitutive model for the evolution of particle concentration in a flowing suspension proposed by Phillips et al., but extended to shear-thinning, suspending fluids. The fluid constitutive equation is Carreau-like in its shear-thinning behavior but also varies with the local particle concentration. The model captures many of the trends found in the experimental data, but does not yet agree quantitatively. In fact, quantitative agreement with a diffusive flux constitutive equation would be impossible without the addition of another fitting parameter that may depend on the shear-thinning nature of the suspending fluid. Because of this, we feel that the Phillips model may be fundamentally inadequate for simulating flows of particles in non-Newtonian suspending fluids without the introduction of a normal stress correction or other augmenting terms.
Fires in aircraft engine nacelles must be rapidly suppressed to avoid loss of life and property. The design of new and retrofit suppression systems has become significantly more challenging due to the ban on production of Halon 1301 for environmental concerns. Since fire dynamics and the transport of suppressants within the nacelle are both largely determined by the available air flow, efforts to define systems using less effective suppressants greatly benefit from characterization of nacelle air flow fields. A combined experimental and computational study of nacelle air flow therefore has been initiated. Calculations have been performed using both CFD-ACE (a Computational Fluid Dynamics (CFD) model with a body-fitted coordinate grid) and VULCAN (a CFD-based fire field model with a Cartesian "brick" shaped grid). A quarter-scale test fixture was designed and fabricated for the purpose of obtaining spatially-resolved measurements of velocity and turbulence intensity in a smooth nacelle. Numerical calculations have been performed for the conditions of the experiment and comparisons with experimental results obtained from the quarter-scale test fixture are discussed. In addition, numerical simulations were performed to assess the sensitivity of the predictions to the grid size and to the turbulence models with and without wall functions. In general, the velocity predictions show very good agreement with the data in the center of the channel but deviate near the walls. The turbulence intensity results tend to amplify the differences in velocity, although most of the trends are in agreement. In addition, there were some differences between VULCAN and CFD-ACE results in the angled wall regions due to the Cartesian grid structure used by the VULCAN code. Also, the experimental data tended to show poorer resolution near the walls of the transition ducts. The increased uncertainty in the data highlights some of the challenges in getting data near the walls due to the low signal to noise ratio. Overall, this effort provided a benchmark case for both the VULCAN and CFD-ACE codes for the application of interest.
The purpose of this study was to generate the material database for carbon and glass composite panels created by the SCRIMP process. The materials tested were glass/polyester composites, two types of carbon/polyester composites, and carbon and glass hybrid composites. The differences between the two types of carbon/polyester, which we call Type 1 and Type 2, are the ply thickness (.037 inch/ply and .048 inch/ply) and slightly different treatment of polyester resin. The tests that were performed for this study are four-point-bending tests, tension tests, panel warping tests, and beam bend-twist coupling tests. The material properties of interest were basic longitudinal and transverse stiffness and strength, residual stress due to curing, and the effect of bend-twist coupling. The bend-twist coupling is a feature that can be added to the composite laminate or structure, such that when it is bent, it will also twist.
There have been significant efforts to develop cognitively plausible software architectures of human information processing in the last three decades. This report summarizes several architectures that continue to be developed. The specific type of cognitive models developed are known as production system architectures, which refers to the characterization of knowledge in terms of procedural (''how-to'' knowledge) condition-action relationships consisting of declarative (''what'' or factual) knowledge. To illustrate the ability for these models to instantiate human cognitive performance, a simulation using ACT-R (Adaptive Control of Thought - Rational) was implemented for a supervisory control task. Correlations between simulated and human learning of the task were measured and yielded correlations as high as 0.93.
Tethered films of poly n-isopropylacrylamide (PNIPAM) films have been developed as materials that can be used to switch the chemistry of a surface in response to thermal activation. In water, PNIPAM exhibits a thermally-activated phase transition that is accompanied by significant changes in polymer volume, water contact angle, and protein adsorption characteristics. New synthesis routes have been developed to prepare PNIPAM films via in-situ polymerization on self-assembled monolayers. Swelling transitions in tethered films have been characterized using a wide range of techniques including surface plasmon resonance, attenuated total reflectance infrared spectroscopy, interfacial force microscopy, neutron reflectivity, and theoretical modeling. PNIPAM films have been deployed in integrated microfluidic systems. Switchable PNIPAM films have been investigated for a range of fluidic applications including fluid pumping via surface energy switching and switchable protein traps for pre-concentrating and separating proteins on microfluidic chips.
These editing tips contain helpful suggestions to assist writers who are writing, editing, and publishing technical publications in the JNWPS. The suggestions clarify some of the most common writing problems and requirements of two publications used in the JNWPS: ''DOE-DTRA TP 1-1, Joint Nuclear Weapons Publications System Operating Procedures, Specifications, and Standards, and United States Government Printing Office Style Manual''. Topics include requirements for abbreviations, formats for drafts, layouts of illustrations and tables, appropriate wording for interim changes, guidance for creating a list of effective pages, how to insert and delete pages and paragraphs, referencing other technical publications, use of revision bars, requirements for safety precautions, use of hyphens, and how to place warnings, cautions, and notes. Also included are a writer's checklist, samples of draft title pages, and a section of helpful tips for the writers who use the department's desktop publishing software program, Adobe{reg_sign} FrameMaker{reg_sign}.
The primary goal of this portion of the LDRD is to develop a vertical programmable diffraction grating that can be fabricated with Sandia's Ultra-planar Multi-level MEMS Technology, the SUMMiT V{trademark} process. This grating is targeted for use in a chemical detection system dubbed the Polychromator. A secondary goal is to design diffraction grating structures with additional degrees of freedom (DOF). Gratings with 2.5 microns of vertical stroke have been realized. In addition, rotational DOF grating structures have been successfully actuated, and a structure has been developed that minimizes residual stress effects.
This project combined nanocomposite materials with microfabricated optical device structures for the development of microsensor arrays. For the nanocomposite materials we have designed, developed, and characterized self-assembling, organic/inorganic hybrid optical sensor materials that offer highly selective, sensitive, and reversible sensing capability with unique hierarchical nanoarchitecture. Lipid bilayers and micellar polydiacetylene provided selective optical response towards metal ions (Pb(II), Hg(II)), a lectin protein (Concanavalin A), temperature, and organic solvent vapor. These materials formed as composites in silica sol-gels to impart physical protection of the self-assembled structures, provide a means for thin film surface coatings, and allow facile transport of analytes. The microoptical devices were designed and prepared with two- and four-level diffraction gratings coupled with conformal gold coatings on fused silica. The structure created a number of light reflections that illuminated multiple spots along the silica surface. These points of illumination would act as the excitation light for the fluorescence response of the sensor materials. Finally, we demonstrate an integrated device using the two-level diffraction grating coupled with the polydiacetylene/silica material.
This paper presents a high-level overview of the algorithms and supporting functionality provided by SIERRA Framework Version 3 for h-adaptive finite-element mechanics application development. Also presented is a fairly comprehensive description of what is required by the application codes to use the SIERRA h-adaptivity services. In general, the SIERRA framework provides the functionality for hierarchically subdividing elements in a distributed parallel environment, as well as dynamic load balancing. The mechanics application code is required to supply an a posteriori error indicator, prolongation and restriction operators for the field variables, hanging-node constraint handlers, and execution control code. This paper does not describe the Application Programming Interface (API), although references to SIERRA framework classes are given where appropriate.
Vibrational spectra can serve as chemical fingerprints for positive identification of chemical and biological warfare molecules. The required speed and sensitivity might be achieved with surface-enhanced Raman spectroscopy (SERS) using nanotextured metal surfaces. Systematic and reproducible methods for preparing metallic surfaces that maximize sensitivity have not been previously developed. This work sought to develop methods for forming high-efficiency metallic nanostructures that can be integrated with either gas or liquid-phase chem-lab-on-a-chip separation columns to provide a highly sensitive, highly selective microanalytical system for detecting current and future chem/bio agents. In addition, improved protein microchromatographic systems have been made by the creation of acrylate-based porous polymer monoliths that can serve as protein preconcentrators to reduce the optical system sensitivity required to detect and identify a particular protein, such as a bacterial toxin.
This report covers research done as part of US Department of Energy contract DE-PS26-99FT14299 with the Fuel Cell Propulsion Institute on the fuel cell RATLER{trademark} vehicle, Lurch, as well as work done on the fuel cells designed for the vehicle. All work contained within this report was conducted at the Robotic Vehicle Range at Sandia National Laboratories in Albuquerque New Mexico. The research conducted includes characterization of the duty cycle of the robotic vehicle. This covers characterization of its various abilities such as hill climbing and descending, spin-turns, and driving on level ground. This was accomplished with the use of current sensors placed in the vehicle in conjunction with a Data Acquisition System (DAS), which was also created at Sandia Labs. Characterization of the two fuel cells was accomplished using various measuring instruments and techniques that will be discussed later in the report. A Statement of Work for this effort is included in Appendix A. This effort was able to complete characterization of vehicle duty cycle elements using battery power, but problems with the fuel cell control systems prevented completion of the characterization of the fuel cell operation on the benchtop and in the vehicle. Some data was obtained characterizing the fuel cell current-voltage performance and thermal rise rate by bypassing elements of the control system.
RMS Guidelines defines the processes and conventions to manage both records and documents for the ASCI Verification and Validation Program at Sandia National Laboratories, employing the ASCI V&V RMS application. It is the definitive source for all information regarding the creation, submittal, use, maintenance, and disposition of records and documents. This document is also used as evidence of meeting records management requirements as stated in DOE Order 414.1A, Quality Assurance, and Sandia National Laboratories Corporate Technical Business Practice TBP-500, Records Management.
We have conducted surface treatment and alloying experiments with Al, Fe, and Ti-based metals on the RHEPP-1 accelerator (0.8 MV, 20 W, 80 ns FHWM, up to 1 Hz repetition rate) at Sandia National Laboratories. Ions are generated by the MAP gas-breakdown active anode, which can yield a number of different beam species including H, N, and C, depending upon the injected gas. Beams of intense pulsed high-power ion beams have been used to produce surface modification by changes in microstructure caused by rapid heating and cooling of the surface. Increase of beam power leads to ablation of a target surface, and redeposition of ablated material onto a separate substrate. Experiments are described in which ion beams are used in an attempt to increase high-voltage breakdown of a treated surface. Surface alloying of coated Pt and Hf layers is also described. This mixing of a previously deposited thin-film layer into a Ti-alloy substrate leads to significantly enhanced surface wear durability, compared to either untreated Ti-alloy alone, or the Ti alloy alone treated with the ion beam. Thin-film layers have been produced from a number of target materials. Films of fine-grain Pt and Er are described, and are compared to conventionally formed films. First attempts to form high-dielectric constant BaTiO{sub 3} are described.
High-energy ion-irradiated 3.3-nm oxynitride film and 2.2-nm SiO2-film MOS capacitors show premature break-down during subsequent electrical stress. This degradation in breakdown increases with increasing ion linear energy transfer (LET), increasing ion fluence, and decreasing oxide thickness. The reliability degradation due to high-energy ion-induced latent defects is explained by a simple percolation model of conduction through SiO2 layers with irradiation and/or electrical stress-induced defects. Monitoring the gate-leakage current reveals the presence of latent defects in the dielectric films. These results may be significant to future single-event effects and single-event gate rupture tests for MOS devices and ICs with ultrathin gate oxides.
The Electricity Generation Cost Simulation Model (GenSim) is a user-friendly, high-level dynamic simulation model that calculates electricity production costs for variety of electricity generation technologies, including: pulverized coal, gas combustion turbine, gas combined cycle, nuclear, solar (PV and thermal), and wind. The model allows the user to quickly conduct sensitivity analysis on key variables, including: capital, O&M, and fuel costs; interest rates; construction time; heat rates; and capacity factors. The model also includes consideration of a wide range of externality costs and pollution control options for carbon dioxide, nitrogen oxides, sulfur dioxide, and mercury. Two different data sets are included in the model; one from the US. Department of Energy (DOE) and the other from Platt's Research Group. Likely users of this model include executives and staff in the Congress, the Administration and private industry (power plant builders, industrial electricity users and electric utilities). The model seeks to improve understanding of the economic viability of various generating technologies and their emissions trade-offs. The base case results, using the DOE data, indicate that in the absence of externality costs, or renewable tax credits, pulverized coal and gas combined cycle plants are the least cost alternatives at 3.7 and 3.5 cents/kwhr, respectively. A complete sensitivity analysis on fuel, capital, and construction time shows that these results coal and gas are much more sensitive to assumption about fuel prices than they are to capital costs or construction times. The results also show that making nuclear competitive with coal or gas requires significant reductions in capital costs, to the $1000/kW level, if no other changes are made. For renewables, the results indicate that wind is now competitive with the nuclear option and is only competitive with coal and gas for grid connected applications if one includes the federal production tax credit of 1.8cents/kwhr.
This report describes the process that will be used to develop and submit for approval designs for the Sandia Extended Network (SXN). The process follows the steps in the Change Management Process used in the Telecommunication Operations Department's quality management system. Those steps are planning, requirements review, detailed design analysis, implementation, verification, and validation. Two companion reports complete a description of the designs to date: ''Sandia Extended Network: Design Requirements and Sandia Extended Network: Conceptual Design Definition.''
This report contains the design requirements for creating a limited-access Sandia Extended Network (SXN), which would be used to collaborate with Nuclear Weapons Complex Labs personnel, university collaborators, industry, and others who may not be allowed accounts on the Sandia Restricted Network (SRN). This document contains the design requirements for creating a limited-access Sandia Extended Network (SXN), which would be used by non-Sandians to collaborate with NWC Labs personnel and others who are not allowed accounts on the Sandia Restricted Network (SRN). Its main purpose is to articulate the requirements upon which the design options and hardware costs for the Sandia eXtended Network (SXN) can be based and in turn presented to 8900 and 9300 Management. The requirements are further addressed in reports outlining its security architecture and in the five-volume set of network architecture reports: An Architecture for the Sandia Extended Network: Overview; Detailed Description of the Architecture, Design of the Model, and Balanced Protections; Background of the Architecture and its Relevance to Sandia; Terminology and Concepts Relevant to Networks; and Policy-Based Networks and Information Management.
The goal of this project was to increase the power of vertical cavity surface emitting lasers and to convert their wavelength into the blue/ultraviolet and the infrared for sensing applications. We have increased the power to the multi-watt level and have generated several milliwatts of blue light using optical pumping. Electrical pump has been less successful, but we have identified the problems and begun work to overcome them using a bottom emitting design.
This report describes the various methods and circuits that have been developed to detect an islanding condition for photovoltaic applications and presents three methods that have been developed to test those methods and circuits. Passive methods for detecting an islanding condition basically monitor parameters such as voltage and frequency and/or their characteristics and cause the inverter to cease converting power when there is sufficient transition from normal specified conditions. Active methods for detecting the island introduce deliberate changes or disturbances to the connected circuit and then monitor the response to determine if the utility grid with its stable frequency, voltage and impedance is still connected. If the small perturbation is able to affect the parameters of the load connection within prescribed requirements, the active circuit causes the inverter to cease power conversion and delivery of power to the loads. The methods not resident in the inverter are generally controlled by the utility or have communications between the inverter and the utility to affect an inverter shut down when necessary. This report also describes several test methods that may be used for determining whether the anti-islanding method is effective. The test circuits and methodologies used in the U.S. have been chosen to limit the number of tests by measuring the reaction of a single or small number of inverters under a set of consensus-based worst-case conditions.
A new protocol technology is just starting to emerge from the laboratory environment. Its stated purpose is to provide an additional means in which networks, and the services that reside on them, can be protected from adversarial compromise. This report has a two-fold objective. First is to provide the reader with an overview of this emerging Dynamic Defenses technology using Dynamic Network Address Translation (Dynat). This ''structure overview'' is concentrated in the body of the report, and describes the important attributes of the technology. The second objective is to provide a framework that can be used to help in the classification and assessment of the different types of dynamic defense technologies along with some related capabilities and limitations. This information is primarily contained in the appendices.
An Interferometric Moving Target Indicator radar can be used to measure the tangential velocity component of a moving target. Multiple baselines, along with the conventional radial velocity measurement, allow estimating the true 3-D velocity vector of a target.
We have undertaken the synthesis of a thin film ''All Ceramic Battery'' (ACB) using solution route processes. Based on the literature and experimental results, we selected SnO{sub 2}, LiCoO{sub 2}, and LiLaTiO{sub 3} (LLT) as the anode, cathode, and electrolyte, respectively. Strain induced by lattice mismatch between the cathode and bottom electrode, as estimated by computational calculations, indicate that thin film orientations for batteries when thicknesses are as low as 500 {angstrom} are strongly controlled by surface energies. Therefore, we chose platinized silicon as the basal platform based on our previous experience with this material. The anode thin films were generated by standard spin-cast methods and processing using a solution of [Sn(ONep)]{sub 8} and HOAc which was found to form Sn{sub 6}(O){sub 4}(ONep){sub 4}. Electrochemical evaluation showed that the SnO{sub 2} was converted to Sn{sup o} during the first cycle. The cathode was also prepared by spin coating using the novel [Li(ONep)]{sub 8} and Co(OAc){sub 2}. The films could be electrochemically cycled (i.e., charged/discharged), with all of the associated structural changes being observable by XRD. Computational models indicated that the LLT electrolyte would be the best available ceramic material for use as the electrolyte. The LLT was synthesized from [Li(ONep)]{sub 8}, [Ti(ONep){sub 4}]{sub 2}, and La(DIP){sub 3}(py){sub 3} with RTP processing at 900 C being necessary to form the perovskite phase. Alternatively, a novel route to thin films of the block co-polymer ORMOLYTE was developed. The integration of these components was undertaken with each part of the assembly being identifiably by XRD analysis (this will allow us to follow the progress of the charge/discharge cycles of the battery during use). SEM investigations revealed the films were continuous with minimal mixing. All initial testing of the thin-film cathode/electrolyte/anode ACB devices revealed electrical shorting. Alternative approaches for preparing non-shorted devices (e.g. inverted and side-by-side) are under study.
The SIERRA Framework core services provide essential services for managing the mesh data structure, computational fields, and physics models of an application. An application using these services will supply a set of physics models, define the computational fields that are required by those models, and define the mesh upon which its physics models operate. The SIERRA Framework then manages all of the data for a massively parallel multiphysics application.
It was discovered that MgO or Mg(OH){sub 2} when it reacts with water is a very strong sorbent for arsenic. Distribution constants, or K{sub d} values, are as high as 1 x 10{sup 6} L/mole. In this work, Mg(OH){sub 2} and other compounds have been investigated as sorbents for arsenic and other contaminants. This work has resulted in several major accomplishments including: (1) design, construction, and testing of a pressure sand filter to remove Mg(OH){sub 2} after it has sorbed arsenic from water, (2) stabilization of Mg(OH){sub 2} as a Sorrel's cement against reaction with carbonate that results in MgCO{sub 3} formation decreasing the efficiency of Mg(OH){sub 2} to sorb arsenic, and (3) the development of a new, very promising sorbent for arsenic based on zirconium. Zirconium is an environmentally benign material found in many common products such as toothpaste. It is currently used in water treatment and is very inexpensive. In this work, zirconium has been bonded to activated carbon, zeolites, sand and montmorillonite. Because of its high charge in ionic form (+6), zirconium is a strong sorbent for many anions including arsenic. In equilibrium experiments arsenic concentrations in water were reduced from 200 ppb to less than 1 ppb in less than 1 minute of contact time. Additionally, analytical methods for detecting arsenic in water have also been investigated. Various analytical techniques including HPLC, AA and ICP-MS are used for quantification of arsenic. Due to large matrix interferences HPLC and AA techniques are not very selective and are time consuming. ICP-MS is highly efficient, requires a low sample volume and has a high tolerance for interferences. All these techniques are costly and require trained staff, and with the exception of ICP-MS, these methods cannot be used at low ppb arsenic concentration without using a pre-concentration step. An alternative to these traditional techniques is to use a colorimetric method based on leucocrystal violet dye interaction with iodine. This method has been adapted in our facility for quantifying arsenic concentrations down to 14 ppb.
Aspen, a powerful economic modeling tool that uses agent modeling and genetic algorithms, can accurately simulate the economy. In it, individuals are hired by firms to produce a good that households then purchase. The firms decide what price to charge for this good, and based on that price, the households determine which firm to purchase from. We will attempt to discover the Nash Equilibrium price found in this model under two different methods of determining how many orders each firm receives. To keep it simple, we will assume there are only two firms in our model, and that these firms compete for the sale of one identical good.
This report outlines our work on the integration of high efficiency photonic lattice structures with MEMS (MicroElectroMechanical Systems). The simplest of these structures were based on 1-D mirror structures. These were integrated into a variety of devices, movable mirrors, switchable cavities and finally into Bragg fiber structures which enable the control of light in at least 2 dimensions. Of these devices, the most complex were the Bragg fibers. Bragg fibers consist of hollow tubes in which light is guided in a low index media (air) and confined by surrounding Bragg mirror stacks. In this work, structures with internal diameters from 5 to 30 microns have been fabricated and much larger structures should also be possible. We have demonstrated the fabrication of these structures with short wavelength band edges ranging from 400 to 1600nm. There may be potential applications for such structures in the fields of integrated optics and BioMEMS. We have also looked at the possibility of waveguiding in 3 dimensions by integrating defects into 3-dimensional photonic lattice structures. Eventually it may be possible to tune such structures by mechanically modulating the defects.
The purpose of this review was to provide insights and information to Sandia National Laboratories' (SNL) Education Council on the state of technical education and training at SNL in order to address the concern that a change in philosophy surrounding education had occurred. To accomplish this, the status of current and past technical training and education programs were compared, and significant changes at SNL were assessed for their impact on education and training. Major changes in education and training are in the advertisement of course offerings, the course delivery methods, and the funding mechanisms for student and instructor time as well as course costs. The significant changes in SNL which influenced technical training and education are the considerable increase in mandatory or compliance training, a fundamental shift in SNL's management structure from an institutional structure to a more business-like, project-budgeted structure, and the change in SNL's mission at the end of the Cold War. These changes contributed to less time for technical training, reduction of training funds, elimination of some training, and a Service Center approach to paying for training. Most importantly, the overall combined effect has resulted in a shift from a strategic to a tactical training approach. The Corporate Training Department (CTD) has maneuvered to accommodate these changes and keep abreast of constantly changing needs.
Sandia National Laboratories has been investigating the use of remotely operated weapon platforms in Department of Energy (DOE) facilities. These platforms offer significant force multiplication and enhancement by enabling near instantaneous response to attackers, increasing targeting accuracy, removing personnel from direct weapon fire, providing immunity to suppressive fire, and reducing security force size needed to effectively respond. Test results of the Telepresent Rapid Aiming Platform (TRAP) from Precision Remotes, Inc. have been exceptional and response from DOE sites and the U.S. Air Force is enthusiastic. Although this platform performs comparably to a trained marksman, the target acquisition speeds are up to three times longer. TRAP is currently enslaved to a remote operator's joystick. Tracking moving targets with a joystick is difficult; it dependent upon target range, movement patterns, and operator skill. Even well-trained operators encounter difficulty tracking moving targets. Adding intelligent targeting capabilities on a weapon platform such as TRAP would significantly improve security force response in terms of effectiveness and numbers of responders. The initial goal of this project was to integrate intelligent targeting with TRAP. However, the unavailability of a TRAP for laboratory purposes drove the development of a new platform that simulates TRAP but has a greater operating range and is significantly faster to reposition.
The Global Energy Futures Model (GEFM) is a demand-based, gross domestic product (GDP)-driven, dynamic simulation tool that provides an integrated framework to model key aspects of energy, nuclear-materials storage and disposition, environmental effluents from fossil and non fossil energy and global nuclear-materials management. Based entirely on public source data, it links oil, natural gas, coal, nuclear and renewable energy dynamically to greenhouse-gas emissions and 12 other measures of environmental impact. It includes historical data from 1990 to 2000, is benchmarked to the DOE/EIA/IEO 2001 [5] Reference Case for 2000 to 2020, and extrapolates energy demand through the year 2050. The GEFM is globally integrated, and breaks out five regions of the world: United States of America (USA), the Peoples Republic of China (China), the former Soviet Union (FSU), the Organization for Economic Cooperation and Development (OECD) nations excluding the USA (other industrialized countries), and the rest of the world (ROW) (essentially the developing world). The GEFM allows the user to examine a very wide range of ''what if'' scenarios through 2050 and to view the potential effects across widely dispersed, but interrelated areas. The authors believe that this high-level learning tool will help to stimulate public policy debate on energy, environment, economic and national security issues.
Recent terrorist attacks in the United States have increased concerns about potential national security consequences from energy supply disruptions. The purpose of this Laboratory Directed Research & Development (LDRD) is to develop a high-level dynamic simulation model that would allow policy makers to explore the national security consequences of major US. energy supply disruptions, and to do so in a way that would integrate energy, economic and environmental components. The model allows exploration of potential combinations of demand-driven energy supplies that meet chosen policy objectives, including: Mitigating economic losses, measured in national economic output and employment levels, due to terrorist activity or forced outages of the type seen in California; Control of greenhouse gas levels and growth rates; and Moderating US. energy import requirements. This work has built upon the Sandia US. Energy and greenhouse Gas Model (USEGM) by integrating a macroeconomic input-output framework into the model, adding the capability to assess the potential economic impact of energy supply disruptions and the associated national security issues. The economic impacts of disruptions are measured in terms of lost US. output (e.g., GDP, sectoral output) and lost employment, and are assessed either at a broad sectoral level (3 sectors) or at a disaggregated level (52 sectors). In this version of the model, physical energy disruptions result in quantitative energy shortfalls, and energy prices are not permitted to rise to clear the markets.
This LDRD project has involved the development and application of Sandia's massively parallel materials modeling software to several significant biophysical systems. They have been successful in applying the molecular dynamics code LAMMPS to modeling DNA, unstructured proteins, and lipid membranes. They have developed and applied a coupled transport-molecular theory code (Tramonto) to study ion channel proteins with gramicidin A as a prototype. they have used the Towhee configurational bias Monte-Carlo code to perform rigorous tests of biological force fields. they have also applied the MP-Sala reacting-diffusion code to model cellular systems. Electroporation of cell membranes has also been studied, and detailed quantum mechanical studies of ion solvation have been performed. In addition, new molecular theory algorithms have been developed (in FasTram) that may ultimately make protein solvation calculations feasible on workstations. Finally, they have begun implementation of a combined molecular theory and configurational bias Monte-Carlo code. They note that this LDRD has provided a basis for several new internal (e.g. several new LDRD) and external (e.g. 4 NIH proposals and a DOE/Genomes to Life) proposals.
As demonstrated by the anthrax attack through the United States mail, people infected by the biological agent itself will give the first indication of a bioterror attack. Thus, a distributed information system that can rapidly and efficiently gather and analyze public health data would aid epidemiologists in detecting and characterizing emerging diseases, including bioterror attacks. We propose using clusters of adverse health events in space and time to detect possible bioterror attacks. Space-time clusters can indicate exposure to infectious diseases or localized exposure to toxins. Most space-time clustering approaches require individual patient data. To protect the patient's privacy, we have extended these approaches to aggregated data and have embedded this extension in a sequential probability ratio test (SPRT) framework. The real-time and sequential nature of health data makes the SPRT an ideal candidate. The result of space-time clustering gives the statistical significance of a cluster at every location in the surveillance area and can be thought of as a ''health-index'' of the people living in this area. As a surrogate to bioterrorism data, we have experimented with two flu data sets. For both databases, we show that space-time clustering can detect a flu epidemic up to 21 to 28 days earlier than a conventional periodic regression technique. We have also tested using simulated anthrax attack data on top of a respiratory illness diagnostic category. Results show we do very well at detecting an attack as early as the second or third day after infected people start becoming severely symptomatic.
Thermoluminescent dosimeters (TLDs), particularly CaF{sub 2}:Mn, are often used as photon dosimeters in mixed (n/{gamma}) field environments. In these mixed field environments, it is desirable to separate the photon response of a dosimeter from the neutron response. For passive dosimeters that measure an integral response, such as TLDs, the separation of the two components must be performed by post-experiment analysis because the TLD reading system cannot distinguish between photon and neutron produced response. Using a model of an aluminum-equilibrated TLD-400 chip, a systematic effort has been made to analytically determine the various components that contribute to the neutron response of a TLD reading. The calculations were performed for five measured reactor neutron spectra and one theoretical thermal neutron spectrum. The five measured reactor spectra all have dosimetry quality experimental values for aluminum-equilibrated TLD-400 chips. Calculations were used to determined the percentage of the total TLD response produced by neutron interactions in the TLD and aluminum equilibrator. These calculations will aid the Sandia National Laboratories-Radiation Metrology Laboratory (SNL-RML) in the interpretation of the uncertainty for TLD dosimetry measurements in the mixed field environments produced by SNL reactor facilities.
Buried landmines are often detected through the chemical signature in the air above the soil surface by mine detection dogs. Environmental processes play a significant role in the chemical signature available for detection. Due to the shallow burial depth of landmines, the weather influences the release of chemicals from the landmine, transport through the soil to the surface, and degradation processes in the soil. The effect of weather on the landmine chemical signature from a PMN landmine was evaluated with the T2TNT code for Kabul, Afghanistan. Results for TNT and DNT gas-phase and soil solid-phase concentrations are presented as a function of time of the day and time of the year.
This report describes research and development of methods to couple vastly different subsystems and physical models and to encapsulate these methods in a Java{trademark}-based framework. The work described here focused on developing a capability to enable design engineers and safety analysts to perform multifidelity, multiphysics analyses more simply. In particular this report describes a multifidelity algorithm for thermal radiative heat transfer and illustrates its performance. Additionally, it describes a module-based computer software architecture that facilitates multifidelity, multiphysics simulations. The architecture is currently being used to develop an environment for modeling the effects of radiation on electronic circuits in support of the FY 2003 Hostile Environments Milestone for the Accelerated Strategic Computing Initiative.
The quantitative analysis of microstructure and sequence distribution in polysiloxane copolymers using high-resolution solution {sup 29}Si NMR is reported. Copolymers containing dimethylsiloxane (DMS) and diphenysiloxane (DPS) monomer units prepared with either high vinyl content (HVM) or low vinyl content (LVM) were analyzed. The average run length (R{sub exp}), the number average sequence length (l{sub A}, l{sub B}), along with the various linkage probabilities (p{sub AA}, p{sub AB}, p{sub BA}, and p{sub BB}) were determined for different production lots of the LVM97 and HVM97 samples to address the lot variability of microstructure in these materials.
In an effort to recruit and retain skilled workers in the Manufacturing Science and Technology Center (14000), an innovative and highly diverse team at Sandia National Laboratories and the U.S. Department of Energy joined with concerned community constitutents, such as Albuquerque Technical Vocational Institute and the Albuquerque Public Schools, to offer mentoring and on-the-job training to qualified students in high schools and community colleges. Now, within several years of its inception, the educational program called the Advanced Manufacturing Trades Training Program is a model in the community and the nation, while enabling Sandia to have valuable trained and skilled employees to meet its national mission and workforce demands.
This manual describes the use of the Xyce Parallel Electronic Simulator code for simulating electrical circuits at a variety of abstraction levels. The Xyce Parallel Electronic Simulator has been written to support,in a rigorous manner, the simulation needs of the Sandia National Laboratories electrical designers. As such, the development has focused on improving the capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers. (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) A client-server or multi-tiered operating model wherein the numerical kernel can operate independently of the graphical user interface (GUI). (4) Object-oriented code design and implementation using modern coding-practices that ensure that the Xyce Parallel Electronic Simulator will be maintainable and extensible far into the future. The code is a parallel code in the most general sense of the phrase--a message passing parallel implementation--which allows it to run efficiently on the widest possible number of computing platforms. These include serial, shared-memory and distributed-memory parallel as well as heterogeneous platforms. Furthermore, careful attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved even as the number of processors grows. Another feature required by designers is the ability to add device models, many specific to the needs of Sandia, to the code. To this end, the device package in the Xyce Parallel Electronic Simulator is designed to support a variety of device model inputs. These input formats include standard analytical models, behavioral models and look-up tables. Combined with this flexible interface is an architectural design that greatly simplifies the addition of circuit models. One of the most important contribution Xyce makes to the designers at Sandia National Laboratories is in providing a platform for computational research and development aimed specifically at the needs of the Laboratory. With Xyce, Sandia now has an ''in-house''capability with which both new electrical (e.g., device model development) and algorithmic (e.g., faster time-integration methods) research and development can be performed. Furthermore, these capabilities will then be migrated to the end users.
Fixtures are tools used to hold parts in specific positions and orientations so that certain manufacturing steps can be carried out within required accuracies. Despite the importance of fixtures in the production of expensive devices at Sandia National Laboratories, there is little in-house expertise in mathematical design issues associated with fixtures. As a result, fixtures typically do not work as intended when they are first manufactured. Thus, an inefficient and expensive trial-and-error approach must be utilized. This design methodology adversely impacts important mission duties of Sandia National Laboratories, such as the production of neutron generators. The work performed under the support of this LDRD project took steps toward providing mechanical designers with software tools based on rigorous analytical techniques for dealing with fixture stability and tolerance stack-up.
The constitutive behavior of mechanical joints is largely responsible for the energy dissipation and vibration damping in weapons systems. For reasons arising from the dramatically different length scales associated with those dissipative mechanisms and the length scales characteristic of the overall structure, this physics cannot be captured through direct numerical simulation (DNS) of the contact mechanics within a structural dynamics analysis. The difficulties of DNS manifest themselves either in terms of Courant times that are orders of magnitude smaller than that necessary for structural dynamics analysis or as intractable conditioning problems. The only practical method for accommodating the nonlinear nature of joint mechanisms within structural dynamic analysis is through constitutive models employing degrees of freedom natural to the scale of structural dynamics. In this way, development of constitutive models for joint response is a prerequisite for a predictive structural dynamics capability. A four-parameter model, built on a framework developed by Iwan, is used to reproduce the qualitative and quantitative properties of lap-type joints. In the development presented here, the parameters are deduced by matching experimental values of energy dissipation in harmonic loading and values of the force necessary to initiate macro-slip. (These experiments can be performed on real hardware or virtually via fine-resolution, nonlinear quasi-static finite elements.) The resulting constitutive model can then be used to predict the force/displacement results from arbitrary load histories.
This report presents a detailed multi-methods comparison of the spatial errors associated with finite difference, finite element and finite volume semi-discretizations of the scalar advection-diffusion equation. The errors are reported in terms of non-dimensional phase and group speeds, discrete diffusivity, artificial diffusivity, and grid-induced anisotropy. It is demonstrated that Fourier analysis (aka von Neumann analysis) provides an automatic process for separating the spectral behavior of the discrete advective operator into its symmetric dissipative and skew-symmetric advective components. Further it is demonstrated that streamline upwind Petrov-Galerkin and its control-volume finite element analogue, streamline upwind control-volume, produce both an artificial diffusivity and an artificial phase speed in addition to the usual semi-discrete artifacts observed in the discrete phase speed, group speed and diffusivity. For each of the numerical methods considered, asymptotic truncation error and resolution estimates are presented for the limiting cases of pure advection and pure diffusion. The Galerkin finite element method and its streamline upwind derivatives are shown to exhibit super-convergent behavior in terms of phase and group speed when a consistent mass matrix is used in the formulation. In contrast, the CVFEM method and its streamline upwind derivatives yield strictly second-order behavior. While this work can only be considered a first step in a comprehensive multi-methods analysis and comparison, it serves to identify some of the relative strengths and weaknesses of multiple numerical methods in a common mathematical framework.