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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A Flexible and Generic Functional Mock-up Unit Based Threat Injection Framework for Grid-interactive Efficient Buildings: A Case Study in Modelica

Grid-interactive efficient buildings (GEBs) have been considered as an important asset to support the power grid reliability by utilizing the demand flexibility offered by GEBs. GEBs are enabled by advances in sensors and controls, and the communication between building equipment, whole buildings, and the grid. The integration of different building technologies and network-based communication system makes GEBs vulnerable to passive threats such as equipment failure and active threats such as cyber-attacks. Modeling and simulation is an effective way to evaluate the impact of threats on the system performance. This paper proposes a generic and flexible threat injection framework for commonly-used building energy simulators such as EnergyPlus and Modelica to support threat modeling and evaluation. This framework leverages functional mock-up unit (FMU) to develop a general modeling interface for threat injection and simulation. A numerical case study using Modelica as a building energy simulator is conducted to demonstrate the capability of the framework for supporting single/multiple-order threat modeling and simulation of a GEB. Four threats and their combinations are injected on a Modelica-based threat-free building energy and control system, including operating supply fan at its full speed, remotely cycling the chiller on and off, blocking the chiller from receiving the chilled water supply temperature setpoints, and hijacking the global zone air temperature setpoint. Simulation results show that the cyber-attack that leads to short-term signal blocking has small effects on the system operation due to the "self-healing" feature of the heating, ventilation, and air-conditioning (HVAC) interactive control system. The threat that takes control of resetting the global zone air temperature setpoints has the most adverse impact on the system energy use, peak power demand, thermal comfort and the provision of demand flexibility. The combination of four threats have aggregative effects on the system but the effects are less than the additive effects of the individual threat.

Fu, Yanyang↗
DOE OSTI2024

Evaluation of 304L stainless steel deformation and failure under cyclic loading

The motivation for this work came from the simulation of components subjected to impact. In these problems, repeated loading can be generated by vibrational motion between internal parts of the component or by the occurrence of multiple impacts. Graphs of tensile equivalent plastic strain $(\overline{\varepsilon}_T^\rho)$ versus time from such cases had the character shown in Fig. 1. The equivalent plastic strain grew by alternating periods of accumulation and hold. The periods of hold indicate times when $d\overline{\varepsilon}_T^\rho$ did not accumulate. This can be due to the material unloading into the elastic range, or actually plastically reverse loading with hydrostatic stress σ h < 0 as was demonstrated for an impact simulation. Similar observations were made from the simulation of a 2-degree-of-freedom mass-elastic-plastic-spring system subjected to impact.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessing the performance of solar radiation management geoengineering simulations

Offsetting the global warming caused by anthropogenic increases in atmospheric greenhouse gases by deliberate injection of aerosols into the stratosphere is the most studied of solar radiation management geoengineering schemes. The long-term success or failure of such schemes in achieving their stated goals is assessed by comparing simulated geoengineered temperature, precipitation and tropical cyclones metrics to equivalent fields in the simulated targeted climate simulations. Results using available data sets from three single model stabilized climate target experiments and three multimodel climate change reduction experiments are presented and compared against a measure of internal variability. While all but one experimental scheme is successful in achieving their targeted global mean annual surface temperature, their success at regional scales varies significantly and is often larger than the internal variability metric used here.

climate model evaluation↗
DOE OSTI2021

Burst characteristics of advanced accident-tolerant FeCrAl cladding under temperature transient testing

Here assessment of burst characteristics of accident-tolerant fuel (ATF) claddings is essential to evaluate the safety margins of nuclear reactors, and better understanding can enable accelerated licensing for reactor concepts that include the new materials. Therefore, this study investigated the burst behavior of an ATF candidate of C26M, which is an iron–chromium–aluminum alloy (FeCrAl) under transient testing like the simulated loss-of-coolant accident (LOCA) conditions, except the water-quenching phase, in light-water reactors (LWRs). The effect of LOCA specimen length was assessed in terms of post-burst tube parameters. No critical length effect was determined on the burst pressure, burst location, and burst size measurements. Burst temperature showed larger variation likely due to its measurement approach used in this study. The post-test diametral strain was identified as the critical parameter to reduce the specimen length as compared to tube burst length or width. Postmortem optical metrology and digital image correlation were employed to determine the local strain-state during the simulated LOCA, which showed the loading path was close to equibiaxial conditions in regions away from the burst location, rather than internally pressurized conditions (or “plane-strain” tension). Pst-test microstructural characterizations of the FeCrAl revealed that ductile damage was present at the edge of the outer radial surface while the rest of the material was ruptured via cleavage at LOCA burst temperatures. Furthermore, internal grain boundary cracks were observed at locations away from the tube burst region. This behavior was considered to occur during the cooling down period of the simulated LOCA test. Overall, this study aimed to provide essential data for fuel performance code development by considering micromechanics of the deformation and failure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

FGMS-poster

