SEARCH · Engineering Papers

Results for “Creep modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Design by analysis rules for ASME Section III, Division 5, Class B components

The current rules for elevated temperature ASME Section III, Division 5, Class B components, other than piping, were basically adapted from the Design-by-Rule approach for Section VIII, Division 1 vessels. A goal of the proposed new rules is to explicitly account for the design life and cyclic service while recognizing the less-rigorous requirements commensurate with lesser safety consideration. Further goals are to maximize the use of modern computational technology, for example finite element analysis in conjunction with reference stress concepts, and to avoid the need for stress classification. This paper summarizes the work done to address these issues. A summary of the design-by-analysis for primary loads, strain limit evaluation, and creep-fatigue damage assessment is presented. The proposed design-by-analysis creep-fatigue damage calculation approach uses a new elastic follow-up-based Isochronous Stress Strain Curve stress relaxation procedure. A set of sample problems are selected to validate the proposed design-by-analysis rules for creep-fatigue damage assessment. The proposed Class B rules are evaluated against the Class A elastic design rules and the experimental data obtained from a family of a simplified model test-based key-feature test results. The proposed Class B creep-fatigue damage assessment methodology yields conservative design cycles estimates compared to the experimental results.

36 MATERIALS SCIENCE↗

NASA NTRS1981

Creep crack-growth: A new path-independent T sub o and computational studies

Two path independent integral parameters which show some degree of promise as fracture criteria are the C* and delta T sub c integrals. The mathematical aspects of these parameters are reviewed. This is accomplished by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects and large deformations. Two principal conclusions are that delta T sub c is a valid crack tip parameter during nonsteady as well as steady state creep and that delta T sub c has an energy rate interpretation whereas C* does not. An efficient, small displacement, infinitestimal strain, displacement based finite element model is developed for general elastic/plastic material behavior. For the numerical studies, this model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations.

Stonesifer, R. B.↗

NASA NTRS1984

A nonlinear viscoelastic characterization of graphite/epoxy composites

A nonlinear viscoelastic model based on the procedure of Findley et al., (1948) has been used to accurately represent the creep of several graphite/epoxy composites. Applying this approach to unidirectional 0, 90 and 10 deg off-axis tensile specimens, the viscoelastic response of a lamina can be characterized. The resulting lamina model has been useful for representing the behavior of a lamina in a numerical procedure to predict creep and delayed failures of general laminates. Also, independently, the Findley procedure has been used to characterize the nonlinear viscoelastic behavior of general laminates.

Dillard, D. A.↗

DOE OSTI2025

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis↗

NASA NTRS1991

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components, part 2

The technical effort and computer code enhancements performed during the sixth year of the Probabilistic Structural Analysis Methods program are summarized. Various capabilities are described to probabilistically combine structural response and structural resistance to compute component reliability. A library of structural resistance models is implemented in the Numerical Evaluations of Stochastic Structures Under Stress (NESSUS) code that included fatigue, fracture, creep, multi-factor interaction, and other important effects. In addition, a user interface was developed for user-defined resistance models. An accurate and efficient reliability method was developed and was successfully implemented in the NESSUS code to compute component reliability based on user-selected response and resistance models. A risk module was developed to compute component risk with respect to cost, performance, or user-defined criteria. The new component risk assessment capabilities were validated and demonstrated using several examples. Various supporting methodologies were also developed in support of component risk assessment.

Source record↗

NASA NTRS2021

Micromechanics-based Modeling of Laminated SiC/SiC Ceramic Matrix Composites

The behavior and response of ceramic matrix composites (CMCs), in particular silicon carbide fiber reinforced silicon carbide matrix (SiC/SiC), is affected by many factors such as variation of fiber volume fraction, residual stresses resulting from processing of the composites at high temperature, random microstructures, and the presence of matrix flaws (e.g., voids, pores, cracks etc.) as well as general material nonlinearity and heterogeneity that occurs randomly in a composite. Residual stresses arising from the phase change of constituents are evaluated in this paper and it is shown that they do influence composite strength and need to be properly accounted for. Additionally, the microstructures (location of fiber centers, coating thickness etc.) of advanced CMCs are usually disordered (or random) and fiber diameter and strength typically have a distribution. They rarely resemble the ordered fiber packing (square, rectangular, or hexagonal) that is generally assumed in micromechanics-based models with periodic boundary conditions for computational expediency. These issues raise the question of how should one model such systems effectively? Can an ordered hexagonal packed repeating unit cell (RUC) accurately represent the random microstructure behavior? How many fibers need to be included to enable accurate representation? Clearly, the number of fibers within an RUC must be limited to insure a balance between accuracy and efficiency. NASA’s in-house micromechanics-based code MAC/GMC provides a framework to analyze such RUCs for the overall composite behavior and the FEAMAC computer code provides linkage of MAC/GMC to the commercial FEA code, ABAQUS. The appropriate level of discretization of the RUC as well as the analysis method employed, i.e., Generalized Method of Cells (GMC) or High Fidelity Generalized Method of Cells (HFGMC), is investigated in this paper in the context of a unidirectional as well as a cross-ply laminated CMC. Results including effective composite properties, proportional limit stress (an important design parameter) and fatigue are shown utilizing both GMC as well as HFGMC. Finally, a few multiscale analyses are performed on smooth bar test coupons as well as test coupons with features such as open-hole and double notches using FEAMAC. Best practices and guidance are provided to take these phenomena into account and keep a proper balance between fidelity (accuracy) and efficiency. Following these guidelines can account for important physics of the problem and provide significant advantages when performing large multiscale composite structural analyses.

