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

Stress accommodation in nanoscale dolan bridges designed for superconducting qubits

Josephson junctions are the principal circuit element in numerous superconducting quantum information devices and can be readily integrated into large-scale electronics. However, device integration at the wafer scale necessarily depends on having a reliable, high-fidelity, and high-yield fabrication method for creating Josephson junctions. When creating Al/AlO x based superconducting qubits, the standard Josephson junction fabrication method relies on a sub-micron suspended resist bridge, known as a Dolan bridge, which tends to be particularly fragile and can often times fracture during the resist development process, ultimately resulting in device failure. In this work, we demonstrate a unique Josephson junction lithography mask design that incorporates stress-relief channels. Our simulation results show that the addition of stress-relief channels reduces the lateral stress in the Dolan bridge by more than 70% for all the bridge geometries investigated. In practice, our novel mask design significantly increased the survivability of the bridge during device processing, resulting in 100% yield for over 100 Josephson junctions fabricated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗
DOE OSTI2020

Supporting Bioproducts Industry Growth with a System Dynamics Decision-Support Tool

A bio-based economy requires chemical products as well as fuels to be produced from biomass. Although a variety of universities, government agencies, start-ups and established firms have engaged in bioproduct development, many projects have failed to reach the point of commercialization and commercialized bioproducts struggle to capture and maintain market share. To date, there has been no general research into the factors that contribute to bioproduct failure or success. This work presents the Bioproduct Transition Dynamics (BTD) system dynamics model, a decision-support tool that simulates the bioproduct development process from pre-piloting research through construction of the first commercial-scale plant, as well as the processes of obtaining funding from investors and government agencies. The core of the BTD is a feedback loop between bioproduct developers and funders, which relates development progress measured with indicators such as net present value to funders’ decisions to continue investing. External factors such as feedstock prices, market size and growth, and the existence of bioproduct consumers also influence a development project’s chances of receiving follow-on funding. Virtually any bioproduct can be represented with the BTD: direct replacements, performance-advantaged products, niche and commodity markets can all be modeled. The goal of the BTD project is to inform decisions made by developers in both established firms and start-ups, investors, government agencies, and other stakeholders interested in growing the nascent U.S. bioproducts industry. This talk will cover the general structure and functionality of the BTD, and present results from an analysis performed with the BTD to demonstrate its use as a decision-support tool and the insights it can provide. The goal of the analysis is to identify the most critical factors that lead to direct replacement and performance-advantaged bioproduct projects emerging successfully from the “Valley of Death”, and to determine if these factors differ between the two bioproduct types.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Explicit Dynamic Impact Analysis of the Building 3525 Packages

The safe transport and handling of radioactive material containers is critical to ensure the protection of personnel, facilities, and the environment. Accidental drops during handling pose a significant hazard, necessitating robust design and analysis to mitigate potential consequences. This study presents an explicit dynamic impact analysis of two container designs, the Baby Sugarman and the Core Conduction Cooldown Test Facility (CCCTF) Tall Boy, used in Building 3525 at Oak Ridge National Laboratory. The analysis employed nonlinear finite element modeling in LS-DYNA to simulate free-fall impacts from a height of 16.5 ft across eight critical orientations. The study evaluated deformation, strain failure, and shielding integrity under these extreme conditions. Results demonstrated that both designs maintain structural integrity, with no breaches or loss of shielding under specified administrative controls. Notably, a maximum shielding loss of 10 lbs was observed during a top-corner drop, highlighting the need to limit drop heights to less than 1 ft for certain orientations. These findings underscore the effectiveness of engineering controls and design measures in preventing containment breaches and ensuring compliance with safety requirements. This work provides a comprehensive methodology for impact analysis, contributing valuable insights to the field of radioactive materials packaging and transportation safety.

Martinez, Oscar [ORNL] (ORCID:0000000181814046)↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

Gear Test Assembly - Experimental Testing and Analysis of Gears and Bearings - FY2023

