Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Boundary zonal flows in rapidly rotating turbulent thermal convection

Recently, in Zhang et al. ( Phys. Rev. Lett. , vol. 124, 2020, 084505), it was found that, in rapidly rotating turbulent Rayleigh–Bénard convection in slender cylindrical containers (with diameter-to-height aspect ratio $\varGamma =1/2$ ) filled with a small-Prandtl-number fluid ( ${Pr}\approx 0.8$ ), the large-scale circulation is suppressed and a boundary zonal flow (BZF) develops near the sidewall, characterized by a bimodal probability density function of the temperature, cyclonic fluid motion and anticyclonic drift of the flow pattern (with respect to the rotating frame). This BZF carries a disproportionate amount ( ${>}60\,\%$ ) of the total heat transport for ${Pr} < 1$ , but decreases rather abruptly for larger ${Pr}$ to approximately $35\,\%$ . In this work, we show that the BZF is robust and appears in rapidly rotating turbulent Rayleigh–Bénard convection in containers of different $\varGamma$ and over a broad range of ${Pr}$ and ${Ra}$ . Furthermore, tirect numerical simulations for Prandtl number $0.1 \leq {\textit {Pr}} \leq 12.3$ , Rayleigh number $10^7 \leq {Ra} \leq 5\times 10^{9}$ , inverse Ekman number $10^{5} \leq 1/{\textit {Ek}} \leq 10^{7}$ and $\varGamma = 1/3$ , 1/2, 3/4, 1 and 2 show that the BZF width $\delta _0$ scales with the Rayleigh number ${Ra}$ and Ekman number ${\textit {Ek}}$ as $\delta _0/H \sim \varGamma ^{0} Pr^{\{-1/4, 0\}} {Ra}^{1/4} {\textit {Ek}}^{2/3}$ ( $\{{\textit {Pr}}<1, {\textit {Pr}}>1\}$ ) and with the drift frequency scales as $\omega /\varOmega \sim \varGamma ^{0} Pr^{-4/3} {Ra}\,{\textit {Ek}}^{5/3}$ , where $H$ is the cell height and $\varOmega$ the angular rotation rate. The mode number of the BZF is 1 for $\varGamma \lesssim 1$ and $2 \varGamma$ for $\varGamma = \{1,2\}$ independent of ${Ra}$ and ${Pr}$ . The BZF is quite reminiscent of wall mode states in rotating convection.

42 ENGINEERING↗
DOE OSTI2026

Role of Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes for Enhanced Geothermal Systems

The ability to sustain high energy extraction efficiency from Enhanced Geothermal Systems (EGS) is affected by the ability to measure, characterize, and predict the effect of different perturbations that arise from fluid injection rates, fluid temperature, and shut-in conditions on existing and stimulated fracture systems throughout a subsurface reservoir’s lifecycle. The difficulty arises from limited knowledge of the interaction of fluid and temperature driven fractures with frictional interfaces that govern local deformation and frictional behavior in the surrounding rock. Thermo-poro-mechanical coupled processes tend to dominate this interaction and strongly affect the permeability and local stress distributions that impact efficiency and prevent optimal stimulation conditions over multi-decadal time frames.

