Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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.

699 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High Pressure Synthesis of Rubidium Superhydrides

Through laser-heated diamond anvil cell experiments, we synthesize a series of rubidium superhydrides and explore their properties with synchrotron x-ray powder diffraction and Raman spectroscopy measurements, combined with density functional theory calculations. Upon heating rubidium monohydride embedded in H 2 at a pressure of 18 GPa, we form RbH 9 − I , which is stable upon decompression down to 8.7 GPa, the lowest stability pressure of any known superhydride. At 22 GPa, another polymorph, RbH 9 − II is synthesised at high temperature. Unique to the Rb-H system among binary metal hydrides is that further compression does not promote the formation of polyhydrides with higher hydrogen content. Instead, heating above 87 GPa yields RbH 5 , which exhibits two polymorphs ( RbH 5 − I and RbH 5 − II ). All of the crystal structures comprise a complex network of quasimolecular H 2 units and H − anions, with RbH 5 providing the first experimental evidence of linear H 3 − anions. Published by the American Physical Society 2025

Kuzovnikov, Mikhail A.↗
DOE OSTI2023

High-Temperature Creep Behavior of Thin-Walled FeCrAl Alloy Tubes

This report summarizes measured creep properties for FeCrAl alloy C26M measured using subsized specimens machined from prototypic thin-walled tubes. A representative tube fabricated in a previous large batch production effort was machined into dual-gauge axial specimens using electrical-discharge machining. Creep tests, comprising strain-rate jump tests (SRJ), jump stress (JS) and constant-stress (CS) tests, were conducted at temperatures ranging from 600-900ºC to determine creep mechanisms (elucidated from stress exponents) and activation energies to compare with prior data available in the literature at temperatures up to 650ºC. It was found that the results from the SRJ tests indicated that the highest stress exponent of 5.8 occurred at 600ºC, which decreased to the lowest stress exponent of 2.5 at 700ºC. At higher temperatures, the stress exponents increased slightly to a range of 2.9 to 4.3 between 750ºC and 900ºC. The high stress exponent value of 5.8 suggested dislocation glide creep mechanisms dominated creep deformation at 600ºC since high stress values were measured. The low stress exponent of 2.5 at 700ºC suggested diffusional creep mechanisms start to occur that aid in dislocation cross-slip and climb mechanisms. The slightly higher stress exponents of 2.9 to 3.9 may suggest different contributions of creep deformation mechanisms based on diffusional processes associated with Nabarro-Herring and Coble creep and dislocation creep that involve thermally activated cross-slip and climb mechanisms. These deformation mechanisms could have also been affected by the development of recovery and/or recrystallization processes in the microstructure of the C26M2 specimens. From the SJ and single CS tests, two sets of creep activation energies were calculated. The steady state creep values obtained from the SJ tests showed several abnormalities especially with the data at 850ºC. The comparison of the creep activation energies between the data from the SJ and CS tests showed higher creep rates for the SJ tests than the CS tests, indicating microstructure softening occurred in the C26M2 specimens during the SJ tests that led to the higher creep rates at stresses above the initial stress of 10 MPa. For the CS tests, the results showed steady increases in MCR with increasing temperature at each stress. The best linear fit to the MCR values at each stress was with data obtained at 800ºC, 850ºC and 900ºC. The lowest Q c value of 130.8 KJ/mole was observed with the highest stress of 25 MPa. As stress was lowered, the Q c values increased continuously to the highest value of 188.2 KJ/mole at 10 MPa. This trend suggests that the dominant creep mechanism is diffusional based dislocation climb since the decreasing activation energies with increasing stress promotes easier dislocation climb around obstacles. The comparison between the SJ and CS tests indicates that the best type of creep test performed on FeCrAl alloys such as C26M2 is the single CS test.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure–Activity Relationships for Ethanol Dehydrogenation to Acetaldehyde by Silica-Supported Zinc Oxide Catalysts

