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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

A Framework for Performing Multiscale Stochastic Progressive Failure Analysis of Composite Structures

A framework is presented that enables coupled multiscale analysis of composite structures. The recently developed, free, Finite Element Analysis-Micromechanics Analysis Code (FEAMAC) software couples the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) with ABAQUS to perform micromechanics based FEA such that the nonlinear composite material response at each integration point is modeled at each increment by MAC/GMC. As a result, the stochastic nature of fiber breakage in composites can be simulated through incorporation of an appropriate damage and failure model that operates within MAC/GMC on the level of the fiber. Results are presented for the progressive failure analysis of a titanium matrix composite tensile specimen that illustrate the power and utility of the framework and address the techniques needed to model the statistical nature of the problem properly. In particular, it is shown that incorporating fiber strength randomness on multiple scales improves the quality of the simulation by enabling failure at locations other than those associated with structural level stress risers.

Bednarcyk, Brett A.↗
NASA NTRS1999

SiC/SiC Leading Edge Turbine Airfoil Tested Under Simulated Gas Turbine Conditions

Silicon-based ceramics have been proposed as component materials for use in gas turbine engine hot-sections. A high pressure burner rig was used to expose both a baseline metal airfoil and ceramic matrix composite leading edge airfoil to typical gas turbine conditions to comparatively evaluate the material response at high temperatures. To eliminate many of the concerns related to an entirely ceramic, rotating airfoil, this study has focused on equipping a stationary metal airfoil with a ceramic leading edge insert to demonstrate the feasibility and benefits of such a configuration. Here, the idea was to allow the SiC/SiC composite to be integrated as the airfoil's leading edge, operating in a "free-floating" or unrestrained manner. and provide temperature relief to the metal blade underneath. The test included cycling the airfoils between simulated idle, lift, and cruise flight conditions. In addition, the airfoils were air-cooled, uniquely instrumented, and exposed to the same internal and external conditions, which included gas temperatures in excess of 1370 C (2500 F). Results show the leading edge insert remained structurally intact after 200 simulated flight cycles with only a slightly oxidized surface. The instrumentation clearly suggested a significant reduction (approximately 600 F) in internal metal temperatures as a result of the ceramic leading edge. The object of this testing was to validate the design and analysis done by Materials Research and Design of Rosemont, PA and to determine the feasibility of this design for the intended application.

Robinson, R. Craig↗
NASA NTRS1989

Probalistic Finite Elements (PFEM) structural dynamics and fracture mechanics

The purpose of this work is to develop computationally efficient methodologies for assessing the effects of randomness in loads, material properties, and other aspects of a problem by a finite element analysis. The resulting group of methods is called probabilistic finite elements (PFEM). The overall objective of this work is to develop methodologies whereby the lifetime of a component can be predicted, accounting for the variability in the material and geometry of the component, the loads, and other aspects of the environment; and the range of response expected in a particular scenario can be presented to the analyst in addition to the response itself. Emphasis has been placed on methods which are not statistical in character; that is, they do not involve Monte Carlo simulations. The reason for this choice of direction is that Monte Carlo simulations of complex nonlinear response require a tremendous amount of computation. The focus of efforts so far has been on nonlinear structural dynamics. However, in the continuation of this project, emphasis will be shifted to probabilistic fracture mechanics so that the effect of randomness in crack geometry and material properties can be studied interactively with the effect of random load and environment.

Liu, Wing-Kam↗
DOE OSTI2018

Prediction of shock initiation thresholds and ignition probability of polymer-bonded explosives using mesoscale simulations

