SEARCH · Engineering Papers

Results for “fracturing fluid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Advances and trends in structures and dynamics; Proceedings of the Symposium, Washington, DC, October 22-25, 1984

Among the topics discussed are developments in structural engineering hardware and software, computation for fracture mechanics, trends in numerical analysis and parallel algorithms, mechanics of materials, advances in finite element methods, composite materials and structures, determinations of random motion and dynamic response, optimization theory, automotive tire modeling methods and contact problems, the damping and control of aircraft structures, and advanced structural applications. Specific topics covered include structural design expert systems, the evaluation of finite element system architectures, systolic arrays for finite element analyses, nonlinear finite element computations, hierarchical boundary elements, adaptive substructuring techniques in elastoplastic finite element analyses, automatic tracking of crack propagation, a theory of rate-dependent plasticity, the torsional stability of nonlinear eccentric structures, a computation method for fluid-structure interaction, the seismic analysis of three-dimensional soil-structure interaction, a stress analysis for a composite sandwich panel, toughness criterion identification for unidirectional composite laminates, the modeling of submerged cable dynamics, and damping synthesis for flexible spacecraft structures.

Noor, A. K.↗

NASA NTRS2003

MN Carbonates in the Martian Meteorite Nakhla: Possible Evidence of Brine Evaporation

The importance of secondary phases in martian meteorites lies in their potential to provide clues about the martian environments responsible for their formation. During this study, we analyzed a number of carbonate-bearing fracture surfaces from the Nakhla meteorite. Here we describe the physical and chemical properties of several manganese-calcium-rich siderites. Additionally, we describe a potential model for the formation and alteration of these carbonates, and we suggest constraints on the conditions responsible for their precipitation. Nakhla is an olivine-bearing clinopyroxenite with minor amounts of feldspar, FeS, and Fe oxides. Secondary mineral assemblages include vein filling clay with embedded iron oxides, a calcium sulfate, amorphous silica, chlorapatite, halite and carbonates. Bridges and Grady suggested that the carbonates in Nakhla formed from brine evaporation. Isotope studies of the Mn rich siderite are also consistent with formation from hydrothermal fluids with an upper T constraint of ~170 C.

Bailey, J. V.↗

DOE OSTI2021

Advances and challenges in CO2 foam technologies for enhanced oil recovery in carbonate reservoirs

Utilization of surfactants for generating foam has the potential for fluid-mobility control during CO2 flooding processes, leading to improved oil recovery. This manuscript is a review of the most important aspects for the design of CO2-foams with mobility control in carbonate reservoirs, including recent advances in novel surfactants and analytical techniques for analysis of surfactant adsorption on carbonate minerals and their thermal stability. Several key parameters and properties regarding foam transport in porous media are reviewed such as, the minimum pressure gradient required for foam generation, the effects of the partition coefficient of a CO2-soluble surfactant on foam transport in rock cores, visualization of foam flow by either microfluidic or PET/CT imaging processes in heterogeneous or fractured systems and the effect of oil on foam transport under reservoir conditions. Also are discussed the most recent advances of two major foam modeling methods, including the semi-empirical STARS foam model and the population balance model, as well as CO2 foam pilot tests and the factors governing their successes.

Jian, Guoqing↗

NASA NTRS2004

Determination of the Three-Dimensional Morphology of ALH84001 and Biogenic MV-1 Magnetite: Comparison of Results from Electron Tomography and Classical Transmission Electron Microscopy