Ceramic Matrix Composites↗

NASA NTRS1988

The viscoplasticity theory based on overstress applied to the modeling of a nickel base superalloy at 815 C

Short term strain rate change, creep and relaxation tests were performed in an MTS computer controlled servohydraulic testing machine. Aging and recovery were found to be insignificant for test times not exceeding 30 hrs. The material functions and constants of the theory were identified from results of strain rate change tests. Numerical integration of the theory for relaxation and creep tests showed good predictive capabilities of the viscoplasticity theory based on overstress.

Krempl, E.↗

DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

DOE OSTI2025

Computational materials reliability assessment of hydrogen fueled gas turbine power generation engines

The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.

36 MATERIALS SCIENCE↗

NASA NTRS2024

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles With Interrupted Relaxation Periods of TIMETAL 21S at 650 oC

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

Viscoplaticity↗

NASA NTRS2026

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

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

Viscoplaticity↗

DOE OSTI2021

Draft Rules for Alloy 617 Creep-Fatigue Design Using an EPP+SMT Approach

The report describes an improved new creep-fatigue design method for structural components in high temperature nuclear service. The new method uses an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that greatly simplifies the design process by avoiding the separate evaluation of creep and fatigue damages and eliminating the requirement of stress classification. The report includes draft design rules, presented in a format compatible with an ASME nuclear Code Case, as well as a commentary on the rules and the validation of the rules using pressurized SMT (p-SMT) test data. Areas of potential improvement have been identified for further evaluation. The report also develops EPP+SMT design charts for Alloy 617 at temperatures between 800°C and 950°C by extrapolating high strain range, short hold test data to low strain range and long hold times, typical for structural components in high temperature nuclear service. The design charts are currently limited to 10,000 cycles to avoid excessive extrapolation outside the test database. More robust alternate extrapolation procedures will be considered in future work. The report includes several verification problems comparing the new EPP+SMT design method with current ASME creep-fatigue design methods. The comparison demonstrates that the EPP+SMT method significantly reduces the over-conservatism in current methods. Finally, the report includes a fully-documented sample problem with multiple load cases detailing the application of proposed EPP+SMT design rules.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Quantifying trapped magnetic vortex losses in niobium resonators at mK temperatures

Trapped magnetic vortices in niobium introduce microwave losses that degrade the performance of superconducting resonators. While such losses have been extensively studied above 1 K, we report here their direct quantification in the millikelvin and low-photon regime relevant to quantum devices. Using a high-quality factor 3D niobium cavity cooled through its superconducting transition in controlled magnetic fields, we isolate vortex-induced losses and find the resistive component of the sensitivity to trapped flux S to be approximately 2 n Ω/mG at 10 mK and 6 GHz. The decay rate is initially dominated by two-level system (TLS) losses from the native niobium pentoxide, with vortex-induced degradation of T 1 occurring above B trap ∼ 50 mG. In the absence of the oxide, even 10 mG of trapped flux limits performance, Q 0 ∼ 10 10 , or T 1 ∼ 350 ms, underscoring the need for stringent magnetic shielding. The resistive sensitivity, S, decreases with temperature and remains largely field-independent, whereas the reactive component, S′, exhibits a maximum near 0.8 K. These behaviors are well modeled within the Coffey–Clem framework in the zero-creep limit, under the assumption that vortex pinning is enhanced by thermally activated processes. Our results suggest that niobium-based transmon qubits can tolerate vortex-induced dissipation at trapped field levels up to several hundred mG, but achieving long coherence times still requires careful magnetic shielding to suppress lower-field losses from other mechanisms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Advancements in NEAMS Tool Capabilities for Multiphysics Simulation of Fast Reactor Core Bowing and Identification of Validation Test Data

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves Multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Last year, thermal fluids and mechanics were coupled on a multi-assembly benchmark problem based on ABR-1000 design. This year, the reactor physics code Griffin was assessed for readiness of core bowing calculations. Preliminary integration of Griffin’s ring-heterogeneous model with thermal fluids and thermal mechanics solvers was performed. Specifically, thermal-mechanics and reactor physics were coupled for single- and multi-assembly problems, and reactor physics and subchannel methods were coupled for a single assembly model. Finally, the workflow of all three physics was preliminarily demonstrated on a single assembly model. Caveats and future development needed have been identified. To supplement the multiphysics demonstration, verification and assessment efforts of thermos-mechanical capabilities for modeling thermo-mechanical core bowing behavior were continued by analyzing IAEA Verification Problem 5 which includes radiation swelling and creep. Additionally, a small core reactor physics benchmark defined by Japan Atomic Energy Agency (JAEA) was performed to assess neutronics models for estimating reactivity feedback. Finally, Fast Flux Test Facility (FFTF) validation test data for core bowing phenomena has been identified and summarized, with a recommended path forward for validation once this capability is mature.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Quantifying Trapped Magnetic Vortex Losses in Niobium Resonators at mK Temperatures

