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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2017

New Method to Estimate the Permittivity of Dielectric Materials Using MSPS

A new method to estimate the relative dielectric permittivity of material’s sheet is presented. In this method, a conventional microstrip patch sensor (MSPS) is utilized to resonate at 2.45 GHz. An unknown material, considered as a superstrate layer of 5 mm thickness, is located in the near-field region of the MSPS, at a distance of 5 mm above the sensor. From the changes in the MSPS’s parameters, the dielectric property of the unknown material can be deduced. Based on the conformal map theory, theoretical analyses were performed, and hence, simulation analyses were done using the CST-MW software package. This proposed method was validated through measurements. From the theoretical, simulation, and measurement results, it was shown that depending on the frequency response, the relative permittivity can be estimated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

TDEP: Temperature Dependent Effective Potentials

The Temperature Dependent Effective Potential (TDEP) method is a versatile and efficient approach to include temperature in ab initio materials simulations based on phonon theory. TDEP can be used to describe thermodynamic properties in classical and quantum ensembles, and several response properties ranging from thermal transport to Neutron and Raman spectroscopy. A stable and fast reference implementation is given in the software package of the same name described here. The underlying theoretical framework and foundation is briefly sketched with an emphasis on discerning the conceptual difference between bare and effective phonon theory, in both self-consistent and non-self-consistent formulations. References to numerous applications and more in-depth discussions of the theory are given.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Development of a Three-Dimensional, Unstructured Material Response Design Tool

A preliminary verification and validation of a new material response model is presented. This model, Icarus, is intended to serve as a design tool for the thermal protection systems of re-entry vehicles. Currently, the capability of the model is limited to simulating the pyrolysis of a material as a result of the radiative and convective surface heating imposed on the material from the surrounding high enthalpy gas. Since the major focus behind the development of Icarus has been model extensibility, the hope is that additional physics can be quickly added. This extensibility is critical since thermal protection systems are becoming increasing complex, e.g. woven carbon polymers. Additionally, as a three-dimensional, unstructured, finite-volume model, Icarus is capable of modeling complex geometries. In this paper, the mathematical and numerical formulation is presented followed by a discussion of the software architecture and some preliminary verification and validation studies.

Thermal Protection System↗
DOE OSTI2018

The Lorenzo–Hartley’s function for fractional calculus and its applications pertaining to fractional order modelling of anomalous relaxation in dielectrics

Considering the applications and a great significance of the study pertaining to generalized functions in applied sciences, the author investigated the Lorenzo–Hartley’s function which is a generalization of classical functions widely used in fractional calculus. The Laplace transform pairs are derived and the requirements for the Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) to be completely monotonic (for t>0) are investigated. The author has shown the applications of this generalized function to describe the relaxation models, particularly in dielectrics. The Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) of a real variable t is considered to investigate a computable mathematical framework for standard Debye and non-Debye relaxation processes in dielectric materials. The non-negative spectral distribution function is obtained for the corresponding response function. It is also demonstrated that the classical models like Cole–Cole (C–C), Davidson–Cole (D–C), and Havriliak–Negami (H–N) for non-Debye relaxation and standard Debye relaxation are the particular cases of the Lorenzo–Hartley’s function. Some of the study cases are also worked out to visualize the effects of variations of parameters on the response function and corresponding spectral distribution function. A generalized and unified fractional relaxation differential equation which governs response functions for classical dielectrics models pertaining to non-Debye relaxation (C–C, D–C, and H–N) is also established in the present investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dynamic Manipulation in Piezoresponse Force Microscopy: Creating Nonequilibrium Phases with Large Electromechanical Response

