Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Inelastic deformation of metal matrix composites

A theoretical model capable of predicting the thermomechanical response of continuously reinforced metal matrix composite laminates subjected to multiaxial loading was developed. A micromechanical model is used in conjunction with nonlinear lamination theory to determine inelastic laminae response. Matrix viscoplasticity, residual stresses, and damage to the fiber/matrix interfacial zone are explicitly included in the model. The representative cell of the micromechanical model is considered to be in a state of generalized plane strain, enabling a quasi two-dimensional analysis to be performed. Constant strain finite elements are formulated with elastic-viscoplastic constitutive equations. Interfacial debonding is incorporated into the model through interface elements based on the interfacial debonding theory originally presented by Needleman, and modified by Tvergaard. Nonlinear interfacial constitutive equations relate interfacial tractions to displacement discontinuities at the interface. Theoretical predictions are compared with the results of an experimental program conducted on silicon carbide/titanium (SiC/Ti) unidirectional, (O4), and angle-ply, (+34)(sub s), tubular specimens. Multiaxial loading included increments of axial tension, compression, torque, and internal pressure. Loadings were chosen in an effort to distinguish inelastic deformation due to damage from matrix plasticity and separate time-dependent effects from time-independent effects. Results show that fiber/matrix debonding is nonuniform throughout the composite and is a major factor in the effective response. Also, significant creep behavior occurs at relatively low applied stress levels at room temperature.

Lissenden, C. J.↗
DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonlinear modeling and performance analysis of cracked beam microgyroscopes

Cracks constitute a common structural defect in microelectromechanical systems that may arise during manufacturing, mechanical fatigue, or shock loading. Areas weakened by cracks increase the flexibility of the microstructure. Furthermore, depending on the crack severity, the increased flexibility can cause dramatic changes in the static and dynamic behaviors of electrically actuated MEMS. In this study, a numerical investigation of the depth and location of multiple edge cracks on the performance of beam microgyroscopes is conducted. Applying the Griffith strain energy release theorem for an edge crack, the vibrational characteristics of microgyroscopes with multiple cracks including the static deflection, natural frequencies, and mode shapes are studied analytically. Then, the differential quadrature method is employed to model the nonlinear dynamic behaviors of the cracked gyroscope system. Cracks initiated along the driving or sensing directions of the microgyroscope are found to have differing effects on its nonlinear dynamic response. This numerical study reveals that the onset of severe single cracks or a series of multiple cracks can significantly affect the sensitivity of the microgyroscope to base rotations and therefore leads to a degradation in the performance of the damaged MEMS sensor. This study also shows that a potential design approach for broadband microgyroscope systems manufactured with cracks can be proposed.

42 ENGINEERING↗
DOE OSTI2023

Data-driven prediction of α IIb β 3 integrin activation paths using manifold learning and deep generative modeling

The integrin heterodimer is a transmembrane protein critical for driving cellular process and is a therapeutic target in the treatment of multiple diseases linked to its malfunction. Activation of integrin involves conformational transitions between bent and extended states. Some of the conformations that are intermediate between bent and extended states of the heterodimer have been experimentally characterized, but the full activation pathways remain unresolved both experimentally due to their transient nature and computationally due to the challenges in simulating rare barrier crossing events in these large molecular systems. An understanding of the activation pathways can provide new fundamental understanding of the biophysical processes associated with the dynamic interconversions between bent and extended states and unveil new putative therapeutic targets. In this work, we apply nonlinear manifold learning to coarse-grained molecular dynamics simulations of bent, extended, and two intermediate states of αI I b β3 integrin to learn a low-dimensional embedding of the configurational phase space. We then train deep generative models to learn an inverse mapping between the low-dimensional embedding and high-dimensional molecular space and use these models to interpolate the molecular configurations constituting the activation pathways between the experimentally characterized states. Furthermore, this work furnishes plausible predictions of integrin activation pathways and reports a generic and transferable multi-scale technique to predict transition pathways for biomolecular systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗
DOE OSTI2024

Physics informed neural network can retrieve rate and state friction parameters from acoustic monitoring of laboratory stick-slip experiments