Pyrak-Nolte, Laura [Purdue Univ., West Lafayette, ↗
DOE OSTI2022

FSVPy: A python-based package for fluorescent streak velocimetry (FSV)

Predictive constitutive equations that connect easy-to-measure transport properties (e.g., viscosity and conductivity) with system performance variables (e.g., power consumption and efficiency) are needed to design advanced thermal and electrical systems. In this work, we explore the use of fluorescent particle-streak analysis to directly measure the local velocity field of a pressure-driven flow, introducing a new Python package (FSVPy) to perform the analysis. Fluorescent streak velocimetry combines high-speed imaging with highly fluorescent particles to produce images that contain fluorescent streaks, whose length and intensity can be related to the local flow velocity. By capturing images throughout the sample volume, the three-dimensional velocity field can be quantified and reconstructed. We demonstrate this technique by characterizing the channel flow profiles of several non-Newtonian fluids: micellar Cetylpyridinium Chloride solution, Carbopol 940, and Polyethylene Glycol. We then explore more complex flows, where significant acceleration is created due to microscale features encountered within the flow. We demonstrate the ability of FSVPy to process streaks of various shapes and use the variable intensity along the streak to extract position-specific velocity measurements from individual images. Thus, we demonstrate that FSVPy is a flexible tool that can be used to extract local velocimetry measurements from a wide variety of fluids and flow conditions.

Mechanics↗
DOE OSTI2023

Microseismic Monitoring of Laboratory Hydraulic Fracturing Experiments in Granitic Rocks for Different Fracture Propagation Regimes

While hydraulic fracturing (HF) is a widely employed process, the underlying fracturing processes are still heavily contested. The attributes of the HF generated fracture network can exhibit substantial variation when dealing with specific HF propagation regimes encountered in the field. In this study, HF experiments were performed on true-triaxially loaded Barre granite cubes, with microseismic monitoring, to identify and characterize the fracturing mechanisms associated with different viscosity injection fluids. Utilizing fluids with high (oil/1450 cP) and low (water/1 cP) viscosity represented two key HF propagation regimes: viscosity- and toughness-dominated. The experiments conducted with oil involved higher breakdown pressures, larger fluid volumes, and slower fracture propagation speeds. Furthermore, the frequency–magnitude distribution (b value) for all experiments (1.9–2.3) was similar to those encountered for large-scale operations. Slightly larger b values were encountered during the initiation phase (2.4–2.7) relative to the fracture propagation and post-fracturing phases (1.9–2.2). Techniques such as polarity and moment-tensor inversion were utilized to characterize the source mechanisms. For the HF experiments with oil, tensile fractures were most dominant (92%) in the initiation phase compared to fracture propagation and post-fracturing phases (70–75%). Similar tensile fracturing dominancy was not observed with water, attributable to fluid permeation and leak-off. Regardless of the injection fluid or classification criteria employed, tensile fractures were the dominant type consistently, with fewer occurring in water experiments but the specific ratio of crack types varied with different source mechanism criteria employed.

58 GEOSCIENCES↗
DOE OSTI2024

Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide

The project Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. Sandia Technologies, LLC, and co-investigator Conrad Geoscience Corporation, examined the potential for large-scale, permanent CO2 storage in sedimentary strata within the Newark Rift Basin. The Newark Rift Basin underlies an industrialized, developed region comprising parts of New York, New Jersey, and Pennsylvania. The project characterized and investigated the suitability of Triassic age sedimentary formations for potential geologic CO2 storage. The project team drilled and cored two test wells to define the sedimentary geologic formations underlying the basin and to document or reach basement rock. With this geologic characterization phase, an integration of seismic, geologic, borehole, and formation core results provided a higher resolution assessment of CO2 storage potential. The Stockton Formation is known to be a potentially favorable geologic storage formation in the basin. In 2011, the 1-NYSTA Tandem Lot stratigraphic test well was drilled to a depth of 6,855 feet in the northern portion of the Newark Basin in southern New York State. Approximately 9 miles south-southeast on the Lamont Doherty Campus, TW-4 was drilled and cored in 2013 to a depth of 1,802 feet and contacted apparent igneous basement at a depth of 1,712 feet. Both wells penetrated the Palisades Sill ranging from 800 feet thick in the eastern well to approximately 1,800 feet in thickness at the 1-NYSTA Tandem Lot deep drill site. A diabase sill can provide an excellent seal and dense confining layer for potential CO2 storage reservoirs and flow layers that are situated beneath it within the Stockton Sandstone. The Stockton Sandstone was encountered beneath the sill in the TW-4 well on the Lamont campus, and data integration suggests that it was likely observed near total depth in the deep 1-NYSTA Tandem Lot well. The test wells confirm and define reservoirs are present beneath the sill and offer CO2 storage potential. The integration of geologic and reservoir characterization of well logs, formation cores, and formation fluids indicated Triassic age lacustrine playa lake and mudbank shales of the Upper Passaic Group can provide an effective seal for the porous and permeable underlying sandstone reservoir layers. This project acquired seismic data, drilled borehole well logs, acquired core samples, and integrated these findings to provide a better understanding of the subsurface geologic formations in the Newark Rift Basin. These findings have contributed to a higher degree of accuracy in predicting potential geologic storage opportunities, while refining geologic storage capacity estimates for the indicated reservoirs and flow units.

.las↗
DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗
DOE OSTI2025

Transportability of exogenous microbial community correlates with interwell connectivity in deep aquifers

Subsurface resource engineering operations often utilize continuous injection of externally-sourced water into geological reservoirs for formation pressure maintenance, resource recovery or energy/waste storage. Such injected water generally contains naturally occurring microbes. Little is known, however, about how the injectate microbes transport through geological media as a community, how such transportability is affected by injector-producer connectivity, and whether such knowledge can be utilized for flowpath characterization. In this study, we analyzed daily-to-weekly timeseries microbial community data from the injected- and produced-fluids of a ten-month flow test at a deep, well-characterized engineered aquifer. We found that the injectate microbial community was distinct from the indigenous community at the amplicon sequence variant (ASV) level, and that the transportability of injectate community towards a given producer, quantified by an “nASV-Overlap” metric we propose, had strong and significant positive correlation with known injector-producer connectivities at our site. This suggests that the better the connectivity, the higher the probability for more injectate species to flow through the interwell region and arrive at a producer. Because interwell connectivity is an important yet usually unknown parameter in subsurface resource engineering, such correlation in turn points to nASV-Overlap as a useful indicator of interwell connectivity for aquifer characterization and long-term monitoring. Based on our findings, an nASV-Overlap-based microbial tracing approach was developed for characterizing and monitoring the relative connectivities across multiple producers with a given injector. A side-by-side comparison between the new nASV-Overlap approach and traditional artificial tracer methods is presented, and their respective strengths and limitations are discussed.

Deep biosphere↗
DOE OSTI2026

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Broadband Rotational Spectroscopy in Uniform Supersonic Flows: Chirped Pulse/Uniform Flow for Reaction Dynamics and Low Temperature Kinetics

ConspectusThe study of gas-phase chemical reactions at very low temperatures first became possible with the development and implementation of the CRESU (French acronym for Reaction Kinetics in Uniform Supersonic Flows) technique. CRESU relies on a uniform supersonic flow produced by expansion of a gas through a Laval (convergent-divergent) nozzle to produce a wall-less reactor at temperatures from 10 to 200 K and densities of 1016-1018 cm-3 for the study of low temperature kinetics, with particular application to astrochemistry. In recent years, we have combined uniform flows with revolutionary advances in broadband rotational spectroscopy to yield an instrument that affords near-universal detection for novel applications in photodissociation, reaction dynamics, and kinetics. This combination of uniform supersonic flows with chirped-pulse Fourier-transform microwave spectroscopy (Chirped-Pulse/Uniform Flow, CPUF) permits detection of any species with a modest dipole moment, thermalized to the uniform temperature of the gas flow, with isomer, conformer, and vibrational state specificity. In addition, the use of broadband, high-resolution, and time-dependent (microsecond time scale) micro- and mm-wave spectroscopy makes it an ideal tool for characterizing both transient and stable molecules, as well as studying their spectroscopy and dynamics.In this Account, we review recent advances made using the CPUF technique, including studies of photodissociation, radical-radical reaction dynamics, and low temperature kinetics. These studies highlight both the strength of universal and multiplexed detection and the challenges of coupling it to a high-density collisional environment. Product branching and product evolution as a function of time have been measured for astrochemically relevant systems, relying on the detailed characterization of these flow conditions via experiments and fluid dynamics simulations. In the photodissociation of isoxazole, an unusual heterocyclic molecule with a very low-energy conical intersection, we have identified 7 products in 5 reaction channels and determined the product branching, pointing to both direct and indirect pathways. We have also approached the same system from separated NO and C3H3 reactants to explore a broader range of the potential energy surface, demonstrating the power of multichannel branching measurements for complex radical-radical reactions. We determined the product branching in the C3H2 isomers in the photodissociation of the propargyl radical and identified the importance of a hydrogen atom catalyzed isomerization to the lowest energy cyclic form. This then motivated a study of direct D-H exchange reaction in radicals, in which we demonstrate that it is an important and overlooked pathway for deuterium fractionation in astrochemical environments. Recently, we have shown the measurement of low temperature kinetics inside an extended Laval nozzle, after which a shock-free secondary expansion to low temperature and density affords an ideal environment for detection by rotational spectroscopy. These results highlight the power and potential of the CPUF approach, and future prospects will also be discussed in light of these developments.

Dias, Nureshan↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Design of Microwave Resonant Cavity Transducer. Development of sensor performance model of microwave cavity flow meter for advanced reactor high temperature fluids

High-temperature fluid reactors, such as molten salt cooled reactors (MSCR) and sodium fast reactors (SFR), are a promising advanced reactor option. Measurement of high-temperature fluid process variables, in particular the flow inside the pressure vessel, is a challenging task because of harsh environment, which includes high radiation, high temperature, and contact with highly corrosive coolant fluid. We are investigating a microwave cavity-based transducer for high-temperature fluid flow sensing. This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. A cavity is characterized by its resonant frequencies. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the flow sensor is evaluated with signal sensitivity using COMSOL computer simulations. A right cylinder geometry stainless steel cavity with dimeter of 0.8in was investigated. We choose membrane thickness of 10mil, so that corrosion anticipated to proceed at the rate of 1mil/year in liquid sodium would affect no more than 10% of the membrane. Using the properties of liquid sodium fluid, and stainless-steel material property values at 500oC, we calculate frequency shift for a range of values of fluid velocity from 0.5m/s to 2m/s. Results of computer simulations indicate measurable sensitivity to flow for this cavity design. Following these simulations, we have developed a preliminary design for fabrication of a transducer operating in microwave K-band for proof-of-principle tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Thermal Hydraulic Experimental Test Article: Second Year of Testing with Secondary Sodium System (Fiscal Year 2025 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is currently installed in the Mechanisms Engineering Test Loop (METL) 28” test vessel #4. Both the primary and secondary sodium systems remain online to facilitate continued testing. This fiscal year, work was performed using a COMSOL Multiphysics magnetohydrodynamic model to characterize flow more accurately in the secondary electromagnetic flowmeters. Experimental campaigns were then performed to study the thermal hydraulic differences between sodium and water as a surrogate fluid in the THETA geometry as well as a study to better characterize and understand temperature oscillations that exist at the outlet of the core to the hot pool. A peer-reviewed article was published in the ASME Journal of Nuclear Engineering and Radiation Science detailing the THETA facility and providing an overview of a test that was performed with the primary and secondary system online [1]. Work continues to develop a database to house experimental THETA data to better facilitate collaboration with industry and laboratory partners for their use of the data for code benchmarking/validation. THETA remains fully operational and is positioned for continued testing in fiscal year 2026.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Predicting Fracture Porosity Evolution in Sandstone (Final Report)

To better understand porosity, strength, chemical reactivity, and patterns of fractures, we developed methods using mineral deposits to unravel how fracture growth and diagenesis interact to create and destroy fracture porosity. Quartz and dolomite cement textures and associated fluid inclusion assemblage sequences and thermal histories provided data to infer the duration and rates at which fractures open. We developed and tested a model that accounts for how fracture porosity and rock properties change as a function of thermal exposure and how fracture size and arrangement evolve. Effects of cement deposits on growth were investigated with new geomechanical models and new methods to quantify fracture spatial arrangement. We tested these concepts in a wide range of rock types using outcrop and core-based datasets. We documented opening histories of fractures and fracture arrays in contrasting tectonic and thermal regimes. Our methods set a new standard for fracture characterization and our diagenetic-geomechanical model accurately predicts fluid flow characteristics of fractured rocks based on coupled effects of diagenesis and deformation. Results represent a significant advance in fundamental understanding of fractures in the Earth and have application to improving subsurface engineering operations in a wide range of contexts including geothermal and sequestration applications. The breakthrough of successfully reconstructing the timing and opening rates of parts of single fractures and portions of fracture arrays to document the growth and linkage of fracture arrays was an accomplishment of this grant that is being followed up in our current BES grant DE-SC0022968, ‘Reconstructing and Predicting Fracture Pattern Evolution‘

58 GEOSCIENCES↗
DOE OSTI2023

Observation of hydrodynamic plasmons and energy waves in graphene

Thermally excited electrons and holes form a quantum-critical Dirac fluid in ultraclean graphene and their electrodynamic responses are described by a universal hydrodynamic theory. The hydrodynamic Dirac fluid can host intriguing collective excitations distinctively different from those in a Fermi liquid. Here we report the observation of the hydrodynamic plasmon and energy wave in ultraclean graphene. In this study, we use the on-chip terahertz (THz) spectroscopy technique to measure the THz absorption spectra of a graphene microribbon as well as the propagation of the energy wave in graphene close to charge neutrality. We observe a prominent high-frequency hydrodynamic bipolar-plasmon resonance and a weaker low-frequency energy-wave resonance of the Dirac fluid in ultraclean graphene. The hydrodynamic bipolar plasmon is characterized by the antiphase oscillation of massless electrons and holes in graphene. The hydrodynamic energy wave is an electron-hole sound mode with both charge carriers oscillating in phase and moving together. The spatial–temporal imaging technique shows that the energy wave propagates at a characteristic speed of $V_f$ / $ \sqrt2$ near the charge neutrality. Our observations open new opportunities to explore collective hydrodynamic excitations in graphene systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Relative permeability of gas and water flow in hydrate-bearing porous media: A micro-scale study by lattice Boltzmann simulation

The water-gas relative permeability is an important parameter to characterize multiphase flow in sediments. To study the water-gas relative permeability of hydrate-bearing porous media, multiphase flow simulations were carried out at the pore scale using the lattice Boltzmann method. In this work, the effects of hydrate saturation and hydrate-growth habits on the water-gas relative permeability, which is scaled by the relative permeability considering the hydrate only, were evaluated in a two-dimensional porous medium. Results show that the increase of hydrate saturation causes the decrease of water-gas effective permeability as expected. However, the effect of hydrate saturation on the water-gas relative permeability is different from that of hydrate saturation on the water-gas effective permeability. The water-gas relative permeability increases with the increase of hydrate saturation in the pore-filling case. The water-gas relative permeability decreases with the increase of hydrate saturation in the grain-coating case. The wettability of solid phase has a different effect on the relative permeability of wetting phase and nonwetting phase. The Jamin effect (phase blocking) was observed and may exist in the production of gas from natural gas hydrate reservoirs. This seriously affects the multiphase flow characteristics. The changes of microscale fluid distribution effect the changes of water-gas relative permeability. The relationship between the water-gas relative permeability and the characterization parameters of microscale fluid distribution was analyzed.

58 GEOSCIENCES↗
DOE OSTI2021

A chemical kinetic mechanism for combustion and flame propagation of CH 2 F 2 /O 2 /N 2 mixtures

Here in this work, we compiled and evaluated rate expressions for reactions relevant to the decomposition and combustion of CH 2 F 2 (difluoromethane, refrigerant R-32) in CH 2 F 2 /O 2 /N 2 flames. The recommended values have been used in premixed flame calculations, reported elsewhere, to model experimentally derived burning velocities determined using a constant volume spherical flame method. In this work, we also provide a detailed description of the reaction pathways for decomposition and combustion of CH 2 F 2 . This work is part of a larger effort at NIST to characterize and predict the flammability of new refrigerant working fluids and their blends for consideration as replacements of current refrigerants with high global warming potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