Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “failure simulation”

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

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

At least 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Analyzing the Efficacy of Flexible Execution, Replanning, and Plan Optimization for a Planetary Lander

Plan execution in unknown environments poses a number of challenges: uncertainty in domain modeling, stochasticity at execution time, and the presence of exogenous events. These challenges motivate an integrated approach to planning and execution that is able to respond intelligently to variation. We examine this problem in the context of the Europa Lander mission concept, and propose a planning and execution framework that responds to feedback and task failure using two techniques: flexible execution and replanning with plan optimization. We develop a theoretical framework to predict the value of each of these techniques, and we compare these predictions to empirical results generated in simulation. We demonstrate that an integrated approach to planning and execution that is grounded in flexible execution, replanning, and utility maximization will be an enabling technology for future tightly-constrained planetary surface missions.

Basich, Connor↗
DOE OSTI2024

Fuel Performance Evaluation of THOR-C Experiments

The Temperature Heatsink Overpower Response Commissioning (THOR-C) and THOR-Metal (THOR-M) experiments will be performed as part of an ongoing project for testing sodium fast reactor fuels with the Japan Atomic Energy Agency (JAEA). The THOR-C experiments consist of fresh metallic fuel pins and have been analyzed using the ABAQUS, Ansys codes and the BISON fuel performance code. THOR-M-Loss of Flow-1 (THOR-M-LOF-1) is designed to test an EBR-II irradiated fuel pin under LOF conditions. Simulation of the THOR-MLOF-1 experiment required first simulating the base irradiation of the fuel pin in EBR-II. MFUEL module of SAS4A/SASSYS-1 [1] is a physics-based metallic fuel performance model applicable to the normal operation, transient scenarios and fuel failure modeling including scenarios with bulk fuel melting. The model has been validated using EBR-II normal operation, separate effect transient tests as well as TREAT M-Series transient tests [2]. In this study, MFUEL models has been utilized together with a new capsule heat transfer model developed in this project. The new heat transfer model was necessary due to (1) significant amount of heat losses that required 2D heat transfer, (2) the presence of a titanium heat sink, rejecting a significant amount of heat, and (3) stagnant coolant conditions, which are inconsistent with SAS4A/SASSYS-1 (SAS) heat transfer model. Updates to SAS4A/SASSYS-1 and MFUEL has been described below, followed by a preliminary validation effort using the results from THOR-C-2 fresh fuel capsule experiment. A previous study for THOR-C-2 analysis using BISON code is also utilized in this study to model this test [3]. [1] D. O’Grady, A. J. Brunett, L. Ibarra, A. Karahan, T. Kim, T. S. Sumner, R. Thomas, T. H. Fanning, “The SAS4A/SASSYS-2 Version 5.7 Safety Analysis Code System,” Argonne National Laboratory,ANL/NSE-SAS/5.7, (2023). [2] A. Karahan, T. Kim, T. Fanning, D. O’Grady, “Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1,” Argonne National Laboratory, ANL/NSE-23/11, (2023). [3] M. Mihelish, A. Zabriskie, K. Paaren, P. Medvedev, C. Jensen, “Fuel Performance Predictions for the TREAT THOR-C Experiments,” Idaho National Laboratory, INL/RPT-23-73397, Revision 0, (2023)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Resiliency in numerical algorithm design for extreme scale simulations

