Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary finite element”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Computationally Efficient AC Resistance Model for Stator Winding With Rectangular Conductors

The resistance of the electric machine’s windings is an important parameter when analyzing the machine’s performance, as conduction losses in the winding are often the largest loss component. However, this resistance changes with frequency due to the skin and proximity effects, making it challenging to accurately estimate these losses. In this article, a computationally efficient, general method of calculating the ac resistance of an electric machine winding with rectangular conductors is presented. The method uses a 1-D magnetoquasi-static model of the stator slot region to develop boundary transfer relations between individual turns. These relations are then used to create an impedance matrix from which the ac resistance can be calculated. This technique is more flexible than previous solutions since it can accommodate different winding configurations such as fractional-pitch windings. It is shown that a fractional-pitch winding has a lower ac resistance than a full-pitch winding. Finally, finite-element analysis simulations and experimental results are used to validate the proposed method.

42 ENGINEERING↗
DOE OSTI2025

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid AI-ML and FE-based Digital Twin Predictive Modeling Framework for a PWR Coolant System Components: Updates on Multi-Time-Series-3D-Location Dependent Usages Factor Prediction

The long-term operation (LTO) of nuclear power plants (NPP) beyond their original design life of 40 years can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition, the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Requiring frequent NDE inspections may lead to frequent NPP shutdowns which can lead to power outages and additional NDE inspection cost-related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing a more digital approach such as using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. The DT framework is based on limited experimental data, Artificial-intelligence (AI)-Machine-Learning (ML) and multiphysics-computationalmechanics such as finite element- (FE) based models. Given the real-time thermal-hydraulic process measurements from several existing plant sensors, the overall goal of the DT framework is to predict the cumulative usages factors or equivalent fatigue lives in real time and at any random 3D location of the components. This includes inaccessible locations such as inside the thickness or location of a component. This prediction can be at thousands to millions of 3D point clouds or locations like conventional FE-based models, but without running an FE model in real time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Macro-micro multiscale modeling to assist the design of HPDC Al castings microstructure and alloys for EV super-large body structures (Phase 1)

Implementation of High Pressure Die Casting (HPDC) Aluminum (Al) body structures for high volume electrified vehicles (EV) to improve electric efficiency remains a key strategy within many original equipment manufacturer (OEM)s. In addition to high strength for safety requirements, superior Self-Piercing Riveting (SPR) performance is demanded for HPDC Al alloys to be compatible with high volume SPR joining. In this work, it is proposed to extend and validate an existing Contractor finite element multiscale macro-micro modeling approach to quantify the influence of the microstructure of HPDC alloys on the fracture strain/displacement under 3-point bend and clinch testing. The success of this work will allow to replace solution treatment stage with low energy consumption heat treatment (HT) processes, or to design new non heat treatable (NHT) HPDC Al alloys to eliminate HT requirements. Ultimately, this project will facilitate the application of HPDC Al alloys for super-large vehicle structures to significantly reduce vehicle weight, and thus improving energy efficiency. The purpose of this project is to extend and validate an existing finite element code, which is based on the Contractor developed macro-micro multi-scale modeling approach, to numerically simulate the three-point bending and clinch test and study the influences of material microstructural characteristics and phase properties on the rivetability. The macro-micro modeling approach begins with a sample scale model and identify the location, which is mostly prone to failure, the deformation history of the boundaries of that location calculated will be used to drive a microstructure-based sub-models where the material microstructure and microscale properties are considered. Using this approach, the wrap-bending failure for two Al alloys are correctly predicted for the first time. This will start with phase I effort of building a framework of macro-micro three point bending test and clinch test of Al10SiMgMn HPDC alloy in the as-cast and T7 heat treated conditions. Those results will then be validated with experimental test results. The phase II effort will involve the utilization of the knowledge learned in phase I to establish the quantitative correlation between the microstructure characteristics and the riveting performance, which will be further used to guide the optimization of HPDC Al alloy microstructure using heat treatment process to achieve sufficient rivetability to join large thin-wall HPDC alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Combining four-point bending and corrosion to map stress-dependent corrosion susceptibility of 316H stainless steel in FLiNaK

Understanding the effect of stress on the corrosion of steels in molten salts is important for material screening and safety assessment of molten salt reactors. We present a method that combines four-point bending with corrosion testing, enabling the mapping of corrosion susceptibility as a function of local stress from a single specimen. Guided by finite element analysis, a four-point bending setup was used to bend a 316H stainless-steel bar under controlled elastic stress during a 100-hour exposure in FLiNaK at 700 °C. Post-exposure characterization using scanning electron microscopy and energy-dispersive X-ray spectroscopy revealed a pronounced correlation between stress and corrosion. Regions under tensile stress exhibited significantly deeper attack depths and greater chromium depletion along grain boundaries, whereas compressive zones showed shallower corrosion cracks. The stress effect is attributed to stress concentration under tensile loading, which promotes chromium dissolution and crack propagation. This technique can be readily applied to other structural alloys, enabling quantitative measurement of corrosion susceptibility as a function of local stress.

