Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Tailoring Mid–Gap States of Chalcogenide Glass by Pressure–Induced Hypervalent Bonding Towards the Design of Electrical Switching Materials

Phase change memory (PCM) and ovonic threshold switching (OTS) materials using chalcogenide glass are essential elements in advanced 3D memory chips. The mid–gap states, induced by the disorder and defects in the glass, are the physical mechanisms of the electrical switching behavior, while the origin of these trap states is still under debate and the medium–range clusters that break the global octet rule, such as over–coordinated atoms, are known to be responsible in various glass. Here, it is discovered that a large fraction of over–coordinated clusters fails to generate mid–gap states, which are probably caused by hypervalent bonding, a multi–centered covalent bond participated by delocalized lone–pair electrons. This is confirmed by the pressure–driven simulations of amorphous GeSe models, in which it is found that octahedral motifs and hypervalent bonds prevent the over–coordinated medium–range clusters from providing excessive electrons. In practical applications, compatible dopants can be used to change the number of hypervalent bonds, thus controlling the number of mid–gap states and consequently the performance of PCM and OTS materials. Finally, these results reveal the origin of mid–gap states in chalcogenide glasses, enabling extensive control in the development of pioneering electrical switching materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tailoring material microstructure and property in wire-laser directed energy deposition through a wiggle deposition strategy

Developing effective strategies to directly control material microstructure, property, and anisotropy is an active research area in metal additive manufacturing. This work develops a wiggle deposition pattern for wire-laser directed energy deposition (DED) of 316L stainless steel (SS) to modify the solidification texture, particularly in the building direction, in as-deposited samples. Through multi-physics simulation, operando near-infrared imaging, and synchrotron x-ray characterization, it is found that the wiggle deposition strategy induces highly dynamic melt flow and oscillating thermal gradient in the melt pool, which is responsible for the variation of preferable grain growth direction and crystallographic texture in the sample. The specific texture reduces the anisotropy in the tensile strength of as-printed 316L SS samples cut along different directions. Also, it largely increases the ductility along the build direction. Crystal plasticity simulation is performed to correlate the sample texture with mechanical property. In conclusion, this work offers a unique approach for tailoring local properties through the control of melt pool instability by applying different tool paths.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings

The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.

58 GEOSCIENCES↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantifying the Impacts of Grain-scale Heterogeneity on Mechanical Response

Driven by the exceedingly high computational demands of simulating mechanical response in complex engineered systems with finely resolved finite element models, there is a critical need to optimally reduce the fidelity of such simulations. The minimum required fidelity is constrained by error tolerances on the simulation results, but error bounds are often impossible to obtain a priori. One such source of error is the variability of material properties within a body due to spatially non-uniform processing conditions and inherent stochasticity in material microstructure. This study seeks to quantify the effects of microstructural heterogeneity on component- and system-scale performance to aid in the choice of an appropriate material model and spatial resolution for finite element analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of interlayer confinement and hydration on capacitive charge storage in birnessite

Nanostructured birnessite exhibits high specific capacitance and nearly ideal capacitive behaviour in aqueous electrolytes, rendering it an important electrode material for low-cost, high-power energy storage devices. The mechanism of electrochemical capacitance in birnessite has been described as both Faradaic (involving redox) and non-Faradaic (involving only electrostatic interactions). To clarify the capacitive mechanism, we characterized birnessite’s response to applied potential using ex situ X-ray diffraction, electrochemical quartz crystal microbalance, in situ Raman spectroscopy and operando atomic force microscope dilatometry to provide a holistic understanding of its structural, gravimetric and mechanical responses. These observations are supported by atomic-scale simulations using density functional theory for the cation-intercalated structure of birnessite, ReaxFF reactive force field-based molecular dynamics and ReaxFF-based grand canonical Monte Carlo simulations on the dynamics at the birnessite–water–electrolyte interface. We show that capacitive charge storage in birnessite is governed by interlayer cation intercalation. We conclude that the intercalation appears capacitive due to the presence of nanoconfined interlayer structural water, which mediates the interaction between the intercalated cation and the birnessite host and leads to minimal structural changes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Detailed Characterization of CZT Detector Response for Improved Coded-Aperture Imaging Performance