Here this work is based on the seminar titled ‘Resiliency in Numerical Algorithm Design for Extreme Scale Simulations’ held March 1–6, 2020, at Schloss Dagstuhl, that was attended by all the authors. Advanced supercomputing is characterized by very high computation speeds at the cost of involving an enormous amount of resources and costs. A typical large-scale computation running for 48 h on a system consuming 20 MW, as predicted for exascale systems, would consume a million kWh, corresponding to about 100k Euro in energy cost for executing 10 23 floating-point operations. It is clearly unacceptable to lose the whole computation if any of the several million parallel processes fails during the execution. Moreover, if a single operation suffers from a bit-flip error, should the whole computation be declared invalid? What about the notion of reproducibility itself: should this core paradigm of science be revised and refined for results that are obtained by large-scale simulation? Naive versions of conventional resilience techniques will not scale to the exascale regime: with a main memory footprint of tens of Petabytes, synchronously writing checkpoint data all the way to background storage at frequent intervals will create intolerable overheads in runtime and energy consumption. Forecasts show that the mean time between failures could be lower than the time to recover from such a checkpoint, so that large calculations at scale might not make any progress if robust alternatives are not investigated. More advanced resilience techniques must be devised. The key may lie in exploiting both advanced system features as well as specific application knowledge. Research will face two essential questions: (1) what are the reliability requirements for a particular computation and (2) how do we best design the algorithms and software to meet these requirements? While the analysis of use cases can help understand the particular reliability requirements, the construction of remedies is currently wide open. One avenue would be to refine and improve on system- or application-level checkpointing and rollback strategies in the case an error is detected. Developers might use fault notification interfaces and flexible runtime systems to respond to node failures in an application-dependent fashion. Novel numerical algorithms or more stochastic computational approaches may be required to meet accuracy requirements in the face of undetectable soft errors. These ideas constituted an essential topic of the seminar. The goal of this Dagstuhl Seminar was to bring together a diverse group of scientists with expertise in exascale computing to discuss novel ways to make applications resilient against detected and undetected faults. In particular, participants explored the role that algorithms and applications play in the holistic approach needed to tackle this challenge. This article gathers a broad range of perspectives on the role of algorithms, applications and systems in achieving resilience for extreme scale simulations. The ultimate goal is to spark novel ideas and encourage the development of concrete solutions for achieving such resilience holistically.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Mitigation of distortion of Al/steel part under simulated paint baking condition: Experiment and numerical model studies

Multi-material joining of lightweight structures is essential to reduce vehicle weight for more energy savings and less greenhouse gas emission. However, mismatch of thermal expansion coefficient for dissimilar materials during the paint baking process can induce part distortion and joint failure for adhesive bonding. Here, in the present work, a thermomechanical model based on contact mechanics and large deformation theory was developed for dissimilar high-strength Al alloy and steel components to study the distortion mechanism and influential factors of the residual gap. The established model was used to optimize joint conditions, such as pitch distance and part geometry. When a weld pitch is shorter than 100 mm, the maximum gap between Al and steel part can be greatly reduced to 0.1 mm, and the local stress and plastic strain around the joint during the oven heating and cooling cycle are also substantially reduced compared with the long pitch case (900 mm). The numerical modeling results revealed that a comparable bending stiffness ratio between the steel and Al cross sections is critical to the minimization of gap and distortion under paint baking condition. Digital image correlation technique was used to measure the overall part distortion and local strain distribution that were used to validate the model prediction. Weld bonding (adhesive bonding with friction bit joining) process was successfully employed to join Al to steel component without gap opening in adhesive after the paint baking and cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stress Controls Rupture Extent and Maximum Magnitude of Induced Earthquakes

Abstract Seismic hazard forecasts of induced seismicity often require estimates of the maximum possible magnitude ( M max ). Empirical models suggest that maximum magnitudes, or expected number of earthquakes, are related to the volume of injected fluid. We perform a suite of 3D physics‐based earthquake simulations with rate‐ and state‐dependent friction, systematically varying the area of the pressurized region and the amplitude of the initial homogeneous or heterogeneous shear stress. Using the resulting catalog we explore the conditions that result in pressure‐controlled versus runaway ruptures that extend outside the pressurized zone. We find that proposed empirical scaling laws correctly predict M max when shear stresses are further from failure (≤90% of maximum shear stress) and for high amplitude stress fields. Runaway ruptures are observed for higher initial shear stresses and smoother stress fields. In these cases, runaway ruptures occur early after the onset of injection and rarely preceded by foreshock activity.

58 GEOSCIENCES↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2020

Using optical fibers to examine thermal mixing of liquid sodium in a pool-type geometry

Sodium fast reactors (SFR) are poised to be a leading candidate for the next generation of commercial nuclear reactor deployment. Companies are commissioning designs for SFR increasing the need for experimental and computational analysis to improve the economics and safety of these reactors. Analysis of SFR behavior during reactor transients directly informs the understanding of reactor safety. Particularly, investigating the thermal transients associated with the reactor loss of flow transient allows for informed design decisions to be made regarding reactor safety. During loss of flow transients the flow rate of the coolant is significantly reduced, but not necessarily stopped. The reduction in flow rate results in coolant temperatures exiting the reactor core that can be drastically different than the coolant in upper pool. When these different temperature fluids interact, thermally stratified layers can form causing cyclic thermal fatigue on the reactor vessel potentially leading to pre-material failure. Past experimental studies have lacked high resolution temperature measurements. This research focuses on demonstrating the reliability of novel techniques for acquiring high resolution temperature measurements in a scaled liquid sodium facility. Tests are conducted that simulate postulated loss of flow transients and provide unprecedented spatial temperature distribution measurements.

