Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Predicting long-term stress relaxation on Alloy 709 using the crystal plasticity finite element method

This report summarizes the results of physics-based, crystal plasticity simulations for the long-term stress relaxation behavior of Alloy 709. The purpose of the study was to provide insight into five key questions related to long-term behavior in high temperatures materials which are difficult or impossible to answer experimentally: (1) is there a threshold stress for long-term relaxation? (2) is there strain threshold for relaxation damage, below which significant damage does not accumulate? (3) does damage continue to accumulate as the material relaxes or will damage accumulation plateau under some loading conditions? (4) does stress relaxation loading inevitably lead to failure? and (5) which, if any, engineering models for relaxation damage accumulation reasonably match the simulation results? The report summarizes the numerical simulations used to address these five questions and provides at least partial answers to each question.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗
DOE OSTI2024

Improving and Automating Building Model Data Exchange

There are many instances throughout a project’s lifecycle where there arises a need for quick and accurate risk assessment of building designs. For example, an unexpected design change during construction may necessitate structural engineers to perform a seismic risk assessment on analytical models of the updated building design using high fidelity structural analysis software, such as ANSYS or Abaqus. However, the efficiency of such workflows often depends upon the interoperability of architectural design software and structural analysis software. When the quality of this interoperability is lacking or even non-existent, the efficiency of virtual engineering workflows is hampered, which increases project costs. A McGraw Hill industry survey of professional users of Building Information Modeling (BIM) technologies found that there is high demand for BIM interoperability for structural analysis, but that the value/difficulty ratio is currently too low for practical use. There have been efforts by the academic community to facilitate model data exchange between the architectural design and structural analysis domains, but such solutions have not been widely adopted by industry, face technical challenges, and oftentimes are limited in applicability for users of various BIM software. Therefore, INL is developing capabilities to improve, automate, and generalize model data exchange between architectural BIM software (e.g., Revit) and structural analysis software (e.g., SAP2000, ANSYS). The goal is to help expedite and automate as much of the pre-processing step for creating analytical models in finite element analysis software as reasonably as possible. Such a "BIM-to-FEA" conversion tool should provide direct benefit to end-users through accuracy, automation, quick turn-around, and wide applicability. To generalize the application of this BIM-to-FEA conversion tool and increase its useability among the many different commercial BIM software currently used by industry, the program is being developed with the concept of openBIM. OpenBIM is the application of non-proprietary, open data standards that allow for BIM model data exchange in a format that is accessible, retainable, and useable for all users. The most widely used open, non-proprietary data exchange format for BIM is the Industry Foundation Classes (IFC) schema. IFC is developed by buildingSMART international and is ISO certified (ISO 16739-1:2018). The BIM-to-FEA conversion tool is being developed for compatibility with typical commercial building designs of steel framed structures. The tool is currently capable of importing architectural BIM data of framed building structures, recognizing and extracting the aspects of the model that are required for structural analysis, adjusting the connectivity of frame members, and finally exporting to an analytical model stored in the IFC format. The exported IFC analytical model can then be imported into various openBIM compliant software, such as SAP2000. Such capabilities have already been tested on commercial software, as shown above, and continue to be improved. Work is underway to test the conversion on various commercial BIM software, develop a user-friendly interface, incorporate the program into the broader DeepLynx data warehouse project being developed by INL, and to eventually open-source the tool for the benefit of the community. Future development of the tool envisions the ability for efficient iterative risk assessment of generative building designs, all within a workflow utilizing open-source tools. One such open-source tool will be MOOSE, an advanced finite element analysis tool developed at INL. The conversion tool will also branch out from typical commercial building designs and will aim to incorporate nuclear construction. The aim will be to convert both structural and non-structural components of nuclear facilities, such as curved concrete containment structures and piping systems, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Environmental influence on cracking and debonding of electrically conductive adhesives

Electrically conductive adhesives (ECAs) are starting to replace metallic solders in recent designs of photovoltaic (PV) modules. This transition represents a significant material change, and a proper understanding of the durability and reliability of the new interconnect needs to be established. This paper presents our continued work on developing a degradation model for ECA interconnects in PV modules. Here, we characterize the fracture mechanics properties of an epoxy-based ECA, for both critical and subcritical, mode I and mode II loading conditions. Emphasis is on the influence of different environmental conditions such as temperature and humidity. We use the Finite Element Method to account for residual stresses, induced by temperature changes and moisture absorption, and correct the apparent fracture toughness. We found that high moisture levels not only can weaken the fracture resistance of the ECA interconnect, but can also promote subcritical debonding at significantly lower driving forces than in dry environments.