Various machine learning (ML) and deep learning (DL) techniques have been recently applied to the forecasting of laboratory earthquakes from friction experiments. The magnitude and timing of shear failures in stick-slip cycles are predicted using features extracted from the recorded ultrasonic or acoustic emission (AE) signals. In addition, the Rate and State Friction (RSF) constitutive laws are extensively used to model the frictional behavior of faults. In this work, we use data from shear experiments coupled with passive acoustic (variance, kurtosis, and AE rate) interleaved with active source ultrasonic monitoring (transmitted wave amplitude) to develop physics-informed neural network (PINN) models incorporating the RSF law and AE rate generation equation with wave amplitude serving as a proxy for friction state variable. This PINN framework allows learning RSF parameters from stick-slip experiments rather than measuring them through a series of velocity step experiments. We observe that when the stick-slip cycles are irregular, the PINN models outperform the data-driven DL models. Transfer learning (TL) PINN models are also developed by pre-training on data collected at one normal stress level followed by forecasting shear failures and retrieving RSF parameters at other stress levels (i.e., with different recurrence intervals) after retraining on a limited amount of new data. Our findings suggest that TL models perform better compared to standalone models. Both standalone and TL PINN-estimated RSF parameters and their ground truth values show excellent agreements thus demonstrating that RSF parameters can be retrieved from laboratory stick-slip experiments using the corresponding acoustic data and that the transmitted wave amplitude provides a good representation of the evolving frictional state during stick-slips.

58 GEOSCIENCES↗
DOE OSTI2023

Stress-strain and work hardening relationships of 304L AM alloy

A new approach to analytically derive constitutive stress-strain relationships from modeling the work hardening behavior of alloys was developed for assessing the strength and ductility of the Ti-6Al-4V alloy. This new approach is now successfully applied for assessing the quasi-static stress-strain behavior of an additively manufactured 304L sample. A predictive capability of this modelling approach may then be extended to model material stress-strain behavior at higher strain rates of loading.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Radial deformation of the earth by oceanic tidal loading

A high-degree spherical harmonic series is used to compute the radial deformation of the Earth by oceanic tidal loading. By exploiting fast numerical transforms, this approach is found to be much more efficient, but no less accurate, than the traditional Green's function approach. The method is used to derive an atlas of load tide maps for 10 constitutents of the NSWC ocean tide model.

Ray, R. D.↗
NASA NTRS1925

Wire Suspensions in Wind Tunnel Experiments

The elimination of the rigid supports for models and their replacement by wires constitute a great improvement by rendering negligible the interactions of support. There are disadvantages to wire, namely, the aerodynamic resistance is very large and their use is rather difficult because the whole suspension lacks rigidity and easily becomes distorted. We will here investigate the nature of these distortions, evaluate the errors they entail and describe the methods for taking account of or avoiding them.

WIRE SUSPENSIONS - WIND TUNNEL TESTS↗
NASA NTRS1996

First Cycle Heterogeneous Deformation Behavior and Cyclic Shakedown Phenomena of Nitinol Near A(sub f) Temperatures

Experimental observations on the cyclic behavior of a NiTi alloy (Nitinol) at temperatures in the neighborhood of the A(sub f) (austenite finish) temperature are presented. The strongly heterogeneous nature of the deformation behavior of this material at temperatures within this regime during the first cycle is examined with emphasis placed on the difficulties that the existence of such phenomena pose on the formulation of realistic constitutive relations. It is further demonstrated that this heterogeneity of deformation persists on subsequent cycles with the result that the hysteretic cyclic behavior of these alloys can exhibit a point to point variation in an otherwise uniform geometry. The experimental observations on the deformation behavior of this alloy show that it is strongly dependent on temperature and prior deformation history of the sample, thus resulting in an almost intractable problem with respect to capturing an adequate constitutive description from either experiment or modeling.

Jones, H. N.↗
NASA NTRS2004

A Real-Time Satellite-Based Icing Detection System

Aircraft icing is one of the most dangerous weather conditions for general aviation. Currently, model forecasts and pilot reports (PIREPS) constitute much of the database available to pilots for assessing the icing conditions in a particular area. Such data are often uncertain or sparsely available. Improvements in the temporal and areal coverage of icing diagnoses and prognoses would mark a substantial enhancement of aircraft safety in regions susceptible to heavy supercooled liquid water clouds. The use of 3.9 microns data from meteorological satellite imagers for diagnosing icing conditions has long been recognized (e.g., Ellrod and Nelson, 1996) but to date, no explicit physically based methods have been implemented. Recent advances in cloud detection and cloud property retrievals using operational satellite imagery open the door for real-time objective applications of those satellite datasets for a variety of weather phenomena. Because aircraft icing is related to cloud macro- and microphysical properties (e.g., Cober et al. 1995), it is logical that the cloud properties from satellite data would be useful for diagnosing icing conditions. This paper describes the a prototype realtime system for detecting aircraft icing from space.

