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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Stochastic-Strength-Based Damage Simulation of Ceramic Matrix Composite Laminates

The Finite Element Analysis-Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to characterize and predict the progressive damage response of silicon-carbide-fiber-reinforced reaction-bonded silicon nitride matrix (SiC/RBSN) composite laminate tensile specimens. Studied were unidirectional laminates [0] (sub 8), [10] (sub 8), [45] (sub 8), and [90] (sub 8); cross-ply laminates [0 (sub 2) divided by 90 (sub 2),]s; angled-ply laminates [plus 45 (sub 2) divided by -45 (sub 2), ]s; doubled-edge-notched [0] (sub 8), laminates; and central-hole laminates. Results correlated well with the experimental data. This work was performed as a validation and benchmarking exercise of the FEAMAC/CARES program. FEAMAC/CARES simulates stochastic-based discrete-event progressive damage of ceramic matrix composite and polymer matrix composite material structures. It couples three software programs: (1) the Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC), (2) the Ceramics Analysis and Reliability Evaluation of Structures Life Prediction Program (CARES/Life), and (3) the Abaqus finite element analysis program. MAC/GMC contributes multiscale modeling capabilities and micromechanics relations to determine stresses and deformations at the microscale of the composite material repeating-unit-cell (RUC). CARES/Life contributes statistical multiaxial failure criteria that can be applied to the individual brittle-material constituents of the RUC, and Abaqus is used to model the overall composite structure. For each FEAMAC/CARES simulation trial, the stochastic nature of brittle material strength results in random, discrete damage events that incrementally progress until ultimate structural failure.

composite structures↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1985

GPS high dynamic receiver tracking demonstration results

Demonstration results are presented for a high dynamic GPS receiver. The receiver tested is a breadboard unit capable of tracking one simulated satellite signal in pseudorange and range rate. The receiver makes approximate maximum likelihood estimates of pseudorange and range rate each 20 ms, and tracks these observables using a third order filter with a time constant of 0.14 s. Carrier phase is not tracked, which eliminates the typical failure mode of loss of carrier lock associated with PLLs at high dynamics. The receiver tracks with pseudorange lag errors of under 0.06 m when subjected to simulated 50 g turns with 40 g/s peak jerk. Pseudorange errors due to receiver noise alone are approximately 0.6 m rms at a carrier power to noise spectral density ratio of 34 dB-Hz. The tracking threshold SNR is approximately 28 dB-Hz, which provides 12 dB margin relative to the 40 dB-Hz that occurs with minimum specified satellite signal strength, 3.5 dB system noise figure, and 0 dBi antenna gain.

Hurd, W. J.↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗
NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗
NASA NTRS2017

Development of a Numerical Model of Hypervelocity Impact into a Pressurized Composite Overwrapped Pressure Vessel

