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

How extreme rainfall and failing dams unleashed the Derna flood disaster

On September 11, 2023, Storm Daniel unleashed unprecedented rainfall over the Wadi Derna watershed, triggering one of the most devastating floods in modern history, striking Derna, a coastal city in Libya. This study reconstructs the disaster using an integrated modeling approach that combines satellite imagery, hydrologic, hydraulic, and geotechnical simulations, machine learning, eyewitness accounts, and digital elevation data to assess the impact of cascading dam failures. Our findings reveal that the region’s dams, even if structurally sound, would have provided minimal protection against the extreme runoff. However, their failure unleashed a destructive surge wave, amplifying the disaster’s magnitude and devastation. Here, we show that the collapse of aging flood control infrastructures, compounded by inadequate risk assessment and emergency preparedness, dramatically escalated the disaster’s impact. Our findings underscore the urgent need for systematic dam safety evaluations, enhanced flood forecasting, and adaptive risk management strategies that address climate extremes and infrastructure vulnerabilities.

Hydrology↗
DOE OSTI2020

Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI

The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λ L ), transversal diffusivity (λ T ), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Simulation-Based Recovery Action Analysis Using the EMRALD Dynamic Risk Assessment Tool

A recovery action is defined as the action that prevents deviant conditions from producing unwanted effects. It generally indicates a kind of countermeasure performed in response to a failure of human action. The recovery actions especially play an important role in complex systems like nuclear power plants (NPPs), which consist of highly sophisticated controllers to ensure that desired performance and safety must be achieved and maintained. This is because a combination of human error and its recovery failure may be able to cause a catastrophic effect on a system. Analyzing recovery actions has been a critical part of HRA, which is a technique to evaluate human errors and provide human error probabilities (HEPs) for application in probabilistic safety assessment (PSA). If recovery actions are not adequately analyzed and applied to PSA models, the PSA results may be under-estimated or be not able to reasonably account for the failure of human actions in the context of PSA. For this reason, some regulatory documents such as ASME/ANS RA-Sb-2013 by the American Society for Mechanical Engineers and the American Nuclear Society and NUREG-1792 by U.S. Nuclear Regulatory Commission have emphasized the importance of recovery analysis within the HRA. A couple of existing HRA methods, such as the Technique for Human Error-Rate Prediction (THERP), the Cause-Based Decision Tree (CBDT), and the Korean Standard HRA (K-HRA), have respectively suggested their own approaches to the HRA recovery analysis. However, there are a couple of limitations to treating recovery actions using only the current HRA methods available. The biggest limitation is that the existing recovery analysis does not explicitly consider a variety of recovery action types and recovery sequences as they occur in actual NPPs. To handle the limitations of existing recovery analysis, this study proposes a simulation-based recovery analysis method using the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software. The EMRALD software is a dynamic simulation tool for PSA. It supports realistic and dynamic modeling of human actions as they would be performed at NPPs. It is also favorable to simultaneously model the specific moment at which an action is performed, the time it takes to perform the action, and the failure probability of that action. In this paper, a detailed methodology for modeling recovery actions in the simulation platform is proposed with a couple of examples. Then, outputs from the simulation are discussed as reviewing if this novel approach can complement the challenges of existing recovery analyses.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Development of mechanistic, microstructure-informed BISON models for fission product-induced failure mechanisms in advanced nuclear fuels

The US continues to prioritize commercial deployment of advanced reactors—particularly those utilizing U-Zr metallic and TRISO particle fuels. Fission product-induced failure mechanisms for these systems include fuel–cladding chemical interaction in metallic fuel rods, which threatens cladding integrity, and Pd penetration in TRISO particles, which degrades SiC layer properties. Empirical models can be applied within the bounds of existing irradiation databases, but their utility is limited when considering new designs or investigating fundamental material behaviors. The Nuclear Energy Advanced Modeling and Simulation Program is therefore developing mechanistic models for these behaviors, which are expected to aid in fuel design, assist in development of failure mitigation strategies, and provide support for qualification and licensing. These efforts leverage multiscale capabilities to capture the microstructural features and processes that govern these behaviors. This talk summarizes past and ongoing efforts to develop and validate these models for the BISON fuel performance code.

