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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗

DOE OSTI2021

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗

NASA NTRS2004

Beyond Mechanism Design

The field of mechanism design is concerned with setting (incentives superimposed on) the utility functions of a group of players so as to induce desirable joint behavior of those players. It arose in the context of traditional equilibrium game theory applied to games involving human players. This has led it to have many implicit restrictions, which strongly limits its scope. In particular, it ignores many issues that are crucial for systems that are large (and therefore far off-equilibrium in general) and/or composed of non-human players (e.g., computer-based agents). This also means it has concentrated on issues that are often irrelevant in those broader domains (e.g., incentive compatibility). This paper illustrates these shortcomings by reviewing some of the recent theoretical work on the design of collectives, a body of work that constitutes a substantial broadening of mechanism design. It then presents computer experiments based on a recently suggested nanotechnology testbed that demonstrates the power of that extended version of mechanism design.

Wolpert, David H.↗

DOE OSTI2020

Initial High Temperature Inelastic Constitutive Model for Alloy 617

This report describes progress on a Alloy 617 inelastic constitutive model slated for inclusion in the ASME Boiler & Pressure Vessel Code in Code Case N-898 “Use of Alloy 617 (UNS N06617) for Class A Elevated Temperature Construction, Section III, Division 5” as part of an appendix on inelastic modeling. This appendix will provide a reference material model for A617 for use in conjunction with the ASME design by inelastic analysis criteria. This report focuses on the implementation and evaluation of an existing Alloy 617 model developed to represent the data collected at Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) as part of the Alloy 617 Code Case experimental campaign. While the model accurately captures the cyclic stress/strain behavior of Alloy 617, it fails to capture the long-term creep and lower temperature tensile data. In addition, the model form is complicated and has over 50 configurable parameters. As such, we recommend the development of a simplified model, incorporating key features of the ORNL model but reducing the model complexity.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

The Role of Unit-Cell Topology in Modulating the Compaction Response of Additively Manufactured Cellular Materials using Simulations and Validation Experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗

NASA NTRS1994

Attention focusing and anomaly detection in systems monitoring

Any attempt to introduce automation into the monitoring of complex physical systems must start from a robust anomaly detection capability. This task is far from straightforward, for a single definition of what constitutes an anomaly is difficult to come by. In addition, to make the monitoring process efficient, and to avoid the potential for information overload on human operators, attention focusing must also be addressed. When an anomaly occurs, more often than not several sensors are affected, and the partially redundant information they provide can be confusing, particularly in a crisis situation where a response is needed quickly. The focus of this paper is a new technique for attention focusing. The technique involves reasoning about the distance between two frequency distributions, and is used to detect both anomalous system parameters and 'broken' causal dependencies. These two forms of information together isolate the locus of anomalous behavior in the system being monitored.

Doyle, Richard J.↗

DOE OSTI2024

The role of unit cell topology in modulating the compaction response of additively manufactured cellular materials using simulations and validation experiments

36 MATERIALS SCIENCE↗

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.↗

DOE OSTI2025

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗

NASA NTRS2001

Spatiotemporal relationships between growth and microtubule orientation as revealed in living root cells of Arabidopsis thaliana transformed with green-fluorescent-protein gene construct GFP-MBD

Arabidopsis thaliana plants were transformed with GFP-MBD (J. Marc et al., Plant Cell 10: 1927-1939, 1998) under the control of a constitutive (35S) or copper-inducible promoter. GFP-specific fluorescence distributions, levels, and persistence were determined and found to vary with age, tissue type, transgenic line, and individual plant. With the exception of an increased frequency of abnormal roots of 35S GFP-MBD plants grown on kanamycin-containing media, expression of GFP-MBD does not appear to affect plant phenotype. The number of leaves, branches, bolts, and siliques as well as overall height, leaf size, and seed set are similar between wild-type and transgenic plants as is the rate of root growth. Thus, we conclude that the transgenic plants can serve as a living model system in which the dynamic behavior of microtubules can be visualized. Confocal microscopy was used to simultaneously monitor growth and microtubule behavior within individual cells as they passed through the elongation zone of the Arabidopsis root. Generally, microtubules reoriented from transverse to oblique or longitudinal orientations as growth declined. Microtubule reorientation initiated at the ends of the cell did not necessarily occur simultaneously in adjacent neighboring cells and did not involve complete disintegration and repolymerization of microtubule arrays. Although growth rates correlated with microtubule reorientation, the two processes were not tightly coupled in terms of their temporal relationships, suggesting that other factor(s) may be involved in regulating both events. Additionally, microtubule orientation was more defined in cells whose growth was accelerating and less stringent in cells whose growth was decelerating, indicating that microtubule-orienting factor(s) may be sensitive to growth acceleration, rather than growth per se.

