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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dark Lock-in Thermography Identifies Solder Bond Failure as the Root Cause of Series Resistance Increase in Fielded Solar Modules

Correlating the electrical performance of photovoltaic modules to the spatially resolved photoluminescence, electroluminescence, and dark lock-in thermography (DLIT) images is an important long-term goal for developing solar cell technology. These images offer highly sophisticated and detailed information about the spatial distribution and (if imaged at successive time intervals) temporal degradation of local series and shunt resistances. There have been extensive studies to correlate these imaging techniques to local characteristics at the cell level. However, it has been difficult to extract and quantify module-level information from the techniques. In this study, we interpret module current-voltage (I-V) measurements along with corresponding DLIT images within a module-scale simulation framework, demonstrating that series resistance degradation of the module I-V characteristics, in this case, can be attributed to solder bond failures. Our simulations highlight how current crowding associated with a failed solder joint (or a section of the solder pad) translates to the characteristic point-like (asymmetric doublet) heating pattern in neighboring solder joints (or the neighboring regions of the solder pad). Overall, the correlation between series resistance and solder joint degradation could inform expedited diagnosis of field degradation of solar modules.

42 ENGINEERING↗
DOE OSTI2024

Characterizing and improving the performance of molten-salt-steam heat exchangers in concentrating solar power plants

Shell-and-tube heat exchangers (HXs) for steam generation from molten salts in concentrating solar power (CSP) plants experience thermal fatigue due to significant temperature gradients and inherent transient operation. Molten salt-steam HX design lifespans exceed actual lifespans, and, as a consequence, designers overpredict plant profitability and operators neglect appropriate prescriptions to optimize these lifetimes. Here, this study refines HX lifespan estimates with data benchmarked against thermal-fluid mechanical modeling of stress and accumulated fatigue. Reduced-order thermal models of the molten salt-steam, shell-and-tube evaporator and superheater predict transient temperature profiles along the two HXs salt-steam flow paths. The modeled evaporator and superheater temperature profiles enable assessment of cyclic stresses within the HX tubesheets, where molten-salt HX failures are most common. Evaporator and superheater performance data from a current 110 MW elec commercial CSP plant provide a basis for validating the reduced-order HX models. HX life predictions derived from stochastic failure distributions serve as inputs for simulating and optimizing existing plant operations. The impact of the updated lifespans on overall plant revenue depends on operating scenarios. This study suggests that typical ramping rates for a CSP plant with a high-temperature Rankine cycle result in an evaporator and superheater life of approximately 10 and 25 years, respectively, compared to the design target of 30 years. Reduced HX lifespans decrease operational plant revenue on average by 4.6-5.1%. Furthermore, there may be as many as four HX replacements over the 30-year lifetime of the plant; and, purchase agreement loss due to failure to meet contractual production requirements can have ramifications that include the risk of bankruptcy.

14 SOLAR ENERGY↗
DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydropower Role in a Power Grid with a High Proportion of Variable Renewables

The importance of hydropower increases as the power grid evolves with the higher variable renewable contribution. As conventional thermal power plants are retired, the importance of hydropower contribution increases to balance the variability of solar and wind generation. However, reservoir water resources are constrained by multiple constraints, and variability of water inflow to the reservoirs creates limitations to dam water releases for power grid needs. Coordination of multiple tools, including water resources modeling, ecological modeling, and technical and economic power grid modeling, informs dam water releases. Better representation of ecological constraints and other water use details in the suit of power grid models are required to identify the hydropower flexibility to provide multiple power grid services. Impacts of water planning and hydropower operation changes to a power grid are better understood by production cost model simulations with technical and economic data, including unserved energy, reserve failures, transmission congestion to serve load, and local marginal price increases. Similarly, PCM simulation of current and future projected power grid information, including LMP, and carbon emission rates, informs short-term water release plans and long-term investments for power plant upgrades for technical and ecological needs. Hydropower operations of several case studies are simulated to understand operation patterns, revenue, water release decisions in multiple time scales (month, day, minutes), hydropower contracts scheduling, and water variability impacts.

environmental↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automated Detection of Instability-Inducing Channel Geometry Transitions in Saint-Venant Simulation of Large-Scale River Networks

A new sweep-search algorithm (SSA) is developed and tested to identify the channel geometry transitions responsible for numerical convergence failure in a Saint-Venant equation (SVE) simulation of a large-scale open-channel network. Numerical instabilities are known to occur at “sharp” transitions in discrete geometry, but the identification of problem locations has been a matter of modeler’s art and a roadblock to implementing large-scale SVE simulations. The new method implements techniques from graph theory applied to a steady-state 1D shallow-water equation solver to recursively examine the numerical stability of each flowpath through the channel network. The SSA is validated with a short river reach and tested by the simulation of ten complete river systems of the Texas–Gulf Coast region by using the extreme hydrological conditions recorded during hurricane Harvey. The SSA successfully identified the problematic channel sections in all tested river systems. Subsequent modification of the problem sections allowed stable solution by an unsteady SVE numerical solver. The new SSA approach permits automated and consistent identification of problem channel geometry in large open-channel network data sets, which is necessary to effectively apply the fully dynamic Saint-Venant equations to large-scale river networks or for city-wide stormwater networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