Gamma-ray imaging is a powerful method for locating and quantifying sources of radiation. The coded-aperture technique demonstrates superior angular resolution in comparison to other methods (e.g., Compton reconstruction). In this method, a mask constructed of highly attenuating material encodes the scene as a shadow pattern on a position-sensitive detector; this pattern can then be used to recreate the origin(s) of incident radiation. This is typically done through convolution of the mask and shadow patterns. Iterative methods which attempt to reconstruct the observed shadow pattern using a weighted combination of simulated patterns may also be employed. In either case, errors in event position reconstruction due to detector imperfections alter the shadow pattern and will therefore degrade system performance and may introduce imaging artifacts. These effects can be mitigated with a detailed understanding of such errors – allowing for the generation of representative simulations that include the errors and/or correction of raw imager data to remove the errors. We present a calibration process for a commercially available cadmium zinc telluride (CZT) gamma imager which provides a comprehensive characterization of the spatial and energy dependence of event reconstruction. By illuminating a mask featuring a regular grid of pinholes with a calibration source, the localized response of the detector can be measured with fine granularity. These local responses are combined to generate a full detector response map which can be used to distort simulations in a manner that is representative of the observed detector data. Details of the calibration procedure and an assessment of the impact of its end products on the performance of iterative imaging methods will be presented.

Ziock, Klaus-Peter↗
DOE OSTI2023

Autonomous Polycrystalline Material Decomposition For Hyperspectral Neutron Tomography

Hyperspectral neutron tomography is an effective method for analyzing crystalline material samples with complex compositions in a non-destructive manner. Since the counts in the hyperspectral neutron radiographs directly depend on the neutron cross-sections, materials may exhibit contrasting neutron responses across wavelengths. Therefore, it is possible to extract the unique signatures associated with each material and use them to separate the crystalline phases simultaneously.We introduce an autonomous material decomposition (AMD) algorithm to automatically characterize and localize polycrystalline structures using Bragg edges with contrasting neutron responses from hyperspectral data. The algorithm estimates the linear attenuation coefficient spectra from the measured radiographs and then uses these spectra to perform polycrystalline material decomposition and reconstructs 3D material volumes to localize materials in the spatial domain. Our results demonstrate that the method can accurately estimate both the linear attenuation coefficient spectra and associated reconstructions on both simulated and experimental neutron data.

Samin nur chowdhury, Mohammad↗
DOE OSTI2021

Dynamic Strength of Copper at High Pressures Using Pressure Shear Plate Experiments

Measurement of strength at high pressures is critical to understanding the behavior of materials subjected to extreme loading conditions. Recent developments have extended the pressure shear plate impact (PSPI) technique to the high-pressure regime. Such experiments provide a unique opportunity to extract the complete stress-strain behavior of materials at relatively high pressures. The modified technique includes a new fiber-optic heterodyne transverse velocity interferometer system and new analysis methods to account for the inelastic response of the anvil plates. In this study, PSPI experiments are conducted on oxygen-free high conductivity (OFHC) copper at pressures ranging from 10 to 43 GPa and at strain rates of ~10 5 s -1 . Complete stress-strain curves of copper are obtained using a hybrid methodology that relies on numerical simulations and particle velocity records gathered in these experiments. Experiments are conducted at similar pressures and strain rates using two types of anvils that undergo different levels of plastic strains: tool steel and tungsten carbide, to check the robustness of the hybrid method. The behavior of copper at these high pressures and strain rates is compared with previous literature data, and the effect of pressure and strain rate on the strength behavior of copper are discussed. Here, it was observed that the rate of change in yield strength change in pressure was 2.5 times the rate of change of shear modulus with pressure. The strength models that predict the response of copper at these extreme conditions are also examined and validated.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Mechanisms of local stress sensing in multifunctional polymer films using fluorescent tetrapod nanocrystals

