Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Creep modeling”

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 991 records · Page 55

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

An Investigation of SiC/SiC Woven Composite Under Monotonic and Cyclic Loading

The desirable properties in ceramic matrix composites (CMCs), such as high temperature strength, corrosion resistance, high toughness, low density, or good creep resistance have led to increased use of CMCs in high-speed engine structural components and structures that operate in extreme temperature and hostile aero-thermo-chemical environments. Ceramic matrix composites have been chosen for turbine material in the design of 21 st-century civil propulsion systems to achieve high fuel economy, improved reliability, extended life, and reduced cost. Most commercial CMCs are manufactured using a chemical vapor infiltration (CVI) process. However, a lower cost fabrication known as melt-infiltration process is also providing CMCs marked for use in hot sections of high-speed civil transports. The scope of this paper is to report on the material and mechanical characterization of the CMCs subjected to this process and to predict the behavior through an analytical model. An investigation of the SiC/SiC 8-harness woven composite is ongoing and its tensile strength and fatigue behavior is being characterized for room and elevated temperatures. The investigation is being conducted at below and above the matrix cracking stress once these parameters are identified. Fractography and light microscopy results are being studied to characterize the failure modes resulting from pure uniaxial loading. A numerical model is also being developed to predict the laminate properties by using the constituent material properties and tow undulation.

Lang, J.↗

MURR LEU structural and thermal hydraulics analyses: Part II – Impacts of irradiation thermo-mechanical behavior on thermal hydraulics safety analyses

A series of structural analyses have been performed to support the conversion of the University of Missouri Research Reactor (MURR) from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, thermal expansion, as well as thermal resistance from the oxide layer growth for the MURR LEU element in prototypic thermal and irradiation conditions as presented in Part I of this article. Overall, this irradiation thermo-mechanical analysis predicts smaller gap thickness reductions in previously limiting regions, and larger reductions in the middle of the outermost end channels where power density is not typically a maximum. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This article presents the results from the steady-state safety analysis for normal operation as well as the two most limiting accident scenarios. The calculation models were revised in order to account for the spatial and temporal variation of the channel gap thicknesses. The results show that sufficient safety margins are still maintained for normal operation as well as during the postulated accident transients. This work provides a methodology of incorporating the irradiation thermo-mechanical behavior of plate-type fuel into the thermal hydraulics safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Synergistic Effects of Physical Aging and Damage on Long-Term Behavior of Polymer Matrix Composites

The research consisted of two major parts, first modeling and simulation of the combined effects of aging and damage on polymer composites and secondly an experimental phase examining composite response at elevated temperatures, again activating both aging and damage. For the simulation, a damage model for polymeric composite laminates operating at elevated temperatures was developed. Viscoelastic behavior of the material is accounted for via the correspondence principle and a variational approach is adopted to compute the temporal stresses within the laminate. Also, the effect of physical aging on ply level stress and on overall laminate behavior is included. An important feature of the model is that damage evolution predictions for viscoelastic laminates can be made. This allows us to track the mechanical response of the laminate up to large load levels though within the confines of linear viscoelastic constitutive behavior. An experimental investigation of microcracking and physical aging effects in polymer matrix composites was also pursued. The goal of the study was to assess the impact of aging on damage accumulation, in ten-ns of microcracking, and the impact of damage on aging and viscoelastic behavior. The testing was performed both at room and elevated temperatures on [+/- 45/903](sub s) and [02/903](sub s) laminates, both containing a set of 90 deg plies centrally located to facilitate investigation of microcracking. Edge replication and X-ray-radiography were utilized to quantify damage. Sequenced creep tests were performed to characterize viscoelastic and aging parameters. Results indicate that while the aging times studied have limited ]Influence on damage evolution, elevated temperature and viscoelastic effects have a profound effect on the damage mode seen. Some results are counterintuitive, including the lower strain to failure for elevated temperature tests and the catastrophic failure mode observed for the [+/- 45/9O3](sub s), specimens. The fracture toughness for transverse cracks increases with increasing temperature for both systems: transverse cracking was completely absent prior to failure in [+/- 45/903](sub s), and was suppressed for [02/903](sub s). No significant effect of damage on aging or viscoelastic parameters was observed.

Brinson, L. Cate↗

