Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constitutive Modeling”

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.

796 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Development, Characterization, and Validation of Elevated-Temperature Constitutive Models: Deformation and Damage

This report provides a brief review of experimentally observed hereditary and nonhereditary material behavior along with background information on standard as well as advanced internal state variable constitutive modeling at elevated temperature. A description of exploratory, characterization, and validation testing is presented along with a detailed outline of what constitutes “sufficient” data content (i.e., quality and quantity) for developing or enhancing, characterizing, and validating a particular sophisticated nonlinear time- and history-dependent (hereditary) class of constitutive models known as GVIPS (generalized viscoplasticity with potential structure). The tests described are necessary to reveal a material’s behavior in the reversible (or viscoelastic) and irreversible (or viscoplastic) regimes, for the identification of both deformation and damage model parameters. Results presented are primarily for metallic materials. In all cases, both uniaxial and multiaxial tests are described, and the linkage between the specific tests and the parameters within the model that can be characterized from the results of these tests are also discussed. Discussion is also provided relative to the role information management must play relative to material data collection, analysis, maintenance, and dissemination. The need for such an information system is particularly important as both analyst and designer move toward utilizations of sophisticated, nonlinear time- and history dependent (hereditary) constitutive models. Lastly, the concept of state space and its utility in understanding and establishing constitutive models is addressed throughout. The intent behind this document is to help both the modeler and experimentalist understand each other’s specific points of view and provide guidance for both model development and characterization. Emphasis has been placed on providing the mechanician with information regarding how tests are performed and what issues to be aware of when interpreting results. It is hoped that experimentalists will take away a new perspective on the types of information that modelers and analysts are looking for from them.

Constitutive Modeling↗
NASA NTRS2026

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles with Interrupted Relaxation Periods of TIMETAL 21S at 650 °C

The predictive performance of a multimechanism GVIPS model with saturating hardening function is assessed against a set of high temperature cyclic experiments under both strain-control and stress-control at varying rates of loading as well as several complex interrupted cyclic responses. The material specimens were composed of the titanium alloy, TIMETAL 21S, and tested at 650 °C, see Lissenden et al., 2007. The model did a very good job of predicting the variety of tests, particularly given the fact that only monotonic loading tests (tensile, creep, relaxation) and a single fully reversed cycle was used for characterization. It was particularly interesting that the model was capable of reasonably capturing the rate of accumulated strain during ratchetting under tensile mean stress. Given complex interrupted cyclic response, the model was able to both qualitatively and even quantitatively predict reasonably well both the cyclic and associated relaxation periods (in all quadrants of the stress-strain space) thus confirming the validity of the chosen functional forms for both hardening and thermal recovery. The quantitative inaccuracy is associated with the fact that the material specimens tested by Lissenden et al., 2007 was significantly “softer” than those tested by Castelli in Saleeb et al., 1994 which were used for characterization of the GVIPS model parameters.

Viscoplaticity↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Calibration of the Nonlinear Viscoelastic SPECtacular Model for 828-CTBN/DEA/GMB and 828/Z/AlOx

Nonlinear viscoelastic model calibrations are developed for two filled epoxies: an adduct of Epon 828 DGEBA and CTBN cross-linked with diethanolamine (DEA) filled with glass microballoons, called 828-CTBN/DEA/GMB, and Epon 828 DGEBA cross-linked with Ancamine Z and filled with alumina, called 828/Z/AlOx. The thermo-mechanical responses of these materials are represented using the SPECtacular constitutive model, an extension of the Simplified Potential Energy Clock (SPEC) model. The model couples