Nanoscale stress-sensing can be used across fields ranging from detection of incipient cracks in structural mechanics to monitoring forces in biological tissues. We demonstrate how tetrapod quantum dots (tQDs) embedded in block-copolymers act as sensors of tensile/compressive stress. Remarkably, tQDs can detect their own composite dispersion and mechanical properties, with a switch in optomechanical response when tQDs are in direct contact. Using experimental characterizations, atomistic simulations and finite-element analyses, we show that under tensile stress, densely-packed tQDs exhibit a photoluminescence peak shifted to higher energies (“blue-shift”) due to volumetric compressive stress in their core; loosely-packed tQDs exhibit a peak shifted to lower energies (“red-shift”) from tensile stress in the core. The stress-shifts result from the tQD's unique branched morphology in which the CdS arms act as antennas that amplify the stress in the CdSe core. Our nanocomposites exhibit excellent cyclability and scalability with no degraded properties of the host polymer. Colloidal tQDs allow sensing in many materials to potentially enable auto-responsive, smart structural nanocomposites that self-predict impending fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Approach and rotation of reconnecting topological defect lines in liquid crystal

Topological defects are a universal concept across many disciplines, such as crystallography, liquid-crystalline physics, low-temperature physics, cosmology, and even biology. In nematic liquid crystals, topological defects called disclinations have been widely studied. For their three-dimensional (3D) dynamics, however, only recently have theoretical approaches dealing with fully 3D configurations been reported. Further, recent experiments have observed 3D disclination line reconnections, a phenomenon characteristic of defect line dynamics, but detailed discussions were limited to the case of approximately parallel defects. In this paper, we focus on the case of two disclination lines that approach at finite angles and lie in separate planes, a more fundamentally 3D reconnection configuration. Observing and analyzing such reconnection events, we find the square-root law of the distance between disclinations and the decrease of the interdisclination angle over time. We compare the experimental results with theory and find qualitative agreement on the scaling of distance and angle with time, but quantitative disagreement on distance and angle relative mobilities. To probe this disagreement, we derive the equations of motion for systems with reduced twist constant and also carry out simulations for this case. These, together with the experimental results, suggest that deformations of disclinations may be responsible for the disagreement. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2019

Summary of Grizzly Development for Advanced Reactor Structural Materials

This summary report contains an overview of work performed during Fiscal Year 2019 under the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) Fuels Product Line work package entitled "MS-19IN020106 - GRIZZLY - INL". Grizzly is a code based on the MOOSE framework that has been under development for several years funded by the US Department of Energy’s Light Water Reactor Sustainability (LWRS) program. Prior to 2019, Grizzly development was focused on development of capabilities needed to simulate the effects of aging mechanisms in light water reactor (LWR) components. The goal of this LWR-focused work has been to develop tools that can be used to assess the ability of existing LWRs to perform safety in long term operation (LTO) scenarios. This LWR-focused work has primarily focused on the effects of reactor pressure vessel (RPV) embrittlement and aging of reinforced concrete structures. Grizzly development to support LWR degradation needs is still ongoing, but starting in 2019, Grizzly’s capabilities are being expanded to also meet growing needs for a tool to address structural component integrity and degradation in advanced reactors. A simulation capability for advanced reactor components has many of the same basic requirements that one for LWRs does, such as the ability to solve coupled systems of partial differential equations that arise from the fact that many of the problems of interest involve multiple coupled physics. The materials issues also inherently involve multiple length scales, with the response at the engineering scale being strongly affected by mechanisms occurring at atomistic and mesoscopic scales.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Reliability Barriers to Enhanced Solar PV Deployment: Selected Research Findings - State Outreach Technical Sessions

This presentation provides summaries of three reports selected from a body of work that examines potential reliability challenges related to planning future power systems with higher solar photovoltaic (PV) penetrations in the western United States. The three topics covered are 1) Power system flexibility requirements and supply; 2) Resource adequacy and the capacity credit of solar; 3) Simulating distributed energy resource responses to transmission system-level faults considering IEEE 1547 performance categories on three major WECC transmission paths. The full reports and related resources are available at the project website, https://www.westernenergyboard.org/western-interstate-energy-board/barrier-mitigation-to-enhanced-distributed-solar-photovoltaic/. This presentation is targeted toward a technical audience. A companion policy presentation provides higher-level summaries of the same material.

distributed energy resource↗