316H↗
DOE OSTI2022

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Stress relaxation cracking susceptibility evaluation in 347H stainless steel welds

Abstract Stress relaxation cracking (SRC) has been reported as a dominant failure mechanism for 347H stainless steel welds at elevated service temperatures, occurring in either the heat-affected zone (HAZ) or fusion zone (FZ). In this study, SRC susceptibility of physically simulated HAZ and cross-welded E347-347H stainless steel welds was studied using a four-step accelerated SRC test with a thermomechanical physical simulator. The stress and temperature range studied represents weld-induced residual stress at 150–600 MPa and post-weld heat treatment temperatures between 800 and 1050 °C. For all temperature conditions, the cross-welded samples failed at a lower critical stress threshold than the simulated HAZ. A finite element model of a 50.8-mm-thick single-V groove weld was used to generate the residual stress maps induced by a 40-pass welding procedure, which in combination with the experimental threshold stress map predicted potential SRC failure locations. Postmortem microstructural evaluations were performed to identify contributing characteristics to SRC. It was found that HAZ samples exhibited a combination of creep void development, secondary cracks with intergranular fracture near grain boundary carbides, and eutectic phases. FZ samples showed brittle fracture with cracks propagating along straight, non-tortuous interdendritic grain boundaries.

42 ENGINEERING↗
DOE OSTI2025

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗
DOE OSTI2020

Validation of finite-element models using full-field experimental data: Levelling finite-element analysis data through a digital image correlation engine

Full-field data from digital image correlation (DIC) provide rich information for finite-element analysis (FEA) validation. However, there are several inherent inconsistencies between FEA and DIC data that must be rectified before meaningful, quantitative comparisons can be made, including strain formulations, coordinate systems, data locations, strain calculation algorithms, spatial resolutions and data filtering. As such, in this paper, we investigate two full-field validation approaches: (1) the direct interpolation approach, which addresses the first three inconsistencies by interpolating the quantity of interest from one mesh to the other, and (2) the proposed DIC-levelling approach, which addresses all six inconsistencies simultaneously by processing the FEA data through a stereo-DIC simulator to ‘level’ the FEA data to the DIC data in a regularisation sense. Synthetic ‘experimental’ DIC data were generated based on a reference FEA of an exemplar test specimen. The direct interpolation approach was applied, and significant strain errorswere computed, even though therewas no model form error, because the filtering effect of theDIC enginewas neglected. In contrast, the levelling approach provided accurate validation results, with no strain error when no model form error was present. Next, model form error was purposefully introduced via a mismatch of boundary conditions. With the direct interpolation approach, the mismatch in boundary conditions was completely obfuscated, while with the levelling approach, it was clearly observed. Finally, the ‘experimental’ DIC datawere purposefully misaligned slightly fromthe FEA data. Both validation techniques suffered from the misalignment, thus motivating continued efforts to develop a robust alignment process. In conclusion, direct interpolation is insufficient, and the proposed levelling approach is required to ensure that the FEA and the DIC data have the same spatial resolution and data filtering. Only after the FEA data have been ‘levelled’ to the DIC data can meaningful, quantitative error maps be computed.

42 ENGINEERING↗
DOE OSTI2022

Measurement of grain boundary strength of Inconel X-750 superalloy using in-situ micro-tensile testing techniques in FIB/SEM system

Grain boundaries (GBs), known as two-dimensional defects, are omnipresent in polycrystalline metallic alloys and thus influence a wide range of mechanical properties under different environmental conditions like irradiation and corrosion. Therefore, quantifying the strength of individual GBs is critical for understanding the degradation of mechanical properties of materials under different conditions. In this study we developed an efficient approach for the fabrication of micro-tensile specimens with a GB almost perpendicular to the tensile direction, which is expected to advance the development of individual GB tensile testing at micro or nanoscale in a wide scope of materials. Further, an in-situ cantilever micro-tensile testing method was developed and used to quantify the strength of a Σ3 GB in Inconel X-750 with the combination of finite element modeling. The average ultimate tensile strength (UTS) of a non-irradiated Σ3 GB is estimated at around 1.4 GPa, comparable to that of a neutron-irradiated Σ3 GB with a dose of ~1.5 dpa (1.3 GPa). Moreover, the in-situ push-to-pull micro-tensile testing technique developed in this work provides valuable insights into the high-angle GB deformation and fracture behavior. This method generates qualitatively similar ductility behavior before and after neutron irradiation as the bulk material testing. However, the ductility and UTS values obtained from this method are different from bulk measurements due to vastly different specimen dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Response of 11 B enriched ZrB 2 ultra-high temperature ceramic to neutron irradiation at elevated temperatures

