Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive behavior”

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

Observation of a massive phason in a charge-density-wave insulator

The lowest-lying fundamental excitation of an incommensurate charge density wave (CDW) material is widely believed to be a massless phason -- a collective modulation of the phase of the CDW order parameter. However, as first pointed out by Lee and Fukuyama, long-range Coulomb interactions should push the phason energy up to the plasma energy of the CDW condensate, resulting in a massive phason and a fully gapped spectrum. Whether such behavior occurs in a CDW system has been unresolved for more than four decades. Using time-domain THz emission spectroscopy, we investigate this issue in the material (TaSe 4 ) 2 I, a classical example of a quasi-one-dimensional CDW insulator. Upon transient photoexcitation at low temperatures, we find the material strikingly emits coherent, narrow-band THz radiation. The frequency, polarization and temperature-dependence of the emitted radiation imply the existence of a phason that acquires mass by coupling to long-range Coulomb interaction. Our observations constitute the first direct evidence of the massive "Lee-Fukuyama" phason and highlight the potential applicability of fundamental collective modes of correlated materials as compact and robust sources of THz radiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A comparison of plastic and visco-plastic constitutive models for accurate predictions in PM-HIP

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that can efficiently produce near net shape parts with high material utilization and uniform microstructures. While the PM-HIP process is used frequently to produce small-scale components, its application to large-scale components is still limited due to inadequate understanding of its complex processes that cause unpredictable post-HIP shape distortions. A computational model can be a cost-effective alternative to exhaustive experimentation required to understand the PM-HIP process fully. Therefore, we present a comparison of plastic and visco-plastic models that are frequently used to model the PM-HIP process. We show that the predictions of powder densification behavior and shape distortions obtained from both the models are similar under most conditions and agree well with experimental observations.

Sarkar, Subrato [ORNL]↗
DOE OSTI2024

Dissimilar Material Joining via Interlocking Metasurfaces

Background The integration of dissimilar materials poses a significant challenge in engineering, necessitating innovative solutions for robust and reliable joining. Interlocking metasurfaces (ILMs) are a new joining technology comprising arrays of autogenous features patterned across two surfaces that interlock to form robust structural joints. Objective Here, this study elucidates the factors influencing the tensile performance of ILM joints formed between dissimilar materials. Methods We employed parametric optimization to identify optimal unit cell geometries for maximal yield strength based on the hypothesis that the elastic tensile properties of the materials are the primary determinants of tensile performance. Experimental validation was performed by mechanically testing the theorized optimal ILM geometry and a range of ILM geometries to capture the overall behavior trends of joints between two additively manufactured polymers, VeroPureWhite (VW) and RGDA8430-DM (8430). Results Experimental validation of optimized designs revealed that additional factors, e.g. flexural strength and localized plasticity, also strongly influenced the tensile performance of T-slot ILMs joining dissimilar materials. The proposed optimal design remained the best performer. Conclusions This study demonstrates the viability of ILMs as a joining method for dissimilar materials. ILMs can join dissimilar materials with no loss in joint yield strength compared to joints composed solely of the weaker of the two constitutive materials. ILMs demonstrated their potential as a versatile and effective joining technology in diverse engineering applications.

Elbrecht, Benjamin James [Sandia National Laborato↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗
DOE OSTI2024

Exploring structure-function relationships governing transport by and regulation of the major cyanobacterial bicarbonate transporters SbtA and BicA

Cyanobacteria play a large and essential role in the global carbon cycle, while also representing promising microbial platforms for the direct bioproduction of fuels and chemicals from solar energy and carbon dioxide (CO 2 ). In both cases, a critical step in the overall photosynthetic process is the transport of inorganic carbon (Ci) across the cell membrane, which occurs via different mechanisms and associated transporter proteins that target either CO 2 or bicarbonate (HCO 3 - ). In the latter case, BicA (constitutively-expressed, high-flux/low-affinity Na + /HCO 3 - symporter) and SbtA (low-flux/high-affinity Na + /HCO 3 - symporter, expressed only under low CO 2 conditions) have been identified as the two main cyanobacterial HCO 3 - transporters. This project entailed combining detailed functional and high-resolution structural studies to help fill key knowledge gaps regarding the mechanisms of HCO 3 - transport by BicA and SbtA. Specifically, major goals of this project included: i) identifying key residues and structural features of BicA and SbtA that are responsible for modulating their uptake flux and affinity behaviors, and ii) obtaining high-resolution structures of full-length BicA and SbtA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Nanotubes from the Misfit Layered Compound (SmS) 1.19 TaS 2 : Atomic Structure, Charge Transfer, and Electrical Properties