We discuss how domain walls and topological defects in ferroelectric materials have emerged as a powerful tool for functional electronic devices including memory and logic. Similarly, wall interactions and dynamics underpin a broad range of mesoscale phenomena ranging from giant electromechanical responses to memory effects. Exploring the functionalities of individual domain walls, their interactions, and controlled modifications of the domain structures is crucial for applications and fundamental physical studies. However, the dynamic nature of these features severely limits studies of their local physics since application of local biases or pressures in piezoresponse force microscopy induce wall displacement as a primary response. Here, we introduce an approach for the control and modification of domain structures based on automated experimentation, whereby real-space image-based feedback is used to control the tip bias during ferroelectric switching, allowing for modification routes conditioned on domain states under the tip. This automated experiment approach is demonstrated for the exploration of domain wall dynamics and creation of metastable phases with large electromechanical response.

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

Optimization of electrical treatment strategy for surface roughness reduction in conducting thin films

The surface roughness of deposited conducting thin films is responsible for various materials reliability problems in nanoelectronics and nanofabrication technologies. In this work, we report a modeling and simulation study that aims at optimizing the electrical surface treatment of deposited conducting thin films as a physical processing strategy for their surface roughness reduction. Our research is based on a continuum model of film surface morphological evolution that accounts for the residual stress in the deposited conducting thin film, the film’s wetting of the substrate layer that it is deposited on, film texture and surface diffusional anisotropy, and surface electromigration. Through systematic linear stability analysis and dynamical simulation protocols, we examine in detail the effects of film surface crystallographic orientation and applied electric field direction toward minimizing the electric field strength required for film surface smoothening. We believe that the critical electric field strength requirement for surface roughness reduction on {110}, {100}, and {111} surfaces of face-centered cubic crystalline conducting thin films exhibits a very strong dependence on the applied electric field direction, expressed as the electric field misalignment with respect to the principal residual stress directions in the film and the fast surface diffusion directions. Due to these findings, we optimize the electrical treatment strategy for surface roughness reduction of conducting thin films with respect to all relevant processing and material parameters.

42 ENGINEERING↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2018

A Methodology for Determining the Concentration of Naturally Occurring Radioactive Materials in an Urban Environment

The detection of radioactive sources in an urban setting is greatly complicated by natural background radiation, which emanates from various materials including roadways, sidewalks, soil, and building exteriors. The method presented and demonstrated here represents an effort to characterize the concentration of naturally occurring radioactive material (NORM) in these types of materials. The location surveyed in this work was the Fort Indiantown Gap Combined Arms Collective Training Facility in Lebanon County, Pennsylvania. Over 70 measurements with a high-purity germanium detector were performed to ascertain the NORM concentrations present in the soil, asphalt, gravel, concrete, and walls found throughout the site. Monte Carlo radiation transport calculations were used to obtain detector responses for these various geometries and materials to convert these measurements into NORM concentration estimates. Lastly, synthetic spectra were simulated using the predicted source terms and compared to actual measurements, showing acceptable agreement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Photorefractive Correction of the Coupling Ratio of an Integrated Optical Directional X-Coupler on a Lithium Niobate Substrate

The possibility of correcting the power coupling ratio of an integrated optical directional X-coupler upon excitation of a local photorefractive response in the substrate material (LiNbO{sub 3}) has been analyzed. The X-coupler operation has been numerically simulated, and the regions of its maximum sensitivity to photorefraction have been determined. The correspondence between the calculated and experimental data has been obtained. The maximum correction value of the coupling ratio is 1–2%. The photorefractive correction of the coupling ratio has been used to increase the extinction ratio of a Mach–Zehnder modulator to 47 dB.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
DOE OSTI2020

Multiscale Thermomechanical Damage Model with Internal State Variables for Ceramic Matrix Composites

A length scale-dependent matrix damage model has been developed to simulate the response of ceramic matrix composites (CMCs). This model captures damage initiation and propagation in the brittle matrix material by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress-strain constitutive relationships at each length scale. Material degradation due to matrix cracking is captured by a damage variable determined using fracture mechanics and the self-consistent scheme. An additional state variable is introduced to capture the nucleation and growth of matrix porosity, a key mechanism contributing to matrix nonlinearity. Porosity occurs as a result of material diffusion around grain boundaries and the growth of the porosity state variable is related to the material entropy dissipation and the volumetric strain. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide, unidirectional (UD) carbon fiber/silicon carbide matrix (C/SiC) CMC, and five-harness satin (5HS) woven C/SiC CMC.