MLEC-Sim: A Simulator for Evaluating Multi-Level Erasure Coding

We present MLEC-Sim, a sophisticated simulator for Multi-Level Erasure Coding (MLEC), developed in approximately 13 KLOC. The simulator is engineered to analyze the impact of various system configurations and erasure coding policies on system durability and network overhead. It supports a comprehensive range of parameters including disk capacity, disk I/O bandwidth, failure rates, network bandwidth, and system scale, accommodating various erasure coding approaches such as Single-Level Erasure Coding (SLEC), Multi-Level Erasure Coding (MLEC), and Local Reconstruction Codes (LRC). MLEC-Sim provides support for multiple chunk placement policies, including clustered parity and declustered parity, and encompasses a variety of repair methods like Repair-ALL, Repair-FCO, Repair-HYB, and Repair-MIN. It is capable of simulating disk failures through a variety of means, including distribution-based or trace-based mechanisms, and can handle complex multi-level (de)clustered placements and repair processes. A key feature of MLEC-Sim is its adoption of the splitting simulation method for evaluating system durabilities at extremely high levels, which are challenging to assess with traditional simulation approaches. This feature allows for a detailed evaluation of system resilience under a range of conditions, aiding in the selection of appropriate erasure coding solutions for enhancing system durability. MLEC-Sim contributes to the field of data storage and reliability by providing a tool for the detailed evaluation of the durability and efficiency of erasure coding configurations, intended for use by researchers and practitioners in the design and optimization of storage systems.

Wang, Meng↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Metamodels for Rapid Analysis of Large Sets of Building Designs for Robotic Constructability: Technology Demonstration Using the NASA 3D Printed Mars Habitat Challenge

Disruptive robotic construction technologies such as additive deposition of cementitious materials like concrete (or "3D concrete printing") require the synchronous operation of multiple pieces of equipment in the production setup. In such an environment, it is crucial to simulate the robotic motions (for toolpath clashes) and the cementitious material behavior (for toolpath failures) to ensure fail-proof constructability of the envisioned building geometry. However, toolpath clash detection requires 4D simulations of the production setup, which are computationally graphics intensive, whereas toolpath failure detection requires actual 3D printing of test parts from the geometry to identify areas prone to failure while 3D printing, which is physically tedious. Both these processes, being computationally and physically intensive, have largely curtailed designers from simulating and exploring large sets of design options with varying geometries and toolpath configurations. To overcome this and allow designers to explore large sets of design possibilities, this paper proposes two novel computational metamodels capable of performing robotic toolpath clash detection and failure detection with significantly reduced times than the earlier approaches. The developed metamodels were used to rapidly simulate large sets of building design options for robotic constructability in the NASA 3D-Printed Mars Habitat Challenge.

clash detection↗
DOE OSTI2020

Bicontinuous nanoporous design induced homogenization of strain localization in metallic glasses

Bicontinuous nanoporous metallic glasses (MG) synergize the outstanding properties of MGs and open-cell nanoporous materials. The low-density and high-specific-surface-area of bicontinuous nanoporous structures have the potential to enhance the applicability of MGs in catalysis, sensors, and lightweight structural designs. In this work, we report molecular dynamics simulations of tensile loading deformation and failure of bicontinuous nanoporous Cu 64 Zr 36 MG with 55% porosity and 4.4 nm ligament size. Results indicate an anomalous mechanical behavior featuring delocalized plastic deformation preceding ductile failure. The deformation follows two mechanisms: i) Necking of ligaments aligned with the loading direction and ii) progressive alignment of randomly oriented ligaments. Failure occurs at 0.16 strain, following massive rupture of ligaments. This work indicates that a bicontinuous nanoporous design is able to effectively delocalize strain localization in a MG due to a combination of size effect on the ductility of MGs resulting in nano ligaments necking and progressive asynchronous alignment of ligaments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reliability Prediction using FMEA, FTA, and Related Techniques [Slides]

Summary: We presented a reliability analysis framework. We made point estimates of reliability using reliability block diagrams, fault trees, success trees and estimates with uncertainty using expert elicitation, Monte Carlo simulation, Bayesian analysis. Expert elicitation of failure modes and probabilities is labor-intensive, but critical. Bayesian analysis updates information from expert elicitation with data from reliability and aging tests (aging/compatibility data are needed to estimate lower-bound reliabilities at end of life). Estimation by more than one method helps insure consistency and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Hanford Double Shell Waste Tank Corrosion Studies- Final Report FY2024