Factors Influencing Dwell Fatigue Cracking in Notches of Powder Metallurgy Superalloys

The influences of heat treatment and cyclic dwells on the notch fatigue resistance of powder metallurgy disk superalloys were investigated for low solvus high refractory (LSHR) and ME3 disk alloys. Disks were processed to produce material conditions with varied microstructures and associated mechanical properties. Notched specimens were first subjected to baseline dwell fatigue cycles having a dwell at maximum load, as well as tensile, stress relaxation, creep rupture, and dwell fatigue crack growth tests at 704 C. Several material heat treatments displayed a bimodal distribution of fatigue life with the lives varying by two orders-of-magnitude, while others had more consistent fatigue lives. This response was compared to other mechanical properties, in search of correlations. The wide scatter in baseline dwell fatigue life was observed only for material conditions resistant to stress relaxation. For selected materials and conditions, additional tests were then performed with the dwells shifted in part or in total to minimum tensile load. The tests performed with dwells at minimum load exhibited lower fatigue lives than max dwell tests, and also exhibited early crack initiation and a substantial increase in the number of initiation sites. These results could be explained in part by modeling evolution of peak stresses in the notch with continued dwell fatigue cycling. Fatigue-environment interactions were determined to limit life for the fatigue cycles with dwells.

Gabb, T. P.↗

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗

Sample problems for Section III, Division 5 design by inelastic analysis of Grade 91 components

This report works through two sample design analyses of representative Grade 91 components using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Appendix HBB-T rules for the strain and deformation and creep-fatigue design criteria using a design by inelastic analysis approach. The component geometries and loadings were selected to span a wide variety of potential applications. The sample components are an endplate geometry, which is fairly representative of typical pressure vessel construction, and a core block, representing a critical component in a heat pipe microreactor, which has a complicated geometry and complex thermal stress history. These sample problems serve several functions. First, they are complete worked examples of how to apply the current design by inelastic analysis rules. Second, they serve as test cases to evaluate the current rules and provide feedback on how they could be optimized to deal with complicated component geometries or to be better automated in modern finite element analysis software. Finally, the analysis here uses an inelastic constitutive model for Grade 91 developed by Argonne National Laboratory and proposed for inclusion in the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules as a reference constitutive model. Using this constitutive model in the sample problem analyses tests the model in the context of the ASME Code rules. The final conclusions of this report are that the material model and the current Code rules are satisfactory, but that the Code rules could be optimized to simplify automation and reduce over conservatism. Future work should address the specific optimization topics discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DGPS measurements of solifluction and hillslopes creep at the Teller 47 field site from 2017-2022, Seward Peninsula, AK

DGPS measurements of benchmarks locations were collected during the 2017, 2018, 2019, 2021 and 2022 summers at the NGEE Arctic Teller-47 solifluction site on the Seward Peninsula, Alaska. Benchmarks were installed in 2017 on several solifluction lobes and hillslope locations and re surveyed in 2018 and 2019 with a Trimble R10 DGPS instrument. In 2018, coincident with a UAS-based laser altimetry collection campaign, more ground control targets were installed and subsequently remeasured in 2019. Additional ground control targets were installed in 2021. Points were corrected using NOAA’s Online Positioning User Software (OPUS) and Trimble Business Center (TBC) software. An update to this dataset was made in July 2024, adding new data from 2021 and 2022. The dataset now contains ten total *.csv files: five are data files plus one file-level metadata and five data dictionaries, two *.jpg, four shape file folders, two *.kml, one archive folder of raw data, and one *.pdf. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Utilizing a Dynamic Segmentation Convolutional Neural Network for Microstructure Analysis of Additively Manufactured Superalloy 718

Additive manufacturing (AM) is revolutionizing almost all industries through the production of intricate geometries previously prohibited by cost or machinability. Ni-based superalloys form a primary alloy class for high temperature applications in the petrochemical, aerospace, and nuclear industries because of their intrinsic resistance to creep and the ability to heat treat the superalloy for high strength. Despite these attractive properties, the extreme work hardening of Ni-based superalloys makes traditional manufacturing of complex shapes difficult and thus, these alloys are an attractive target for AM. Superalloy 718 was chosen as an example superalloy because of the wide variety of precipitates that can form within its composition space from the repetitive heating and cooling cycles of the AM process. The precipitates and other microstructure features, such as grain boundaries and dislocations, will dictate the mechanical properties and thus, there is an extensive challenge to characterize the size, number density, composition, and volume fraction of each microstructural feature from AM fabrication using analytical electron microscopy. This work focused on the application of a pixel-wise classification machine learning (ML) model called a dynamic segmentation convolutional neural network (DSCNN) to identify the microstructural features of an as-fabricated additively manufactured superalloy 718.