Idaho National Laboratory (INL) performs post irradiation examination (PIE) of tri-structural isotropic (TRISO)-coated particle fuel to help qualify it for high temperature gas cooled reactors. TRISO fuel compacts are re-irradiated in the Neutron Radiography Reactor (NRAD) to generate the short lived fission products needed for fission product release testing. The Fuel Accident Condition Simulator (FACS) furnace and newly added Screen Neutron Irradiated Fuel for Failure (SNIFF) furnace heat the compacts in helium to temperatures of up to 2,000°C, prompting fission product release—predominantly gaseous xenon and krypton isotopes and condensable products such as cesium—from failed particles. These released isotopes are transported to a fission gas monitoring system (FGMS 1 or FGMS 3), where they accumulate in cryogenic cold traps and are quantified using high-purity germanium (HPGe) detectors. The addition of SNIFF and FGMS 3 increases throughput by enabling simultaneous testing of multiple compacts. Furthermore, automated INL developed software provides continuous, near-real time monitoring of fission product inventories and manages the liquid nitrogen cooling of the traps. These system enhancements improve the efficiency, data quality, and testing capacity of TRISO fuel performance evaluations.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2022

Chromium-coated cladding analysis under simulated LOCA burst conditions

Of the near-term accident-tolerant fuel concepts identified, chromium coatings have been shown to slow cladding oxidation without altering fuel system geometry or greatly affecting neutronic performance. To identify how coated cladding tubes perform under high-temperature accident conditions, pressurized-tube burst tests were conducted in the Severe Accident Test Station at Oak Ridge National Laboratory. To analyze how the coating affects cladding behavior, coated and uncoated cladding burst tests were simulated by using the BISON fuel performance code. Thermocouple data from these tests were fit into axial and azimuthal profiles and combined to generate 3D cladding surface temperatures, and pressure transducer data were compared until cladding failure. Additionally, the cladding temperatures at failure and the internal gas pressure evolution show relatively good agreement between the simulation and experiment results. Simulations were then extrapolated to demonstrate the coatings effectiveness to increasing the cladding burst temperature by using a parametric evaluation of the initial tube gas pressure. This work demonstrates the possibility of an increased cladding failure margin under transient conditions due to the addition of chromium coating, and more pertinently, an increased cladding failure margin from more accurate experimental characterization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanistic modeling of CO 2 well leakage in a generic abandoned well through a bridge plug cement-casing gap

Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.

58 GEOSCIENCES↗
DOE OSTI2022

Role of tilt grain boundaries on the structural integrity of WSe 2 monolayers

Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe 2 ) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe 2 , of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe 2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe 2 monolayer with tilted GBs. WSe 2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures.

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

Phase-Field Modeling of Thermally-Grown Oxide and the Induced Damage Evolution in Environmental Barrier Coatings

The advent of next-generation hydrogen-based engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the critical role of thermally grown oxides (TGO) in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and microstructural evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

Cheng, Tianle↗
DOE OSTI2025

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗
DOE OSTI2023

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Report on local data recovery approaches suitable for weather and climate prediction (Deliverable 1.3) (V.1.0)

Numerical weather and climate prediction rates as one of the scientific applications whose accuracy improvements greatly depend on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increasing average frequency of hardware and software failures forces users to review their algorithms and systems in order to protect simulations from breakdown. This report surveys approaches for fault-tolerance in numerical algorithms and system resilience in parallel simulations from the perspective of numerical weather and climate prediction systems. A selection of existing strategies is analyzed, featuring interpolation-restart and compressed checkpointing for the numerics, in-memory checkpointing, user-level failure mitigation-based and backup-based methods for the systems. Numerical examples showcase the performance of the techniques in addressing faults, with particular emphasis on iterative solvers for linear systems, a staple of atmospheric fluid flow solvers. The potential impact of these strategies is discussed in relation to current development of numerical weather prediction algorithms and systems towards the exascale. Trade-offs between performance, efficiency and effectiveness of resiliency strategies are analyzed and some recommendations outlined for future developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Mechanical properties of heterogeneous metallic glasses: Insights from brick-and-mortar designs

The brick-and-mortar architectural paradigm is commonly applied in structural design to achieve an exceptional combination of strength and plasticity due to its highly tunable mechanical attributes. This study applies molecular dynamics simulations of tensile loading to investigate the mechanical properties and failure mechanisms of brick-and-mortar metallic glasses (BMMGs). The focus is primarily on the implications of the aspect ratio of the bricks and interlayer thickness on the strength, plasticity, and deformation mechanisms. Results indicate that the failure mode is typified by multiple shear bands localized within the softer mortar regions, creating a staggered network under smaller aspect ratios that significantly enhances plastic deformation. However, a transition to a single dominant SB occurs when the aspect ratio exceeds 4.1. In addition, results indicate that a concurrent improvement in strength and plasticity can be achieved by modulating the distance between brick layers in BMMGs. This enhancement originates from the extensive generation and interaction of shear transformation zones. Furthermore, this study highlights that a desirable balance between strength and plasticity can be obtained in BMMGs with appropriate brick aspect ratios and interlayer thicknesses, providing a potential design strategy for the advancement of novel metallic glasses with superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