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 217 records · Page 12

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗

5-2428: Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior [Slides]

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameter uncertainty and better integrate friction constitutive laws for faults to improve our ability to predict and mitigate seismic hazard at Utah FORGE and future EGS sites. We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, using lab measurements of friction parameters, herein referred to as “THMC+E” models. These simulations will enable exploration of various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗

Thermo-hydro-mechanical analysis of subsurface ice-based thermal energy storage

Ice-based thermal energy storage systems are widely utilized for cooling and managing peak electrical demand globally, offering daily or weekly storage capabilities for both individual homes and larger office buildings. However, scaling these systems for district-level cooling or integrating them with renewable energy sources presents challenges, especially in accommodating larger volumes and addressing seasonal storage requirements in densely populated urban areas. This paper proposes a novel solution by evaluating subsurface ice-based thermal energy storage, in which the underground is subjected to seasonal freeze/thaw cycles. However, these cycles may influence ground behavior, affecting pore pressure and inducing ground movement. To systematically investigate these challenges, we enhance the TOUGH-FLAC simulator by integrating water/ice phase change capabilities and updating the effective stress–strain constitutive relation. Both modifications are validated against analytical solutions or experimental data. Through numerical simulations spanning a decade with ten seasonal freeze/thaw cycles, we evaluate the performance and long-term stability of a generic subsurface ice-based thermal energy storage system, considering factors such as ground permeability, freezing pipe spacing, freeze/thaw damage, and glycol solution temperature. The simulations indicate that ice formation induces pore pressure variations that drive seasonal surface heave and settlement, controlled by ground permeability, pipe spacing, and glycol solution temperature, along with tensile and localized shear deformation around freeze pipes. This highlights the need for accurate ground property characterization and geomechanical analysis for subsurface ice-based thermal energy storage.

Ground freezing↗

Near-constant retreat rate of a terrestrial margin of the Laurentide Ice Sheet during the last deglaciation

The Laurentide Ice Sheet (LIS) was the largest ice sheet during the last glacial period. An accurate representation of its behavior during the last deglaciation is critical to understanding its influence on and response to a changing climate. Here we use 10 Be dating and Bayesian modeling to track the recession of the southwest sector of the Labrador Dome of the LIS along an ~500-km-long transect west of Lake Superior during the last deglaciation. This transect reflects terrestrial ice-margin retreat and crosses multiple moraine sets, with the southwestern part of the transect deglaciated by ca. 19 ka and the northeastern part deglaciated by ca. 10 ka. The predominant behavior of the ice margin during this interval is near-constant retreat with retreat rates varying between ~59 m/a and 38 m/a. The moraine sets mark standstills and/or readvances that in total constitute only ~17% of the retreat interval. The spatial and temporal pattern of ice-margin retreat tracked here differs from existing reconstructions that are based on using isochrons to define ice-margin positions. Acknowledging the uncertainties associated with the modeled ages of ice-margin retreat, we suggest that the overall retreat pattern is consistent with forcing by a gradual increase in Northern Hemisphere, high-latitude summer insolation. The pattern of ice-margin retreat is inconsistent with Greenland ice-core temperature records, and thus these records may not be suitable to drive models of the LIS.

58 GEOSCIENCES↗

Uncertainty Quantification for Component Modeling Using the Discrete-Direct Approach

Threaded fastener behavior can be an important aspect of complex component and system behavior, but there is no one-size-fits-all finite element analysis technique. Proper modeling of threaded fastener joints requires careful consideration of many details, from test setup and data acquisition to constitutive modeling and uncertainty quantification approaches. This report details analysis of a “mini-radax” bolted-joint exemplar where a Discrete-Direct uncertainty quantification approach is employed to evaluate margin of the component. The mini-radax geometry is tested to failure on a drop table, and single-coupon tests of individual fasteners serve as foundational data for the analysis. Analysis predictions complement the test data well and provide additional context for engineering decision-making.

42 ENGINEERING↗

Online and Offline Data Quality Monitoring for the Mu2e Calorimeter