Minnis, Patrick↗
NASA NTRS2011

Kalman Orbit Optimized Loop Tracking

Under certain conditions of low signal power and/or high noise, there is insufficient signal to noise ratio (SNR) to close tracking loops with individual signals on orbiting Global Navigation Satellite System (GNSS) receivers. In addition, the processing power available from flight computers is not great enough to implement a conventional ultra-tight coupling tracking loop. This work provides a method to track GNSS signals at very low SNR without the penalty of requiring very high processor throughput to calculate the loop parameters. The Kalman Orbit-Optimized Loop (KOOL) tracking approach constitutes a filter with a dynamic model and using the aggregate of information from all tracked GNSS signals to close the tracking loop for each signal. For applications where there is not a good dynamic model, such as very low orbits where atmospheric drag models may not be adequate to achieve the required accuracy, aiding from an IMU (inertial measurement unit) or other sensor will be added. The KOOL approach is based on research JPL has done to allow signal recovery from weak and scintillating signals observed during the use of GPS signals for limb sounding of the Earth s atmosphere. That approach uses the onboard PVT (position, velocity, time) solution to generate predictions for the range, range rate, and acceleration of the low-SNR signal. The low- SNR signal data are captured by a directed open loop. KOOL builds on the previous open loop tracking by including feedback and observable generation from the weak-signal channels so that the MSR receiver will continue to track and provide PVT, range, and Doppler data, even when all channels have low SNR.

Young, Lawrence E.↗
DOE OSTI2024

Model Formulations: Integrating Distributed Energy Resources (DER) using Advanced Unit Commitment Models and DER Aggregation Methodologies

A distribution energy resource aggregator (DERA) constitutes a group of distribution energy resources with small generation capacities which meet the threshold to participate in the electricity wholesale market. This document provides the proposed DERA model formulation that will be implemented in the SCUC simulation’s architecture for the SCUC-DER project. Different economical assessment methodologies have been developed to incorporated bids for individual distributed resources, which include solar cost dispatch and cost model, BESS opportunity cost offer algorithm, and price sensitive demand response model. Detailed methods are proposed to aggregate individual cost offers to a DERA cost curve to bid in SCUC market while three methods are proposed to simulate DER actual dispatch. Based on the DERA models in this document, the SCUC-DER project will be able to assess the impacts of DERA on the distribution system’s operation and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1999

Third and Final Shuttle Mission of the Isothermal Dendritic Growth Experiment Conducted: Highest Supercooling Ever Recorded Achieved

Dendrites describe the treelike crystal morphology commonly assumed in metals and alloys that freeze from supercooled or supersaturated melts. There remains a high level of engineering interest in dendritic solidification because the size, shape, and orientation of the dendrites determine the final microstructure of a material. It is the microstructure that then determines the physical properties of cast or welded products. Although it is well known that dendritic growth is controlled by the transport of latent heat from the moving solid-liquid interface, an accurate and predictive model has not yet been developed. The effects of gravity-induced convection on the transfer of heat from the interface have prevented adequate testing, under terrestrial conditions, of solidification models. The Isothermal Dendritic Growth Experiment (IDGE) constituted a series of three microgravity experiments flown aboard the Space Shuttle Columbia. The apparatus was used to grow and record dendrite solidification in the absence of gravity-induced convective heat transfer, thereby producing a wealth of benchmark-quality data for testing solidification models and theories.

Glicksman, Martin E.↗
NASA NTRSDate not supplied

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗
NASA NTRS2017

DPW-VI Results Using FUN3D with Focus on k-kL-MEAH2015 (k-kL) Turbulence Model

The Common Research Model wing-body configuration is investigated with the k-kL-MEAH2015 turbulence model implemented in FUN3D. This includes results presented at the Sixth Drag Prediction Workshop and additional results generated after the workshop with a nonlinear Quadratic Constitutive Relation (QCR) variant of the same turbulence model. The workshop provided grids are used, and a uniform grid refinement study is performed at the design condition. A large variation between results with and without a reconstruction limiter is exhibited on "medium" grid sizes, indicating that the medium grid size is too coarse for drawing conclusions in comparison with experiment. This variation is reduced with grid refinement. At a fixed angle of attack near design conditions, the QCR variant yielded decreased lift and drag compared with the linear eddy-viscosity model by an amount that was approximately constant with grid refinement. The k-kL-MEAH2015 turbulence model produced wing root junction flow behavior consistent with wind tunnel observations.

Abdol-Hamid, K. S.↗