Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deformation estimation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the correlation between plastic strain and misorientation in polycrystalline body-centered-cubic microstructures with an emphasis on the grain size, loading history, and crystallographic orientation

In this work, the correlation between plastic strain and crystallographic misorientation, grain size, grain orientation, distance from grain boundary, and loading history were investigated experimentally and numerically for body-centered-cubic (BCC) polycrystalline microstructures based on a misorientation deviation (MD) approach. Nine monotonic tensile experiments were performed on two BCC titanium alloys inside a scanning electron microscope (SEM). The influence of reference orientation was explored both at the grain scale and at the mesoscale using electron backscattered diffraction (EBSD). The correlation between global plastic strain and the MD was quantified. The tendency for orientation change was quantified as a function of plastic strain and grain orientation for three crystallographic orientations (i.e., [100], [110], and [111]) with respect to tensile direction. The subpopulation of small grains exhibited a lower level of misorientation dispersion compared with larger grains, and this discrepancy became more pronounced at higher strains. An empirical equation was proposed to estimate the level of misorientation dispersion for individual grains as a function of grain size and global plastic strain level. Two interrupted in-situ SEM experiments were performed, and this resulted in a significantly increased misorientation compared with uninterrupted tests performed to similar plastic strain levels.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intern Poster

Digital Image Correlation (DIC) is an optical technique that combines image registration and tracking methods for accurate two-dimensional and three-dimensional changes in images. DIC software can be used to track the contour, deformation, and strain of a sample. In the Advanced Test Reactor (ATR) at INL (Idaho National Laboratory) there exists a small working window of samples that can become irradiated. Hundreds of graphite disks called piggybacks have undergone this irradiation as part of the Advanced Reactor Technologies (ART) program. After irradiation, it is desirable to understand the change in tensile strength as a function of dose. Due to the limited space in the ATR, typical dog bone tensile tests reduce the number of graphite samples from hundreds to tens. However, there does exist an ASTM standard, D8289, which uses disc compression of graphite to estimate the tensile strength of the specimen with the Brazilian Disk test fixture. While only used as an estimate, which is typically off by a third, it is the purpose of this study to identify how to amend D8289 to remove the word "estimate" with the help of DIC.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural, Criticality, and Radiation Dose Calculations to Support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m-long and 45.7-cm-diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at INL’s CPP-603 facility. The ISS includes a 316L stainless-steel basket that holds 12 PB2 SNF rods. The PB2 basket rests on top of six individual, 78.7-cm-long, 316L stainless-steel columns that are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister above the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture Version 6.2, including ENDF/B-V continuous-energy cross-section libraries, considering intact SNF in a single storage overpack or two multistorage overpack configurations and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and a 27-group neutron-source spectrum were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of 1 m from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Furthermore, this study provides confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Structural, Criticality, and Radiation Dose Calculations to support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m long and 45.7-cm diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at the CPP-603 facility of the INL. The ISS includes a 316L stainless steel basket that holds twelve PB2 SNF rods. Further, six individual, 78.7-cm long, 316L stainless steel columns are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister on top of the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture version 6.2 including ENDF/B-V continuous energy cross-section libraries, considering intact SNF in a single storage overpack or two multi-storage overpack configurations, and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and 27-group neutron-source spectra were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Further, they increase confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Alpha-Quartz Plastic Strength Investigation Via Diffraction Experiments on Novaculite Using A D-Dia and Elastic Plastic Self- Consistent Interpretation [Thesis}

X-Ray Tomography. Porosity is commonly measured using mercury injection (MI) or water immersion porosimeter (WIP). Both MI and WIP utilize pressure to fill open voids with mercury or water respectively. A measure of the volume change of the sample is then used to estimate the percentage of open voids. However, for the purpose of rheological study it is important to get an accurate measure of open and closed voids. X-ray tomography has been selected as it offers a three-dimensional view into the sample that can quantify all pores limited only by the voxel size which is in the micron range for this study. Radiographs for tomography were collected at the Material Science and Technology Division at Los Alamos National Lab using the Carl Zeiss Xradio 520 instrument and Scout-and-Scan version 16.1 operating software. Two samples were imaged, the starting material and the deformed sample SiO2_65. 3001 radiographs were taken of the starting material with a 6 second exposure time using a 4x objective lens. The x-ray beam was set to 60 kilovoltage peak (kVp) and 5 watts. 1901 radiographs were taken of SiO2_65 with a 25 second exposure time using a 10x objective lens. The x-ray beam was set to 80 kVp and 7 watts. Radiograph files were analyzed by Brian Patterson using Avizo. Void and inclusion volumes were output by voxel sized (1.03 μm) slices, used to calculate a total percentage volume for the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗
DOE OSTI2026

Statistical Correlation of Heliostat Pointing Deviation With Wind

This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.

