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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Improving dynamic collision frequencies: Impacts on dynamic structure factors and stopping powers in warm dense matter

Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, in this work, we systematically investigate variations in dynamic electron–ion collision frequencies ν(ω) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν(ω). These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 6061-T651 Aluminum Alloy from Plate

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Sand Impact Tests of a Half-Scale Crew Module Boilerplate Test Article

Although the Orion Multi-Purpose Crew Vehicle (MPCV) is being designed primarily for water landings, a further investigation of launch abort scenarios reveals the possibility of an onshore landing at Kennedy Space Center (KSC). To gather data for correlation against simulations of beach landing impacts, a series of sand impact tests were conducted at NASA Langley Research Center (LaRC). Both vertical drop tests and swing tests with combined vertical and horizontal velocity were performed onto beds of common construction-grade sand using a geometrically scaled crew module boilerplate test article. The tests were simulated using the explicit, nonlinear, transient dynamic finite element code LS-DYNA. The material models for the sand utilized in the simulations were based on tests of sand specimens. Although the LSDYNA models provided reasonable predictions for peak accelerations, they were not always able to track the response through the duration of the impact. Further improvements to the material model used for the sand were identified based on results from the sand specimen tests.

Vassilakos, Gregory J.↗
DOE OSTI2016

Anomalous thermal conductivity of monolayer boron nitride

In this paper, we use nonequilibrium molecular dynamics modeling to investigate the thermal properties of monolayer hexagonal boron nitride nanoribbons under uniaxial strain along their longitudinal axis. Our simulations predict that hexagonal boron nitride shows an anomalous thermal response to the applied uniaxial strain. Contrary to three dimensional materials, under uniaxial stretching, the thermal conductivity of boron nitride nanoribbons first increases rather than decreasing until it reaches its peak value and then starts decreasing. Under compressive strain, the thermal conductivity of monolayer boron nitride ribbons monolithically reduces rather than increasing. We use phonon spectrum and dispersion curves to investigate the mechanism responsible for the unexpected behavior. Our molecular dynamics modeling and density functional theory results show that application of longitudinal tensile strain leads to the reduction of the group velocities of longitudinal and transverse acoustic modes. Such a phonon softening mechanism acts to reduce the thermal conductivity of the nanoribbons. On the other hand, a significant increase in the group velocity (stiffening) of the flexural acoustic modes is observed, which counteracts the phonon softening effects of the longitudinal and transverse modes. The total thermal conductivity of the ribbons is a result of competition between these two mechanisms. At low tensile strain, the stiffening mechanism overcomes the softening mechanism which leads to an increase in the thermal conductivity. At higher tensile strain, the softening mechanism supersedes the stiffening and the thermal conductivity slightly reduces. Our simulations show that the decrease in the thermal conductivity under compressive strain is attributed to the formation of buckling defects which reduces the phonon mean free path.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

SuperScreen: An open-source package for simulating the magnetic response of two-dimensional superconducting devices

Quantitative understanding of the spatial distribution of magnetic fields and Meissner screening currents in two-dimensional (2D) superconductors and mesoscopic thin film superconducting devices is critical to interpreting the results of magnetic measurements of such systems. Here, we introduce SuperScreen, an open-source Python package for simulating the response of 2D superconductors to trapped flux and applied time-independent or quasi-DC magnetic fields for any value of the effective magnetic penetration depth, Λ. Given an applied magnetic field, SuperScreen solves the 2D London equation using an efficient matrix inversion method to obtain the Meissner currents and magnetic fields in and around structures composed of one or more superconducting thin films of arbitrary geometry with spatially nonuniform magnetic penetration depth. SuperScreen can be used to model screening effects and calculate self- and mutual-inductance in superconducting devices, and simulate the magnetic response of inhomogeneous 2D superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The cosmic muon and detector simulation framework of the extreme energy events (EEE) experiment

Abstract This paper describes the simulation framework of the extreme energy events (EEE) experiment. EEE is a network of cosmic muon trackers, each made of three multi-gap resistive plate chambers (MRPC), able to precisely measure the absolute muon crossing time and the muon integrated angular flux at the ground level. The response of a single MRPC and the combination of three chambers have been implemented in a GEANT4-based framework (GEMC) to study the telescope response. The detector geometry, as well as details about the surrounding materials and the location of the telescopes have been included in the simulations in order to realistically reproduce the experimental set-up of each telescope. A model based on the latest parametrization of the cosmic muon flux has been used to generate single muon events. After validating the framework by comparing simulations to selected EEE telescope data, it has been used to determine detector parameters not accessible by analysing experimental data only, such as detection efficiency, angular and spatial resolution.