42 ENGINEERING↗
DOE OSTI2022

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗
DOE OSTI2024

R&D to Ensure a Scientific Basis for Qualification Tests and Standards (Final Report)

Project return on investment in a photovoltaic (PV) system depends increasingly on maintaining high energy yields, and the system lifetime is a major factor in levelized cost of electricity (LCOE). Thus, the rate of PV deployment and the success of these assets depends upon reliable long-term power generation. The overarching objective of this program is to improve photovoltaic (PV) module reliability via development of tests and standards. Where reliability problems or risk are discovered, we can design tests to ensure that these liabilities don't affect future generations of products. Customers can use these tests to understand which products are susceptible to certain degradation mechanisms, and manufacturers can use the tests to design unwanted characteristics out of their products. The work under this program identifies PV reliability needs, performs characterization that provides scientific understanding of targeted degradation mechanisms, and translates those data into practical and predictive test protocols and standards. Major accomplishments include: A model for polarization-type potential induced degradation (PID-p) was developed and validated against experimental data. NREL is currently leading a new edition of IEC 62804-1 for PID detection. PID-p can cause large losses in current and voltage for some module designs on cloudy days. Finite element modeling (FEM) and experiment was used to determine when cells crack in a module. It was shown that cells in landscape orientation are much more likely to crack than those on portrait orientation. Shortly thereafter, the first products with portrait-oriented cells were introduced. Studies of how to test for light and elevated temperature degradation (LeTID) culminated with the publication of IEC TS 63342. Software to predict the progression of LeTID was developed, validated, and made publicly available. Field validated tests and international standards for durability of PV module coatings abrasion, backsheets, and encapsulants were developed. Examples are IEC 62788-1-1, IEC 62788-2 ED2, IEC TS 62788-7-2, IEC 62788-7-3 ED1, IEC 63209-2. NREL led the development a high-temperature testing technical specification, and published guidelines that enable installers to determine whether higher-temperature testing is needed, simply based on location and mounting configuration. In a number of our case studies, variations in the bills of materials or workmanship have been associated with variations in reliability. These observations emphasize the importance of quality assurance to reliability. A framework for criticality (i.e. Pareto) analysis was developed and published. The framework helps us and other researchers determine what problems should be addressed for reliability research to have the biggest industry impact. NREL continues to participate actively in international standards development and stakeholder engagement activities, including organizing an annual PV Reliability Workshop. These activities are important for ensuring we address issues that are relevant and timely, and that we convey our results to those who may benefit.

14 SOLAR ENERGY↗
DOE OSTI2021

The Role of Interface in Additively Manufactured Interpenetrating Composites

Additively Manufactured Interpenetrating Composites (AMIPCs) are a relatively new metal-metal chain composite in development for use in high energy absorption systems. In this system, reinforcing phase of additively manufactured continuous lattice configurations 316L austenitic stainless-steel is in melt infiltrated with a matrix phase of A356 aluminum-silicon casting alloy. Measurements and observations of this material system have shown that weakly bonded or open/porous interface between the reinforcement and matrix phases exhibits dramatically different mechanical properties of AMIPCs, which is not currently well understood. In this work, Finite Element Models (FEM) are used to model the effects of interfaces between the composite phases. Mechanical tensile tests measurements of various composite volume fractions and varying degrees of casting infiltration are also examined and used to show consistency with the FEM results. The outcome provides insight into material design criteria and performance predictions for new hybrid material systems with exceptional damage tolerance.

Allen, Jason↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2021

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Cladding Profilometry Analysis of Experimental Breeder Reactor-II Metallic Fuel Pins with HT9, D9, and SS316 Cladding

BISON finite element method fuel performance simulations were conducted using an existing automated process that couples the Fuels Irradiation & Physics Database (FIPD) and the Integral Fast Reactor Materials Information System database by writing input files and comparing the BISON output to post-irradiation fuel pin profilometry measurements contained within the databases. The importance of this work is to demonstrate the ability to benchmark fuel performance metallic fuel models within BISON using Experimental Breeder Reactor-II fuel pin data for a number of similar pins, while building off previous modeling efforts. Changes to the generic BISON input file include implementing pin specific axial power and flux profiles, pin specific fluences, frictional contact, and irradiation-induced volumetric swelling models for cladding. A statistical analysis of irradiation-induced volumetric swelling models for HT9, D9, and SS316 was performed for experiments X421/X421A, X441/X441A, and X486. Between these three experiments, there were 174 post-irradiation examination (PIE) profilometries used for validating the swelling models presented using a standard error of the estimate (SEE) method. Implementation of the volumetric swelling models for D9 and SS316 claddings was found to have a significant impact on the BISON profilometry simulated, where HT9 clad pins had an insignificant change due to low fluence values. BISON profilometry simulated for HT9, D9, and SS316 fuel pins agreed with PIE profilometry measurements, with assembly SEE values being 4.4 × 10−3 for X421A, 2.0 × 10−3 for X441A, and 2.8 × 10−3 for X486. D9 clad pins in X421/X421A had the highest SEE values, which is due to the BISON simulated profilometry being shifted axially. While this work accomplished its purpose to demonstrate the modeling of multiple fuel pins from the databases to help validate models, the results suggest that the continued development of metallic fuel models is necessary for qualifying new metallic fuel systems to better capture some physical performance phenomena, such as the hot pressing of U-Pu-Zr and the fuel cladding chemical interaction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The 3D printing and modeling of functionally graded Kelvin foams for controlling crushing performance

