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 523 records · Page 29

New Physics in the Hot Electron Furnace (Final Report)

The primary purpose of this award was experiments and simulations related to the idea of a ‘Hot Electron Furnace’ where nonlocal heat transport would allow unique new experimental opportunities at x-ray free electron lasers. The period of this award included beamtime at the MEC endstation of the LCLS, and the execution of those experiments was a primary goal. Unfortunately, the beamtime had limited success due to (1) inability of the LCLS to deliver the expected two-color operations mode, and (2) failure of the MEC endstation’s own software to save the cavity timing data that would have allowed inference of time delays for the few hours when two-color mode briefly operated. The beamtime was then switched to operate in single-color (one pulse) mode, and studies were performed of x-ray heating and probe, simultaneously. Two interesting manuscripts are finally near submission, having been delayed by student graduations, need for collaborator simulations of x-ray heating processes, and then COVID. I am hopeful that they will soon be submitted (especially the paper on x-ray heating of MgO) because one of the principal students involved has returned from industry to take a postdoc position at SLAC. Below, we first list all publications and other work products (including the two that are nearing submission) and then summarize several outcomes from this effort. During the term of this award three UW physics graduate students were at least partially supported, and all three received their Ph.D.’s with significant components of their supported work in the dissertation.

47 OTHER INSTRUMENTATION↗

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗

Temperature Rise on Liner Surfaces of Fuel Cell Electric Vehicle Tanks during Fueling Process

Herein, the temperature of the inner tank walls, or plastic liners, of composite pressure tanks in fuel cell electric vehicles during fueling using three‐dimensional computational fluid dynamics (CFD) models is evaluated. The liner materials are limited to 85.0 °C to prevent thermal stress causing material failure that would result in a hydrogen leak. Therefore, the temperature of hydrogen gas must be limited to below 85.0 °C during the fueling process as dictated by the current fueling protocol. However, there are limited experimental or simulation data confirming that the temperature changes do not exceed the threshold. Herein, the liner temperatures with CFD tank models for two sizes of type IV tanks representative of the upper and lower system bounds that are close to the SAE J2601 fueling protocol (36.0 and 244.0 L) are evaluated. First, each model's reliability is validated with experimental data and then analyzed, and the data are used to evaluate the maximum hydrogen and liner temperatures under real‐world fueling conditions. The evaluation shows that the maximum liner surface temperature of each tank model is at least 7 °C lower than that of the hydrogen. Additionally, there is at least 12 °C difference found between the upper limit and actual liner temperatures.

33 ADVANCED PROPULSION SYSTEMS↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗

The influence of increasing atmospheric CO 2 , temperature, and vapor pressure deficit on seawater-induced tree mortality

We report increasing seawater exposure is killing coastal trees globally, with expectations of accelerating mortality with rising sea levels. However, the impact of concomitant changes in atmospheric CO 2 concentration, temperature, and vapor pressure deficit (VPD) on seawater-induced tree mortality is uncertain. We examined the mechanisms of seawater-induced mortality under varying climate scenarios using a photosynthetic gain and hydraulic cost optimization model validated against observations in a mature stand of Sitka-spruce (Picea sitchensis) trees in the Pacific Northwest, USA, that were dying from recent seawater exposure. The simulations matched well with observations of photosynthesis, transpiration, nonstructural carbohydrates concentrations, leaf water potential, the percentage loss of xylem conductivity, and stand-level mortality rates. The simulations suggest that seawater-induced mortality could decrease by ~16.7% with increasing atmospheric CO 2 levels due to reduced risk of carbon starvation. Conversely, rising VPD could increase mortality by ~5.6% because of increasing risk of hydraulic failure. Across all scenarios, seawater-induced mortality was driven by hydraulic failure in the first two years after seawater exposure began, with carbon starvation becoming more important in subsequent years. Changing CO 2 and climate appear unlikely to have a significant impact on coastal tree mortality under rising sea levels.

54 ENVIRONMENTAL SCIENCES↗

Modeling gas, hydrates, and slope stability on the U.S. Atlantic margin during Pleistocene glacial cycles

Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.

Hydrate dissociation↗

Laser ablation spectrometry for studies of uranium plasmas, reactor monitoring, and spent fuel safety

Nuclear security is one of the defining challenges of our time. Nuclear threats range from deliberate dispersal of radioactive material to contaminate the vital infrastructure to diversion and smuggling of special nuclear material for clandestine nuclear programs and nuclear terrorism, respectively. There is an associated need to develop and sustain the nuclear forensics capabilities, which requires the understanding of complex processes that occur in plasmas of nuclear materials. The area of nuclear safety has seen a resurgence of public interest, and there is a concomitant need to safely store used nuclear fuel and detect structural material failure in nuclear power systems, especially in innovative reactor designs envisioned for future adoption. Laser-produced plasmas are complicated extreme environments that can generate intense and rich, highly specific signatures of nuclear and radiological materials, which can then be explored in a wide range of applications. They include interdiction and rapid detection of nuclear materials, including their isotopic composition, detection over long distances, laboratory simulation of weapons effects, monitoring the condition of structural materials in dry cask storage containers, and novel instrumentation for nuclear power systems. We present a compilation of recent representative examples of the application of laser spectroscopy, and laser-induced breakdown spectroscopy in particular, to nuclear safety and security problems. A case is made that spectroscopic techniques based on laser-produced plasmas offer complementary, and sometimes unique, capabilities that motivate the continued exploration of their efficient production and understanding of the signatures they produce