The design of new materials requires establishment of macroscopic measures of material performance as functions of microstructure. Traditionally, this process has been an empirical endeavor. In this paper, an approach to computationally predict the probabilistic ignition thresholds of polymer-bonded explosives (PBXs) using mesoscale simulations is developed. The simulations explicitly account for microstructure, constituent properties, and interfacial responses and capture processes responsible for the development of hotspots and damage. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, and heat conduction. The probabilistic analysis uses sets of statistically similar microstructure samples to directly mimic relevant experiments for quantification of statistical variations of material behavior due to inherent material heterogeneities. The particular thresholds and ignition probabilities predicted are expressed in James type and Walker–Wasley type relations, leading to the establishment of explicit analytical expressions for the ignition probability as function of loading. Specifically, the ignition thresholds corresponding to any given level of ignition probability and ignition probability maps are predicted for PBX 9404 for the loading regime of U p = 200–1200 m/s where U p is the particle speed. The predicted results are in good agreement with available experimental measurements. A parametric study also shows that binder properties can significantly affect the macroscopic ignition behavior of PBXs. Finally, the capability to computationally predict the macroscopic engineering material response relations out of material microstructures and basic constituent and interfacial properties lends itself to the design of new materials as well as the analysis of existing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Assessment of Multi-Step CADIS Parameter Sensitivity to Neutron Flux Spectrum in Shutdown Dose Rate Calculations

Shutdown dose rate (SDDR) analysis requires (a) a neutron transport calculation to estimate neutron flux fields, (b) an activation calculation to compute radionuclide inventories and associated photon sources, and (c) a photon transport calculation to estimate final dose rate. In some applications, accurate full-scale Monte Carlo (MC) SDDR simulations are needed for very large systems with massive amounts of shielding materials. For example, SDDR assessments are required throughout the biological shield (bio-shield) of the ITER experimental facility to evaluate the required waiting period after the shutdown of ITER and to identify the locations where human accessibility should be prohibited. The bio-shield is a large cylindrical concrete structure (30 m tall and 34 m in diameter) surrounding the very complex tokamak machine. Determining the effects on SDDR of important factors, such as the cross talk (interactions) between the different ports of ITER, is possible only through full-scale simulations that involve all the complex inner details of the ITER tokamak machine. These simulations are impractical because calculation of space- and energy-dependent neutron fluxes throughout the structural materials is needed to estimate the distribution of radioisotopes causing the SDDR. Determining a good importance function that can be used in biasing a neutron MC calculation for SDDR is complicated because it is difficult to explicitly express the SDDR response function, which depends on subsequent computational steps. The Multi-Step Consistent Adjoint Driven Importance Sampling (MS-CADIS) hybrid MC/deterministic method was developed to speed up the SDDR MC neutron transport calculation using an importance function that represents the neutron importance to the final SDDR. The MS-CADIS method uses the CADIS method, which has been Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for the United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan). successfully used for more than a decade in radiation shielding applications to generate consistent source biasing and weight window (WW) variance reduction parameters to accelerate continuous-energy MC simulations. However, because MS-CADIS focuses on multi-step shielding calculations such as SDDR calculations, it develops an importance function for the initial radiation transport calculation (e.g., the neutron calculations in SDDR simulations) that represents the importance of particles to the final response of the overall simulation. To develop this importance function, the MS-CADIS method uses an adjoint calculation with a specific neutron adjoint source that represents the neutron's contribution to the final SDDR. Brute force calculation of the MS-CADIS neutron adjoint source requires a deterministic neutron transport calculation coupled with activation calculations in each element of the deterministic mesh. These calculations are needed to determine the MS-CADIS parameter for the neutron interactions that produce radioisotopes contributing to the SDDR. The MS-CADIS parameter is the microscopic cross section of the radioisotope-producing interaction, multiplied by the mass of this radioisotope existing at the end of the irradiation and decay scenario, and divided by the steady-state rate of the interaction. Because the activation calculations are performed for each radionuclide in each mesh element, implementation of the MS-CADIS method can be computationally expensive. However, these neutron transport and activation calculations are not needed if we assume that the radioisotopes are created due to single neutron interactions and their concentrations do not change as a result of burnup by neutron irradiation. The use of this approximation, which renders the MS-CADIS parameter independent of the neutron flux, seems reasonable in calculating the MS-CADIS importance function because this function is needed only to speed up the MC calculation, not to calculate the final SDDR. In this work, we investigate the validity of this approximation by investigating the behavior of the MS-CADIS parameter as a function of the neutron flux spectra in the Fusion Evaluated Nuclear Data Library (FENDL) Neutronics and Shielding Benchmark problem that contains representative isotopes and flux environments of fusion energy systems. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2007

Ballistic Impact Response of Kevlar 49 and Zylon under Conditions Representing Jet Engine Fan Containment