Trapped magnetic vortices in niobium introduce microwave losses that degrade the performance of superconducting resonators. While such losses have been extensively studied above 1 K, we report here their direct quantification in the millikelvin and low-photon regime relevant to quantum devices. Using a high-quality factor 3D niobium cavity cooled through its superconducting transition in controlled magnetic fields, we isolate vortex-induced losses and find the resistive component of the sensitivity to trapped flux S to be approximately 2 nOhm/mG at 10 mK and 6 GHz. The decay rate is initially dominated by two-level system (TLS) losses from the native niobium pentoxide, with vortex-induced degradation of T1 occurring above Btrap~50 mG. In the absence of the oxide, even 10 mG of trapped flux limits performance, Q0 ~ 1010, or T1 ~350 ms, underscoring the need for stringent magnetic shielding. The resistive sensitivity, S, decreases with temperature and remains largely field-independent, whereas the reactive component, S , exhibits a maximum near 0.8 K. These behaviors are well modeled within the Coffey Clem framework in the zero-creep limit, under the assumption that vortex pinning is enhanced by thermally activated processes. Our results suggest that niobium-based transmon qubits can tolerate vortex-induced dissipation at trapped field levels up to several hundred mG, but achieving long coherence times still requires careful magnetic shielding to suppress lower-field losses from other mechanisms

Bafia, Daniel P. [Fermilab]↗

NASA NTRS1970

Generation of long time creep data on refractory alloys at elevated temperatures

Creep tests were conducted on two tantalum alloys (ASTAR 811C and T-111 alloy), on a molybdenum alloy (TZM), and on CVD tungsten. The T-111 alloy 1% creep life data have been subjected to Manson's station function analysis, and the progress on this analysis is described. In another test program, the behavior of T-111 alloy with continuously varying temperatures and stresses has been studied. The results indicated that the previously described analysis predicts the observed creep behavior with reasonable accuracy. In addition to the T-111 test program, conventional 1% creep life data have been obtained for ASTAR 811C alloy. Previously observed effects of heat treatment on the creep strength of this material have been discussed and a model involving carbide strengthening primarily at the grain boundaries, rather than in a classical dispersion hardening mechanism, has been proposed to explain the observed results.

Sheffler, K. D.↗

NASA NTRS1975

Supercritical wing sections 2, volume 108

A mathematical theory for the design and analysis of supercritical wing sections was previously presented. Examples and computer programs showing how this method works were included. The work on transonics is presented in a more definitive form. For design, a better model of the trailing edge is introduced which should eliminate a loss of fifteen or twenty percent in lift experienced with previous heavily aft loaded models, which is attributed to boundary layer separation. How drag creep can be reduced at off-design conditions is indicated. A rotated finite difference scheme is presented that enables the application of Murman's method of analysis in more or less arbitrary curvilinear coordinate systems. This allows the use of supersonic as well as subsonic free stream Mach numbers and to capture shock waves as far back on an airfoil as desired. Moreover, it leads to an effective three dimensional program for the computation of transonic flow past an oblique wing. In the case of two dimensional flow, the method is extended to take into account the displacement thickness computed by a semi-empirical turbulent boundary layer correction.

Bauer, F.↗

NASA NTRS1984

Low cycle fatigue of MAR-M 200 single crystals at 760 and 870 deg C

Fully reversed low cycle fatigue tests were conducted on single crystals of the nickel-base superalloys Mar-M 200 at 760 C and 870 C. At 760 C, planar slip (octahedral) lead to orientation-dependent strain hardening and cyclic lives. Multiple slip crystals strain hardened the most, resulting in relatively high stress ranges and low lives. Single slip crystals strain hardened the least, resulting in relatively low stress ranges and higher lives. A preferential crack initiation site which was related to slip plane geometry was observed in single slip orientated crystals. At 870 C, the trends were quite different, and the slip character was much more homogeneous. As the tensile axis orientation deviated from 001 , the stress ranges increased and the cyclic lives decreased. Two possible mechanisms were proposed to explain the behavior: one is based on Takeuchi and Kuramoto's cube cross-slip model, and the other is based on orientation-dependent creep rates.

Milligan, W. W.↗