The Gear Test Assembly (GTA) is an experimental test apparatus built and installed in the Mechanisms Engineering Test Loop (METL) at Argonne National Laboratory (ANL). The focus of GTA has been to test larger radial spur gears sets machined from Inconel 718 along with a variety of bearing components that are to be used in the fuel handling machines built for use in pool-type sodium fast reactors. The testing in GTA also informs the design and material choices for the full-scale Gripper Test Assembly (GrTA) being developed by Argonne. To date, the Inconel 718 gears used in all six GTA experimental campaigns have completed over 23.2 million revolutions under various loads equivalent to approximately 37,156 simulated fuel assembly maneuvers (removal from and reinsertion to core) using a conventional height of a core assembly. While the heat treated and untreated tapered roller bearings and cylindrical pin thrust bearings used in the first five experimental campaigns have experienced mechanical or material failure after completing a range of fuel assembly maneuvers ranging from as early as 575 operations in Campaign #5 to 9800 operations in Campaign #1, Campaign #6 was the most successful campaign, achieving 22,563 maneuvers and was concluded before any catastrophic failure occurred. In this report we will present the results from the most recent campaigns, #5 and #6. The longevity of Campaign #6 was likely aided by a variety of adjustments made to the materials and components used and to the experimental testing procedures employed. In terms of materials and components, new gearboxes were fabricated, as opposed to using the old gear boxes which had been impacted by the previous bearing failures, and ball bearings were replaced with tapered roller bearings. The adjustments to the experimental procedure focused on a lowering of the maximum torque applied to the gears and bearings and the increase of the duration of each simulated fuel assembly maneuver conducted. This report will present and discuss the changes made to the materials, components, and procedures and the subsequent results for Campaign #5 and Campaign #6.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

An Integrated High-performance Computing and Digital Real-time Simulation Testbed to Benchmark Closed-loop Load Shedding Algorithms in Power Systems

An integrated testbed using digital real-time simulator (DRTS) and a high-performance computing (HPC) cluster is presented here to compare speed and performance of computational schemes to mitigate time-critical issues in electric power systems. The first approach in this testbed validation is taken by running a set of closed-loop load shedding algorithms to compare and contrast two paradigms of arresting cascading failure propagation. Two algorithms involve solving DC and AC power flow model-based optimization problems to compute load shedding at different buses, while a model-based stochastic search using parallel computing provides a viable alternative. The algorithms are implemented in the DRTS-HPC testbed for the IEEE 14-bus benchmark transmission system. As a proof of the concept, simulation results are presented for implementation of closed-loop load-shedding algorithms for cascading failures in the DRTS-HPC testbed

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Towards Polymer-Free, Femto-Second Laser-Welded Glass/Glass Solar Modules

This project explores the use of femto-second (fs) lasers to form glass-to-glass welds for hermetically sealed, polymer-free solar modules. Low iron solar glass coupons were welded together without the use of glass filler using a fs laser with dedicated optics to elongate the focal plane parallel to the incident beam. The resulting welds were then stress tested to failure to reveal the critical stress intensity factor, KIc. These values were used in a structural mechanics model of a 1 m by 2 m glass/glass module under a simulated static load test. The results show that the fs laser welds are strong enough for a suitably framed module to pass the IEC 61215 static load test with a load of 5400 Pa. Key to this finding is that the module must be framed and braced, and the glass must be ribbed to allow pockets for the cells and welds inside the border of the module. The result is a module design that is completely polymer-free, hermetically sealed, has improved thermal properties, and is easily recycled.

femto-second laser welded↗
DOE OSTI2026

Gear Test Assembly - Experimental Testing and Analysis of Gears and Bearings - FY2026

The Gear Test Assembly (GTA) is an experimental test apparatus built and installed in the Mechanisms and Engineering Test Loop (METL) at Argonne National Laboratory (ANL). While designed to accommodate a variety of components intended for use in advanced compact in-vessel transfer machine, the focus of GTA been the testing of large radial spur gears. The performance of the components used in GTA motivates the choices made in the design, components, and material choices for the forthcoming METL test articles; the Gripper Test Assembly (GrTA), and the Bearing Test Assembly (BTA). To date, GTA has completed nine experimental testing campaigns and achieved over 21 million revolutions under prototypic sodium fast reactor loads and operating conditions. This report provides an overview and results for the two most recent experiments Campaign #8 and Campaign #9. `Campaign #8 was the first campaign to make use of a set of nickel alloy radial spur gears, made from Hastelloy C-22HS, and combined with heat-treated tapered roller bearings. The bearings and torque profile used allowed Campaign #8 to be directly compared with Campaign #1 where Inconel 718 gears were used instead. Campaign #8 attained 1,058,880 shaft revolutions within 23% of the 1,314,855 revolutions achieved in Campaign #1. These results demonstrate that the heat-treated tapered roller bearings show consistent increased lifetimes, approximately three times longer than for non heat-treated taper rollers, and that the particular material of the gear does not play a substantial role in performance. Similar to Campaign #7 where bearing failure resulted in fragmentation and transport of bearing material into the gear teeth led to tooth damage in the Inconel gears, a similar process led to the damage of the Hastelloy gears. Subsequent NDE analysis showed that the damage was only limited to the tooth surface and does not extend into the gear interiors. Based on the NDE reports, the Hastelloy gears can be reconditioned and put back into service for future campaigns as was done for the Inconel gears. Reconditioned Inconel 718 gears damaged in Campaign #7 were returned to service in Campaign #9 and paired with non heat-treated tapered roller bearings. In Campaign #9 a novel set of operating conditions were implemented whereby the gears and bearings were cycled between hot, submerged refueling conditions where fuel handling operations occurred and cold, dry, and inert conditions where the sodium was drained and the vessel was allowed to cool simulating long, ex-vessel storage of the fuel handling machine. A total of eight of these operating cycles were achieved before a bearing failure occurred after the accumulation of 475,692 shaft revolutions were achieved. During Campaign #9, one of the drive motors was damaged and over the course of a 240 day standby period, the entire GTA drive system was upgraded from a 480VAC system to a 240VAC system that provided the opportunity to consolidate wall-space in the METL facility and provide more space for additional test articles being installed in the coming year. The over 475,000 shaft revolutions achieved in Campaign #9 eclipses the previous record for lifetime of non heat-treated tapered roller bearings held by Campaign #3 of 392,000 shaft revolutions. Further analysis and future campaigns hope to shed light on the cause of this approximately 17% increase in lifetime.