Taller, Stephen↗

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗

Influence of High Cycle Thermal Loads on Thermal Fatigue Behavior of Thick Thermal Barrier Coatings

Thick thermal barrier coating systems in a diesel engine experience severe thermal Low Cycle Fatigue (LCF) and High Cycle Fatigue (HCF) during engine operation. In the present study, the mechanisms of fatigue crack initiation and propagation, as well as of coating failure, under thermal loads which simulate engine conditions, are investigated using a high power CO2 laser. In general, surface vertical cracks initiate early and grow continuously under LCF and HCF cyclic stresses. It is found that in the absence of interfacial oxidation, the failure associated with LCF is closely related to coating sintering and creep at high temperatures, which induce tensile stresses in the coating after cooling. Experiments show that the HCF cycles are very damaging to the coating systems. The combined LCF and HCF tests produced more severe coating surface cracking, microspallation and accelerated crack growth, as compared to the pure LCF test. It is suggested that the HCF component cannot only accelerate the surface crack initiation, but also interact with the LCF by contributing to the crack growth at high temperatures. The increased LCF stress intensity at the crack tip due to the HCF component enhances the subsequent LCF crack growth. Conversely, since a faster HCF crack growth rate will be expected with lower effective compressive stresses in the coating, the LCF cycles also facilitate the HCF crack growth at high temperatures by stress relaxation process. A surface wedging model has been proposed to account for the HCF crack growth in the coating system. This mechanism predicts that HCF damage effect increases with increasing temperature swing, the thermal expansion coefficient and the elastic modulus of the ceramic coating, as well as the HCF interacting depth. A good agreement has been found between the analysis and experimental evidence.

Zhu, Dongming↗

Ab-Initio Study of Stacking Fault Segregation Behavior in Ni-based Superalloys at High Temperatures

Ni-based superalloys demonstrate extraordinary creep properties compared to traditional metallic systems due to the formation of L12 (γ’) precipitates. Despite experimental evidence of alloying element segregation to the stacking faults, prior ab-initio efforts indicated repulsion of Nb from superlattice intrinsic stacking faults (SISFs) . Due to the importance of this element in key strengthening mechanisms, we developed a new technique to account for multi-component segregation behavior at temperature combining ab-initio molecular dynamics (AIMD) and Monte-Carlo (MC) simulations. Using this technique, we replicated experimentally observed segregation behavior of Nb, Co, and Ti to the SISF in L12 (Figure 1). Furthermore, we investigated the co-segregation behaviors of individual pairs of alloying elements to determine the interaction effects which lead to this segregation behavior in the full alloy system. By simulating the reduced-chemistry models using our approach, we reconciled previous ab-initio calculations demonstrating repulsion of Nb from the SISF and the experimentally observed segregation behavior.

Molecular Dynamics↗

Concurrent Precipitation of Nb(C,N) and Metastable M 23 C 6 in Alloy 347H at 700°C and 750°C: Computer Simulations and Comparison to Experiment

Here, we present our results for the concurrent precipitation of metastable M 23 C 6 , Nb(C,N) secondary precipitates, and the Nb(C,N) primary crystals in 347H austenitic stainless steel. For precipitation modeling, we have accounted for the elastic contribution to interfacial energy, and for the Fe-spin-polarization for NbC/Fe and M 23 C 6 /Fe interfacial energy values: for NbC/Fe ~ 0.63 J/m 2 . For M 23 C 6 precipitates, an error function was used to describe the interfacial energy growth with particle size. In precipitation simulations, the average size of the primary Nb(C,N) particles remained ~ 1 μm at 700°C and ~ 0.3 μm at 750°C. The M 23 C 6 precipitates at 750°C dissolved after 120 h (our simulations) compared to 300 h (experiments). The Nb(C,N)/Fe interfacial energy was not affected by the nitrogen additions. With these modifications, reasonable agreement with the available experimental data was obtained, which allows using them in the development of the 2nd-phase particle-informed creep theory.