Abbrescia, M.↗
DOE OSTI2021

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗
NASA NTRS2011

Two-Dimensional Finite Element Ablative Thermal Response Analysis of an Arcjet Stagnation Test

The finite element ablation and thermal response (FEAtR, hence forth called FEAR) design and analysis program simulates the one, two, or three-dimensional ablation, internal heat conduction, thermal decomposition, and pyrolysis gas flow of thermal protection system materials. As part of a code validation study, two-dimensional axisymmetric results from FEAR are compared to thermal response data obtained from an arc-jet stagnation test in this paper. The results from FEAR are also compared to the two-dimensional axisymmetric computations from the two-dimensional implicit thermal response and ablation program under the same arcjet conditions. The ablating material being used in this arcjet test is phenolic impregnated carbon ablator with an LI-2200 insulator as backup material. The test is performed at the NASA, Ames Research Center Interaction Heating Facility. Spatially distributed computational fluid dynamics solutions for the flow field around the test article are used for the surface boundary conditions.

Dec, John A.↗
DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
NASA NTRSDate not supplied

Phase transition energetics-based mass loss modeling of chondritic Near Earth Objects

Mitigating potentially hazardous asteroids relies on accurate knowledge of their composition and critical physical characteristics such as shape and mass. Ground-based observatories and sample return missions in recent years have focused on such characterization to track Near Earth Objects. Historically, entry speeds of asteroids or mass change due to break-up/airburst have been considered in risk assessment studies as critical parameters. To this effect, most hydrocode simulations that consider energy deposition techniques rely on gross meteor modeling with simplistic spherical shapes while neglecting heat of ablation considerations. This is the general approach for fragment-cloud, pancake and hybrid fragmentation models. Alternatively, aerothermal fragmentation modeling that accounts for meteoritic shape change have typically focused on environmental modeling to understand break-up points primarily due to radiative heating. A physics-based risk assessment could benefit by modeling the environment and material response of typical meteoroidal materials as precursors to fragmentation. A greater focus on the material response to applied environmental heating will be emphasized here. This work presents a forward analysis of thermal and material response of stony meteorites by ac-counting for mass loss and shape change effects through a boundary layer formalism. This allows for simulations starting from simplistic shapes that could eventually lead to alterations in contours due to ablation. These effects would be characterized through phase transitions of the surface, such as melting and vaporization, in response to aerodynamic heating. Coupled ablation and radiative heating effects have been found to lead to reduced heat transfer coefficients in the computation of meteor mass loss rates. This has the potential to increase ground damage footprints for hazardous impacts. The effort here quantifies the heat transfer coefficient by accounting for ablation assisted by phase transition. The melt and vaporization mediated ablation of chondritic meteors provides essential recession estimates based upon surface environments. The altered mass due to such phenomena would be presented for a H-chondrite meteoritic sample.

Pratibha Raghunandan↗
DOE OSTI2018

GeV ion irradiation of NiFe and NiCo: Insights from MD simulations and experiments

Concentrated solid solution alloys have attracted rapidly increasing attention due to their potential for designing materials with high tolerance to radiation damage. To tackle the effects of chemical complexity in defect dynamics and radiation response, we present in this paper a computational study on swift heavy ion induced effects in Ni and equiatomic Ni -based alloys (Ni 50 Fe 50 , Ni 50 Co 50 ) using two-temperature molecular dynamics simulations (2T-MD). The electronic heat conductivity in the two-temperature equations is parameterized from the results of first principles electronic structure calculations. A bismuth ion (1.542 GeV) is selected and single impact simulations performed in each target. We study the heat flow in the electronic subsystem and show that alloying Ni with Co or Fe reduces the heat dissipation from the impact by the electronic subsystem. Simulation results suggest no melting or residual damage in pure Ni while a cylindrical region melts along the ion propagation path in the alloys. In Ni 50 Co 50 the damage consists of a dislocation loop structure (d = 2 nm) and isolated point defects, while in Ni 50 Fe 50 , a defect cluster (d = 4 nm) along the ion path is, in addition, formed. The simulation results are supported by atomic-level structural and defect characterizations in bismuth-irradiated Ni and Ni 50 Fe 50 . Finally, the significance of the 2T-MD model is demonstrated by comparing the results to those obtained with an instantaneous energy deposition model without consideration of e-ph interactions in pure Ni and by showing that it leads to a different qualitative behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
NASA NTRS1987