A ballistic impact test program was conducted to provide validation data for the development of numerical models of blade out events in fabric containment systems. The impact response of two different fiber materials - Kevlar 49 (E.I. DuPont Nemours and Company) and Zylon AS (Toyobo Co., Ltd.) was studied by firing metal projectiles into dry woven fabric specimens using a gas gun. The shape, mass, orientation and velocity of the projectile were varied and recorded. In most cases the tests were designed such that the projectile would perforate the specimen, allowing measurement of the energy absorbed by the fabric. The results for both Zylon and Kevlar presented here represent a useful set of data for the purposes of establishing and validating numerical models for predicting the response of fabrics under conditions simulating those of a jet engine blade release situation. In addition some useful empirical observations were made regarding the effects of projectile orientation and the relative performance of the different materials.

Pereira, J. Michael↗
NASA NTRS2008

FAA Development of Reliable Modeling Methodologies for Fan Blade Out Containment Analysis: Ballistic Impact Testing - Part 2

This report summarizes the ballistic impact testing that was conducted to provide validation data for the development of numerical models of blade out events in fabric containment systems. The ballistic impact response of two different fiber materials - Kevlar 49 (E.I. DuPont Nemours and Company) and Zylon AS (Toyobo Co., Ltd.) was studied by firing metal projectiles into dry woven fabric specimens using a gas gun. The shape, mass, orientation and velocity of the projectile were varied and recorded. In most cases the tests were designed such that the projectile would perforate the specimen, allowing measurement of the energy absorbed by the fabric. The results for both Zylon and Kevlar presented here represent a useful set of data for the purposes of establishing and validating numerical models for predicting the response of fabrics under conditions simulating those of a jet engine blade release situations. In addition some useful empirical observations were made regarding the effects of projectile orientation and the relative performance of the different materials.

Revilock, Duane M.↗
DOE OSTI2025

Property optimized energy absorber for automotive bumpers utilizing multi-material and structural design strategies

This study proposes a novel design for automotive bumper using optimized lattice structures and multi-materials to balance low-speed collision and high-speed pedestrian impact performance. Different blends of 20 % carbon fiber-reinforced acrylonitrile butadiene styrene with thermoplastic polyurethane were used to tailor material properties. The energy absorber features lattice structures with customized mechanical responses, created by varying the incline angle θ from 0 to 180°. We conducted 576 finite element simulations on a half-scale model to optimize energy absorption and stiffness, leading to 66 optimized designs that met both low-speed and high-speed impact criteria. Two sub-scale optimized energy absorbers with different peak forces—both meeting low-speed impact requirements—were 3D printed and validated through drop-weight testing. The one with lower peak stress demonstrated a more compliant response, exhibiting approximately 90 % lower initial peak force and an increase in energy absorption of around 33 % (from 24 J to 32 J). Compared to the baseline triangular lattice, the optimized absorber increased energy absorption by 68 % from (19 J to 32 J) and reduced peak stress by 70 %. It also showed near-complete recovery with minimal fractures, making it suitable for repeated use. This design improves safety while offering a lightweight, durable, and cost-effective bumper system.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Intercalation-Induced Amorphization Boosts Aqueous Magnesium-Ion Storage

The design of aqueous battery cathode materials that can store divalent ions with high capacity and satisfactory reversibility is of great technical importance and challenge. Here, we report that divalent Mg 2+ storage is facilitated by an intercalation-induced amorphization of vanadate electrode materials. Electrokinetic analyses and in situ synchrotron X-ray diffraction and absorption spectroscopy collectively demonstrate that vanadate layered materials (Li–V 3 O 8 ) undergo a structural transformation to amorphization induced by Mg 2+ intercalation, and a reversible restoration of crystalline structure upon Mg 2+ deintercalation. Debye scattering simulations suggest that intercalation-induced turbostratic disorder, especially random rotations, translational shifts, oscillatory motions, or varied interlayer spacing of adjacent V–O molecular layers, could be responsible for the observed amorphization. The highly distorted local structure, in turn, facilitates Mg 2+ intercalation across the vanadate electrode materials, responsible for nearly 3/7 of the total Mg 2+ ions intercalated. The study presented reveals an intriguing relationship between ion transport and the reversible amorphization-to-crystallization dynamics it induces, opening a paradigm for designing advanced aqueous battery electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Generalized Backward Euler algorithm for the numerical integration of a viscous breakage model