36 MATERIALS SCIENCE↗

Geometric Acoustics for Aircraft Noise Scattering

This paper discusses aircraft noise scattering by geometric acoustics, which consists of the basic features of sound propagation and reflection in rays or ray tubes, diffraction by smooth geometry in terms of surface creeping waves, and diffraction by abrupt geometry features such as the wing trailing edges. Based on the classic theories of these features, prediction methodologies can be constructed for aircraft noise scattering. The methodologies, however, require important modifications and extensions to account for unique features in aircraft noise applications, for both current conventional aircraft and future unconventional designs. These modifications and extensions include the derivation of a general reflection coefficient that contains the effects of surface geometry curvature, surface impedance, and mean flow. Corrections to the basic formulation of diffraction, both smooth geometry and sharp edges, are formulated to account for the finite dimensions of practical applications. For smooth geometry diffraction, a second order correction is used to continuously transition from insolified to shadow zones and an approach is presented to compute the properties of the geodesic path of surface wave propagation. A model is developed to make use of the ray propagation formulation for incoherent and partially coherent propagation and scattering, which is an important phenomenon in aircraft noise. Mean flow effect is also included in the methodology for low Mach number flows. Examples of calculations based on these theoretical developments are presented to illustrate the unique features in aircraft noise applications.

Yueping Guo↗

Determination of convergence rates across the Ventura Basin, Southern California, using GPS and historical triangulation

Comparison of angles from historical triangulation observations dating as far back as 1932 with Global Positions System (GPS) measurements taken in 1987 indicates that rapid convergence may be taking place on decade timescales in the central and eastern part of the Ventura basin, an east-west trending trough bounded by thrust faults. Changes in angles over this time were analyzed using Prescott's modified Frank's method and in terms of a model which assumes that the regions to the north and south of the basin are rigid blocks undergoing relative motion. For the two block model, inversion of the observed angle changes over the last 28 years for the relative motion vector leads to north-south convergence across the basin of 30 + or - 5 mm/yr, with a left lateral component of 10 + or - 1 mm/yr in the Fillmore-Santa Paula area in the central part of the basin. The modified Frank's method yields strain rates of approximately 2 microrad/yr in both the east and central parts of the basin for measurements spanning the 1971 San Fernando earthquake. Assuming no east-west strain yeilds north-south compression of approximately 3.5 + or - .2 cm/yr. Comparison of triangulation data prior to the earthquake shows no strain outside the margin of error. The convergence rates determined by geodetic techniques are consistent with geologic observations in the area. Such large geodetic deformation rates, with no apparent near-surface creep on the major thrust, can be understood if these faults become subhorizontal at relatively shallow depths and if the subhorizontal portions of the faults are creeping. An alternative explanation of the large displacement rates might be that the pumping of oil in the vicinity of the benchmarks caused large horizontal motions, although it is unlikely that meter scale horizontal motions are due to oil withdrawal. These and other hypotheses are evaluated to better constrain the tectonics of this active region.

Donnellan, Andrea↗

Evolution of Pore Structure and Permeability of Rocks Under Hydrothermal Conditions (Final Report)

The physical and transport properties of porous rocks can be altered by a variety of diagenetic, metamorphic, and tectonic processes, and the changes that result are of critical importance to such industrial applications as resource recovery, carbon dioxide sequestration, and waste isolation in geologic formations. The associated interrelationships between rocks, pore fluids and deformation are also key to understanding many natural processes, including dynamic metamorphism, fault mechanics, fault stability and pressure solution creep. In this project, we investigated the changes of permeability and pore geometry owing to inelastic deformation by solution-transfer, brittle fracturing and dislocation creep. In particular, we studied the coupling between pore fluids and deformation in fluid-filled quartz, calcite, mudstones and ultramafic rocks and examine the effects of coupled mechanical and chemical processes on the evolution of porosity and permeability under hydrothermal conditions. The investigations used a combination of techniques, that included triaxial laboratory experiments; uniaxial compressive loading with in situ observations of microstructure; permeability under hydrostatic compression; observations of microstructure using optical and electron microscopes, and micro-CT imaging; and numerical calculations. Laboratory experiments were designed to provide mechanical and transport data under conditions that isolate the particular mechanisms responsible for the changes. The data obtained were used to quantify changes in surface roughness, porosity, pore dimensions and their spatial fluctuations. The results of the experiments and data from image analyses were compared to the results of network, finite-difference and other numerical models to verify the validity of experimentally established relations between permeability and other rock properties. A bibliography for the period from 2015-2019 is appended below. New results obtained in the final year and the unfunded extension period are detailed below.