Mechanical impact protection is an important consideration in many applications, ranging from product transportation to sports. Cellular materials are typically used due to their desirable energy absorption properties and light weight. However, their large deformation and rate dependent responses (especially of polymer foams) are challenging to consider in design. Additionally, the use of foams with uniform properties, such as uniform density and uniform stiffness, often restricts the designed foams to only be suitable for a narrow range of mechanical impact conditions whereas real applications commonly face unpredictable situations. 3D printing offers fabrication flexibility and thus opens the door to create foams with tailored properties. In this work, we investigate the feasibility of using 3D printing for functionally graded foams (FGFs) that are optimal over a broad range of mechanical environments. The foams are fabricated by the recently developed grayscale digital light processing (g-DLP) method which can print parts with locally designed properties. These foams are tested under drop test conditions and with slower displacement control. We also model the large deformation behavior of FGFs using finite element analysis in which we account for the different viscoelastic behaviors of the distinct grayscale regions. We then use the model to examine the impact mitigation capabilities of FGFs in different loading scenarios. Finally, we show how FGFs can be used to satisfy real-world design goals using the case study of a motorcycle helmet. In contrast to prior work, we investigate continuous, functionally graded foams of a single density that differ in their viscoelastic responses. This work provides further insight into the benefits of viscoelastic properties and modulus graded foams and presents a manufacturing approach that can be used to produce the next generation of flexible lattice foams as mechanical absorbers.

3D printing↗
DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗
DOE OSTI2023

2023 Annual Report of Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve

DOE maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns at the U.S. Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR’s ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. This report covers impacts on drawdown availability due to SPR operations during Calendar Year 2022. A cavern has an available drawdown if, after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of available drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent red flags for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a pillar-to-diameter ratio of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for oil sales, purchases, or exchanges authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern’s drawdown capacity may be continued. A methodology for assessing and tracking the available drawdowns for each cavern is reiterated. This report is the latest in a series of annual reports, and it includes the baseline available drawdowns for each cavern, and the most recent assessment of the evolution of drawdown expenditures. A total of 222 million barrels of oil were released in calendar-year 2022. A nearly-equal amount of raw water was injected, resulting in an estimated 34 million barrels of cavern leaching. Twenty caverns have now expended a full drawdown. Cavern BC 18 has expended all its baseline available drawdowns, and has no drawdowns remaining. Cavern BM 103 has expended one of its two baseline drawdowns, and is now a single-drawdown cavern. All other caverns with an expenditure went from at-least-5 to at-least-4 remaining drawdowns.

02 PETROLEUM↗
DOE OSTI2020

Reification of latent microstructures: On supervised unsupervised and semi-supervised deep learning applications for microstructures in materials informatics

Machine learning (ML), including deep learning (DL), has become increasingly popular in the last few years due to its continually outstanding performance. In this context, we apply machine learning techniques to "learn" the microstructure using both supervised and unsupervised DL techniques. In particular, we focus (1) on the localization problem bridging (micro)structure (localized) property using supervised DL and (2) on the microstructure reconstruction problem in latent space using unsupervised DL. The goal of supervised and semi-supervised DL is to replace crystal plasticity finite element model (CPFEM) that maps from (micro)structure (localized) property, and implicitly the (micro)structure (homogenized) property relationships, while the goal of unsupervised DL is (1) to represent high-dimensional microstructure images in a non-linear low-dimensional manifold, and (2) to discover a way to interpolate microstructures via latent space associating with latent microstructure variables. At the heart of this report is the applications of several common DL architectures, including convolutional neural networks (CNN), autoencoder (AE), and generative adversarial network (GAN), to multiple microstructure datasets, and the quest of neural architecture search for optimal DL architectures.

36 MATERIALS SCIENCE↗