Skinner, Travis D.↗
NASA NTRS2010

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗
DOE OSTI2021

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method (DEM) has been intensively used to study the constitutive behaviour of granular materials. However, to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear. This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials. Three groups of models with spherical and clumped particles are investigated from four perspectives: (i) deviatoric stress and volumetric behaviour; (ii) critical state behaviour; (iii) stress-dilatancy relationship; and (iv) the evolution of principal stress ratio against axial strain. The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses, volumetric behaviour, and critical state behaviour observed in experiments for granular materials. On the other hand, some qualitative deviations between experiments and the investigated DEM simulations are also observed, in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain, which are proven to be critical for constitutive modelling. The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective. It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Development of Computational Materials Workflows for Additively Manufactured Metallic Materials to Enable Accelerated Prediction of Fatigue Performance

The maturation of computational materials approaches for fatigue performance prediction in a qualification and certification process is stifled by the ability to validate complex, microstructure-based simulations. Such a validation strategy bears immediate challenges including generating accurate virtual microstructures, efficiently solving physics-based mechanical simulations over relevant spatial and temporal scales, and acquiring high-fidelity calibration and validation data at the appropriate length scale. This presentation will overview these common challenges and present a case study to demonstrate a computational materials workflow for additively manufactured metallic materials. In this study, process-specific defects are characterized using segmented X-Ray micro-computed tomography measurements and overlaid on virtual microstructures. Accelerated crystal plasticity-based fatigue simulations are performed to demonstrate cyclic evolution and localization of mechanical fields in the vicinity of defects in response to their precise spatial configuration. An example of how this computational materials workflow may support next-generation qualification is discussed.

computational materials↗
DOE OSTI2021

Spectral Effects in Albedo and Rearside Irradiance Measurement for Bifacial Performance Estimation

We investigate the impact of spectral dependence of ground surface reflectivity on albedo and rearside irradiance measurements necessary for bifacial photovoltaic (PV) module performance estimation and monitoring. Because PV modules are spectrally selective, albedo and irradiance measurements performed with common irradiance sensors may require spectral mismatch corrections when used for performance prediction. We investigate via simulation the differences in spectrally responsive albedo measured with thermopile pyranometers and crystalline silicon PV reference cells in comparison to a typical crystalline-silicon bifacial PV module. Simulations are performed for nine different representative ground surface materials using simulated solar spectra together with spectral reflectivity data distributed with the SMARTS simulation software. For the materials considered, the results show that albedo spectral mismatch relative to the bifacial module is distributed over a range of ±9.2% for thermopile pyranometers versus only ±3.7% for a typical PV reference cell. We consider the impact of this spectrally-responsive albedo mismatch on bifacial PV module rearside irradiance measurements. Using synthesized rearside spectral irradiance distributions, we find that for the nine different ground surface materials the predicted rearside irradiance measurement deviates from the effective irradiance observed by the PV module by on the order of 16.5 W/m 2 for the pyranometer and 3.6 W/m 2 for the PV reference cell. We discuss the implications for bifacial albedo and irradiance measurement.

41 EE - Solar Energy Technologies Office (EE-4S)↗
NASA NTRS2005

Primal and Dual Integrated Force Methods Used for Stochastic Analysis

At the NASA Glenn Research Center, the primal and dual integrated force methods are being extended for the stochastic analysis of structures. The stochastic simulation can be used to quantify the consequence of scatter in stress and displacement response because of a specified variation in input parameters such as load (mechanical, thermal, and support settling loads), material properties (strength, modulus, density, etc.), and sizing design variables (depth, thickness, etc.). All the parameters are modeled as random variables with given probability distributions, means, and covariances. The stochastic response is formulated through a quadratic perturbation theory, and it is verified through a Monte Carlo simulation.

Patnaik, Surya N.↗