Silica-supported ZnO efficiently catalyzes the nonoxidative dehydrogenation of ethanol to acetaldehyde, which is relevant for production of 1,3-butadiene from bioethanol. Characterization with in situ spectroscopies under dehydrated conditions (high sensitivity-low energy ion scattering (HS-LEIS), diffuse reflectance (DR) UV–vis, X-ray absorption spectroscopy (XAS), diffuse reflectance Fourier transform infrared spectroscopy (DRIFTS), inelastic neutron scattering (INS), and UV Raman), and ammonia adsorption probed by temperature-programmed desorption followed by DRIFTS and mass spectrometry (DRIFTS-MS NH 3 -TPD), and DFT calculations revealed that the supported ZnO x phase was present as isolated surface ZnO x sites on SiO 2 , with the vast majority coordinated by two siloxane bonds and one silicon atom with two nonbridging oxygens ((≡SiO) 2 Zn 2+ O 2 Si=), anchored at 4-, 5-, and 6-membered siloxane rings. A minor fraction of surface ZnO x sites possessed Lewis acidity, and even fewer sites possessed a Bro̷nsted acidic Zn(OH) + Si moiety. Ethanol temperature-programmed surface reaction-mass spectrometry (TPSR-MS) with various oxidative or ethanol reaction pretreatments indicated that only sites with Lewis and Bro̷nsted acidic character (Zn(OH) + Si) were active for ethanol dehydrogenation, while the majority surface (≡SiO) 2 Zn 2+ O 2 Si= sites were inactive. Greater heterogeneity among all surface ZnO x sites, as assessed by in situ DR UV–vis spectroscopy, was associated with a greater number of ZnO x sites that were active for ethanol dehydrogenation as well as lower enthalpic barriers for acetaldehyde production among the most active surface ZnO x sites. Turnover frequencies and the apparent activation energy for ethanol dehydrogenation were determined from steady-state kinetics. Together, these findings suggested that anchoring inactive surface (≡SiO) 2 Zn 2+ O 2 Si= sites on the silica support caused a greater number of active surface ZnO x sites to adopt a more strained configuration, promoting ethanol dehydrogenation catalysis. Pretreatments and catalysts that promoted desorption of ethanol during TPSR, taken as a marker of surface dehydroxylation, were associated with an increased number of the most active surface (Zn(OH) + Si) sites. Such findings suggested that inactive surface ZnO x sites were activated for ethanol dehydrogenation by dehydroxylation of the support and/or decreased coordination to hemilabile siloxane ligands.

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

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

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

Loads Response That is Due to Wake Steering on a Pair of Utility-Scale Wind Turbines

The goal of this report is to present the resulting turbine loads of a full-scale wind turbine wake steering field campaign. In this experiment, wake steering controls were applied to an upwind turbine (T2) and the downwind turbine (T3) was instrumented to measure mechanical loads. The test subjects (T2 & T3) were GE 1.5SLE CWE turbines located on a wind farm with a prevailing wind direction coming from the northwest. The data collection strategy involved toggling the wake steering controls of T2 between on and off. This resulted in two databases (baseline and wake steered) with similar turbulence intensities. Valid data was extracted and processed. The analysis involved scaling data to engineering units, applying coordinate transformations where applicable, calculating 10-minutes statistics, and calculating damage equivalent loads (DELs) to assess fatigue. Figures are shown as statistics and DELs binned by wind speed with supplemental scatter plots provided in the appendices. Results between the baseline and wake steered databases were compared. Overall, the binned statistics showed minimal differences in loading between the two cases. However, the DELs of the wake steered case were observed to be consistently smaller than the baseline for most of the turbine components. A potential for future work was recognized based on the findings of this report. Some areas of further study include an analysis for more granular loads sensitivity to different yaw offset angles and rotor wake overlap, an in-depth comparison of loads between the upwind and downwind turbines, and further fatigue analysis to quantify differences observed in this experiment.

17 WIND ENERGY↗
DOE OSTI2025

Helium Release During Fracture and Granular Fragmentation of Rocks