NASA Discipline Plant Biology↗

NASA NTRS2010

Oxidation Characterization of Hafnium-Based Ceramics Fabricated by Hot Pressing and Electric Field-Assisted Sintering

Ceramic borides, such as hafnium diboride (HfB2) and zirconium diboride (ZrB2), are members of a family of materials with extremely high melting temperatures referred to as Ultra High Temperature Ceramics (UHTCs). UHTCs constitute a class of promising materials for use in high temperature applications, such as sharp leading edges on future-generation hypersonic flight vehicles, because of their high melting points. The controlled development of microstructure has become important to the processing of UHTCs, with the prospect of improving their mechanical and thermal properties. The improved oxidation resistance of HfB2 has also become important if this material is to be successfully used at temperatures above 2000 C. Furthermore, the use of UHTCs on the leading edges of vehicles traveling at hypersonic speeds will mean exposure to a mixed oxidation environment comprised of both molecular and atomic oxygen. The current study has investigated the high-temperature oxidation behavior of HfB2-based materials in a pure O2 environment, as well as in environments containing different levels of dissociated oxygen (O/O2). Materials were processed by two techniques: conventional hot pressing (HP) and electric field-assisted sintering (FAS). Their oxidation behavior was evaluated in both a tube furnace at 1250 C for 3 hours and in a simulated re-entry environment in the Advanced Heating Facility (AHF) arcjet at NASA Ames Research Center, during a 10-minute exposure to a cold wall heat flux of 250W/sq cm and stagnation pressure of 0.1-0.2 atm. The microstructure of the different materials was characterized before and after oxidation using scanning electron microscopy (SEM).

Gasch, Matt↗

NASA NTRS2001

Monitoring Programs Using Rewriting

We present a rewriting algorithm for efficiently testing future time Linear Temporal Logic (LTL) formulae on finite execution traces, The standard models of LTL are infinite traces, reflecting the behavior of reactive and concurrent systems which conceptually may be continuously alive in most past applications of LTL, theorem provers and model checkers have been used to formally prove that down-scaled models satisfy such LTL specifications. Our goal is instead to use LTL for up-scaled testing of real software applications, corresponding to analyzing the conformance of finite traces against LTL formulae. We first describe what it means for a finite trace to satisfy an LTL property end then suggest an optimized algorithm based on transforming LTL formulae. We use the Maude rewriting logic, which turns out to be a good notation and being supported by an efficient rewriting engine for performing these experiments. The work constitutes part of the Java PathExplorer (JPAX) project, the purpose of which is to develop a flexible tool for monitoring Java program executions.

Havelund, Klaus↗

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↗

DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗

DOE OSTI2024

CIRCLES: Congestion Impacts Reduction via CAV-in-the-loop Lagrangian Energy Smoothing

The energy efficiency of today’s vehicular mobility relies on the un-integrated combination of i) control via static assets (traffic lights, metering, variable speed limits, etc.); and ii) onboard vehicle automation (adaptive cruise control (ACC), ecodriving, etc.). These two families of control were not co-designed and are not engineered to work in coordination. Recent studies have shown i) limitations of controls, and even ii) negative impacts of ACC. This project focused on the technology development, implementation and prototyping, and validation of Mobile Traffic Control (MTC). MTC can be viewed as an extension of classical traffic control (in which static infrastructure actuates traffic flow). In the MTC paradigm, automated vehicles actuate the entire flow via their behavior, offering enhanced possibilities to optimize the energy footprint of traffic, if designed correctly. We set out to demonstrate for the first time that considerably reduced fuel consumption of all vehicles in traffic can be achieved via distributed control of a small proportion of CAVs. Compared to baseline vehicular technologies, our work offers a significant design departure: control algorithms for the CAVs consider the impact one vehicle can have on overall traffic, improving resulting overall fuel consumption. We focus on using a few vehicles (as traffic controllers via CAV technology) to improve the energy efficiency of traffic flow to further optimize energy efficiency. A live-traffic demonstration in November 2022 featured the deployment of 100 specially-equipped CAVs, constituting an approximate local penetration rate upwards of 2%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