42 ENGINEERING↗
DOE OSTI2021

Modelling yield cap evolution in sandstone based on brittle creep experiments

The Zenifim Formation is being considered as a potential disposal formation for a deep borehole nuclear repository concept in Israel. Site selection and repository construction are intended to ensure that waste is separated from circulating groundwater, but long-term deformation of the wellbore could potentially create fluid flow pathways. To understand how time-dependent rock strength could affect wellbore stability, here we conducted creep tests under low to moderate confining pressures on retrieved core from the Zenifim formation. During creep, samples strain slowly as gradual damage accumulation progressively weakens the samples. Failure eventually occurred through the near-instantaneous formation of a shear fracture. Experimental results were used to calibrate a continuum damage poro-elastic model for sandstones. The calibrated damage-poro-elastic model successfully simulates different types of loading experiments including quasi-static and creep. The state of strain in experiments is close to yield during loading as the yield cap continuously evolves with damage accumulation. For creep tests, most damage occurs during triaxial loading. Minor damage accumulation occurs under constant load until the final stage of creep, where damage accelerates and promotes unstable fracturing.

58 GEOSCIENCES↗
DOE OSTI2020

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Gyro and accelerometer failure detection and identification in redundant sensor systems.

An algorithm for failure detection and identification for redundant non-orthogonal arrays of single degree of freedom gyros and accelerometers is described which is optimum in the sense that detection occurs as soon as it is no longer possible to account for the instrument outputs as the outputs of good instruments operating within their noise tolerances, and identification occurs as soon as it is true that only a particular instrument failure could account for the actual instrument outputs within the noise tolerance of good instruments. An estimation algorithm is described which minimizes the maximum possible error magnitude for the given set of instrument outputs. Monte Carlo simulation results are presented for the application of the algorithms to an inertial reference unit consisting of six gyros and six accelerometers in two alternate configurations.

Potter, J. E.↗
NASA NTRS1973

Fused slurry silicide coatings for columbium alloys reentry heat shields. Volume 1: Evaluation analysis

The R-512E (Si-20Cr-20Fe) fused slurry silicide coating process was optimized to coat full size (20in x 20in) single face rib and corrugation stiffened panels fabricated from FS-85 columbium alloy for 100 mission space shuttle heat shield applications. Structural life under simulated space shuttle lift-off stresses and reentry conditions demonstrated reuse capability well beyond 100 flights for R-512E coated FS-85 columbium heat shield panels. Demonstrated coating damage tolerance showed no immediate structural failure on exposure. The FS-85 columbium alloy was selected from five candidate alloys (Cb-752, C-129Y, WC-3015, B-66 and FS-85) based on the evaluation tests which have designed to determine: (1) change in material properties due to coating and reuse; (2) alloy tolerance to coating damage; (3) coating emittance characteristics under reuse conditions; and (4) new coating chemistries for improved coating life.

Fitzgerald, B.↗
NASA NTRS1989

A scheme for fault tolerance in earth sensors

A system is presented that uses dual-redundant earth sensors to measure pitch and roll errors of a three-axis stabilized spacecraft, with provision for (1) autonomously detecting and identifying a faulty earth sensor, and (2) automatically selecting the outputs of the fault-free sensor for closed-loop attitude control, before failures cause major problems. A brief description is given of the system, and various failure modes of earth sensors and their effects are discussed. Novel techniques and algorithms for automatic fault detection, identification, and reconfiguration (FDIR) of dual-redundant earth sensors are developed. The algorithms are validated through computer simulations, and the results are presented. The proposed scheme can easily be implemented without much penalty on hardware, power consumption, and processing time.

Murugesan, S.↗
NASA NTRS1991

High-frequency data observations from space shuttle main engine low pressure fuel turbopump discharge duct flex joint tripod failure investigation

Observations made by Marshall Space Flight Center (MSFC) engineers during their participation in the Space Shuttle Main Engine (SSME) low pressure fuel turbopump discharge duct flex joint tripod failure investigation are summarized. New signal processing techniques used by the Component Assessment Branch and the Induced Environments Branch during the failure investigation are described in detail. Moreover, nonlinear correlations between frequently encountered anomalous frequencies found in SSME dynamic data are discussed. A recommendation is made to continue low pressure fuel (LPF) duct testing through laboratory flow simulations and MSFC-managed technology test bed SSME testing.