This thesis presents the design, implementation, and validation of a calorimeter Data Quality Monitoring (DQM) toolchain for the Mu2e experiment at Fermilab. Mu2e searches for charged lepton flavor violation via coherent muon-to-electron conversion in the field of an aluminum nucleus, $\mu^- Al \rightarrow e^-Al$, a process whose observation would constitute clear evidence of physics beyond the Standard Model. Achieving target sensitivity requires stringent control of detector performance and data integrity during acquisition, as subtle issues in readout configuration, data formatting, or electronics behavior can compromise reconstruction and bias downstream analyzes. To address these challenges, this work develops a multi-layer DQM approach spanning both raw data validation and reconstructed digi-level diagnostics. At the low level, a fragment analysis component performs word- and bit-field decoding of calorimeter readout blocks, enabling sanity checks of the expected structure and producing detailed error and integrity statistics useful for commissioning and troubleshooting. At the digi level, the CaloDigiDQM analyzer is implemented within the art framework and transforms each CaloDigiCollection into a structured hierarchy of ROOT histograms designed for fast drill-down diagnostics. The module generates coherent monitoring views at global, disk, board, and channel granularity, including occupancy, waveform-derived features (baseline, RMS, peak amplitude and position), and left-right sensor consistency metrics. Detector-aware channel-to-electronics mapping is performed through the conditions system (CaloDAQMap), ensuring that diagnostics remain aligned with hardware identifiers used in operations. For end-to-end testing without reliance on live DAQ data, a synthetic CaloDigi producer is developed to generate realistic waveforms with controlled noise and pulse shapes. The resulting system supports both offline ROOT-file production and online operation, including optional histogram streaming through otsdaq via ots::HistoSender. This toolchain provides a practical and scalable foundation for calorimeter commissioning and stable data collection, enabling early detection of anomalies and reducing operational risk for Mu2e.

Vakulenko, Mark [Drew U.] (ORCID:0009000276197818)↗

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

On the peridynamic effective force state and multiphase constitutive correspondence principle

This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Finally, numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.

42 ENGINEERING↗

Supercritical fluids behave as complex networks

Abstract Supercritical fluids play a key role in environmental, geological, and celestial processes, and are of great importance to many scientific and engineering applications. They exhibit strong variations in thermodynamic response functions, which has been hypothesized to stem from the microstructural behavior. However, a direct connection between thermodynamic conditions and the microstructural behavior, as described by molecular clusters, remains an outstanding issue. By utilizing a first-principles-based criterion and self-similarity analysis, we identify energetically localized molecular clusters whose size distribution and connectivity exhibit self-similarity in the extended supercritical phase space. We show that the structural response of these clusters follows a complex network behavior whose dynamics arises from the energetics of isotropic molecular interactions. Furthermore, we demonstrate that a hidden variable network model can accurately describe the structural and dynamical response of supercritical fluids. These results highlight the need for constitutive models and provide a basis to relate the fluid microstructure to thermodynamic response functions.

42 ENGINEERING↗

Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

We develop a flexible, relativistically covariant parametrization of the dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code smash, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second-order cumulant ratio, and constitute a clear signature of the crossover region above the critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

User-based I/O Profiling for Leadership Scale HPC Workloads

I/O constitutes a significant portion of most of the application run-time. Spawning many such applications concurrently on an HPC system leads to severe I/O contention. Thus, understanding and subsequently reducing I/O contention induced by such multi-tenancy is critical for the efficient and reliable performance of the HPC system. In this study, we demonstrate that an application’s performance is influenced by the command line arguments passed to the job submission. We model an application’s I/O behavior based on two factors: past I/O behavior within a time window and user-configured I/O settings via command-line arguments. We conclude that I/O patterns for well-known HPC applications like E3SM and LAMMP are predictable, with an average uncertainty below 0.25 (A probability of 80%) and near zero (A probability of 100%) within a day. However, I/O pattern variance increases as the study time window lengthens. Additionally, we show that for 38 users and at least 50 applications constituting approximately 93000 job submissions, there is a high correlation between a submitted command line and the past command lines made within 1 to 10 days submitted by the user. We claim the length of this time window is unique per user.

Yazdani, Ahmad Hossein [Virginia Tech, Blacksburg]↗

Accelerated Creep Test Qualification of Creep-Resistance Using the Wilshire–Cano–Stewart Constitutive Model and Stepped Isostress Method

In this paper, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire–Cano–Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing is a long-term continuous process; in fact, the ASME B&PV III requires that 10,000+ h of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The SSM has been developed, which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750 °C and predictions are compared to conventional creep testing data. The WCS model has proven to make long-term predictions for stress-rupture, MCSR, creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are postaudit validated using the National Institute of Material Science database. Similarly, the MCSR predictions are compared from previous studies. Finally, the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and industrial gas turbines.

42 ENGINEERING↗

Tonoplast Sucrose Trafficking Modulates Starch Utilization and Water Deficit Behavior in Poplar Leaves