58 GEOSCIENCES↗

University of Missouri Research Reactor (MURR) LEU Fuel Fluid-Structure Interaction Analysis

The University of Missouri Research Reactor (MURR ® ) is one of five U.S. high performance research reactors (USHPRR), plus one critical facility, that are actively collaborating with the National Nuclear Security Administration (NNSA) Material Management and Minimization (M 3 ) Reactor Conversion Program to convert from highly enriched uranium (HEU, ≥ 20 wt.% U-235) to low-enriched uranium (LEU, < 20 wt.% U-235) fuel. A new type of LEU fuel with very high density, based on an alloy of uranium and 10 weight percent molybdenum (U-10Mo), is expected to allow the conversion to LEU of USHPRR that have been found unable to be converted with previously qualified uranium silicide-aluminum (U 3 Si 2 -Al) dispersion fuel. MURR has been working with the USHPRR Reactor Conversion (RC) Pillar at Argonne National Laboratory to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analyses, and accident safety analyses in preparation for the conversion of MURR and to support a preliminary Safety Analysis Report (SAR) for conversion to LEU fuel. In this work, Fluid-structure interaction (FSI) analysis at the fuel element level, as compared to the fuel plate level of the previous work. is performed which models all components of the MURR LEU fuel element, including fuel plates and the supporting structures e.g., side plates, end fittings, and combs. Therefore, the effect of supporting structures on the coolant flow distribution, the fuel plate deflection, and the resulting coolant channel gap reduction can be evaluated. In addition to the nominal element geometry and flow rate, the tolerances in the geometry dimensions of coolant channel and plate thickness, the effect of a comb on plate deflection, and the uncertainty of the flow rate per element are considered in this work. The effect of comb on plate displacement is quantified through two bounding cases: the case assuming a perfect bond between the comb and plates and the case neglecting the comb effect. Note that in this analysis, the FSI has been decoupled from the other structural effects caused by irradiation (e.g., swelling and irradiation creep). Assessment of combined effects is planned for a later stage of this project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microbial Curing of Cement for Energy Applications

Rutgers University, Lawrence Livermore National Laboratory, and the University of Arizona executed this program over 36 months with a 4-Task Program: (M1)-Microbial Engineering (M2)-Microstructure Modelling, (M3) Cement and Concrete Formulation, (M4)-Techno-economic analysis (TEA) and lifecycle assessment (LCA). This program developed a new carbonate cement concrete manufacturing process called microbial curing (MBC). MBC utilizes in-situ microbial production of CO 2 that dissolves into pore-bound water to carbonate a cementitious material creating a bonding matrix of CaCO 3 and SiO 2 that hardens and densifies the material. This is the first work of this type where calcium silicate was used in microbial studies. All other work reported in the literature always used a soluble form of calcium. This final report describes the work done in the final quarter, the best procedures and results, and the final techno-economic and lifecycle analyses. This quarter, we cast and cured twenty-five (4”x8”) cylinder samples. The compressive strength, split tensile, Young’s modulus, chloride permeability, and creep measurements were performed. These measurements demonstrated that MBC of calcium silicate concrete exhibits scaling problems due to the escape of ammonia gas, a product of the microbial reaction to generate CO 2 . As curing proceeds, the ability for the gas to escape from within the cylinder is restricted by the outer cylindrical portions densifying, making the outer portion highly impermeable to gas flow. In contrast, the small samples cured are uniformly cured, forming materials whose mechanical properties are 7x better than the large samples. This problem needs to be solved before this technology can be commercialized. The techno-economic and lifecycle analyses indicate that the technology developed in this program exhibits a significant opportunity to reduce the cost of cement and CO 2 emissions associated with concrete, provided the curing issue associated with larger samples can be addressed.

36 MATERIALS SCIENCE↗