Dated at approximately 3.9 billion years of age, carbonate disks, found within fractures of the host rock of Martian meteorite ALH84001, have been interpreted as secondary minerals that formed at low temperature in an aqueous medium. Heterogeneously distributed within these disks are magnetite nanocrystals that are of Martian origin. Approximately one quarter of these magnetites have morphological and chemical similarities to magnetite particles produced by magnetotactic bacteria strain MV-1, which are ubiquitous in aquatic habitats on Earth. Moreover, these types of magnetite particles are not known or expected to be produced by abiotic means either through geological processes or synthetically in the laboratory. The remaining three quarters of the ALH84001 magnetites are likely products of multiple processes including, but not limited to, precipitation from a hydrothermal fluid, thermal decomposition of the carbonate matrix in which they are embedded, and extracellular formation by dissimilatory Fe-reducing bacteria. We have proposed that the origins of magnetites in ALH84001 can be best explained as the products of multiple processes, one of which is biological. Recently the three-dimensional (3-D) external morphology of the purported biogenic fraction of the ALH84001 magnetites has been the subject of considerable debate. We report here the 3-D geometry of biogenic magnetite crystals extracted from MV-1 and of those extracted from ALH84001 carbonate disks using a combination of high resolution classical and tomographic transmission electron microscopy (TEM). We focus on answering the following questions: (1) which technique provides adequate information to deduce the 3-D external crystal morphology?; and, (2) what is the precise 3-D geometry of the ALH84001 and MV-1 magnetites?

Thomas-Keprta, Kathie L.↗

DOE OSTI2020

Direct Injection of Biomineralizing Agents to Restore Injectivity and Wellbore Integrity

Summary In this manuscript, we describe the second of two field demonstrations of microbially induced calcium carbonate precipitation (MICP) performed in a failed waterflood injection well in Indiana. In 2012, fracture-related flow pathways developed in the wellbore cement, causing injection water to bypass the oil-bearing formation and enter a high-permeability sandstone thief zone, thereby substantially decreasing injection pressure. In the first field demonstration, our study team characterized the well's mode of failure and successfully applied MICP to decrease flow through the defective cement. However, because the MICP treatment was conducted using a bailer delivery system, the degree of permeability reduction achievable was not adequate to fully restore the historic injection pressure of 1,400 psi at 1 gal/min. For the second field demonstration (reported herein), a direct injection system was developed that substantially increased the injection volume of MICP-promoting fluids. Two strategies were implemented to produce more ureolytic microbes: resuspending concentrated frozen cells immediately before injection and scaling up the bioreactor growth capacity. Multiple pulses of microbes and urea-calcium media were pumped into a string of 1-in.-diameter tubing separated by brine spacers and injected continuously at a flow rate of 3.4 to 1.4 gal/min. During the third day of injection, an injection pressure of 1,384 psi at a flow rate of 1.4 gal/min was achieved, and the experiment was terminated. This study demonstrates that MICP can be successfully used in large-volume applications where the time frame for the delivery of reactants is limited. This finding has significant relevance for commercialization of the MICP biotechnology in the oil and gas industry.

Engineering↗

DOE OSTI2021

Thermal Stimulation and Injectivity Testing at Raft River, ID EGS Site

The Raft River geothermal field is the site of an innovative Department of Energy Enhanced Geothermal System (EGS) project to determine the viability of using combined thermal and hydraulic stimulation techniques to improve energy production. Well RRG-9 is currently undergoing a stimulation program using injectate from the US Geothermal Raft River Power Plant and cold water from a cooling tower make-up water well. The stimulation began on 13 June 2013 with injection from the power plant at a temperature of about 39 °C and a pressure of 275 psig. Next, two positive displacement plunger type pumps were used to increase the injection pressure and flow rate for about one month. The highest rate achieved was 258 gpm at a pressure of 741 psig. During this time, fluid from the cooler water well was injected for about 2 weeks at various pressures. Then, the pumps were removed and plant injection resumed on 25 September. Plant injection will continue until the spring of 2014, when a high pressure hydraulic stimulation will be conducted. A series of seismic monitoring stations deployed around the well are providing data on seismic events occurring at the site. Over the past year, 51 microseismic events have been recorded, all less than Magnitude 1. During injection, several diagnostic tests were conducted to gain a better understanding of the well and reservoir. A step-rate test was performed on 22 August to measure the in-situ stress and aid in modeling in-situ fractures. A tracer was injected into the well on 9 September. No tracer was detected in adjacent production wells after several months. A second borehole televiewer survey was conducted for comparison to pre-stimulation images. Here, a third borehole televiewer survey is planned after the high pressure stimulation. Injection test data is evaluated in real time. A modified Hall plot analysis indicates the effective permeability is increasing. The injectivity index supports the results of the modified Hall plot analysis. As the thermal stimulation has continued, the injectivity index has consistently followed an upward trend from 0.1 gpm/psi to 0.53 gpm/psi.