(020.3260) Isotope shifts, (140.3440) Laser-induce↗

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗

Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites

A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.

Al-SiC↗

Failure Probability Constrained AC Optimal Power Flow

Despite cascading failures being the central cause of blackouts in power transmission systems, existing operational and planning decisions are made largely by ignoring their underlying cascade potential. This paper posits a reliability-aware AC Optimal Power Flow formulation that seeks to design a dispatch point which has a low operator-specified likelihood of triggering a cascade starting from any single component outage. By exploiting a recently developed analytical model of the probability of component failure, our Failure Probability-constrained ACOPF (FP-ACOPF) utilizes the system's expected first failure time as a smoothly tunable and interpretable signature of cascade risk. Here, we use techniques from bilevel optimization and numerical linear algebra to efficiently formulate and solve the FP-ACOPF using off-the-shelf solvers. Extensive simulations on the IEEE 118-bus case show that, when compared to the unconstrained and N-1 security-constrained ACOPF, our probability-constrained dispatch points can significantly lower the probabilities of long severe cascades and of large demand losses, while incurring only minor increases in total generation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗

Advanced Monitoring and Control in the ANL METL Facility Using an Engineering Digital Twin

The potential benefits of using an engineering digital twin to achieve greater autonomy for monitoring and control functions in advanced reactors was investigated for the Mechanisms Engineering Test Loop (METL) liquid sodium facility at Argonne National Laboratory. The METL sodium purification system served as a representative system as it requires significant human in-the-loop interaction to accomplish its design function. The objective was to demonstrate how real-time operation could be automated while preserving oversight of the operator for ensuring that the system design functions are met. A digital twin model of the purification system was developed for both the cold trap purification loop and plugging meter diagnostic loop using information from the METL piping and instrumentation diagram (P&ID). Automated monitoring and diagnosis of component degradation in the METL facility was demonstrated in tests using the PRO-AID health monitoring software with the digital twin model incorporated in the library of components. Component failures were introduced and were successfully diagnosed in real time. These tests serve to demonstrate an advanced monitoring capability able to differentiate sensor degradation from component degradation, to generate a rank ordering of probabilities of different failure mechanisms that serves to circumvent the false alarm problem with current anomaly detection methods, and how facility monitoring can be transformed from anomaly detection to identification of a specific fault. Automated control of the purification system was demonstrated through simulations that exercised a model predictive controller designed using the digital twin model. Results of these simulations compared favorably with experimental data showing very good reference tracking response with negligible overshoot. In conclusion, these pilot tests and simulations successfully demonstrated the use of a digital twin for improved automation of monitoring and control. It was shown how the digital twin enables switching between control modes from cold trap operation where impurities are removed to plugging meter operation where impurity concentrations are measured.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Wildfire-Power Grid Interactions: Feedback, Impacts, Monitoring, Modeling, and Mitigation Strategies

Wildfires are increasingly interacting with electric power systems through a two-way hazard chain: fires damage grid assets and trigger cascading outages, while grid faults can ignite new fires under hot, dry, and windy conditions. This review synthesizes the state of knowledge across five domains: (i) physical impacts of flames, heat, and smoke on lines, towers, insulators, and substations; (ii) power-infrastructure-initiated ignitions via conductor clash, high-impedance faults, and corona discharge; (iii) widespread blackouts and disproportionate societal impacts; (iv) multi-scale monitoring spanning laboratory tests, in-situ and grid-integrated sensors, and Earth observation; (v) coupled modeling that links fire behavior with grid operations; and (vi) technological and strategic mitigation pathways spanning prevention, response, and recovery. We integrate these domains into a novel 'feedback-aware' socio-technical framework. Through a longitudinal analysis (2005-2025) of global incidents, we identify that while vegetation contact remains the most frequent ignition source, aging infrastructure failure has emerged as a critical driver of catastrophic 'mega-fires'. We further identify persistent gaps, including limited interoperability of high-frequency grid and environmental data, scarce real-time data assimilation, and under-developed equity metrics for outage management. We conclude by outlining a research agenda to (1) deploy interoperable sensing architectures, (2) advance feedback-coupled fire-grid simulations, and (3) evaluate mitigation portfolios through techno-economic and fairness lenses. Recognizing wildfire-grid interactions as coupled socio-technical systems is essential for protecting infrastructure and communities and for ensuring reliable, sustainable electricity in a changing world.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Virtual Pressure Sensor for Electronic Expansion Valve Control in a Vapor Compression Refrigeration System

Virtual sensor technology, which uses simulation models/algorithms to calculate a value to represent an unmeasured variable or replace a directly measured reading, has attracted many studies in the heating, ventilation, air conditioning and refrigeration (HVAC&R) industry. However, most virtual sensor technologies are developed for fault detection and diagnostics (FDD) purposes, which generally compare the virtual sensor values with actual measured values to detect if any fault occurred and identify the causes that led to the fault. It is rare to see studies focus on control performance of virtual sensors after substituting an actual sensor. This is particularly important for the system with no redundant sensor since a virtual sensor is the most effective way to operate the system in the desirable region when any sensor failure occurs. To address this gap, this paper develops a new virtual pressure sensor technology to substitute the actual pressure measurement for electronic expansion valve (EXV) control in a vapor compression refrigeration system by integrating compressor and valve characteristics. The control performance of this proposed virtual pressure sensor technology under various operating conditions is validated with experimental data. Closed loop EXV control simulations with the proposed virtual pressure sensor are conducted, and the results are analyzed.

47 OTHER INSTRUMENTATION↗

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