42 ENGINEERING↗
DOE OSTI2021

Thermophysical Properties of Mixtures of Thorium and Uranium Nitride

Abstract The miscibility, lattice parameter, and thermophysical properties of (Th 0.2 U 0.8 )N and (Th 0.5 U 0.5 )N have been investigated. It is shown that additions of thorium nitride (ThN) to uranium nitride (UN) increases the thermophysical performance of the mixed nitride fuel form in comparison to reference UN. In the more dilute limit, additions of ThN serve as a burnable neutronic poison and reduces the change in k eff over the lifecycle of the fuel. At higher concentrations, additions of ThN serve as a significant fertile source of 233 U. Where appropriate, comparisons to previous work on UN + PuN mixtures are made, as this is a comparable fuel form for potential fast reactor concepts, and a suitable point of contrast in the possible design space afforded by mixed (Th x U 1 − x )N fuel forms. The data from this work are the input parameters for finite element modeling of the temperature distribution in a compact reactor. The results of modeling and simulation of this core design are shown for the case of steady-state operation and during double, adjacent heat pipe failure.

Parker, S. S. (ORCID:000000018597440X)↗
DOE OSTI2021

Application of spectroscopic evidence to diffuse layer model (DLM) parameter estimation for cation adsorption onto ferrihydrite in single- and bi-solute systems

The diffuse layer model (DLM) is one of the most frequently used surface complexation models (SCM) for predicting cation adsorption onto ferrihydrite. However, the DLM with its existing database sometimes fails to predict cation adsorption, especially in multi-solute systems. Notably, the failure is often attributed to inaccurate selection of the surface complexes and/or underestimation of the site density for DLM formulation. DLM simulations of batch adsorption data conducted with cadmium, copper, lead, and zinc on ferrihydrite in CO2 free systems demonstrates the ability of the DLM to fit single-solute data reasonably well but highlights its limitations with respect to predicting bi-solute data. Modifications to the DLM were made in two steps. First, spectroscopic evidence was used to guide cation complexation reactions at the ferrihydrite surface. Second, surface characteristics (reactive site density, specific surface area, and acidity constants) of ferrihydrite were updated according to recent literature. The results from this study suggest that when spectroscopic verification is employed and surface site density is increased, a simple two-site DLM can predict cation adsorption on ferrihydrite in single-solute systems for pH data above 3.5 and in most of the bi-solute systems tested.

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

Comparing modelling approaches for a generic nuclear waste repository in salt

This paper contains a comparison of five modelling approaches for a simplified nuclear waste repository in a domal salt formation. It is the result of a four-year collaboration between five international teams on Task F of the DECOVALEX-2023 project on performance assessment modelling. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. This work demonstrates how these objectives are accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Participating teams made a wide range of model assumptions, ranging from compartmentalized networks to full 3D models of the salt formation and repository. Despite differences in the modelling strategies, all models indicate that salt compaction and diffusion of radionuclides in brine are key processes in the repository. For the isothermal spent nuclear fuel and vitrified waste scenario with multiple early failures considered, all models indicate little of the disposed radionuclides will migrate beyond the repository seal over the 100,000-year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. Disparities between the models are believed to be due to differing simplifications from the conceptual model.

DECOVALEX↗
DOE OSTI2025

The effect of grain boundary yield and disconnections on grain boundary fracture strength

During high-temperature interfacial strain, e.g., creep or sintering, grain boundaries will sample non-equilibrium states associated with the non-equilibrium disconnections that mediate strain. This work tests the hypothesis that those non-equilibrium states should affect the mechanical properties of the grain boundary. High-temperature bicrystal tensile strain experiments performed on Sc-doped ZrO 2 as a function of strain rate enable the observation of fracture before and after grain boundary yield. Yielded boundaries exhibit a median fracture strength ≈ 50 % lower than those that did not yield. Molecular dynamics simulations performed on boundaries with and without interstitial-type disconnection loops indicate that the presence of those loops reduce the failure stress.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fracture Failure Analysis for MQXFA Magnet Aluminum Shells