17 WIND ENERGY↗
DOE OSTI2021

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λ l , for several high-order multipoles (l = 5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M ⊙ NS, R 1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, R M . We find that there is an optimal value of M for every M such that the uncertainty in the estimate of R M is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effect of grain dissolution on sloping ground

The static and dynamic stability of natural or constructed slopes can be affected by dissolution or dissolution-like phenomena. Their underlying mechanisms, however, remain unclear. New experimental results and discrete element simulations provide particle-level and macroscale information on the consequences of mineral dissolution on slope behavior. At the microscale, load-carrying grain arches develop around dissolving particles, the porosity increases, and contact force chains evolve to form a honeycomb topology. At the macroscale, while vertical settlements are the prevailing deformation pattern, lateral granular movements that create mass wasting are prominent in sloping ground, even under the quasi-static granular loss. Horizontal grain displacement is maximum at the surface and decreases linearly with the distance from the slope surface to become zero at the bottom boundaries, much like vertical granular displacement along the depth. Sediments with smaller friction angles and steeper slopes experience greater displacement, both vertically and horizontally. Slopes become flatter after dissolution, with the reduction in slope angle directly related to the loss in ground elevation, ΔH/H o . Yet, because of the porous fabric that results from dissolution, vertical shortening is less than the upper bound, estimated from the loss in the solid mass fraction, ΔH/H o ≈SF. Under water-saturated conditions, the post-dissolution fabric may lead to sudden undrained shear and slope slide.

42 ENGINEERING↗
DOE OSTI2026

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗
DOE OSTI2022

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanical behavior analysis using small-size graphite disc compression testing with digital image correlation methods

The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Safety issues of defective lithium-ion batteries: identification and risk evaluation

Lithium-ion batteries inevitably suffer minor damage or defects caused by external mechanical abusive loading, e.g. , penetration, deformation, and scratch without triggering a hard/major short circuit. The replacement of cells becomes a dilemma if the safety risk of the defective batteries remains unknown. Herein, by introducing a representative defect form, i.e. , screw indentation, we demonstrate the safety characteristics of defective batteries. We prove that defective batteries have a significantly increased thermal risk and deteriorated mechanical integrity, but can go undetected due to prompt voltage recovery and insignificant local temperature increase. We discover that the voltage curve within the first few cycles contains sufficient information to identify defective batteries from otherwise good ones and propose methodologies to monitor the cells. Capacity loss and current leakage are two characteristics that can be estimated using the voltage curve. According to the defect size and position, the capacity loss could be 1 to 10 2 mA h and the leakage current could be 5–50 mA. Results remove the barriers for defective battery safety risk evaluation, enabling identification, monitoring, and early warning of minor damaged batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Unraveling the Mechanics: Computational Modeling of Residual Stresses in U10Mo Fuel with Aluminum Cladding

LEU (Low Enriched Uranium) plate-type fuel elements consist of a high-density, low-enrichment U–Mo alloy-based fuel foil encapsulated in an aluminum cladding and is fabricated through the Hot Isostatic Pressing (HIP) technique. Understanding the mechanism of possible failure modes of nuclear fuel is critical to mitigate potential consequences. One of the major contributing factors in various failure modes is stress and it greatly affects the integrity under temperature, pressure, and irradiation. These stresses could originate from various sources, including manufacturing, and operation, and can lead to deformation, cracking, or even complete failure of fuel elements. This study focuses on the development of a computational model that accurately predicts the residual stresses generated in the U10Mo fuel during fabrication. It has been observed that cladding creep has a substantial impact on the residual stresses in the U10Mo fuel post the HIP fabrication process. Furthermore, during the HIP bonding process the fuel plate system is heated to a temperature of 560 oC and as a result, the aluminum cladding transitions from (Al 6061-T6 to -O). This presents a challenge in capturing the change in material properties accurately. Therefore, a temperature dependent creep model such as hyperbolic sine creep model is considered to estimate the creep properties of the Aluminum cladding. The proposed calibrated creep model accurately predicts the residual stresses in the U10Mo fuel foil and agrees well with the experiments.

42 ENGINEERING↗
DOE OSTI2020

Transport of coexisting Ni-Cu sulfide liquid and silicate melt in partially molten peridotite

Transport of coexisting sulfide and silicate melts in partially molten peridotite contributes to the redistribution of chalcophile elements within the upper mantle as well as the genesis of magmatic Ni-Cu sulfide deposits, but has not been investigated systematically. Using laboratory experiments, theoretical calculations, and X-ray synchrotron microtomography, this study documents the topology and considers controls on the extraction of two immiscible liquids during partial melting of mantle peridotite. Under hydrostatic conditions, the measured dihedral angle at silicate melt-mineral-mineral junctions is 13.7–21.3°. Silicate melt is distributed along grain edges forming incompletely interconnected melt channels that disconnected by some dead ends at melt fractions ~7–9 vol%. Application of theoretically predicted permeability (k ~10 -14 -10 -16 m 2 ) permits estimation of the extraction velocity of silicate melt of 0.7-11.1 μm/day within a single interconnected melt channel. In the absence of silicate melt, isolated sulfide droplets (3.77 vol%) show a sulfide-olivine-olivine dihedral angle of 91.5–101.3°. However, in the presence of silicate melt, sulfide droplets (average size ~2.53±2.14μm, 1σ) are partially surrounded by silicate melt and stranded in triple junctions or melt pockets due to the limitation of the smallest dimension (0.3μm) of melt channels. Thus, the extraction of sulfide liquid is highly restricted by these dead ends and the smallest dimension of melt channels during porous flow of silicate melt. In contrast, during large-strain shear deformation (strain ~1.6–2.5), initially stranded sulfide droplets were elongated and extracted with silicate melt into liquid-rich sheets with a length of several hundred microns, constantly oriented at 14.3±4.5°to the shear plane and antithetic to the shear direction. The angle is lower than that (18-30°) of those sheets containing sulfide liquid only, indicating that silicate melt dominates liquid-rich sheets. Driven by stress, silicate melt-dominated liquid-rich sheets open the appropriately oriented grain boundaries between silicate minerals, thus providing an efficient pathway for the extraction of sulfide liquid during deformation. When such sheet-like channels remain open, sulfide droplets (> millimeter-scale) can be potentially mobile through high strain domains of the upper mantle, contributing to the addition of chalcophile elements and the fertilization of the lithospheric mantle.

high P-T experiment↗
DOE OSTI2021

Towards stacking fault energy engineering in FCC high entropy alloys

Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the explicit use of SFE as an alloy design parameter. In this work, we combine DFT calculations (with necessary configurational averaging), machine-learning (ML) and physics-based models to predict the SFE in the fcc CoCrFeMnNiV-Al high-entropy alloy space. The best-performing ML model is capable of accurately predicting the SFE of arbitrary compositions within this 7-element system. Finally, this efficient model along with a recently developed model to estimate intrinsic strength of fcc HEAs is used to explore the strength–SFE Pareto front, predicting new-candidate alloys with particularly interesting mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sedimentary Geothermal Resources in Nevada, Utah, Colorado, and Texas

The objectives of this project were to (1) perform a literature review of sedimentary geothermal resources, (2) identify data sources and develop data-collection methodologies that characterize selected resources, (3) screen sedimentary basins and formations for sedimentary geothermal potential, and (4) evaluate the technical feasibility of one or more selected locations. Numerous publications have characterized geothermal resources within sedimentary basins. A literature search reviewed publications describing resources located in Colorado, Louisiana, Nevada, Texas, Utah, and Wyoming. The most attractive resources have high temperature gradients, low drilling costs, and reservoir permeabilities greater than 10 millidarcies (mD). Prospects in Colorado, Nevada, Texas, and Utah exhibit attractive characteristics and were chosen for further analysis. Sedimentary resources in Nevada and Utah are most attractive, followed by tested resources in Texas and untested resources in Colorado. The identified resources in Wyoming and Louisiana had lower geothermal gradients and were not evaluated. Reservoir modeling and techno-economic analysis were performed at Marys River Basin–North in Nevada. Geothermal energy production at this location is expected to have a levelized cost of energy (LCOE) ranging between 10 and 20 cents/kWh. Additional work may result in lower LCOE estimates at this location and at other attractive prospects in these three regions. The Great Basin carbonate and alluvial aquifer system of eastern Nevada and western Utah includes a lower carbonate aquifer unit, which has the potential for hosting both conduction- and convection-dominated geothermal systems. Mapping has identified lateral thickness variability expressed as a roughly 120–150-km-wide central corridor, which hosts the thickest and most continuous formations and extends from near Las Vegas north to the Idaho border. Purchase and analysis of privately held legacy seismic data could potentially compensate for the lack of sufficient data documenting measured depth to the lower carbonate aquifer unit. Multiple orogenies, extension episodes, and intrusive events have deformed and displaced the target formations of the lower carbonate aquifer unit. This structural complexity and potential for dipping reservoirs emphasizes the need for detailed geological, geophysical, and reservoir modeling. Heat flow within three Colorado sedimentary basins reviewed as part of this study was calculated in targeted studies by the Colorado Geologic Survey and Colorado School of Mines. These calculations are based on bottom-hole temperature data sets with significant limitations and some variability but produce values consistently higher than the global continental average of 65 mW/m2 for all three basins. Heat flow in the Raton Basin is the highest; however, permeability measurements from specific sedimentary formations with high heat flow have not been obtained. Promising formations for sedimentary geothermal systems were found in all three regions studied—Nevada-Utah, Colorado, and Texas. The next steps for developing sedimentary geothermal resources vary due to differences in available data and resource uncertainties. Additional scopes of work are recommended for identified basins in these three regions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Microstructural Model for Creep-Fatigue Interaction in Grade 91 Steel

This report describes the completion of a crystal plasticity finite element simulation framework for creep-fatigue deformation and damage in Grade 91 steel. In addition to a description of the framework and a set of verification calculations, the report details an initial set of simulations applying the framework to examine low cycle fatigue and creep-fatigue interaction in Grade 91 at 600° C. Grade 91 is a likely material for future advanced microreactors and creep-fatigue will be a limiting structural damage mechanism, particularly for mobile reactors or load following systems. The framework described here could provide a better understanding of creep-fatigue damage mechanisms and directly examine creep, fatigue, and creep-fatigue interaction for Grade 91 at realistic component service conditions. This in turn could lead to improved design methods and component life estimation approaches, leading to more efficient and economical component designs.

36 MATERIALS SCIENCE↗