Geogenic Helium-4 ( 4 He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of 4 He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total 4 He released during fragmentation of the rocks. We observed that the highest 4 He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total 4 He released during crushing. The range of in-place 4 He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped 4 He was not released. However, the 4 He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that 4 He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for 4 He enrichment in a lithology are met.

Deformation signals↗
DOE OSTI2023

Friction Experiment Data Analysis

The Spent Fuel Waste Disposition (SFWD) program under the U.S. Department of Energy (DOE) is planning a seismic shake table test of full-scale dry storage systems of spent nuclear fuel (SNF) to close the gap related to the seismic loads on the fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly during representative earthquakes. A concrete layer will be installed on the shake table before the test to simulate conditions representative of an ISFSI pad. In the shake table tests with the vertical cask, the cask will be free-standing because this is representative of all, except two, ISFSIs in the U.S. with vertical dry storage casks. The static and dynamic friction coefficients between the steel bottom of the cask and the concrete layer on the shake table are important parameters that will affect cask behavior during the test. These parameters must be known for the pre- and post-test modelling, data analysis, and model validation. The friction experiment was performed at the Engineering Department of the University of New Mexico (UNM) to determine the friction coefficients between a steel plate with the same finish as the bottom of the vertical cask manufactured for the test and different concrete surfaces. In this experiment the steel plate was fixed and the concrete sample was pulled over the plate with a constant displacement rate using an MTS machine. This allowed for collecting continuous horizontal force data over the length of the steel plate. Four displacement rates and three vertical loads were used. The tests were performed with four concrete blocks with different degrees of the surface roughness - light sandblast, light to medium sandblast, medium bush hammer, and heavy sandblast. The total number of tests was 48. The data were used to calculate static and dynamic friction coefficients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Spherical Congeners of Polyaromatic Compounds Approaching C 20 - and C 60 -Fullerene-Type Structures

A series of three symmetric, hollow spherical, and shape-persistent molecular organic cages analogous to C 20 and C 60 were examined by computational modeling, analyzing structural elements, strain indicators, and physical properties relevant for potential applications. The compounds are covalent aromatic cages based on 1,3,5-substituted benzene nodes linked by paraphenylene or para-pyrenylene-connectors, with diameters varying from 2.3 to 4.2 nm. The apertures in the cage interior are varied by virtue of the cage type (C 20 - or C 60 -type cage) and the linear connectors placed between the C 6 H 3 -units. NBO and MESP analyses indicate the presence of electrophilic and nucleophilic sites in the molecular skeleton. In the cages with the phenylene-connectors, the HOMO−LUMO gaps are close to 4.0 eV. In the cage coated with an enlarged polyaromatic spacer (pyrene-unit), the gap is reduced by approximately 0.4 eV.

Aromatic compounds↗
DOE OSTI2025

Interplay Between Metastability and Mechanically Induced Structural Instability in CsPbBr 3 Photovoltaic Perovskite

An atomic-level understanding of the underlying structural metastability is still absent in halide perovskite photovoltaic systems. Focusing on the model material CsPbBr 3 , the impact of mechanically induced atomic structure alteration is elucidated through structural modeling and X-ray diffraction measurements. Sudden Cs–Br bond breaking drives the system metastability, where the first-order transition arises from cation-halide bond strain relief by severing the corner connectivity of the PbBr 6 units. The pressure–volume and entropy terms govern the Gibbs free-energy landscape in halide perovskites. In conclusion, metastability is revealed as a critical factor limiting the performance of lead halide perovskites under extreme but natural environments.

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

Final Technical Report for DE-FG02-86ER45268

This work addresses several key issues for understanding the unusual superconductivity in IBS, already shown to be technologically useful materials. 1a). Are reactive gasses like fluorine, chlorine, or water vapor behaving simply as electron donors to cause the superconductivity (at T conset =21 K) in undoped IBS like BaFe 2 As 2 ? 1b). Does exposure to fluorine increase T c in electron-doped IBS materials that are underdoped, e. g. Ba(Fe 0.96 Co 0.04 ) 2 As 2 ? 1c). Can fluorine as a reactive electron-donor, with its relatively greater ease of introduction into the lattice, be used as a rapid tool to search for superconductivity in new materials under investigation? Method of investigation: Using electron probe microanalysis, x-ray photoemission spectroscopy, and dynamical mean field theory calculations, we found that fluorine substitutes for As in the BaFe 2 As 2 lattice, and should therefore behave as an electron donor. We will F-dope Ba(Fe 0.96 Co 0.04 ) 2 As 2 and several new compounds under investigation for superconductivity upon doping to determine fluorine’s effect. Potential impact: If – like for fluorine - a brief exposure to hydrogen indeed causes superconductivity to a depth of 10 microns, this technique could be useful in creating long lengths of superconducting tape without the chemical aggressiveness of fluorine. Using the F ion, which substitutes for As - away from the Fe planes, in underdoped IBS like Ba(Fe 0.96 Co 0.04 ) 2 As 2 may form optimally doped IBS with less impurity scattering of the superconducting Fe 3d electrons resulting in improved superconducting properties. Another potential impact is the speed of checking for superconductivity in new materials under investigation for superconductivity upon doping. 2. What are the effects of Be substitution into IBS, beginning with FeSe and its derivatives (FeSe 1-x S x and FeSe 1-x Te x ) and BaFe 2-x Co x As 2 ? Method of investigation: We recently succeeded in adding Be to the FeSe lattice and found, surprisingly, a factor of four improvement of the low temperature residual resistivity ratio, consistent with improved electronic ordering. Using single crystal x-ray diffraction and x-ray dispersive spectroscopy, determine the concentration and location of the Be dopant. Using resistivity, determine if there is improvement in the residual resistivity achieved by Be doping in the FeSe derivatives (using AlCl 3 :KCl flux growth) and Co-doped BaFe 2 As 2 (using self-flux single crystal growth). Measuring resistivity under strain, determine the changes of the nematic susceptibility at the structural ordering transition causes by the increase in electronic order caused by Be doping in FeSe. Potential impact: Further understanding of IBS, including FeSe which - under physical pressure - has a T c increase to 40 K and -prepared as a monolayer - has a T c reported to be over 100 K. 3. In collaboration, prepare light atom compounds involving Be and look for high pressure-induced superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Quantifying the dislocation content of atomically resolved grain boundary line defects using the Nye tensor