For fiscal year (FY) 2024, the Savannah River National Laboratory (SRNL) focused on two experimental tasks related to Hanford Double Shell Tank (DST) chemistry and integrity. The first task focused on understanding risk of corrosion due to formation of either continuous layers or discrete patches of solids on the tanks’ inner sidewalls and bottoms. The objective of this task was to determine the effect of solid deposits on the corrosion risk to the tank bottom, and whether a combination of scale, saltcake and loose solids lead to under deposit corrosion. Electrochemical testing in a single cell arrangement was employed to determine a threshold inhibitor level for the interstitial liquid in the solids, above which carbon steel was not susceptible to corrosion beneath the solids. The inhibitor levels tested were based on the previously determined probability of failure for carbon steel in a simulated waste (i.e., Pitting Factor).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Sodium Pool

Mechanistic source term (MST) analyses are likely to be an important part of advanced reactor licensing applications. For the purpose, an MST analysis code applicable to newly introduced advanced reactors, called SRT (Simplified Radionuclide Transport) code, has been developed by Argonne National Laboratory. SRT can track overall behaviors of radionuclides especially in metal fuel-based sodium fast reactors (SFRs) and microreactors. Throughout the simulation, migration inside fuel pins before failure, interaction with coolant (for SFR), removal/leakage in cover gas and containment (or confinement), and environmental dose impacts are considered alongside radioactive decay for short-lived nuclides. Among the postulated process, pool scrubbing phenomenon, especially under sodium pool condition, has been identified as high importance with limited supportive data. The phenomenon plays a crucial role in assessing the degree of radiological impacts as radioactive aerosols or vapors are efficiently and effectively removed during the process. To provide validation basis for SRT in assessing pool scrubbing performance inside sodium pools, the University of Wisconsin-Madison performed tests including extensive parametric effects. Separate effect tests were conducted to directly evaluate the SRT models and to estimate degree of contribution by each contributing factor. Specifically, bubble size, aerosol size, aerosol density, aerosol concentration, pool depth, system temperature, and bubble swarm effects were considered. According to the parametric effects, decontamination performance enhances with decreasing bubble size, large density, and deeper pool height. Aerosol concentration provides no effect for the whole range of interest, and pool temperature variation shows minor effects under the considered temperature condition. When multiple bubbles are injected generating a bubble swarm condition, DF performance further enhances by bubble interactions and turbulence characteristics. The measurement shows the exceptional importance of aerosol size range considered, with the lowest decontamination, where most radionuclides are expected to escape.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Reliability estimation of an advanced nuclear fuel using coupled active learning, multifidelity modeling, and subset simulation

Tristructural isotropic (TRISO)-coated particle fuel is a robust nuclear fuel and determining its reliability is critical for the success of advanced nuclear technologies. However, TRISO failure probabilities are small and the associated computational models are expensive. Here we used coupled active learning, multifidelity modeling, and subset simulation to estimate the failure probabilities of TRISO fuels using several 1D and 2D models. With multifidelity modeling, we replaced expensive high-fidelity (HF) model evaluations with information fusion from two low-fidelity (LF) models. For the 1D TRISO models, we considered three multifidelity modeling strategies: only Kriging, Kriging LF prediction plus Kriging correction, and deep neural network (DNN) LF prediction plus Kriging correction. While the results across these multifidelity modeling strategies compared satisfactorily, strategies employing information fusion from two LF models called the HF model least often. Next, for the 2D TRISO model, we considered two multifidelity modeling strategies: DNN LF prediction plus Kriging correction (data-driven) and 1D TRISO LF prediction plus Kriging correction (physics-based). The physics-based strategy, as expected, consistently required the fewest calls to the HF model. However, the data-driven strategy had a lower overall simulation time since the DNN predictions are instantaneous, and the 1D TRISO model requires a non-negligible simulation time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Power Electronics Materials and Bonded Interfaces - Reliability and Lifetime

The emergence of wide-bandgap (WBG) devices in power electronics packages has brought significant attention to the reliability and thermal performance of the bonded interfaces. These interfaces, such as die-attach and substrate-attach, need to withstand the severe stress induced on them by the high-temperature (> 200°C) operation of the WBG devices. Sintered silver and transient liquid phase alloys are key bonded interfaces that can enable the high-temperature operation of WBG devices, however it is essential to evaluate their reliability and study the possible failure mechanisms before they can be used in power electronics packages. Sintered silver is already being adopted in commercial packages as a die-attach, but large-area substrate-attach offers additional reliability challenges. In this project, we investigate the reliability of these interface materials under accelerated thermal cycling and correlate these results with finite element simulations to develop lifetime prediction models. The failure mechanisms of these bonded interfaces are also studied.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Modeling Time to Failure in Potential-Induced Degradation of Silicon Solar Modules based on Quantitative Sodium Kinetics

We present a physics-based simulation environment for the prediction of potential-induced-degradation (PID) failure time in silicon solar modules. Our model is based on experimental sodium migration kinetics obtained from our recently developed trap-corrected bias-stress-temperature method and from secondary ion mass spectrometry measurements. Furthermore, we model sodium migration to stacking faults and subsequent formation of conductive pathways through the p-n junction, which causes a decrease of the shunt resistance and loss of efficiency.

42 ENGINEERING↗