ZrB 2 , an ultra-high temperature ceramic (UHTC) is being considered for use in fusion reactor first-wall structures, yet its response to irradiation remains poorly understood. This study employed scanning/transmission electron microscopy (S/TEM), synchrotron X-ray diffraction (XRD), finite element calculations, and thermal property measurements to thoroughly investigate the neutron-irradiation effects on 11 B-enriched ZrB 2 . Neutron irradiations were conducted at 220 °C and 620 °C, with a neutron fluence of 2.2 × 10 25 neutron/m 2 (energy > 0.1 MeV), resulting in 3.9 dpa and 4200 appm He. The study revealed the unusual prevalence of prism loops and a > c anisotropic lattice swelling, likely linked to the low c/a ratio of ZrB 2 , leading to grain boundary microcracking. Reducing the grain sizes was effective in reducing intergranular cracking and macroscopic swelling. The observation of cavities in ZrB 2 irradiated at 620 °C, as opposed to 220 °C, prompts questions about the temperature at which vacancies in ZrB 2 become mobile, and the role of neutron absorption by 10 B in elevating irradiation temperatures. Isotopic enrichment in 11 B proves to be a viable strategy for mitigating helium production in transition-metal diborides, which is a critical consideration for nuclear applications. Irradiation-induced defects reduce the thermal diffusivity and conductivity of ZrB 2 by a factor of 4–9, which has important implications for its role as a plasma-facing material in fusion reactors that drive high heat fluxes through first-wall materials. Here, this comprehensive study lays the foundation for understanding ZrB 2 behavior under neutron irradiation and highlights important phenomena to consider for various material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗
DOE OSTI2026

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗
DOE OSTI2022

Tunable phononic bandgap materials designed via topology optimization

Topology optimization is used to design phononic bandgap materials that are tunable by mechanical deformation. A periodic media is considered, which due to the assumption of length scale separation, allows the dispersion relations to be obtained by analyzing a single unit cell subjected to Floquet–Bloch boundary conditions. A finite macroscopic deformation is applied to the unit cell to affect its geometry and hence dispersion. We tune the dispersion–deformation relation to our liking by solving a topology optimization problem using nonlinear programming. The adjoint method is employed to compute the sensitivities, and the non-differentiability of degenerate eigenvalues is avoided using symmetric polynomials. Several tunable phononic crystal designs are presented. Also, a verification analysis is performed, wherein the optimized design is interpreted and analyzed using a conforming finite element mesh.

42 ENGINEERING↗
DOE OSTI2024

Nonperturbative running of the tensor operator for N f = 3 QCD from the chirally rotated Schrödinger functional

We study the renormalization group (RG) running of the nonsinglet tensor operator, for N f = 3 QCD with Wilson fermions in a mixed action setup, with standard Schrödinger functional (SF) boundary conditions for sea quarks and chirally rotated Schrödinger functional ( χ SF ) boundary conditions for valence quarks. Based on a recursive finite-size scaling technique we compute nonperturbatively the tensor step-scaling function for an energy range between a hadronic scale and an electroweak scale, above which perturbation theory may be safely applied. Our result is expressed as the RG-running factor T RGI / [ T ( μ had ) ] R , where the numerator is the scale independent (renormalization group invariant—RGI) tensor operator and the denominator is its renormalized counterpart at a hadronic scale μ had = 233 ( 8 ) MeV in a given scheme. We determine the step-scaling function in four distinct renormalization schemes. We also compute the renormalization parameters of these schemes at μ had which, combined with the RG-running factor, gives the scheme-independent quantity Z T RGI ( g 0 2 ) in four schemes and for a range of bare gauge couplings in which large volume hadronic matrix element simulations are performed by the CLS consortium in N f = 2 + 1 QCD. All four results are compatible and also agree with a recent determination based on a unitary setup for Wilson quarks with Schrödinger functional boundary conditions [arXiv:2309.04314]. This provides a strong universality test. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
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↗
DOE OSTI2021

Identification of High-Frequency Transverse Acoustic Modes in Multi-Nozzle Can Combustors

Abstract High frequency thermoacoustic instabilities are problematic for lean-premixed gas turbines. Identifying which acoustic mode is being excited is important, in that it provides insight into potential mitigation measures and mechanical stress/life. However, the frequency spacing between modes becomes quite close for high frequency instabilities in a can combustor. This makes it difficult to distinguish between the modes (e.g., the first transverse mode vs. a higher order axial/mixed mode) based upon frequency calculations alone, which inevitably have uncertainties in boundary conditions, temperature profiles, and combustion response. This paper presents a methodology to simultaneously identify the acoustic mode shapes in the axial and azimuthal directions from acoustic pressure measurements. Multiple high temperature pressure transducers, located at distinct axial and azimuthal positions, are flush mounted in the combustor wall. The measured pressure oscillations from each sensor are then used to reconstruct the pressure distributions by using a least squares method in conjunction with a solution of a three dimensional wave equation. In order to validate the methodology, finite element method (FEM) calculations with estimated post-flame temperature is used to provide the candidate frequencies and corresponding mode shapes. The results demonstrate the reconstructed mode shapes and standing/spinning character of transverse waves, as well as the associated frequencies, both of which are consistent with the FEM predictions. Nodal line location was also extracted from the experimental data during the instabilities in the pressure data.

Kim, J.↗