Leaf osmotic adjustment by the active accrual of compatible organic solutes (e.g. sucrose) contributes to drought tolerance throughout the plant kingdom. In Populus tremula x alba, PtaSUT4 encodes a tonoplast sucrose–proton symporter, whose downregulation by chronic mild drought or transgenic manipulation is known to increase leaf sucrose and turgor. While this may constitute a single drought tolerance mechanism, we now report that other adjustments which can occur during a worsening water deficit are damped when PtaSUT4 is constitutively downregulated. Specifically, we report that starch use and leaf relative water content (RWC) dynamics were compromised when plants with constitutively downregulated PtaSUT4 were subjected to a water deficit. Leaf RWC decreased more in wild-type and vector control lines than in transgenic PtaSUT4-RNAi (RNA-interference) or CRISPR (clustered regularly interspersed short palindromic repeats) knockout (KO) lines. The control line RWC decrease was accompanied by increased PtaSUT4 transcript levels and a mobilization of sucrose from the mesophyll-enriched leaf lamina into the midvein. The findings suggest that changes in SUT4 expression can increase turgor or decrease RWC as different tolerance mechanisms to reduced water availability. Evidence is presented that PtaSUT4-mediated sucrose partitioning between the vacuole and the cytosol is important not only for overall sucrose abundance and turgor, but also for reactive oxygen species (ROS) and antioxidant dynamics. Interestingly, the reduced capacity for accelerated starch breakdown under worsening water-deficit conditions was correlated with reduced ROS in the RNAi and KO lines. A role for PtaSUT4 in the orchestration of ROS, antioxidant, starch utilization and RWC dynamics during water stress and its importance in trees especially, with their high hydraulic resistances, is considered.

59 BASIC BIOLOGICAL SCIENCES↗

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗

Polycrystal thermo-elasticity revisited: theory and applications

The self-consistent (SC) theory is the most commonly used mean-field homogenization method to estimate the mechanical response behavior of polycrystals based on the knowledge of the properties and orientation distribution of constituent single-crystal grains. The original elastic SC method can be extended to thermo-elasticity by adding a stress-free strain to an elastic constitutive relation that expresses stress as a linear function of strain. With the addition of this independent term, the problem remains linear. Although the thermo-elastic self-consistent (TESC) model has important theoretical implications for the development of self-consistent homogenization of non-linear polycrystals, in this paper, we focus on TESC applications to actual thermo-elastic problems involving non-cubic (i.e. thermally anisotropic) materials. To achieve this aim, here we provide a thorough description of the TESC theory, which is followed by illustrative examples involving cooling of polycrystalline non-cubic metals. The TESC model allows studying the effect of crystallographic texture and single-crystal elastic and thermal anisotropy on the effective thermo-elastic response of the aggregate and on the internal stresses that develop at the local level.

36 MATERIALS SCIENCE↗

ACCELERATED CREEP TEST (ACT) QUALIFICATION OF CREEP RESISTANCE USING THE WCS CONSTITUTIVE MODEL AND STEPPED ISOSTRESS METHOD (SSM)

In this study, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire-Cano-Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing (CCT) is a long-term continuous process, in fact, the ASME B&PV III requires that 10,000+ hours of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle (TTSSP) that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The stepped isostress method (SSM) has been developed which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750°C and predictions are compared to CCT data. The WCS model has proven to make long-term predictions for stress-rupture, minimum-creep-strain-rate (MCSR), creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are post-audit validated using the National Institute of Material Science (NIMS) database. Similarly, the MCSR predictions are compared from previous studies. Finally the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and Industrial Gas Turbines (IGT).

36 MATERIALS SCIENCE↗

A Sobering Examination of the Feasibility of Aqueous Aluminum Batteries

Recent reports on aqueous aluminum (Al) batteries boast encouraging metrics including high energy density, nonflammability, and low raw material cost. However, the poorly understood solvation environment around the trivalent aluminum ions (Al3+) raises uncertainty regarding their feasibility. Furthermore, the corrosiveness of the acidic solutions and degree of hydration of aluminum trifluoromethanesulfonate (Al(OTf)3) salt is routinely mischaracterized. This study provides a rigorous examination of aqueous Al electrolytes, with the first compelling evidence for a dynamic octahedral solvation structure around Al3+ dominated by labile water and OH-, without Al-OTf contact ion pairs, at high concentrations. This solvation behavior and the concomitant electrostatic hydrolysis differs greatly from previously reported water-in-salt electrolytes, and occurs due to the high charge density of the Lewis acidic Al3+. Nuclear magnetic resonance spectroscopy and other physicochemical measurements quantitatively reveal how species activity evolves with concentration. The high proton activity observed in transport and electrochemical measurements constitutes a prohibitive barrier to the feasibility of rechargeable aluminum batteries with aqueous electrolyte.

Pastel, Glenn R.↗