Zoladz, T. F.↗
NASA NTRS1997

The Design of Model-Based Training Programs

This paper proposes a model-based training program for the skills necessary to operate advance avionics systems that incorporate advanced autopilots and fight management systems. The training model is based on a formalism, the operational procedure model, that represents the mission model, the rules, and the functions of a modem avionics system. This formalism has been defined such that it can be understood and shared by pilots, the avionics software, and design engineers. Each element of the software is defined in terms of its intent (What?), the rationale (Why?), and the resulting behavior (How?). The Advanced Computer Tutoring project at Carnegie Mellon University has developed a type of model-based, computer aided instructional technology called cognitive tutors. They summarize numerous studies showing that training times to a specified level of competence can be achieved in one third the time of conventional class room instruction. We are developing a similar model-based training program for the skills necessary to operation the avionics. The model underlying the instructional program and that simulates the effects of pilots entries and the behavior of the avionics is based on the operational procedure model. Pilots are given a series of vertical flightpath management problems. Entries that result in violations, such as failure to make a crossing restriction or violating the speed limits, result in error messages with instruction. At any time, the flightcrew can request suggestions on the appropriate set of actions. A similar and successful training program for basic skills for the FMS on the Boeing 737-300 was developed and evaluated. The results strongly support the claim that the training methodology can be adapted to the cockpit.

Polson, Peter↗
NASA NTRS2003

Hydrostatic Stress Effect on the Yield Behavior of Inconel 100

Classical metal plasticity theory assumes that hydrostatic stress has negligible effect on the yield and postyield behavior of metals. Recent reexaminations of classical theory have revealed a significant effect of hydrostatic stress on the yield behavior of various geometries. Fatigue tests and nonlinear finite element analyses (FEA) of Inconel 100 (IN100) equal-arm bend specimens and new monotonic tests and nonlinear finite element analyses of IN100 smooth tension, smooth compression, and double-edge notch tension (DENT) test specimens have revealed the effect of internal hydrostatic tensile stresses on yielding. Nonlinear FEA using the von Mises (yielding is independent of hydrostatic stress) and the Drucker-Prager (yielding is linearly dependent on hydrostatic stress) yield functions were performed. A new FEA constitutive model was developed that incorporates a pressure-dependent yield function with combined multilinear kinematic and multilinear isotropic hardening using the ABAQUS user subroutine (UMAT) utility. In all monotonic tensile test cases, the von Mises constitutive model, overestimated the load for a given displacement or strain. Considering the failure displacements or strains for the DENT specimen, the Drucker-Prager FEM s predicted loads that were approximately 3% lower than the von Mises values. For the failure loads, the Drucker Prager FEM s predicted strains that were up to 35% greater than the von Mises values. Both the Drucker-Prager model and the von Mises model performed equally-well in simulating the equal-arm bend fatigue test.

Allen, Phillip A.↗
NASA NTRS2003

Probabilistic Evaluation of Advanced Ceramic Matrix Composite Structures

The objective of this report is to summarize the deterministic and probabilistic structural evaluation results of two structures made with advanced ceramic composites (CMC): internally pressurized tube and uniformly loaded flange. The deterministic structural evaluation includes stress, displacement, and buckling analyses. It is carried out using the finite element code MHOST, developed for the 3-D inelastic analysis of structures that are made with advanced materials. The probabilistic evaluation is performed using the integrated probabilistic assessment of composite structures computer code IPACS. The affects of uncertainties in primitive variables related to the material, fabrication process, and loadings on the material property and structural response behavior are quantified. The primitive variables considered are: thermo-mechanical properties of fiber and matrix, fiber and void volume ratios, use temperature, and pressure. The probabilistic structural analysis and probabilistic strength results are used by IPACS to perform reliability and risk evaluation of the two structures. The results will show that the sensitivity information obtained for the two composite structures from the computational simulation can be used to alter the design process to meet desired service requirements. In addition to detailed probabilistic analysis of the two structures, the following were performed specifically on the CMC tube: (1) predicted the failure load and the buckling load, (2) performed coupled non-deterministic multi-disciplinary structural analysis, and (3) demonstrated that probabilistic sensitivities can be used to select a reduced set of design variables for optimization.

Abumeri, Galib H.↗