Search NASA⌕ Search

SEARCH · Search NASA

Results for “FRACTURE PROPERTIES”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗
DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

U 3 Si 2 and UO 2 composites densified by spark plasma sintering for accident-tolerant fuels

This work reports the synthesis and characterization of the U 3 Si 2 and UO 2 composites sintered by spark plasma sintering (SPS) with controlled microstructures for accident-tolerant fuel application. The U 3 Si 2 and UO 2 composites with different silicide and oxide ratios were synthesized by SPS at temperatures from 1000 to 1300 °C for 5 minutes. The microstructure and phase composition of the SPS densified composite fuels were characterized with scanning electron microscopy, X-ray diffraction (XRD), and energy dispersed spectroscopy (EDS). A systematic study of the thermal and mechanical properties was conducted using microhardness testing and laser flash apparatus, along with oxidation resistance measurements using thermogravimetric analysis (TGA). The results show that the synthesis of composite fuels can be achieved with a 90% theoretical density (TD) at 1000 °C and over 95% TD when sintered at 1300 °C. XRD and EDS results confirmed that the dominant phases in the composites are U 3 Si 2 and UO 2 . Improved physical density generally leads to improved hardness, fracture toughness, thermal diffusivity, and onset temperature during the oxidation process. U 3 Si 2 was found to play a dominant role in determining the mechanical and oxidation properties of the composite fuels, whereas UO 2 had a more important impact on controlling the thermal diffusivity of the composites. The composite with 50 wt% UO 2 sintered at 1300 °C displayed the onset oxidation temperature of 500 °C by dynamic oxidation testing using TGA at a ramp degree of 10 °C/min. The composite also achieved a high fracture toughness of ~3.5 MPa m ½ . Finally, these results highlight the potential of composite fuel forms densified by SPS with simultaneously enhanced fissile element density, fracture toughness, thermal transport properties, and oxidation resistance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The effect of Ni doping on the mechanical behavior of tungsten

The effect of nickel (Ni) on activated sintering of tungsten has been known and investigated for decades. However, there are few reports on the mechanical properties of tungsten when Ni is added as a low alloying additive with less than 1 wt%. In this work, the mechanical behaviors of tungsten doped with 0.2–1.0 wt% nickel from room temperature up to 400 °C were investigated. The results show that the compression ductility of as-sintered Ni doped tungsten is drastically higher than that of as-sintered pure tungsten at all tested temperatures including room temperature, and the fracture toughness at room temperature is increased by 92%. The results suggest that the ductile-to-brittle transition temperature (DBTT) of tungsten can be reduced substantially by doping with Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Creep Behavior of a High-$\gamma^{\prime}$ Ni-Based Superalloy Fabricated via Electron Beam Melting

Additive manufacturing enables the fabrication of complex engineering components previously inaccessible through traditional processes. Nickel-base superalloys with large γ' volume fraction are typically considered non-weldable and therefore exhibit a propensity for cracking during the fusion process. These crack-prone materials, however, are of great importance in gas turbine engines due to their excellent high temperature creep resistance. In this study we investigate the creep behavior of IN738LC produced by the electron beam melting process. Here, we find that with appropriate post-build heat treatment the creep response of material oriented in the build direction exhibits deformation and rupture behavior comparable to that of conventionally cast IN738 & IN738LC. In the transverse direction properties fall below the expected cast behavior, however, we argue this is likely due to differences in grain scale and crystallographic texture. It may be possible to coarsen the grain morphology with appropriate process-parameter optimization in order to reduce the severity of intergranular fracture in the transverse direction. These results illustrate that high temperature properties exhibited by additively manufactured IN738LC are suitable for the hot section of gas turbine engines.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of Composition on Mechanical Properties and Conductivity of the Dual-Ion Conductor Na 1+ x Mn x /2 Zr 2– x /2 (PO 4 ) 3 for Solid-State Batteries

Here we report the dependence of mechanical properties and ionic conductivity of the Na 1+ x Mn x /2 Zr 2– x /2 (PO 4 ) 3 (designated as NMZP) ionic conductor on its composition (with x = 0.5, 1.0, 1.5, and 2.0). Local mechanics was studied by instrumented nanoindentation, and the properties of interest included the elastic modulus, hardness, and fracture toughness. Overall, the material becomes more compliant and soft with the increase in the Na content. The elastic modulus reduces from 120 to 80 GPa, and nanoindentation hardness reduces from 7.8 to 4.2 GPa when x changes from 0.5 to 2.0. The relationship with fracture toughness ( K C ) is highly nonlinear: fracture toughness first decreases and then increases with the increase in the sodium content, reaching a maximum value of 0.89 MPa × m 1/2 at x = 2.0. Such a relationship was found to correlate with the formation of a glassy phase in NMZP at intermediate Na concentration. The maximum ionic conductivity coincides with the maximum fracture toughness in this material.

25 ENERGY STORAGE↗
DOE OSTI2021

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗
DOE OSTI2023

Using Artificial Microstructures to Understand Microstructure Property Relationship-Toughening Mechanisms in Metallic Glass (Final Report)

