Magnetostrictive CoFe films were investigated for use as magnetoelastic tags or sensors. The ability to electrodeposit these films enables batch fabrication processes to pattern a variety of geometries while controlling the film stoichiometry and crystallography. In current research looking at CoFe, improved magnetostriction was achieved using a co-sputtering, annealing, and quenching method1. Other current research has reported electrodeposited CoFe films using a sulfate based chemistry resulting in film compositions that are Fe rich in the range of Co0.3-0.4Fe0.7-0.6 and have problems of codeposition of undesirables that can have a negative impact on magnetic properties. The research presented here focused on maximizing magnetostriction at the optimal stoichiometry range of Co0.7-0.75Fe0.3-0.25, targeting the (fcc+bcc)/bcc phase boundary, and using a novel chemistry and plating parameters to deposit films without being limited to “line of sight” deposition.
This document contains a description of the system architecture for the IDC Re-Engineering Phase 2 project. This is a draft version that primarily provides background information for understanding delivered Use Case Realizations.
The continued exponential growth of photovoltaic technologies paves a path to a solar-powered world, but requires continued progress toward low-cost, high-reliability, high-performance photovoltaic (PV) systems. High reliability is an essential element in achieving low-cost solar electricity by reducing operation and maintenance (O&M) costs and extending system lifetime and availability, but these attributes are difficult to verify at the time of installation. Utilities, financiers, homeowners, and planners are demanding this information in order to evaluate their financial risk as a prerequisite to large investments. Reliability research and development (R&D) is needed to build market confidence by improving product reliability and by improving predictions of system availability, O&M cost, and lifetime. This project is focused on understanding, predicting, and improving the reliability of PV systems. The two areas being pursued include PV arc-fault and ground fault issues, and inverter reliability.
A moment-of-fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment-of-fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi-Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block-structured adaptive grid. The new method is applied to multiphase problems illustrating contact-line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two-dimensional, three-dimensional axisymmetric (R-Z), and three-dimensional (X-Y-Z) coordinate systems.
We report on experiments demonstrating the transition from thermally-dominated K-shell line emission to non-thermal, hot-electron-driven inner-shell emission for z pinch plasmas on the Z machine. While x-ray yields from thermal K-shell emission decrease rapidly with increasing atomic number Z, we find that non-thermal emission persists with favorable Z scaling, dominating over thermal emission for Z=42 and higher (hn ≥ 17keV). Initial experiments with Mo (Z=42) and Ag (Z=47) have produced kJ-level emission in the 17-keV and 22-keV Kα lines respectively. We will discuss the electron beam properties that could excite these non - thermal lines. We also report on experiments that have attempted to control non - thermal K - shell line emission by modifying the wire array or load hardware setup.
The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the art of local models with the flexibility and accuracy of the nonlocal peridynamic model. In the mixed locality method this coupling occurs across scales, so that the nonlocal model can be used to communicate material heterogeneity at scales inappropriate to local partial differential equation models. Additionally, the computational burden of the weak form of the peridynamic model is reduced dramatically by only requiring that the model be solved on local patches of the simulation domain which may be computed in parallel, taking advantage of the heterogeneous nature of next generation computing platforms. Addition- ally, we present a novel Galerkin framework, the 'Ambulant Galerkin Method', which represents a first step towards a unified mathematical analysis of local and nonlocal multiscale finite element methods, and whose future extension will allow the analysis of multiscale finite element methods that mix models across scales under certain assumptions of the consistency of those models.
One of the more severe environments for a store on an aircraft is during the ejection of the store. During this environment it is not possible to instrument all component responses, and it is also likely that some instruments may fail during the environment testing. This work provides a method for developing these responses from failed gages and uninstrumented locations. First, the forces observed by the store during the environment are reconstructed. A simple sampling method is used to reconstruct these forces given various parameters. Then, these forces are applied to a model to generate the component responses. Validation is performed on this methodology.
Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Hand book, DOE-HDBK-3010, Airborne Release Fractions/Rates and Resp irable Fractions for Nonreactor Nuclear Facilities , to determine source terms. In calcula ting source terms, analysts tend to use the DOE Handbook's bounding values on airbor ne release fractions (ARFs) and respirable fractions (RFs) for various cat egories of insults (representing potential accident release categories). This is typica lly due to both time constraints and the avoidance of regulatory critique. Unfort unately, these bounding ARFs/RFs represent extremely conservative values. Moreover, th ey were derived from very limited small- scale table-top and bench/labo ratory experiments and/or fr om engineered judgment. Thus the basis for the data may not be re presentative to the actual unique accident conditions and configura tions being evaluated. The goal of this res earch is to develop a more ac curate method to identify bounding values for the DOE Handbook using the st ate-of-art multi-physics-based high performance computer codes. This enable s us to better understand the fundamental physics and phenomena associated with the ty pes of accidents for the data described in it. This research has examined two of the DOE Handbook's liquid fire experiments to substantiate the airborne release frac tion data. We found th at additional physical phenomena (i.e., resuspension) need to be included to derive bounding values. For the specific cases of solid powder under pre ssurized condition and mechanical insult conditions the codes demonstrated that we can simulate the phenomena. This work thus provides a low-cost method to establis h physics-justified sa fety bounds by taking into account specific geometri es and conditions that may not have been previously measured and/or are too costly to do so.
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been de- signed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: Capability to solve extremely large circuit problems by supporting large-scale parallel com- puting platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. Device models that are specifically tailored to meet Sandia's needs, including some radiation- aware devices (for Sandia users only). Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase -- a message passing parallel implementation -- which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows. Trademarks The information herein is subject to change without notice. Copyright c 2002-2015 Sandia Corporation. All rights reserved. Xyce TM Electronic Simulator and Xyce TM are trademarks of Sandia Corporation. Portions of the Xyce TM code are: Copyright c 2002, The Regents of the University of California. Produced at the Lawrence Livermore National Laboratory. Written by Alan Hindmarsh, Allan Taylor, Radu Serban. UCRL-CODE-2002-59 All rights reserved. Orcad, Orcad Capture, PSpice and Probe are registered trademarks of Cadence Design Systems, Inc. Microsoft, Windows and Windows 7 are registered trademarks of Microsoft Corporation. Medici, DaVinci and Taurus are registered trademarks of Synopsys Corporation. Amtec and TecPlot are trademarks of Amtec Engineering, Inc. Xyce 's expression library is based on that inside Spice 3F5 developed by the EECS Department at the University of California. The EKV3 MOSFET model was developed by the EKV Team of the Electronics Laboratory-TUC of the Technical University of Crete. All other trademarks are property of their respective owners. Contacts Bug Reports (Sandia only) http://joseki.sandia.gov/bugzilla http://charleston.sandia.gov/bugzilla World Wide Web http://xyce.sandia.gov http://charleston.sandia.gov/xyce (Sandia only) Email xyce@sandia.gov (outside Sandia) xyce-sandia@sandia.gov (Sandia only)
IEEE Transactions on Parallel and Distributed Systems
Shan, Tzu-Ray; Aktulga, Hasan M.; Knight, Chris; Coffman, Paul; Jiang, Wei
Hybrid parallelism allows high performance computing applications to better leverage the increasing on-node parallelism of modern supercomputers. In this paper, we present a hybrid parallel implementation of the widely used LAMMPS/ReaxC package, where the construction of bonded and nonbonded lists and evaluation of complex ReaxFF interactions are implemented efficiently using OpenMP parallelism. Additionally, the performance of the QEq charge equilibration scheme is examined and a dual-solver is implemented. We present the performance of the resulting ReaxC-OMP package on a state-of-the-art multi-core architecture Mira, an IBM BlueGene/Q supercomputer. For system sizes ranging from 32 thousand to 16.6 million particles, speedups in the range of 1.5-4.5x are observed using the new ReaxC-OMP software. Sustained performance improvements have been observed for up to 262,144 cores (1,048,576 processes) of Mira with a weak scaling efficiency of 91.5% in larger simulations containing 16.6 million particles.
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide [1] . The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce . This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide [1] . Trademarks The information herein is subject to change without notice. Copyright c 2002-2015 Sandia Corporation. All rights reserved. Xyce TM Electronic Simulator and Xyce TM are trademarks of Sandia Corporation. Portions of the Xyce TM code are: Copyright c 2002, The Regents of the University of California. Produced at the Lawrence Livermore National Laboratory. Written by Alan Hindmarsh, Allan Taylor, Radu Serban. UCRL-CODE-2002-59 All rights reserved. Orcad, Orcad Capture, PSpice and Probe are registered trademarks of Cadence Design Systems, Inc. Microsoft, Windows and Windows 7 are registered trademarks of Microsoft Corporation. Medici, DaVinci and Taurus are registered trademarks of Synopsys Corporation. Amtec and TecPlot are trademarks of Amtec Engineering, Inc. Xyce 's expression library is based on that inside Spice 3F5 developed by the EECS Department at the University of California. The EKV3 MOSFET model was developed by the EKV Team of the Electronics Laboratory-TUC of the Technical University of Crete. All other trademarks are property of their respective owners. Contacts Bug Reports (Sandia only) http://joseki.sandia.gov/bugzilla http://charleston.sandia.gov/bugzilla World Wide Web http://xyce.sandia.gov http://charleston.sandia.gov/xyce (Sandia only) Email xyce@sandia.gov (outside Sandia) xyce-sandia@sandia.gov (Sandia only)
Subwavelength-thin metasurfaces have shown great promises for the control of optical wavefronts, thus opening new pathways for the development of efficient flat optics. In particular, Huygens’ metasurfaces based on all-dielectric resonant meta-atoms have already shown a huge potential for practical applications with their polarization insensitivity and high transmittance efficiency. Here, we experimentally demonstrate a holographic Huygens’ metasurface based on dielectric resonant meta-atoms capable of complex wavefront control at telecom wavelengths. Our metasurface produces a hologram image in the far-field with 82% transmittance efficiency and 40% imaging efficiency. Such efficient complex wavefront control shows that Huygens’ metasurfaces based on resonant dielectric meta-atoms are a big step towards practical applications of metasurfaces in wavefront design related technologies, including computer-generated holograms, ultra-thin optics, security and data storage devices.
Efforts are being pursued to develop and qualify a system-level model of a reactor core isolation (RCIC) steam-turbine-driven pump. The model is being developed with the intent of employing it to inform the design of experimental configurations for full-scale RCIC testing. The model is expected to be especially valuable in sizing equipment needed in the testing. An additional intent is to use the model in understanding more fully how RCIC apparently managed to operate far removed from its design envelope in the Fukushima Daiichi Unit 2 accident. RCIC modeling is proceeding along two avenues that are expected to complement each other well. The first avenue is the continued development of the system-level RCIC model that will serve in simulating a full reactor system or full experimental configuration of which a RCIC system is part. The model reasonably represents a RCIC system today, especially given design operating conditions, but lacks specifics that are likely important in representing the off-design conditions a RCIC system might experience in an emergency situation such as a loss of all electrical power. A known specific lacking in the system model, for example, is the efficiency at which a flashing slug of water (as opposed to a concentrated jet of steam) could propel the rotating drive wheel of a RCIC turbine. To address this specific, the second avenue is being pursued wherein computational fluid dynamics (CFD) analyses of such a jet are being carried out. The results of the CFD analyses will thus complement and inform the system modeling. The system modeling will, in turn, complement the CFD analysis by providing the system information needed to impose appropriate boundary conditions on the CFD simulations. The system model will be used to inform the selection of configurations and equipment best suitable of supporting planned RCIC experimental testing. Preliminary investigations with the RCIC model indicate that liquid water ingestion by the turbine decreases the developed turbine torque; the RCIC speed then slows, and thus the pump flow rate to the RPV decreases. Subsequently, RPV water level decreases due to continued boiling and the liquid fraction flowing to the RCIC decreases, thereby accelerating the RCIC and refilling the RPV. The feedback cycle then repeats itself and/or reaches a quasi-steady equilibrium condition. In other words, the water carry-over is limited by cyclic RCIC performance degradation, and hence the system becomes self-regulating. The indications achieved to date with the system model are more qualitative than quantitative. The avenues being pursued to increase the fidelity of the model are expected to add quantitative realism. The end product will be generic in the sense that the RCIC model will be incorporable within the larger reactor coolant system model of any nuclear power plant or experimental configuration.
Turbulent fluctuation behavior is approximately modeled using a pulsatile flow model analogy.. This model follows as an extension to the turbulent laminar sublayer model developed by Sternberg (1962) to be valid for a fully turbulent flow domain. Here unsteady turbulent behavior is modeled via a sinusoidal pulsatile approach. While the individual modes of the turbulent flow fluctuation behavior are rather crudely modeled, approximate temporal integration yields plausible estimates for Root Mean Square (RMS) velocity fluctuations. RMS pressure fluctuations and spectra are of particular interest and are estimated via the pressure Poisson expression. Both RMS and Power Spectral Density (PSD), i.e. spectra are developed. Comparison with available measurements suggests reasonable agreement. An additional fluctuating quantity, i.e. RMS wall shear fluctuation is also estimated, yielding reasonable agreement with measurement.
We examine the role of periodic sinusoidal free-stream disturbances on the inner law law-of-the-wall (log-law) for turbulent boundary layers. This model serves a surrogate for the interaction of flight vehicles with atmospheric disturbances. The approximate skin friction expression that is derived suggests that free-stream disturbances can cause enhancement of the mean skin friction. Considering the influence of grid generated free stream turbulence in the laminar sublayer/log law region (small scale/high frequency) the model recovers the well-known shear layer enhancement suggesting an overall validity for the approach. The effect on the wall shear associated with the lower frequency due to the passage of the vehicle through large (vehicle scale) atmospheric disturbances is likely small i.e. on the order 1% increase for turbulence intensities on the order of 2%. The increase in wall pressure fluctuation which is directly proportional to the wall shear stress is correspondingly small.
This initial draft document contains formative data model content for select areas of Re-Engineering Phase 2 IDC System. The purpose of this document is to facilitate discussion among the stakeholders. It is not intended as a definitive proposal.
Flywheel energy storage systems are in use globally in increasing numbers . No codes pertaining specifically to flywheel energy storage exist. A number of industrial incidents have occurred. This protocol recommends a technical basis for safe flywheel de sign and operation for consideration by flywheel developers, users of flywheel systems and standards setting organizations.
The objective of this project was to devolop a plan for modifying the Turbulent Combustion Laboratory (TCL) with the necessary infrastructure to produce a cold (near liquid temperature) hydrogen jet. The necessary infrastructure has been specified and laboratory modifications are currently underway. Once complete, experiments from this platform will be used to develop and validate models that inform codes and standards which specify protection criteria for unintended releases from liquid hydrogen storage, transport, and delivery infrastructure.
This report analyzes the permeation resistance of a novel and proprietary polymer coating for hydrogen isotope resistance that was developed by New Mexico State University. Thermal gravimetric analysis and thermal desoprtion spectroscopy show the polymer is stable thermally to approximately 250 deg C. Deuterium gas-driven permeation experiments were conducted at Sandia to explore early evidence (obtained using Brunauer - Emmett - Teller) of the polymer's strong resistance to hydrogen. With a relatively small amount of the polymer in solution (0.15%), a decrease in diffusion by a factor of 2 is observed at 100 and 150 deg C. While there was very little reduction in permeability, the preliminary findings reported here are meant to demonstrate the sensitivity of Sandia's permeation measurements and are intended to motivate the future exploration of thicker barriers with greater polymer coverage.
With rising grid interconnections of solar photovoltaic (PV) systems, greater attention is being trained on lifecycle performance, reliability, and project economics. Expected to meet production thresholds over a 20-30 year timeframe, PV plants require a steady diet of operations and maintenance (O&M) oversight to meet contractual terms. However, industry best practices are only just beginning to emerge, and O&M budgets—given the arrangement of the solar project value chain—appear to vary widely. Based on insights from in-depth interviews and survey research, this paper presents an overview of the utility-scale PV O&M budgeting process along with guiding rationales, before detailing perspectives on current plant upkeep activities and price points largely in the U.S. It concludes by pondering potential opportunities for improving upon existing O&M budgeting approaches in ways that can benefi t the industry at-large.
In this report, we present a multi-scale computational model to simulate plastic deformation of tantalum and validating experiments. In atomistic/ dislocation level, dislocation kink- pair theory is used to formulate temperature and strain rate dependent constitutive equations. The kink-pair theory is calibrated to available data from single crystal experiments to produce accurate and convenient constitutive laws. The model is then implemented into a BCC crystal plasticity finite element method (CP-FEM) model to predict temperature and strain rate dependent yield stresses of single and polycrystalline tantalum and compared with existing experimental data from the literature. Furthermore, classical continuum constitutive models describing temperature and strain rate dependent flow behaviors are fit to the yield stresses obtained from the CP-FEM polycrystal predictions. The model is then used to conduct hydro- dynamic simulations of Taylor cylinder impact test and compared with experiments. In order to validate the proposed tantalum CP-FEM model with experiments, we introduce a method for quantitative comparison of CP-FEM models with various experimental techniques. To mitigate the effects of unknown subsurface microstructure, tantalum tensile specimens with a pseudo-two-dimensional grain structure and grain sizes on the order of millimeters are used. A technique combining an electron back scatter diffraction (EBSD) and high resolution digital image correlation (HR-DIC) is used to measure the texture and sub-grain strain fields upon uniaxial tensile loading at various applied strains. Deformed specimens are also analyzed with optical profilometry measurements to obtain out-of- plane strain fields. These high resolution measurements are directly compared with large-scale CP-FEM predictions. This computational method directly links fundamental dislocation physics to plastic deformations in the grain-scale and to the engineering-scale applications. Furthermore, direct and quantitative comparisons between experimental measurements and simulation show that the proposed model accurately captures plasticity in deformation of polycrystalline tantalum.
Final report for Cognitive Computing for Security LDRD 165613. It reports on the development of hybrid of general purpose/ne uromorphic computer architecture, with an emphasis on potential implementation with memristors.
Assembled mechanical systems often contain a large number of bolted connections. These bolted connections (joints) are integral aspects of the load path for structural dynamics, and, consequently, are paramount for calculating a structure's stiffness and energy dissipation prop- erties. However, analysts have not found the optimal method to model appropriately these bolted joints. The complexity of the screw geometry causes issues when generating a mesh of the model. This report will explore different approaches to model a screw-substrate connec- tion. Model parameters such as mesh continuity, node alignment, wedge angles, and thread to body element size ratios are examined. The results of this study will give analysts a better understanding of the influences of these parameters and will aide in finding the optimal method to model bolted connections.
This project was intended to enable SNL-CA to produce appropriate specimens of relevant stainless steels for testing and perform baseline testing of weld heat-affected zone and weld fusion zone. One of the key deliverables in this project was to establish a procedure for fracture testing stainless steel weld fusion zone and heat affected zones that were pre-charged with hydrogen. Following the establishment of the procedure, a round robin was planned between SNL-CA and SRNL to ensure testing consistency between laboratories. SNL-CA and SRNL would then develop a comprehensive test plan, which would include tritium exposures of several years at SRNL on samples delivered by SNL-CA. Testing would follow the procedures developed at SNL-CA. SRNL will also purchase tritium charging vessels to perform the tritium exposures. Although comprehensive understanding of isotope-induced fracture in GTS reservoir materials is a several year effort, the FY15 work would enabled us to jump-start the tests and initiate long-term tritium exposures to aid comprehensive future investigations. Development of a procedure and laboratory testing consistency between SNL-CA and SNRL ensures reliability in results as future evaluations are performed on aluminum alloys and potentially additively-manufactured components.
The XVis project brings together the key elements of research to enable scientific discovery at extreme scale. Scientific computing will no longer be purely about how fast computations can be performed. Energy constraints, processor changes, and I/O limitations necessitate significant changes in both the software applications used in scientific computation and the ways in which scientists use them. Components for modeling, simulation, analysis, and visualization must work together in a computational ecosystem, rather than working independently as they have in the past. This project provides the necessary research and infrastructure for scientific discovery in this new computational ecosystem by addressing four interlocking challenges: emerging processor technology, in situ integration, usability, and proxy analysis.
This report documents work done for the ITER International Fusion Energy Organization (Sponsor) under a Funds-In Agreement FI 011140916 with Sandia National Laboratories. The work consists of preparing and analyzing samples for an experiment to measure material erosion and deposition in the EAST Tokamak. Sample preparation consisted of depositing thin films of carbon and aluminum onto molybdenum tiles. Analysis consists of measuring the thickness of films before and after exposure to helium plasma in EAST. From these measurements the net erosion and deposition of material will be quantified. Film thickness measurements are made at the Sandia Ion Beam Laboratory using Rutherford backscattering spectrometry and nuclear reaction analysis, as described in this report. This report describes the film deposition and pre-exposure analysis. Results from analysis after plasma exposure will be given in a subsequent report.
A series of quasi - static indentation experiments are conducted on carbon fiber reinforced polymer laminates with a systematic variation of thicknesses and fixture boundary conditions. Different deformation mechanisms and their resulting damage mechanisms are activated b y changing the thickn ess and boundary conditions. The quasi - static indentation experiments have been shown to achieve damage mechanisms similar to impact and penetration, however without strain rate effects. The low rate allows for the detailed analysis on the load response. Moreover, interrupted tests allow for the incremental analysis of various damage mechanisms and pr ogressions. The experimentally tested specimens are non - destructively evaluated (NDE) with optical imaging, ultrasonics and computed tomography. The load displacement responses and the NDE are then utilized in numerical simulations for the purpose of model validation and vetting. The accompanying numerical simulation work serves two purposes. First, the results further reveal the time sequence of events and the meaning behind load dro ps not clear from NDE . Second, the simulations demonstrate insufficiencies in the code and can then direct future efforts for development.
Ensuring Real Applications perform well on Trinity is key to success. Four components: ASC applications, Sustained System Performance (SSP), Extra-Large MiniApplications problems, and Micro-benchmarks.