enhanced geothermal systems↗

NASA NTRS1991

The Role of Calcium in Osteoporosis

Calcium requirements may vary throughout the lifespan. During the growth years and up to age 25 to 30, it is important to maximize dietary intake of calcium to maintain positive calcium balance and achieve peak bone mass, thereby possibly decreasing the risk of fracture when bone is subsequently lost. Calcium intake need not be greater than 800 mg/day during the relatively short period of time between the end of bone building and the onset of bone loss (30 to 40 years). Starting at age 40 to 50, both men and women lose bone slowly, but women lose bone more rapidly around the menopause and for about 10 years after. Intestinal calcium absorption and the ability to adapt to low calcium diets are impaired in many postmenopausal women and elderly persons owing to a suspected functional or absolute decrease in the ability of the kidney to produce 1,25(OH)2D2. The bones then become more and more a source of calcium to maintain critical extracellular fluid calcium levels. Excessive dietary intake of protein and fiber may induce significant negative calcium balance and thus increase dietary calcium requirements. Generally, the strongest risk factors for osteoporosis are uncontrollable (e.g., sex, age, and race) or less controllable (e.g., disease and medications). However, several factors such as diet, physical activity, cigarette smoking, and alcohol use are lifestyle related and can be modified to help reduce the risk of osteoporosis.

Arnaud, C. D.↗

DOE OSTI2024

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗

NASA NTRSDate not supplied

Indigenous Organic Matter Associated with Unusual Low Temperature Aqueous Secondary Phases in the Mars Meteorite Nakhla

In freshly fractured chips of the Mars meteorite Nakhla, we found carbonaceous phases intimately associated with secondary aqueous alteration phases, both mineral and amorphous, interpreted to have formed through the low-temperature, mineral-water dissolution of the host Martian basalt.We have characterized the nature of carbonaceous matter and spatially associated secondary phases within Nakhla using a multidisciplinary suite of optical and electron beam instrumentation(Thomas-Keprta et al., 2022). Secondary Phases: The preserved secondary mineral phases are best described as ‘juvenile’ hydrated amorphous silica formed when water interacts with mafic rocks. Although amorphous silica on Earth rarely persists through geologic time due to alteration into other phases, e.g.,microcrystalline opal-CT and eventually quartz, the lack of maturation of primarily amorphous silica with embedded nano phase ferrihydrite in Nakhla infers these phases formed in a limited water environment by low temperature, circumneutral pH fluids stabilized by the presence of coexisting amorphous silica.Additionally, we report the presence of the ferrous hydroxycarbonate mineral chukanovite (Fe2(OH)2CO3)within some of the secondary mineral aggregates studied. Organic Matter:Indigenous Nakhla carbonaceous matter is present both in condensed phases and in a dispersed state spatially associated with secondary alteration phases. In the former, discrete refractory micron to submicron assemblages were identified that appear macromolecular in nature and, in several cases, associated with F and, in one case, significant N. While organically bound Fhas yet to be detected remotely on Mars, the identification of organic N within one of the carbonaceous phases is considered a prerequisite for Martian astrobiological potential. Summary: Neither the intimate association of carbonaceous matter with secondary phases nor the identification of chukanovite have been previously reported in any of the Martian meteorites.A variety of habitability-related metastable solid-phase, oxidation-reduction, and organic-molecular sample attributes in Nakhla formed hundreds of millions of years ago near Mars’ surface; these alteration products persisted on Mars until very recentlyand resemble the amorphous materials observed at Gale crater by Mars Science Laboratory CuriosityRover (Rampe et al., 2020).Similar phases may be expected to occur elsewhere on Mars (e.g., Jezerocrater). References: Rampe E. B. et al.(2020) Mineralogy and geochemistry of sedimentary rocks and eolian sediments in Gale crater, Mars: A review after six Earth years of exploration with Curiosity.Geochemistry80, 125605; Thomas-Keprta K.L. et al. (2022) Indigenous carbon-rich matter associated with unusual aqueous alteration features in Nakhla: Complex formation and preservation history. Geochim. et Cosmochim. Acta320, 41-78 (online).