FCCI↗
DOE OSTI2021

Unraveling disadhesion mechanism of epoxy/CSH interface under aggressive conditions

Almost all FRP-reinforced concrete exposed to the marine environment are susceptible to durability deficiencies due to failure of the FRP/epoxy/concrete interface, especially the vulnerable epoxy/concrete interface, which is frequently attacked by water and ions. Herein, molecular dynamics simulations of epoxy/CSH interface under different environmental exposures are utilized to specifically elucidate the degradation principles of water and ions on the bonding properties and mechanical responses of the epoxy/concrete interface. The environmental sequences associated with the interfacial bonding degradation are arranged as Na{sub 2}SO{sub 4} + NaCl>NaCl>Na{sub 2}SO{sub 4} > Water>Dry. Remarkably, water molecules weaken the interaction between the epoxy and CSH by robbing interfacial H- and Ca-O bonds, thereby reducing the energy barrier to interfacial peel and shear failure. The weakening effects are exacerbated by the presence of aggressive ions, as Na{sup +} ions attracted to the CSH surface accumulate more water molecules in the interfacial region by forming hydrated ions clusters, which ultimately accelerates the epoxy/CSH interfacial degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Damage and Failure in a Statistical Crack Model

In this paper I will take a close look at a statistical crack model (SCM) as is used in engineering computer codes to simulate fracture at high strain rates. My general goal is to understand the macroscopic behavior effected by the microphysical processes incorporated into an SCM. More specifically, I will assess the importance of including local interactions between cracks into the growth laws of an SCM. My strategy will be to construct a numerical laboratory that represents a single computational cell containing a realization of a statistical distribution of cracks. The cracks will evolve by the microphysical models of the SCM, leading to quantifiable damage and failure of the computational cell. I will use the numerical data generated by randomly generated ensembles of the fracture process to establish scaling laws that will modify and simplify the implementation of the SCM into large scale engineering codes.

42 ENGINEERING↗
DOE OSTI2022

Predicting the location of shear band initiation in a metallic glass

We report atomistic simulation results which indicate that the location of shear banding in a metallic glass (MG) can be ascertained with reasonably high accuracy solely from the undeformed static structure. Correlation is observed between the location of the initiation of shear bands in a simulated MG and the initial distribution of the density of fertile sites (DFS) for stress-driven shear transformations identified a priori based on a deep learning model devised in our recent work. In addition, we demonstrated that one can judge whether a glass is brittle or ductile solely based upon its initial DFS distribution. These validate that shear bands in MG arise from nonlinear instabilities and that the as-quenched glass structure contains inhomogeneities that influence these instabilities. In this study, we also reveal an important subtlety regarding the nondeterministic nature of athermal quasistatic shear simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mitigating Cascading Outages in Severe Weather Using Simulation-Based Optimization

Severe weather events can trigger cascading power outages and lead to significant losses. In this work, we investigate cascading outage mitigation under severe weather conditions. Given day-ahead weather forecasts and component failure models, we aim to identify a set of power lines that can be hardened to minimize the expected impact of potential cascading outages. Since the expected load shedding cannot be expressed as an explicit function of line hardening decisions and system states, we developed a cascading outage simulator to estimate the expected value of load shedding under various initial weather-related disruption scenarios generated using a weather forecast. To avoid massive enumeration of all possible combinations of line hardening decisions and reduce the simulation efforts, we employed an efficient simulation-based optimization approach that quickly identifies the (near) optimal line hardening decisions in the presence of both large simulation noises due to the highly variable initial disturbances and system states, and significant randomness in the subsequent cascades. Furthermore, the algorithm is also able to utilize parallel computing to dramatically reduce computation time to support decision making in preparation for severe weather conditions. We performed a case study on the Northeast Power Coordinating Council (NPCC) 140-bus system model to demonstrate that our approach can significantly improve power grid resilience to adverse weather events.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

LaueMatching: an approach for rapid and robust indexing of Laue diffraction patterns