Cundiff, Kenneth Noel [Sandia National Laboratorie↗
NASA NTRS2026

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗
DOE OSTI2025

Quantum graph models for transport in filamentary switching

The formation of metallic nanofilaments bridging two electrodes across an insulator is a mechanism for resistive switching. Examples of such phenomena include atomic synapses, which constitute a distinct class of memristive devices the behavior of which is closely tied to the properties of the filament. Until recently, experimental investigation of the low-temperature regime and quantum transport effects has been limited. However, with growing interest in understanding the true impacts of the filament on device conductance, comprehending quantum effects has become crucial for quantum neuromorphic hardware. Here, we discuss quantum transport resulting from filamentary switching in a narrow region where the continuous approximation of the contact is not valid, and only a few atoms are involved. In this scenario, the filament can be represented by a graph depicting the adjacency of atoms and the overlap between atomic orbitals. Using the theory of quantum graphs with locally diffusive node scattering, we calculate the scattering amplitude of charge carriers on this graph and explore the interplay between filamentary formation and quantum transport effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk↗
DOE OSTI2026

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗
DOE OSTI2024

Characterizing the γ-Ray Emission from FR0 Radio Galaxies

Abstract FR0 galaxies constitute the most abundant jet population in the local Universe. With their compact jet structure, they are broadband photon emitters and have been proposed as multimessenger sources. Recently, these sources have been detected for the first time inγrays. Using a revised FR0 catalog, we confirm that the FR0 population as a whole areγ-ray emitters, and we also identify two significant sources. For the first time, we find a correlation between the 5 GHz core radio luminosity andγ-ray luminosity in the 1–800 GeV band, having a 4.8σstatistical significance. This is clear evidence that the jet emission mechanism is similar in nature for FR0s and the well-studied canonical FR (FRI and FRII) radio galaxies. Furthermore, we perform broadband spectral energy distribution modeling for the significantly detected sources as well as the subthreshold source population using a one-zone synchrotron self-Compton model. Within the maximum jet power budget, our modeling shows that the detectedγrays from the jet can be explained as inverse Compton photons. To explain the multiwavelength observations for these galaxies, the modeling results stipulate a low bulk Lorentz factor and a jet composition far from equipartition, with the particle energy density dominating over the magnetic field energy density.

Astronomy & Astrophysics↗
DOE OSTI2025

Quantum mechanical dataset of 836k neutral closed-shell molecules with up to 5 heavy atoms from C, N, O, F, Si, P, S, Cl, Br

Abstract We introduce the Vector-QM24 (VQM24) dataset comprehensively covering all possible neutral closed-shell small organic and inorganic molecules with up to five heavy (p-block) atoms: C, N, O, F, Si, P, S, Cl, Br. All valid stoichiometries, Lewis-rule-consistent graphs, and stable conformers (identified via GFN2-xTB) were enumerated combinatorially, yielding 577k conformational isomers spanning 258k constitutional isomers and 5,599 unique stoichiometries. DFT (ωB97X-D3/cc-pVDZ) optimizations were performed for all, and diffusion quantum Monte Carlo (DMC@PBE0(ccECP/cc-pVQZ)) energies are provided for 10,793 lowest-energy conformers with up to 4 heavy atoms. VQM24 includes structures, vibrational modes, rotational constants, thermodynamic properties (Gibbs free energies, enthalpies, ZPVEs, entropies, heat capacities), and electronic properties such as atomization, electron interaction, exchange-correlation, dispersion energies, multipole moments (dipole to hexadecapole), alchemical potentials, Mulliken charges, and wavefunctions. Machine learning models of atomization energies on this dataset reveal significantly higher complexity than QM9, with none achieving chemical accuracy. VQM24 offers a rigorous, high-fidelity benchmark for evaluating quantum machine learning models.

Science & Technology - Other Topics↗
DOE OSTI2024

Black-hole powered quantum coherent amplifier

Abstract Atoms falling into a black hole (BH) through a cavity are shown to enable coherent amplification of light quanta powered by the BH-gravitational vacuum energy. This process can harness the BH energy towards useful purposes, such as propelling a spaceship trapped by the BH. The process can occur via transient amplification of a signal field by falling atoms that are partly excited by Hawking radiation reflected by an orbiting mirror. In the steady-state regime of thermally equilibrated atoms that weakly couple to the field, this amplifier constitutes a BH-powered quantum heat engine. The envisaged effects substantiate the thermodynamic approach to BH acceleration radiation.

Physics↗
NASA NTRSDate not supplied

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation↗
NASA NTRS1990

Gaseous Optical Contamination of the Spacecraft Environment: A Review

Interactions between the ambient atmosphere and orbiting spacecraft, sounding rockets, and suborbital vehicles, and with their effluents, give rise to optical (extreme UV to LWIR) foreground radiation which constitutes noise that raises the detection threshold for terrestrial and celestial radiations, as well as military targets. Researchers review the current information on the on-orbit optical contamination. Its source species are created in interaction processes that can be grouped into three categories: (1) Reactions in the gas phase between the ambient atmosphere and desorbates and exhaust; (2) Reactions catalyzed by exposed ram surfaces, which occur spontaneously even in the absence of active material releases from the vehicles; and (3) Erosive excitative reactions with exposed bulk (organic) materials, which have recently been identified in the laboratory though not as yet observed on spacecraft. Researchers also assess the effect of optical pumping by earthshine and sunlight of both reaction products and effluents.

Tran, N. H.↗
DOE OSTI2024

Osmotic and phoretic competition explains chemotaxic assembly and sorting

Microscale objects responding to chemical gradients by migrating toward or away from a preferred species is a simple yet constitutive mechanism by which transport occurs in biological organisms. Synthetic chemotaxis provides key physical descriptions of simplified systems that can be used in biological models, or in the creation of advanced responsive material systems. In this article, we provide a quantitative framework for understanding synthetic chemotaxis of microparticles which involves a competition between phoresis and osmosis. We present separate quantitative measurements of phoresis and osmosis acting on individual taxing particles, finding that phoresis follows the long-predicted v ∼ 1 / r 2 scaling while the osmotic contribution depends on the geometry and details of the system, and must be solved on a case-by-case basis. Through this, we are able to develop a more accurate picture of particle transport at the single particle level. Equipped with this approach, we go on to describe how high concentrations of particles in a symmetric chemical gradient grow close-packed hives that reach a steady-state size tunable through light intensity or particle size. Last, we demonstrate that mixed particles experiencing the same chemical gradient will selectively migrate toward or away depending on the nature of the particle surface, thereby locally sorting out a particular species. We anticipate these results will be important in describing both biological and synthetic chemotaxis in phoretic systems and should bring a wealth of studies that take advantage of competing osmotic flows to illicit unexpected dynamic active behavior.

Science & Technology - Other Topics↗
DOE OSTI2025

Nuclear quantum effects in molecular liquids across chemical space

Abstract Nuclear quantum effects (NQEs) influence many physical and chemical phenomena, particularly those involving light atoms or occurring at low temperatures. However, their impact has been carefully quantified in few systems-like water-and is rarely considered more broadly. Here we use path-integral molecular dynamics to systematically investigate NQEs on thermophysical properties of 92 organic liquids at ambient conditions. Depending on chemical constitution, we find substantial impact across thermal expansivity, compressibility, dielectric constant, enthalpy of vaporization, and notably molar volume, which shows consistent, positive quantum-classical differences up to 5%; similar, less pronounced trends manifest as isotope effects from deuteration. Using data-driven analysis, we identify three features-molar mass, classical hydrogen density, and classical thermal expansivity-that accurately predict NQEs and facilitate understanding of how characteristics like branching and heteroatom content influence behavior. This work highlights the broad relevance of NQEs in molecular liquids, while also providing a conceptual and practical framework to anticipate their impact.

Science & Technology - Other Topics↗
NASA NTRS2004

Sample Return Science by Hayabusa Near-Earth Asteroid Mission

Assigning the material species to each asteroid spectral type and finding out the corresponding meteorite category is crucial to make the global material map in the whole asteroid belt and to understand the evolution of the asteroid belt. Recent direct observations by spacecrafts are revealing new intriguing aspects of asteroids which cannot be obtained solely from ground-based observations or meteorite studies. However identification of the real material species constituting asteroids and their corresponding meteorite analogs are still ambiguous. Space weathering makes difficult to identify the true material, and there is still a great gap between the remote sensing data on the global surface and the local microscopic data from meteorites. Sample return from asteroids are inevitable to solve these problems. For this purpose sample return missions to asteroids belonging to various spectral classes are required. The HAYABUSA spacecraft (prelaunch name is MUSESC) launched last year is the first attempt on this concept. This report presents outline of the mission with special stress on its science.

A Fujiwara↗