The High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is approaching the production phase of the US-contributed Q1 and Q3 Interaction Region Quadrupoles (MQXFA). The structures for the MQXFA prototypes were design and inspected by the US-LARP (LHC Accelerator Research Program), AUP developed criteria, which will be used for the pre-series structures. As the first two full-length prototypes with 4.2 m magnetic length, MQXFAP1 and MQXFAP2, were designed and assembled at Lawrence Berkeley National Laboratory (LBNL), and tested at Brookhaven National Laboratory (BNL). The end aluminum short shell of MQXFAP2 was fractured along the shell length during the test, and tests were stopped. Analytical and Finite Element analysis were performed in light of the graded procedure defined in the Structure Design Criteria to investigate the fracture failure for MQXFAP2. In this paper, we report the fracture analysis of the current shell design, including the elasto-plastic simulations with submodel technique, and calculations with Linear Elastic Fracture Mechanics (LEFM). Test material properties are also presented. The results of this analysis explain why the end shell of MQXFAP2 failed, and suggest fillets on the end shell notches to meet the margin specified in the Structural Design Criteria.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Dynamic Disruption Resilience in Intermodal Transport Networks: Integrating Flow Weighting and Centrality Measures

Resilient intermodal freight networks are vital for sustaining supply chains amid increasing threats from natural hazards and cyberattacks. Transportation resilience has been widely studied; understanding how random and targeted disruptions affect structural connectivity and functional performance remains a key challenge. To address this, this study evaluates the robustness of the US intermodal freight network, which consists of rail and water modes, using a simulation-based framework that integrates graph-theoretic metrics with flow-weighted centrality measures. Disruption scenarios are examined, including random failures as well as targeted node and edge removals based on static and dynamically updated degree and betweenness centrality. To reflect more realistic conditions, flow-weighted degree centralities (WDC) and partial node degradation are considered. Two resilience indicators are used: (1) the size of the giant connected component to measure structural connectivity; and (2) flow-weighted network efficiency (NE) to assess freight mobility under disruption. The results show that progressively degrading nodes ranked by WDC to 60% of their original functionality causes a sharper decline in normalized NE, for up to approximately 45 affected nodes, than complete failure (100% loss of functionality) applied to nodes targeted by weighted betweenness centrality or selected at random. This highlights how partial degradation of high-tonnage hubs can produce disproportionately large functional losses. The findings emphasize the need for resilience strategies that go beyond network topology to incorporate freight flow dynamics.

42 ENGINEERING↗
DOE OSTI2024

Dataset for manuscript "Consequences of the failure of equipartition for the p-V behavior of liquid water and the hydration free energy components of a small protein"

Previously, we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step dt that is less than or equal to 0.5 fs to ensure equipartition between translational and rotational modes. We extended that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examined integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small time-step, dt, is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, is seen to be sensitive to dt. Using the mean volume estimated from the NpT simulation, we examined the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to dt. In going from a time-step of 2 fs to a time-step of 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein. The data-set contains the simulation metadata and log files that support the claims noted above.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The impacts of horizontal resolution on the seasonally-dependent biases of the northeastern Pacific ITCZ in coupled climate models.

Double-ITCZ has puzzled climate model community for more than two decades. Here, we found that over the eastern Pacific, precipitation and sea surface temperature (SST) biases in the NCAR CESM1 are seasonally-dependent, with positive (negative) precipitation/SST biases during boreal summer/fall (winter/spring), although the easterly wind bias persists around the year. This seasonally-dependent bias of precipitation/SST is found to be caused by the failure of the model to reproduce the climatological seasonal wind reversal of North American monsoon. During winter/spring, the climatological easterly wind dominates, so the stronger wind speed in the model enhances the evaporation and lowers the SST. It is opposite when the climatological wind turns to westerly during summer/fall. The easterly wind bias also occurs in the atmospheric model when the observation SST is prescribed, suggesting it is of atmospheric origin. Further, the easterly wind bias is mainly evident in the lower troposphere and quite independent of time, implying the effect of complex and narrow-ranged Central American topography, which is not depicted well in the coarse model resolution. When the model resolution is doubled, both SST and precipitation are improved with the reduced easterly wind bias. During boreal spring, when ITCZ bias is most significant, the northern and southern ITCZ can be improved by 29% and 18.8% respectively in the higher resolution. Finally, we suggest that these seasonally-dependent biases over the eastern Pacific are quite universal among 37 CMIP5 coupled models, and the ongoing CMIP6 project will provide a valuable opportunity to examine whether the bias can be reduced by increasing model resolution.

Song, Fengfei↗