Traditional Laue diffraction pattern indexing often struggles with noisy data, weak signals, peak overlap and missing reflections, particularly from complex or deformed microstructures. Here, we introduce LaueMatching, a high-throughput indexing algorithm designed to overcome these limitations. LaueMatching utilizes a fundamentally different approach based on direct pattern correlation: experimentally pre-processed images are compared against a comprehensive pre-computed library of simulated diffraction patterns corresponding to a dense grid of possible orientations. This approach bypasses the need for explicit peak identification and fitting, steps that are often a failure point for traditional methods. The algorithm rapidly and robustly indexes multiple crystallographic orientations and crystal systems simultaneously, even from challenging patterns. LaueMatching's effectiveness and accuracy have been rigorously tested and validated on diverse experimental (Ni, Al, EuAl 2 O 4 ) and simulated diffraction patterns, demonstrating high-fidelity orientation refinement. Code to implement this approach on both CPU and GPU resources can be downloaded from https://github.com/AdvancedPhotonSource/LaueMatching.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rear-Side Irradiance Simulation of Field PV Modules

Assessing the durability of photovoltaic (PV) module backsheet is critical to increasing module lifetime. Historically, laboratory-based accelerated testing has insufficient in predicting large-scale failures of commercial polymeric materials. Additionally, there is growing concern that standard condition tests do not reflect non-uniformities in field exposure, and that certain modules experience more severe degradation due to their location. Anisotropy in field exposures is installation-dependent and reflects different levels of exposure to irradiance due to mounting geometry, ground surface albedo, and climatic zone. "Bificial_Radiance" [1], simulates the amount of reflected irradiance on the Polyethylene Naphthalate (PEN) backsheet of a PV array located in Maryland. In our work, site specific weather data are gathered for the entire length of exposure for the modules. The total full-spectrum dose for the incident irradiance on the backsheet was then determined. The simulation results are integrated with historical field survey data to better understand real orld outdoor degradation.

backsheet↗
DOE OSTI2023

Chapter 6: Abuse Response of Batteries Subjected to Mechanical Impact

Electrochemical and thermal models to simulate nominal performance and abuse response of lithium-ion cells and batteries have been reported widely in the literature. Studies on mechanical failure of cell components and how such events interact with the electrochemical and thermal response are relatively less common. This chapter outlines a framework developed under the Computer Aided Engineering for Batteries program to couple failure modes resulting from external mechanical loading to the onset and propagation of electrochemical and thermal events that follow. Starting with a scalable approach to implement failure criteria based on thermal, mechanical, and electrochemical thresholds, we highlight the practical importance of these models using case studies at the cell and module level. The chapter also highlights a few gaps in our understanding of the comprehensive response of batteries subjected to mechanical crash events, the stochastic nature of some of these failure events, and our approach to build safety maps that help improve robustness of battery design by capturing the sensitivity of some key design parameters to heat generation rates under different mitigation strategies.

abuse simulations↗
DOE OSTI2020

Enhancement of risk informed validation framework for external hazard scenario

In recent years, the U.S. Nuclear Regulatory Commission (USNRC) and the International Atomic Energy Agency (IAEA) have developed methodologies to assess the vulnerabilities of nuclear plants against site specific extreme hazards. In many cases, advanced simulation tools are being considered to simulate multi-physics, multi-scale phenomena and to evaluate vulnerability of nuclear facilities. The credibility of advanced simulation tools is assessed based on a formal verification, validation, and uncertainty quantification procedure. One of the key limitations in validation is the lack of relevant experimental data at system-level. This limitation leads to a decrease in the confidence of system-level risk predictions. Therefore, a robust validation framework is needed to formalize the confidence in predictive capability of advanced simulation results. Additionally, this study enhances the existing risk informed validation methodology, originally proposed by Kwag et al. [1] and Bodda et al. [2], by developing additional attributes and a new set of validation indicies for a complete and wider applicability of the framework. In this manuscript, the methodology to identify the critical path that leads to the system-level failure is illustrated. The overall validation is checked for completeness and consistency by comparing the critical path for both the system-level simulation and experimental models. The applicability of the code for an intended application is represented in terms of various maturity levels and helps in the process of decision making.

42 ENGINEERING↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rear-Side Irradiance Simulation of Field PV Modules: Preprint