This paper discusses the formulation and the numerical performance of a fully implicit algorithm used to integrate a rate-dependent model defined within a breakage mechanics framework. For this purpose, a Generalized Backward Euler (GBE) algorithm has been implemented according to two different linearization strategies: The former is derived by a direct linearization of the constitutive equations, while the latter introduces rate effects through a consistency parameter. Additionally, the accuracy and efficiency of the GBE algorithm have been investigated by (1) performing material point analyses and (2) solving initial boundary value problems. In both cases, the overall performance of the underlying algorithm is inspected for a range of loading rates, thus simulating comminution from slow to fast dynamic problems. As the viscous response of the breakage model can be recast through a viscous nucleus function, the presented algorithm can be considered as a general framework to integrate constitutive equations relying on the overstress approach typical of Perzyna-like viscoplastic models.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. Further, this new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Computational materials assessment of the D/Li-stripping neutron source as a prototypical facility for fusion materials testing

As the US fusion materials community awaits the selection and design of a fusion prototypical neutron source (FPNS), a risk reduction exercise has been conducted to (i) provide an updated materials performance evaluation using state-of-the-art computational materials modeling, (ii) expand on legacy analysis based on pure Fe to other relevant fusion structural materials types, and (iii) ensure that materials response under FPNS operational conditions is consistent with referential fusion reactor conditions. The current paper describes the efforts undertaken to assemble a comprehensive computational methodology that includes neutronics, primary damage calculations, atomistic simulations of displacement cascades, chemical inventory evolution calculations, and a computational thermodynamic analysis of emerging phases during irradiation. Our work extends existing studies in pure Fe to reduced-activation ferritic/martensitic steels, tungsten, silicon carbide, and vanadium alloys. We focus on the single-beam deuteron/lithium-stripping neutron source behind the IFMIF-DONES concept, which we assess against ITER, two DEMO designs, and an ideal pure 14-MeV flux. Our analysis indicates that, within standard uncertainties inherent to the models employed, the DONES concept adequately captures fusion conditions in the four materials analyzed. Our work is intended as a comprehensive irradiation damage analysis of fusion-representative neutron sources, to be used for further neutron source evaluation and fusion facility operation.