Misfit layered compounds (MLCs) MX-TX 2 , where M, T = metal atoms and X = S, Se, or Te, and their nanotubes are of significant interest due to their rich chemistry and unique quasi-1D structure. In particular, LnX-TX 2 (Ln = rare-earth atom) constitute a relatively large family of MLCs, from which nanotubes have been synthesized. The properties of MLCs can be tuned by the chemical and structural interplay between LnX and TX 2 sublayers and alloying of each of the Ln, T, and X elements. In order to engineer them to gain desirable performance, a detailed understanding of their complex structure is indispensable. MLC nanotubes are a relative newcomer and offer new opportunities. In particular, like WS 2 nanotubes before, the confinement of the free carriers in these quasi-1D nanostructures and their chiral nature offer intriguing physical behavior. High-resolution transmission electron microscopy in conjunction with a focused ion beam are engaged to study SmS-TaS 2 nanotubes and their cross-sections at the atomic scale. The atomic resolution images distinctly reveal that Ta is in trigonal prismatic coordination with S atoms in a hexagonal structure. Furthermore, the position of the sulfur atoms in both the SmS and the TaS 2 sublattices is revealed. X-ray photoelectron spectroscopy, electron energy loss spectroscopy, and X-ray absorption spectroscopy are carried out. These analyses conclude that charge transfer from the Sm to the Ta atoms leads to filling of the Ta 5d z 2 level, which is confirmed by density functional theory (DFT) calculations. Transport measurements show that the nanotubes are semimetallic with resistivities in the range of 10 –4 Ω·cm at room temperature, and magnetic susceptibility measurements show a superconducting transition at 4 K.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Regulations for Bat Protection in Mexico's Wind Farms

Wind energy development has expanded the fastest globally among all renewable sources during the last 20 years. However, wind farms have documented adverse impacts on bats, including mortality from collisions with turbine blades and disruptions to habitat and behavior. As the world's sixth most attractive economy for renewables, with 70 operating wind farms, Mexico and its bats now face escalating threats from the country's burgeoning wind industry. Despite this rapid growth, few studies have analyzed Mexico's regulatory framework to prevent, evaluate, and mitigate wind farm effects on bats. In this study, we reviewed Mexican laws and treaties that facilitate wind farm permitting, construction, operation, and decommissioning, and searched for guidelines that specifically address bat conservation. We found eight international pacts that promote wind power adoption along with three relevant articles in Mexico's Constitution. The General Law of Ecological Balance and Environmental Protection proved most pertinent for impact management. Supplementary guidelines from the Ministry of Environment and Natural Resources offer general strategies for evaluating wind farm impacts on bats, but adherence remains voluntary. Given expanding wind power investments across Mexico, we highlight the need for more stringent national standards that require preventative and corrective measures to protect bat populations. Tighter legislation and enforcement offer pathways toward environmentally sustainable wind energy development in Mexico.

17 WIND ENERGY↗
DOE OSTI2020

An integrated computational materials engineering framework to analyze the failure behaviors of carbon fiber reinforced polymer composites for lightweight vehicle applications

A bottom-up multi-scale modeling approach is used to develop an Integrated Computational Materials Engineering (ICME) framework for carbon fiber reinforced polymer (CFRP) composites, which has the potential to reduce development to deployment lead time for structural applications in lightweight vehicles. In this work, we develop and integrate computational models comprising of four size scales to fully describe and characterize three types of CFRP composites. In detail, the properties of the interphase region are determined by an analytical gradient model and molecular dynamics analysis at the nano-scale, which is then incorporated into micro-scale unidirectional (UD) representative volume element (RVE) models to characterize the failure strengths and envelopes of UD CFRP composites. Then, the results are leveraged to propose an elasto-plastic-damage constitutive law for UD composites to study the fiber tows of woven composites as well as the chips of sheet molding compound (SMC) composites. Subsequently, the failure mechanisms and failure strengths of woven and SMC composites are predicted by the meso-scale RVE models. Finally, building upon the models and results from lower scales, we show that a homogenized macro-scale model can capture the mechanical performance of a hat-section-shaped part under four-point bending. Along with the model integration, we will also demonstrate that the computational results are in good agreement with experiments conducted at different scales. The present work illustrates the potential and significance of integrated multi-scale computational modeling tools that can virtually evaluate the performance of CFRP composites and provide design guidance for CFRP composites used in structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Compositional analysis of organosolv poplar lignin by using high-performance liquid chromatography/high-resolution multi-stage tandem mass spectrometry