As the outlook for space exploration becomes more ambitious and spacecraft travel deeper into space than ever before, it is increasingly important that propulsion systems perform reliably within the space environment. The increased reliability compels designers to increase design margin at the expense of system mass, which contrasts with the need to limit vehicle mass to maximize payload. Such are the factors that motivate the integration of high specific strength composite materials in the construction of pressure vessels commonly referred to as composite overwrapped pressure vessels (COPV). The COPV consists of a metallic liner for the inner shell of the COPV that is stiff, negates fluid permeation and serves as the anchor for composite laminates or filaments, but the liner itself cannot contain the stresses from the pressurant it contains. The compo-site-fiber reinforced polymer (CFRP) is wound around the liner using a combination of hoop (circumferential) and helical orientations. Careful consideration of wrap orientation allows the composite to evenly bear structural loading and creates the COPV's characteristic high strength to weight ratio. As the CFRP overwrap carries most of the stresses induced by pressurization, damage to the overwrap can affect mission duration, mission success and potentially cause loss-of-vehicle/loss-of-crew. For this reason, it is critical to establish a fundamental understanding of the mechanisms involved in the failure of a stressed composite such as that of the COPV. One of the greatest external threats to the integrity of a spacecraft's COPV is an impact from the meteoroid and orbital debris environments (MMOD). These impacts, even from submillimeter particles, generate extremely high stress states in the CFRP that can damage numerous fibers. As a result of this possibility, initial assumptions in survivability analysis for some human-rated NASA space-craft have assumed that any alteration of the vessel due to impact is considered a catastrophic failure. This assumption is conservative and made due to lack of knowledge on the level of allow-able damage to the composite overwrap that can be sustained and still allow successful completion of the mission. To quantify the allowable damage level to the composite overwrap involves assessing stress redistribution following damage as well as evaluating possible time-dependent mechanisms involved in the COPV response to an impact event. Limited published work in this subject has shown that COPV can withstand at least some level of damage due to high energy impacts. These observations have been confirmed and expanded upon in recent experimental research performed by NASA. This research has demonstrated that there is not only robustness in a COPV to compensate for CFRP damage, but has also identified two significant failure modes for pressurized COPV. The lowest threshold failure mode involves the perforation of the vessel, and the highest threshold failure mode is the catastrophic rupture. While both of these failure modes mean a loss of the COPV, system robustness affords some tolerance to the venting as opposed to the more catastrophic rupture. As a consequence, it is necessary to understand the conditions that result in the transition between these failure modes. The aforementioned experimental research has been performed in both the unpressurized and pressurized condition to identify the damage level that triggered the failure thresh-old. This COPV test program was sponsored by the NASA Engineering and Safety Center (NESC), and tests were performed at NASA White Sands Test Facility (WSTF). Planning and coordination were provided by NASA JSC Hypervelocity Impact Technology (HVIT) group, and the COPVs were provided by the ISS Program. Unpressurized testing has been conducted at the pressure of the vacuum test chamber, while, the pressurized testing has been conducted at 290 +/- 10 bar (4,200  100 psi) using nitrogen as the pressurizing gas, which corresponds to the design pressure for the target COPV. In this research, spherical aluminum projectiles with varying diameter has been chosen as the impactor. For the unpressurized COPV, the dependence of penetration up to the dependence of hole size in the liner has been obtained as a function of impact conditions. For the pressurized research, the dependence of penetration up to rupture has been obtained as a function of im-pact conditions. Two representative post-test photographs of the failed COPV's from a nor-mal impact into the COPV surface are shown in Fig. 1. These images display the dramatic difference between failure modes, venting (Fig. 1a) and rupture (Fig. 1b). For venting, liner perforation, severed composite fibers/tows and ply delamination are commonly observed damage characteristics of this COPV failure mode. In the case of rupture, the COPV typically experienced a separation of its domed regions and severe break-up of the cylindrical region. Fully understanding the transition from venting to rupture experimentally is costly and potentially unachievable for conditions that cannot be generated in the laboratory. These shortcomings have motivated the performance of three-dimensional numerical simulations to expand the existing experimental database. These simulations have been carried out with the nonlinear-structural-dynamics, analysis-tool, CTH. A typical pressure contour plot from an impact simulation of an entire COPV is shown in Fig. 2. To generate the COPV stress state without initiating a shock wave, the pressure in the simulated COPV is ramped up to the final pressure over a millisecond prior to impact of the projectile with nitrogen gas. Figure 2a shows the system in this initial condition. After one millisecond, a projectile is initiated into the simulation and impacts the COPV. Figure 2b shows the system after this impact. In the figure, the onset of venting is represented as the change in pressure (μbar), red to green, at the perforation site. Also seen in the figure is the eroded projectile that had passed into the COPV vessel with the generated shock wave in the pressurant propagating just ahead of the material. In this paper, pertinent experimental details and the development of the material constitutive models necessary for this work along with the efforts to validate their use are dis-cussed. The simulation results are presented and compared with the NASA experimental observations. While work is on-going from this effort, early observations pertinent to the failure threshold are presented.

Garcia, M. A.↗
NASA NTRS2017

High Speed Thermal Imaging on Ballistic Impact of Triaxially Braided Composites