Assessing the durability of photovoltaic (PV) module backsheet is critical to increasing module lifetime. Historically, laboratory-based accelerated testing has insufficient in predicting large-scale failures of commercial polymeric materials. Additionally, there is growing concern that standard condition tests do not reflect non-uniformities in field exposure, and that certain modules experience more severe degradation due to their location. Anisotropy in field exposures is installation dependent and reflects different levels of exposure to irradiance due to mounting geometry, ground surface albedo, and climatic zone. "Bificial Radiance" [1], simulates the amount of reflected irradiance on the Polyethylene Naphthalate (PEN) backsheet of a PV array located in Maryland. In our work, site specific weather data are gathered for the entire length of exposure for the modules. The total full-spectrum dose for the incident irradiance on the backsheet was then determined. The simulation results are integrated with historical field survey data to better understand realworld outdoor degradation.

backsheet↗
DOE OSTI2021

RAVEN User Manual

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The initial development was aimed to provide dynamic risk analysis capabilities to the Thermo-Hydraulic code RELAP-7, currently under development at the Idaho National Laboratory (INL). Although the initial goal has been fully accomplished, RAVEN is now a multi-purpose probabilistic and uncertainty quantification platform, capable to agnostically communicate with any system code. This agnosticism includes providing Application Programming Interfaces (APIs). These APIs are used to allow RAVEN to interact with any code as long as all the parameters that need to be perturbed are accessible by inputs files or via python interfaces. RAVEN is capable of investigating the system response, and investigating the input space using Monte Carlo, Grid, or Latin Hyper Cube sampling schemes, but its strength is focused to- ward system feature discovery, such as limit surfaces, separating regions of the input space leading to system failure, using dynamic supervised learning techniques. The development of RAVEN has started in 2012, when, within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the need to provide a modern risk evaluation framework became stronger. RAVEN principal assignment is to provide the necessary software and algorithms in order to employ the concept developed by the Risk Informed Safety Margin Characterization (RISMC) program. RISMC is one of the pathways defined within the Light Water Reactor Sustainability (LWRS) program. In the RISMC approach, the goal is not just the individuation of the frequency of an event potentially leading to a system failure, but the closeness (or not) to key safety-related events. Hence, the approach is interested in identifying and increasing the safety margins related to those events. A safety margin is a numerical value quantifying the probability that a safety metric (e.g. for an important process such as peak pressure in a pipe) is exceeded under certain conditions. The initial development of RAVEN has been focused on providing dynamic risk assessment capability to RELAP-7, currently under development at the INL and, likely, future replacement of the RELAP5-3D code. Most the capabilities that have been implemented having RELAP-7 as principal focus are easily deployable for other system codes. For this reason, several side activates are currently ongoing for coupling RAVEN with soft- ware such as RELAP5-3D, etc. The aim of this document is the explanation of the input requirements, focalizing on the input structure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Constitutive model development of aluminum alloy 1100 for elevated temperature forming process

Commercially pure aluminum alloy, AA1100, presents good electrical and thermal conductivity, high formability, and low cost. Those favorable characteristics have the potential to enable bipolar plates with improved economics and enhanced performance compared to current stainless steel bipolar plates for proton exchange membrane fuel cells. An accurate constitutive model is essential to develop and optimize processing parameters and effectively control the forming process. Here, the objective of this work is to develop a constitutive model of AA1100 that is able to simulate stress-strain relation, formed geometry, and predict the onset of fracture strain to avoid forming failure. Initially, a set of tensile tests at temperature between 300 and 500°C and strain rate between 0.005 and 1.0/s were conducted to examine the deformation behavior. Then, a set of damage-based unified visco-plastic constitutive equations is proposed and calibrated based on the results of stress-strain data. A genetic algorithm optimization method is applied to search for best fitting material constants in constitutive equations. The proposed model shows good predictability of both the stress-strain relation and fracture strain at low strain rate and high temperature conditions. The accuracy of proposed model is also evaluated statistically. A comparison of the proposed model with three popular models (Arrhenius-type mode, Johnson-Cook model and Zerilli-Armstrong model) was made. The proposed model shows the best experimental agreement with correlation coefficient of 0.96 in contrast to 0.25, 0.38 and 0.75 for the popular models, respectively. The proposed model can help to optimize the elevated temperature forming process and guide die design to enable optimal geometric features in the formed components.

08 HYDROGEN↗
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↗