Marian, Jaime [Univ. of California, Los Angeles, C↗
DOE OSTI2025

LAROMance Grade 91 Model Integration in NEML2

New reactor designs are targeting higher operating temperatures for increased thermal efficiency when compared to the current fleet of light water reactors. Designing structural components for these high temperature environments with reliable long-term operations requires material models that can accurately capture the deformation mechanisms active in these environments. The LAROMance surrogate material models are based on a database of mechanistic crystal plasticity simulations for high-temperature conditions. Inputs to the LAROMance models reflect the microstructural pedigree of the material, like dislocation densities and precipitate contents. Based on the evolution of these microstructural features, the LAROMance model provides the engineering scale constitutive model response. The LAROMance model was recently parameterized for Grade 91, a high temperature alloy. In the present work, the Grade 91 LAROMance model is implemented in the New Material Model Library, version 2 (NEML2). NEML2 provides a modular way to build material models from smaller blocks and was developed to vectorize the material update to efficiently run on modern computational architectures with graphics processing unit accelerators. NEML2 constitutive models can be used in simulations based on the multiphysics object-oriented simulation environment (MOOSE). This report provides details on the implementation of the Grade 91 LAROMance model in NEML2 and its verification of engineering scale finite element simulations in MOOSE.

42 - ENGINEERING↗
DOE OSTI2016

Large-eddy simulation, fuel rod vibration and grid-to-rod fretting in pressurized water reactors

Grid-to-rod fretting (GTRF) in pressurized water reactors is a flow-induced vibration phenomenon that results in wear and fretting of the cladding material on fuel rods. GTRF is responsible for over 70% of the fuel failures in pressurized water reactors in the United States. Predicting the GTRF wear and concomitant interval between failures is important because of the large costs associated with reactor shutdown and replacement of fuel rod assemblies. The GTRF-induced wear process involves turbulent flow, mechanical vibration, tribology, and time-varying irradiated material properties in complex fuel assembly geometries. This paper presents a new approach for predicting GTRF induced fuel rod wear that uses high-resolution implicit large-eddy simulation to drive nonlinear transient dynamics computations. The GTRF fluid–structure problem is separated into the simulation of the turbulent flow field in the complex-geometry fuel-rod bundles using implicit large-eddy simulation, the calculation of statistics of the resulting fluctuating structural forces, and the nonlinear transient dynamics analysis of the fuel rod. Ultimately, the methods developed here, can be used, in conjunction with operational management, to improve reactor core designs in which fuel rod failures are minimized or potentially eliminated. Furthermore, robustness of the behavior of both the structural forces computed from the turbulent flow simulations and the results from the transient dynamics analyses highlight the progress made towards achieving a predictive simulation capability for the GTRF problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Transient Mechanical Test Development for Accident Tolerant Cladding Candidates to Simulate PCMI-like Conditions

This study investigates the transient mechanical response of accident tolerant fuel (ATF) cladding candidates to provide information on potential PCMI-induced deformation and failure of candidate cladding materials. The mechanical testing is applied t to Gen-I FeCrAl and triaxially-braided SiC fiber SiC matrix composite materials under loading conditions like the PCMI during a postulated design basis reactivity insertion accident. To achieve the specific loading conditions, a modified burst test instrument with accurate mechanical pulse control has been designed and constructed. The mechanical response of FeCrAl, SiC-SiC composite tubes were compared with the pre-hydrided zirconium alloy cladding. The specimens subjected to mechanical pulses were investigated in a strain-controlled manner deformation.

Cinbiz, Nedim↗
DOE OSTI2018

Bayesian calibration of strength parameters using hydrocode simulations of symmetric impact shock experiments of Al-5083

We report that predictive modeling of materials requires accurately parameterized constitutive models. Parameterizing models that describe dynamic strength and plasticity require experimentally probing materials in a variety of strain rate regimes. Some experimental protocols (e.g., plate impact) probe the constitutive response of a material using indirect measures such as free surface velocimetry. Manual efforts to parameterize constitutive models using indirect experimental measures often lead to non-unique optimizations without quantification of parameter uncertainty. This study uses a Bayesian statistical approach to find model parameters and to quantify the uncertainty of the resulting parameters. The technique is demonstrated by parameterizing the Johnson-Cook strength model for aluminum alloy 5083 by coupling hydrocode simulations and velocimetry measurements of a series of plate impact experiments. Simulation inputs and outputs are used to calibrate an emulator that mimics the outputs of the computationally intensive simulations. Varying the amount of experimental data available for emulator calibration showed clear differences in the degree of uncertainty and uniqueness of the resulting optimized Johnson-Cook parameters for Al-5083. The results of the optimization provided a numerical evaluation of the degree of confidence in model parameters and model performance. Lastly, given an understanding of the physical effects of certain model parameters, individual parameter uncertainty can be leveraged to quickly identify gaps in the physical domains covered by completed experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The origin of ultrahigh piezoelectricity in relaxor-ferroelectric solid solution crystals

The discovery of ultrahigh piezoelectricity in relaxor-ferroelectric solid solution single crystals is a breakthrough in ferroelectric materials. A key signature of relaxor-ferroelectric solid solutions is the existence of polar nanoregions, a nanoscale inhomogeneity, that coexist with normal ferroelectric domains. Despite two decades of extensive studies, the contribution of polar nanoregions to the underlying piezoelectric properties of relaxor ferroelectrics has yet to be established. Here we quantitatively characterize the contribution of polar nanoregions to the dielectric/piezoelectric responses of relaxor-ferroelectric crystals using a combination of cryogenic experiments and phase-field simulations. The contribution of polar nanoregions to the room-temperature dielectric and piezoelectric properties is in the range of 50–80%. A mesoscale mechanism is proposed to reveal the origin of the high piezoelectricity in relaxor ferroelectrics, where the polar nanoregions aligned in a ferroelectric matrix can facilitate polarization rotation. This mechanism emphasizes the critical role of local structure on the macroscopic properties of ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