Ballistic impact experiments were performed on triaxially braided polymer matrix composites to study the heat generated in the material due to projectile velocity and penetration damage. Quantifying the heat generation phenomenon is crucial for attaining a better understanding of composite behavior and failure under impact loading. The knowledge gained can also be used to improve physics-based models which can numerically simulate impact of composites. Triaxially braided (0/+60/-60) composite panels were manufactured with T700S standard modulus carbon fiber and two epoxy resins. The PR520 (toughened) and 3502 (untoughened) resin systems were used to make different panels to study the effects of resin properties on temperature rise. Ballistic impact tests were conducted on these composite panels using a gas gun, and different projectile velocities were applied to study the effect on the temperature results. Temperature contours were obtained from the rear surface of the panel during the test through a high speed, infrared (IR) thermal imaging system. The contours show that high temperatures were locally generated and more pronounced along the axial tows for the T700S/PR520 composite specimens; whereas, tests performed on T700S/3502 composite panels using similar impact velocities demonstrated a widespread area of lower temperature rises. Nondestructive, ultrasonic C-scan analyses were performed to observe and verify the failure patterns in the impacted panels. Overall, the impact experimentation showed temperatures exceeding 525 K (485degF) in both composites which is well above the respective glass transition temperatures for the polymer constituents. This expresses the need for further high strain rate testing and measurement of the temperature and deformation fields to fully understand the complex behavior and failure of the material in order to improve the confidence in designing aerospace components with these materials.

thermal imaging↗
DOE OSTI2021

WOMBAT (Windfarm Operations and Maintenance cost-Benefit Analysis Tool) WISDEM® [SWR-21-68]

The windfarm operations and maintenance cost-benefit analysis tool (WOMBAT) is software to simulate the operations and maintenance phase of either onshore or offshore windfarms. WOMBAT is a medium-fidelity tool that is designed to model turbine, cable, and substation failures at the subassembly level, but is flexible enough that individual component or asset-level modeling can also be performed. WOMBAT enables users to model technological improvements, the use of a wide range of service equipment from drone repairs and remote resets to heavy lift vessels and crawler cranes, to analyze cost trends with varying trade-offs in conjunction with energy production, availability and a growing range of metrics. This library provides a tool to simulate the operation and maintenance phase (O&M) of distributed, land-based, and offshore windfarms using a discrete event simulation framework. WOMBAT is written around the SimPy discrete event simulation framework. Additionally, this is supported using a flexible and modular object-oriented code base, which enables the modeling of arbitrarily large (or small) windfarms with as many or as few failure and maintenance tasks that can be encoded. Please note that this is still heavily under development, so you may find some functionality to be incomplete at the current moment, but rest assured the functionality is expanding. With that said, it would be greatly appreciated for issues or PRs to be submitted for any improvements at all, from fixing typos (guaranteed to be a few) to features to testing. Also available at: https://pypi.org/project/wombat/

Cooperman, Aubryn↗
NASA NTRS2022

Predicting the Penetration of a Shielded TPS Tile