Influence of electron radiation and temperature on the cyclic, matrix dominated response of graphite-epoxy

The effects of electron radiation and elevated temperature on the matrix-dominated cyclic response of standard T300/934 and a chemically modified T300/934 graphite-epoxy are characterized. Both materials were subjected to 1.0 x 10 to the 10th rads of 1.0 MeV electron irradiation, under vacuum, to simulate 30 years in geosynchronous orbit. Cyclic tests were performed at room temperature and elevated temperature (121 C) on 4-ply unidirectional laminates to characterize the effects associated with irradiation and elevated temperature. Both materials exhibited energy dissipation in their response at elevated temperature. The irradiated modified material also exhibited energy dissipation at room temperature. The combination of elevated temperature and irradiation resulted in the most severe effects in the form of lower proportional limits, and greater energy dissipation. Dynamic-mechanical analysis demonstrated that the glass transition temperature, T(g), of the standard material was lowered 39 C by irradiation, wereas the T(g) of the modified material was lowered 28 C by irradiation. Thermomechanical analysis showed the occurrence of volatile products generated upon heating of the irradiated materials.

Reed, Susan M.↗
DOE OSTI2019

Viscoplastic Interpretation of Localized Compaction Creep in Porous Rock

Recent laboratory evidence shows that compaction creep in porous rocks may develop through stages of acceleration, especially if the material is susceptible to strain localization. This paper provides a mechanical interpretation of compaction creep based on viscoplasticity and nonlinear dynamics. For this purpose, a constitutive operator describing the evolution of compaction creep is defined to evaluate the spontaneous accumulation of pore collapse within an active compaction band. This strategy enables the determination of eigenvalues associated with the stability of the response which is able to differentiate decelerating from accelerating strain. This mathematical formalism was linked to a constitutive law able to simulate compaction localization. Material point simulations were then used to identify the region of the stress space where unstable compaction creep is expected, showing that accelerating strains correspond to pulses of inelastic strain rate. Such pulses were also found in full-field numerical analyses of delayed compaction, revealing that they correspond to stages of inception and propagation of new bands across the volume of the simulated sample. These results illustrate the intimate relation between the spatial patterns of compaction and their temporal dynamics, showing that while homogeneous compaction develops with decaying rates of accumulation, localized compaction occurs through stages of accelerating deformation caused by the loss of strength taking place during the formation of a band. In addition, they provide a predictive modeling framework to simulate and explain the spatio-temporal dynamics of compaction in porous sedimentary formations.

58 GEOSCIENCES↗
DOE OSTI2024

Role of interfaces on the mechanical response of accumulative roll bonded nanometallic laminates investigated via dislocation dynamics simulations

Unraveling the effects of continuous dislocation interactions with interfaces, particularly at the nanometer length scales, is key to a broader understanding of plasticity, to material design and to material certification. To this end, this work proposes a novel discrete dislocation dynamics-based model for dislocation interface interactions tracking the fate of residual dislocation on interfaces. This new approach is used to predict the impact of dislocation/interface reactions on the overall mechanical behavior of accumulative roll bonded nanometallic laminates. The framework considers the dynamic evolution of the interface concurrent with a large network of dislocations, thus, accounting for the local short and long range effects of the dislocations under the external boundary conditions. Specifically, this study focuses on two-phase Fe/Cu nanometallic laminates, and investigates the role of the underlying elastic and plastic contrast of the Fe and the Cu layers on the composite response of the material. Moreover, the role of initial microstructures, resulting from processing is also investigated. Subsequently, the model is used to examine the effect of layer thickness and interface orientation relationship on the residual stresses of the relaxed microstructure. The associated mechanical response of these laminates are compared when loaded under normal direction compression, as well as shear compression. Finally, this work predicts a dominant effect of the layer thickness, as compared to the interface orientation relationship, on the macroscopic response and on the residual stresses of these nanolaminates, while the local dislocation transmission propensity through the interface is significantly influenced by the corresponding orientation relationship.

36 MATERIALS SCIENCE↗