Organosolv treatment is an efficient and environmentally friendly process to degrade lignin into small compounds. The capability of characterizing the individual compounds in the complex mixtures formed upon organosolv treatment is essential for the optimization of the further lignin conversion processes and for the rational genetic engineering of plants used to produce lignin in order to improve lignin properties. In this study, an organosolv poplar lignin sample was initially analyzed by high-resolution mass spectrometry coupled with negative-ion mode electrospray ionization ((-)ESI HRMS). Lignin monomers and dimers were found to constitute the majority of the compounds in the organosolv lignin sample. Larger lignin oligomers, such as trimers and tetramers, and some not lignin-related compounds, were also detected. A high-performance liquid chromatograph/linear quadrupole ion trap/orbitrap mass spectrometer capable of multi-stage high-resolution tandem mass spectrometry experiments (HRMS n ), equipped with an (-)ESI source (HPLC/(-)ESI HRMS n ), was employed to separate the unknown compounds in the organosolv mixture and to obtain structural information for the deprotonated compounds via collision-activated dissociation (CAD) HRMS n experiments. To improve the understanding of the CAD behavior of deprotonated lignin-related compounds, 16 deprotonated model compounds with different functionalities and linkage types were examined. This approach enabled the assignment of likely structures for several lignin monomers, dimers, trimers, and tetramers, and some not lignin-related compounds, most likely fatty acids. Based on the proposed structures, compounds in the organosolv lignin sample contain β-O-4, 5–5, β-5, and possibly also 4-O-5 linkages. Most compounds contain G- and S-monomeric units although a small amount of H-units were also detected.

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

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2023

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around the {N(CH 3 ) 3 } + and the C=O functional groups of the PMETAC chain and one around the Cl - ion. These first hydration shells around the respective moieties overlap and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms/groups around which they are located. For example, the presence of the methyl groups make the {N(CH 3 ) 3 } + group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH 3 ) 3 } + moiety. These water molecules around the N(CH 3 ) 3 } + group also have enhanced tetrahedrality as compared to the water molecules constituting the hydration layer around the C=O group and the Cl - counterion. Our simulations also identify that there is an intervening water layer between the Cl - ion and {N(CH 3 ) 3 } + group: this layer prevents the Cl - ion from coming very close to the {N(CH 3 ) 3 } + group. As a consequence, there is a significantly large mobility of the Cl - ions inside the PMETAC brush layer. Furthermore, the C=O group of the PE chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. As a result, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site, apolar hydration at another site, and ensures large mobility of the supported Cl- counterions.

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

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2026

Input specific neural networks

Neural networks have emerged as powerful tools for mapping between inputs and outputs. However, their black-box nature limits the ability to encode or impose specific structural relationships between inputs and outputs. Many scientific and engineering problems, such as constitutive modeling in solid mechanics, require networks that can enforce convexity, monotonicity, or other structural constraints to ensure physical consistency. Here, we introduce the Input Specific Neural Network (ISNN), a new architecture that enables multiple, distinct constraints to be imposed on different input subsets for scalar-valued outputs. This framework unifies convex, monotone–convex, monotone, and arbitrary mappings within a single network for the first time. Two ISNN architectures with analytical first- and second-order derivatives are developed. We demonstrate the performance on synthetic toy problems, inverse problems in isotropic hyperelasticity, and finite element simulations. ISNNs achieve improved extrapolation behavior, require fewer invariant inputs than standard input convex networks for polyconvex potentials, and enable significant computational savings via manual differentiation. We also show how ISNNs can be used to learn structural relationships between inputs and outputs via a binary gating mechanism. Particularly, ISNNs are employed to model a homogenized anisotropic free energy potential in a decoupled multiscale setting. The network learns whether or not the potential should be modeled as polyconvex and retains only the relevant layers while using the minimum number of inputs. ISNNs provide a flexible foundation for embedding structural priors into neural networks, enhancing both interpretability and stability. They are broadly applicable across computational mechanics and other scientific domains requiring constrained functional relationships.

Jadoon, Asghar A. [Univ. of Texas, Austin, TX (Uni↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