All spacecraft are subject to the possibility of high-speed particle impacts during their mission life. In low earth orbit, those impacts could be the result of collisions with pieces of orbital debris or with meteorites. Beyond LEO, and especially beyond GEO, those impacts will likely be caused by meteorites. Such high-speed impacts on spacecraft surfaces create debris clouds that travel towards and eventually impact other downstream spacecraft components. In addition to the impulsive load that such debris clouds would impart to the spacecraft elements with which they subsequently collide, the largest fragment in these debris clouds poses a significant threat on its own to those spacecraft elements. In order to be able to assess the severity of the threat posed by such a fragment, it is important to be able to predict the extent of damage sustained by the impacted spacecraft element. In this paper, we present a new process for determining the penetration depth in a shielded TPS tile system. The requirements for the process were that it be applicable across a full spectrum of impact low velocities (i.e. from ~ 7 km/s to above 50 km/s), that it be applicable over a wide range of projectile materials densities (i.e. ranging from water to steel, for example), and that it be applicable over a wide range of trajectory obliquities (i.e. not just normal impacts). Furthermore, the process developed should be sufficiently flexible so that it can be used for shield design parameter and / or system configuration trade studies (i.e. be written in terms of as many shield and TPS tile system parameters as possible). Following the initial impact of the projectile on the outer wall of the dual-wall shield, a debris cloud (i.e. the primary debris cloud) is created that travels towards and impacts the inner wall of the shield. Within this debris cloud is a combination of solid, liquid, and vaporized material, depending on the impact velocity and the impedance mismatch of the projectile and outer wall materials. The impact of this debris cloud on the inner wall creates another debris cloud (i.e. the secondary debris cloud) that then travels towards and impacts the TPS tile. Within the particulate distribution of the material in each of these debris clouds is a so-called largest fragment. The process used to calculate the penetration depth in a shielded TPS tile system developed herein consists of the following three steps: 1. Does the initial projectile perforate the outer wall of the shield? If NO, the TPS tile remains undamaged. If YES, calculate the size and speed of the largest fragment exiting the outer wall of the shield and traveling towards the shield’s inner wall. 2. Does the largest fragment in the debris cloud exiting the outer wall of the shield perforate the inner wall of the shield? If NO, the TPS tile again remains undamaged. If YES, calculate the size and speed of the largest fragment exiting the rear of the shield inner wall. 3. Calculate the penetration depth of the largest fragment in the debris cloud exiting the inner wall into the TPS. Perforation of the inner and outer shield walls is predicted using the Cour-Palais penetration depth equation and an associated failure criterion [1]. The largest debris cloud fragment diameters and associated velocities are calculated using regression equations developed from a mix of hydrocode and experimental diameter and velocity information [2]. Finally, TPS penetration depths are calculated using an empirical penetration depth predictor equation based on high-speed impact test data and are calibrated using a limited number of hydrocode penetration depth predictions for water, nylon, aluminum, Dunite, and steel projectiles The predictions of the penetration depth calculation process are compared against the predictions of more than 60 SPHC hydrocode runs. These comparisons show that the process used to calculate TPS tile penetration depths usually yields values that are within 1 mm of the predictions of hydrocode simulations. If a maximum allowable TPS penetration depth is known, the process developed herein can now be used to develop a failure limit equation that would predict whether or not an impacting particle would penetrate deeply enough into the TPS to result in a “failed” or “not failed” end state. References [1] E.L. Christiansen, “Design and performance equations for advanced meteoroid and debris shields”, International Journal of Impact Engineering, Vol. 14, 1993, pp. 145-156. [2] W.P. Schonberg, “Predicting the size of the largest particle fragment in a debris cloud created by an orbital debris impact and its associated velocity”, 2021 Applied Space Environments Conference, NASA Jet Propulsion Laboratory, Pasadena, California, November, 2021.

Meteroid↗
NASA NTRS2025

HabSim: A Modular Coupled Virtual Testbed for Simulating ExtraTerrestrial Habitat Systems

Extraterrestrial habitats involve a tightly coupled combination of hardware, software, and humans while operating in an unforgiving environment that poses many risks, both anticipated and unanticipated. Traditional approaches with such systems of systems focus on reliability, robustness, and redundancy. These approaches seek to avoid failure rather than reduce overall risk. However, faults are inevitable, and understanding and managing the complex and emergent behavior and cascading events of such a complex system is critical. This study describes the development of HabSim, a computational simulation environment intended to support research to establish the know-how to design and operate resilient and autonomous SmartHabs. HabSim is a modular virtual testbed composed of many of the coupled dynamic systems expected in a typical SmartHab. A heterogeneous set of interconnected physics-based and phenomenological models is used to represent the essential functions of a SmartHab. HabSim further considers disruptions and models damage and repair of certain components. This paper discusses a) system and subsystem requirements of the deep space habitat included in the HabSim platform; b) architectural choices made in response to the requirements; c) technical considerations for developing, verifying, configuring, and executing HabSim; and d) illustrative sample results from a simulation of a representative disruption scenario.

Mohsen Azimi↗
DOE OSTI2025

Development of a Robust Reference Electrode in Aggressive Chemical and Radiation Environments in the Hanford Waste Tanks

The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.