CHIPFoam Modification to Allow for Modelling SWIFT Foams

SWIFT (Silicon-water in Familiar Template) foams are silicone foams that have a continuous porosity without having any obvious spherical or ellipsoidal pore shapes. Like other foamed rubber, SWIFT foams can experience large reversible deformations. Because of their microstructure, SWIFT foams have a stiffering behavior in compression that is very abrupt compared to most other foamed rubbers. Work by Benedikt reported by Miller showed that fitting compressive test data with the CHIPFoam model in its current form is impractical and inaccurate. This report describes an extension made to CHIPFoam to allow the modeling of SWIFT foams. While the extension allows the modeling of SWIFT foams, it also extends the range of porosity of foams that can be modelled with CHIPFoam to include very high porosity (>99%) foamed rubber. The previous model was unstable for porosities greater than 70%. The changes made to CHIPFoam include the introduction of a critical porosity at which the compressible Danielsson strain energy function that couples volumetric and deviatoric mechanical behavior begins to increase with compression. This new feature also required the introduction of a third term in the multiplicative decomposition of the relative volume, J . The organization of this report is described here. The first section, SWIFT Foam Structure and Mechanical Behavior, describes the structure of SWIFT foams as contrasted to other foamed rubbers. It also contrasts the mechanical behavior of SWIFT foams to other more traditional foamed rubbers. The second section, CHIPFoam Modifications, describes the specific modifications done to be able to modify the delayed rapid compressive stiffening seen in SWIFT foams. The third section, Results, demonstrates the use of the modified model to fit compressive loading response data from a SWIFT foam sample. The fourth section, Conclusions and Further Work, puts the work in perspective and presents likely further work.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Electrochemical Investigation of Low-Valent Multiply M≡M Bonded Group VI Dimers: A Standard Chemical Reduction Leads to an Unexpected Product

The synthesis and investigation of multimetallic complexes bearing metal–metal bonds have been significantly advanced over the last 30 years with exciting applications across many fields. While the outlook for synthetically innovative metal–metal species is promising, the reactivity of the classic M 2 X 6 alkyl dimers (where M = W, Mo; X = CH 2 CMe 3 , CH 2 SiMe 3 ) is still not fully understood. In particular, the redox characteristics of these complexes have not been explored by cyclic voltammetry (CV). In this work, W 2 (CH 2 CMe 3 ) 6 and Mo 2 (CH 2 CMe 3 ) 6 (1 and 2, respectively) were characterized using CV to show differing electrochemical behavior between the two species. Stoichiometric reactions, guided by the results of the CV experiments, led to the isolation of an alkylidyne-bridged dimer bearing a rare structural motif, W 2 (μ-CCMe 3 ) 2 (CH 2 CMe 3 ) 4 (3). Single-crystal X-ray diffraction (SCXRD) and 1 H and 13 C nuclear magnetic resonance spectroscopy (NMR) were used to establish the connectivity of the structure in the solid and solution phases. Density functional theory (DFT) calculations on 1–3 and the previously reported W 2 (CSiMe 3 ) 2 (CH 2 SiMe 3 ) 4 (4) were used to rationalize the reactivity of 1 to form 3 and the structural differences between 3 and 4. Complexes 3 and 4 constitute a case in which the stereoelectronic properties of the silyl neopentyl ligand can affect differences in structure and bonding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2007

State-Chart Autocoder

A computer program translates Unified Modeling Language (UML) representations of state charts into source code in the C, C++, and Python computing languages. ( State charts signifies graphical descriptions of states and state transitions of a spacecraft or other complex system.) The UML representations constituting the input to this program are generated by using a UML-compliant graphical design program to draw the state charts. The generated source code is consistent with the "quantum programming" approach, which is so named because it involves discrete states and state transitions that have features in common with states and state transitions in quantum mechanics. Quantum programming enables efficient implementation of state charts, suitable for real-time embedded flight software. In addition to source code, the autocoder program generates a graphical-user-interface (GUI) program that, in turn, generates a display of state transitions in response to events triggered by the user. The GUI program is wrapped around, and can be used to exercise the state-chart behavior of, the generated source code. Once the expected state-chart behavior is confirmed, the generated source code can be augmented with a software interface to the rest of the software with which the source code is required to interact.

Clark, Kenneth↗