Mars meteorite↗

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 OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1992

Large impacts in the Baltic shield with special attention to the Uppland structure

Within the Baltic Shield several very large structures have been identified and are suspected to be of meteorite impact origin. Some of these deeply eroded circular features are presented with special attention to the Uppland structure, where several indications point toward an impact origin in the mid-Proterozoic. The structures exceed 100 km in diameter and the topographic expression is inferior or absent. An arcuate arrangement of lithologies occurs around the margin of the structures and the central regions show conform magnetic and positive gravity anomalies. The Uppland structure is approximately 320 km in diameter as expressed by morphological, geological, and geophysical concentric patterns. The central part is topographically remarkably flat and is characterized by an unusual irregular fracture pattern. A subcircular central tonalite with density of 2.81 Mg(sup -3) gives a positive gravity anomaly of 35 mgal and the gravimetric profile is very similar to that of Manicouagan and Vredefort. The tonalite constitutes a huge antiform, 80 km in diameter, probably representing a 12-km structural uplift of infracrustal rocks. The flancs of the tonalite are characterized by recrystallized pseudotachylitic breccia dykes and breccia zones. Around the central parts amphibolite-grade metamorphic rocks appear as large fragments within a fine-grained granite interpreted as a thermally annealed melt rock. Several occurrences of breccia dykes and breccia-bearing melts have been identified about 100 km from the gravimetric center of the structure. Impact-related ore deposits are located around the margin of the structure and are interpreted as preexisting downfaulted iron formations, and deposits formed from remobilization of these preimpact occurrences. The so-called ball ores are interpreted to have formed by fluid injection similar to the formation of breccia dykes. The extensive hydrothermal alteration along the outer margin of the structure have created extreme soda and K-enriched rocks ('leptites') from preexisting gneiss granites and supracrustal sedimentary gneisses.

Henkel, H.↗

NASA NTRS1999

Olivine in Martian Meteorite Allan Hills 84001: Evidence for a High-Temperature Origin and Implications for Signs of Life

Olivine from Martian meteorite Allan Hills (ALH) 84001 occurs as clusters within orthopyroxene adjacent to fractures containing disrupted carbonate globules and feldspathic shock glass. The inclusions are irregular in shape and range in size from approx. 40 microns to submicrometer. Some of the inclusions are elongate and boudinage-like. The olivine grains are in sharp contact with the enclosing orthopyroxene and often contain small inclusions of chromite The olivine exhibits a very limited range of composition from Fo(sub 65) to Fo(sub 66) (n = 25). The delta(sup 18)O values of the olivine and orthopyroxene analyzed by ion microprobe range from +4.3 to +5.3% and are indistinguishable from each other within analytical uncertainty. The mineral chemistries, O-isotopic data, and textural relationships indicate that the olivine inclusions were produced at a temperature greater than 800 C. It is unlikely that the olivines formed during the same event that gave rise to the carbonates in ALH 84001, which have more elevated and variable delta(sup 18)O values, and were probably formed from fluids that were not in isotopic equilibrium with the orthopyroxene or olivine The reactions most likely instrumental in the formation of olivine could be either the dehydration of hydrous silicates that formed during carbonate precipitation or the reduction of orthopyroxene and spinel If the olivine was formed by either reaction during a postcarbonate beating event, the implications are profound with regards to the interpretations of McKay et al. Due to the low diffusion rates in carbonates, this rapid, high-temperature event would have resulted in the preservation of the fine-scale carbonate zoning' while partially devolatilizing select carbonate compositions on a submicrometer scale. This may have resulted in the formation of the minute magnetite grains that McKay et al attributed to biogenic activity.

Shearer, C. K.↗

DOE OSTI2024

Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations

This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗

DOE OSTI2022

Parameter estimation from spontaneous imbibition into volcanic tuff

Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.

58 GEOSCIENCES↗

DOE OSTI2025

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