Sykes, Kiana [Savannah River National Laboratory (↗
DOE OSTI2022

A systematic approach to resolve high impedance of multifunctional energy storage composites

Previous work has proposed and characterized the structural and electrical performance of Multifunctional Energy Storage Composite (MESC) structures: structural elements with embedded lithium-ion batteries which were developed by the Structures and Composite Laboratory (SACL) at Stanford University. This work conducts a comprehensive parametric study over the design variables in the manufacturing process to determine the root cause of increased DC impedance compared to commercial pouch cell batteries. Several MESC samples were assembled under various controlled environmental temperatures and pressures. The samples underwent Hybrid Peak Pulse Characterization (HPPC) to approximate the DC impedance as a function of state of charge (SOC). Failure limits of the battery with respect to assembly temperature and pressure were tested, and an MESC manufacturing simulation model was constructed in Abaqus to validate the manufacturing limits. Samples assembled without enough temperature or pressure failed due to a lack of structural integrity or hermetic seal. Conversely, samples which underwent elevated temperatures and pressures resulted in high DC impedance due to nonuniform electrolyte wetting. This work resolves the primary concern of MESC by achieving comparable DC impedance to commercial pouch cells of the same chemistry and testing procedure.

25 ENERGY STORAGE↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2022

A resilient network recovery framework against cascading failures with deep graph learning

Because of the increasing importance and dependencies of infrastructure networks and the potential for massive cascading failures in real-world network systems, maintenance optimization to effectively reduce system performance loss caused by diverse disruptions is of significant interest among researchers and practitioners. In this work, a new recovery framework was developed to rapidly identify important system components for maintenance to improve network resilience against cascading failures. Here this work provides distinct advantages to determine an optimal maintenance priority by combining real-time network structure importance with other maintenance prioritization based on customer preference. This approach adopts structural graph embedding and deep reinforcement learning to extract real-time network topology information (such as minimum vertex cover) to update the maintenance priority during the recovery process. Based on the case studies on synthetic networks and a US airport network, the proposed recovery framework with real-time network topology awareness shows better performance than other maintenance prioritization strategies regarding resilience enhancement. This work improves the understanding of how the changing network structure influences maintenance effects. It also provides insights of the practical usefulness of advanced deep learning on helping optimal maintenance prioritization to effectively reduce the intensity and extent of cascading failures.

42 ENGINEERING↗
NASA NTRS1974

A spacecraft integrated power/attitude control system

A study to determine the viability and application of a system capable of performing the dual function of power storage/generation and attitude control has been conducted. Results from the study indicate that an integrated power/attitude control system (IPACS) can satisfy future mission requirements while providing significant savings in weight, volume, and cost over conventional systems. A failure-mode configuration of an IPACS was applied to a shuttle-launched RAM free-flyer and simulated using make-do hardware linked to a hybrid computer. Data from the simulation runs indicate that control interactions resulting from heavy power demands have minimal effect on system control effectiveness. The system was shown to be capable of meeting the stringent pointing requirements of 1 arc-second while operating under the influence of an orbital disturbance environment and during periods of momentum variations imposed by energy transfer requirements.

Keckler, C. R.↗
NASA NTRS1979

Behavior of composite bolted joints at elevated temperature

Experimental results from an investigation which examines the combined effects of temperature, joint geometry and out-of-plane constraint upon the response of mechanically fastened composite joints are presented. Data are presented for simulated mechanically fastened joint conditions in two laminate configurations fabricated from Hercules AS/3501-6 graphite-epoxy. Strength and failure mode results are presented for the test temperatures of 21 C, 121 C and 177 C and for a range of the geometric parameters W/D and e/D from 3.71 to 7.43 and 1.85 to 3.69, respectively. A hole diameter, D of 5.16 mm was utilized for all tests. Pin bearing tests with out-of-plane constraint were conducted at room temperature only. All elevated temperature data were generated for pin bearing conditions. Ultrasonic C scan inspection of the failed specimens was employed to assess the damage region and to determine failure mode. Comparative data are presented for pin bearing and out-of-plane constraint conditions for the above mentioned joint configurations. The joint under pin loading was modeled by two dimensional finite element methods. Predicted net section strain concentrations were compared with experimental results.

Wilson, D. W.↗