The Nye tensor, which quantifies the density of Burgers vector for a given dislocation line direction, can be effectively used to characterize dislocation content in bulk crystals from atomic-resolution transmission electron microscopy images. The Nye tensor can be calculated from these images, in part because the reference state is simply defined by the lattice of the perfect crystal. The application of the Nye tensor to interfacial line defects, for which the natural reference state is the dichromatic pattern of the two grains in their reference orientation, poses additional challenges. In this work, we present a method that employs the Nye tensor to characterize the edge dislocation content of line defects at grain boundaries from atomic-resolution images. This approach enables us to rapidly characterize all edge dislocation content along a grain boundary. Additionally, the Nye tensor provides information about line defect core structure. Finally, we demonstrate this method on two exemplar defects: a twin boundary disconnection and a facet junction in face-centered cubic Au.

Crystallographic defects↗
DOE OSTI2025

Magnetic semiconductors from functionalized Cr 2 ⁢C MXenes

Here, we report an ab initio investigation of functionalized and 3⁢𝑑-electron-doped Cr 2 ⁢C MXenes. Upon functionalization, the Cr 2⁢ C becomes chemically, dynamically, and mechanically stable, and it exhibits magnetic semiconducting behavior. Cr 2 ⁢CF 2 stands out as a wide band gap semiconductor, possessing superexchange interaction mediated by F atoms within the layer; however, the applied strain transforms it from an indirect to a direct band gap semiconductor. Strong spin-phonon coupling found in Cr 2 ⁢CH 2 is supported by the distorted Cr spin density due to the hydrogen environment. Two magnon branches, associated with two sublattice spins, are found in the ferromagnetic Cr 2 ⁢CO 2 and antiferromagnetic Cr 2 ⁢CF 2 . Depending on the types of 3⁢𝑑-electron dopants and functionalization, Cr 2 ⁢C MXenes (except for Cr 2 ⁢CO 2 ) change from an indirect band gap magnetic semiconductor to different states of electronic and magnetic matter, including an exotic direct band gap magnetic bipolar semiconductor. In addition, we reveal a band inversion between the two highest valence bands in the Fe-doped Cr 2 ⁢CCl 2 .

2D materials↗
DOE OSTI2021

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