Metallic glasses are a new class of structural materials which exhibit exciting mechanical properties including high strength and elasticity. In terms of fracture toughness, the material class of metallic glasses spans a wide range; Some metallic glasses are extremely brittle and exhibit near ideal brittle behavior whereas others can be exceptional tough with values comparable to the toughest metals out there. Such large range of observed fracture toughness within the material class of metallic glasses is surprising as they have seemingly a very similar atomic structure. Therefore, we developed “artificial microstructures” which allows to decouple the various contribution of sample geometry, imperfection and structure. Specifically, we decouple variations in the alloys’ chemical composition and the atomic structure and quantified the resulting fracture toughness. Atomic structure of a metallic glass can be modified by the fictive temperature. The fictive temperature of a glass is the temperature at which the liquid metallic glass falls out of equilibrium upon colling and forms a glass. Upon further cooling the structure is maintained only thermal oscillations decrease due to a lower absolute temperature. We found that the effect of fictive temperature (same chemistry, different structural stages of the glass) is comparable to the variations of fracture toughness when the chemistry is varied. Hence, it appears that the subtle differences in the glass structure are responsible for the large range of fracture toughness’s observed. Our results reveal that fracture toughness within the material class of metallic glasses varies significantly and we found some example alloys with exceptional high resistance to fracture and others that are almost ideally brittle. Significant influences on the fracture toughness have the structure of the glass, its chemistry, and some imperfections in the structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High Thermoelectric Performance in 2D Sb 2 Te 3 and Bi 2 Te 3 Nanoplate Composites Enabled by Energy Carrier Filtering and Low Thermal Conductivity

Thermoelectrics are an important class of materials with great potential in alternative energy applications. In this study, two-dimensional (2D) nanoplates of the layered chalcogenides, Sb 2 Te 3 and Bi 2 Te 3 , are synthesized and composites of the two are investigated for their thermoelectric properties. The two materials, Sb 2 Te 3 and Bi 2 Te 3 , were synthesized as hexagonal, 2D nanoplates via a colloidal polyol route. The as-synthesized Sb 2 Te 3 and Bi 2 Te 3 vary drastically from one another in their lateral and vertical dimensions as revealed by scanning electron microscopy and atomic force microscopy. The single crystalline nanoplate nature is deduced by high-resolution transmission electron microscopy and selected area electron diffraction. Nanoplates have well-defined hexagonal facets as seen in the scanning and transmission electron microscopy images. The nanoplates were consolidated as an anisotropic nanostructured pellet via spark plasma sintering. Preferred orientation observed in the powder X-ray diffraction pattern and scanning electron microscopy images of the fractured pellets confirm the anisotropic structure of the nanoplates. Thermoelectric properties in the parallel and perpendicular directions were measured, revealing strong anisotropy with a significant reduction to thermal conductivity in the perpendicular direction due to increased phonon scattering at nanoplate interfaces. All compositions, except that of the 25% Bi 2 Te 3 nanoplate composite, behave as degenerate semiconductors with increasing electrical resistivity as the temperature increases. The Seebeck coefficient is also increased dramatically in the nanocomposites, the highest reaching 210 μV/K for 15% Bi 2 Te 3 . The increase in Seebeck is attributed to energy carrier filtering at the nanoplate interfaces. Overall, these enhanced thermoelectric properties lead to a drastic increase in the thermoelectric performance in the perpendicular direction, with zT ~ 1.26, for the 15% Bi 2 Te 3 nanoplate composite at 450 K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A parametric analysis of oilfield design factors affecting the detectability and characterization of electrically conductive hydrofracks

Electrical responses in the vicinity of energized steel-cased well sources offer significant potential for monitoring induced fractures. However, the high complexity of well-fracture-host models spanning multiple length scales compels analysts to simplify their numerical models due to enormous computational costs. This consequently limits our understanding regarding monitoring capabilities and the limitations of electrical measurements on realistic hydraulically fracturing systems. In this paper, we use the hierarchical finite element approach to construct geoelectric models in which geometrically complex fractures and steel-cased wells are discretely represented in 3D conducting media without sacrificing the model realism and computation efficiency. We have discovered systematic numerical analyses of the electrical responses to evaluate the influences of borehole material conductivity and the source type as well as the effects of well geometry, conductivity contrast, source location, fracture growth, and fracture propagation. Furthermore, the numerical results indicate that the borehole material property has a strong control on the electrical potentials along the production and monitoring wells. The monopole source located at a steel-cased well results in a current density distribution that decays away from the source location throughout the well length, whereas the dipole source produces a current density that dominates mainly along the dipole length. Moreover, the conductivity contrast between the fractures and host does not change the overall pattern of the electrical potentials but varies its amplitude. The fracture models near different well systems indicate that the well geometry controls the entire distribution of potentials, while the characteristics of the voltage difference profiles along the wells before and after fracturing are insensitive to the well geometry and the well in which the source is located. Further, the hydraulic-fracturing models indicate that the voltage differences along the production well before and after fracturing have strong sensitivity to fracture growth and fracture set propagation.

58 GEOSCIENCES↗