Search NASA⌕ Search

SEARCH · Search NASA

Results for “Injection testing”

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 1,477 records · Page 82

Simulations of Multi-Mode Combustion Regimes Realizable In a Gasoline Direct Injection Engine

Lean and dilute gasoline compression ignition (GCI) operation in spark ignition (SI) engines are an attractive strategy to attain high fuel efficiency and low NOx levels. However, this combustion mode is often limited to low-load engine conditions due to the challenges associated with autoignition controllability. In order to overcome this constrain, multi-mode (MM) operating strategies, consisting of advanced compression ignition (ACI) at low load and conventional SI at high load, have been proposed. In this three-dimensional computational fluid dynamics study, the concept of multi-mode combustion using two RON98 gasoline fuel blends (Co-Optima Alkylate and E30) in a gasoline direct injection (GDI) engine were explored. To this end, a new reduced mechanism for simulating the kinetics of E30 fuel blend is introduced in this study. To cover the varying engine load demands for multi-mode engines, primary combustion dynamics observed in ACI and SI combustion modes was characterized and validated against experimental measurements. In order to implement part-load conditions, a strategy of mode transition between SI and ACI combustion (i.e., mixed-mode combustion) was then explored numerically by creating a virtual test condition. The results obtained from the mixed-mode simulations highlight an important feature that deflagrative flame propagation regime coexists with ignition-assisted end-gas autoignition. This study also identifies a role of turbulent flow property adjacent to premixed flame front in characterizing the mixed-mode combustion. The employed hybrid combustion model was verified to perform simulations aiming at suitable range of multi-mode engine operations.

Kim, Sayop↗

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management Facility (EMF) stack monitor location was qualified using a combination of scale model stacks to mitigate the risk of identifying that the sampling location does not meet the qualification criteria on the full-scale stack. The LV-S1 scale model stack was used as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling locations were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. The LV-S1, LB-S1, and LV-S2 scale model stack tests have met the criteria of the ANSI/HPS N13.1-1999 standard to demonstrate the stack sampling locations are well mixed. Verification tests of the EMF stack were performed at normal operating conditions. The minimum 1/6 DV value and the maximum 6 DV value from the scale model testing determine the range of stack flow rates for which the full-scale stack may be operated while remaining in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rate along with the average and maximum qualified stack flow rate based on the LV-S1 scale model DV values. The operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were <5°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, all of which were within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all planned fan operating configurations. This includes each combination of ACV fans with DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Microbiological-enhanced mixing across scales during in-situ bioreduction of metals and radionuclides at Department of Energy Sites

Bioreduction is being actively investigated as an effective strategy for subsurface remediation and long-term management of DOE sites contaminated by metals and radionuclides (i.e. U(VI)). These strategies require manipulation of the subsurface, usually through injection of chemicals (e.g., electron donor) which mix at varying scales with the contaminant to stimulate metal reducing bacteria. There is evidence from DOE field experiments suggesting that mixing limitations of substrates at all scales may affect biological growth and activity for U(VI) reduction. Although current conceptual models hold that biomass growth and reduction activity is limited by physical mixing processes, a growing body of literature suggests that reaction could be enhanced by cell-to-cell interaction occurring over length scales extending tens to thousands of microns. Our project investigated two potential mechanisms of enhanced electron transfer. The first is the formation of single- or multiple-species biofilms that transport electrons via direct electrical connection such as conductive pili (i.e. ‘nanowires’) through biofilms to where the electron acceptor is available. The second is through diffusion of electron carriers from syntrophic bacteria to dissimilatory metal reducing bacteria (DMRB). The specific objectives of this work are (i) to quantify the extent and rate that electrons are transported between microorganisms in physical mixing zones between an electron donor and electron acceptor (e.g. U(IV)), (ii) to quantify the extent that biomass growth and reaction are enhanced by interspecies electron transport, and (iii) to integrate mixing across scales (e.g., microscopic scale of electron transfer and macroscopic scale of diffusion) in an integrated numerical model to quantify these mechanisms on overall U(VI) reduction rates. We tested these hypotheses with five tasks that integrate microbiological experiments, unique micro-fluidics experiments, flow cell experiments, and multi-scale numerical models. Continuous fed-batch reactors were used to derive kinetic parameters for DMRB, and to develop an enrichment culture for elucidation of syntrophic relationships in a complex microbial community. Pore and continuum scale experiments using microfluidic and bench top flow cells were used to evaluate the impact of cell-to-cell and microbial interactions on reaction enhancement in mixing-limited bioactive zones, and the mechanisms of this interaction. Some of the microfluidic experiments were used to develop and test models that considers direct cell-to-cell interactions during metal reduction. Pore scale models were incorporated into a multi-scale hybrid modeling framework that combines pore scale modeling at the reaction interface with continuum scale modeling. New computational frameworks for combining continuum and pore-scale models were also developed

54 ENVIRONMENTAL SCIENCES↗

SU-F-R-34: Quantitative Perfusion Measurement in Rectal Cancer Using Three Different Pharmacokinetic Models: Implications for Prospective Study Design

Purpose: To compare three different pharmacokinetic models for analysis of dynamic-contrast-enhanced (DCE)-CT data with respect to different acquisition times and location of region of interest. Methods: Eight rectal cancer patients with pre-treatment DCE-CTs were included. The dynamic sequence started 4–10seconds(s) after the injection of contrast agent. The scan included a 110s acquisition with intervals of 40×1s+15×3s+4×6s. An experienced oncologist outlined the tumor region. Hotspots with top-5%-enhancement were also identified. Pharmacokinetic analysis was performed using three different models: deconvolution method, Patlak model, and modified Toft’s model. Perfusion parameters as blood flow (BF), blood volume (BV), mean transit time (MTT), permeability-surface-area-product (PS), volume transfer constant (Ktrans), and flux rate constant (Kep), were compared with respect to different acquisition times of 45s, 65s, 85s and 105s. Both hotspot and whole-volume variances were also assessed. The differences were compared using the Wilcoxon matched-pairs test and Bland-Altman plots. Results: Moderate correlation was observed for various perfusion parameters (r=0.56–0.72, p<0.0001) but the Wilcoxon test revealed a significant difference among the three models (P < .001). Significant differences in PS were noted between acquisitions of 45s versus longer time of 85s or 105s (p<0.05) using Patlak but not with the deconvolution method. In addition, measurements varied substantially between whole-volume vs. hotspot analysis. Conclusion: The radiation dose of DCE-CT was on average 1.5 times of an abdomen/pelvic CT, which is not insubstantial. To take the DCE-CT forward as a biomarker in oncology, prospective studies should be carefully designed with the optimal image acquisition and analysis technique. Our study suggested that: (1) different kinetic models are not interchangeable; (2) a 45s acquisition might not be sufficient for reliable permeability measurement in rectal cancer using Patlak model, but might be achievable using deconvolution method; and (3) local variations existed inside the tumor, and both whole-volume-averaged and local-heterogeneity analysis is recommended for future quantitative studies. This work is supported by the National High-tech R&D program for Young Scientists by the Ministry of Science and Technology of China (Grant No. 2015AA020917), Natural Science Foundation of China (NSFC Grant No. 81201091).

60 APPLIED LIFE SCIENCES↗

DRIFT-INDUCED DECELERATION OF SOLAR ENERGETIC PARTICLES

We investigate the deceleration of solar energetic particles (SEPs) during their propagation from the Sun through interplanetary space, in the presence of weak to strong scattering in a Parker spiral configuration, using relativistic full orbit test particle simulations. The calculations retain all three spatial variables describing particles’ trajectories, allowing us to model any transport across the magnetic field. Large energy change is shown to occur for protons, due to the combined effect of standard adiabatic deceleration and a significant contribution from particle drift in the direction opposite to that of the solar wind electric field. The latter drift-induced deceleration is found to have a stronger effect for SEP energies than for galactic cosmic rays. The kinetic energy of protons injected at 1 MeV is found to be reduced by between 35% and 90% after four days, and for protons injected at 100 MeV by between 20% and 55%. The overall degree of deceleration is a weak function of the scattering mean free path, showing that, although adiabatic deceleration plays a role, a large contribution is due to particle drift. Current SEP transport models are found to account for drift-induced deceleration in an approximate way and their accuracy will need to be assessed in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

WEST full tungsten operation with an ITER grade divertor

The mission of WEST (tungsten-W Environment in Steady-state Tokamak) is to explore long pulse operation in a full tungsten (W) environment for preparing next-step fusion devices (ITER and DEMO) with a focus on testing the ITER actively cooled W divertor in tokamak conditions. Following the successful completion of phase 1 (2016-2021), phase 2 started in December 2022 with the lower divertor made entirely of actively cooled ITER-grade tungsten mono-blocks. A boronization prior the first plasma attempt allowed for a smooth startup with the new divertor. Despite the reduced operating window due to tungsten, rapid progress has been made in long pulse operation, resulting in discharges with a pulse length of 100 s and an injected energy of around 300 MJ per discharge. Plasma startup studies were carried out with equatorial boron nitride limiters to compare them with tungsten limiters, while Ion Cyclotron Resonance Heating assisted startup was attempted. High fluence operation in attached regime, which was the main thrust of the first campaigns, already showed the progressive build up of deposits and appearance of dust, impacting the plasma operation as the plasma fluence increased. In total, the cumulated injected energy during the first campaigns reached 43 GJ and the cumulated plasma time exceeded 5 h. Demonstration of controlled X-Point Radiator regime is also reported, opening a promising route for investigating plasma exhaust and plasma-wall interaction issues in more detached regime. This paper summarises the lessons learned from the manufacturing and the first operation of the ITER-grade divertor, describing the progress achieved in optimising operation in a full W environment with a focus on long pulse operation and plasma wall interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First Ever Field Pilot on Alaska's North Slope to Validate the Use of Polymer Floods for Heavy Oil EOR a.k.a Alaska North Slope Field Laboratory (ANSFL)

Alaska’s high viscosity oil resources that range between 20–30+ billion barrels represent about a third of known North Slope original oil in place (OOIP). These resources are primarily concentrated in the Schrader Bluff formation (also called West Sak on the Western North Slope) and Ugnu reservoirs and are categorized as “viscous oils” and “heavy oils” owing to their in-situ viscosities between 5–10,000 cP and up to a million+ cP respectively. The viscous oil deposits are relatively deeper (2,000 – 5,000 ft), whereas the heavy oils are somewhat shallower (2,000 – 4,000 ft). The typically shallow depths and the proximity to the continuous permafrost results in relatively lower formation temperatures and pressures, and consequently higher viscosities. The vertical depth vs. viscosity delineated in Paskvan et al. (2016) differentiates the viscous and heavy oils. As depicted in Paskvan et al. (2016), currently the main focus (referred to as “developing”) is on the viscous oils in the Schrader Bluff formation in the Milne Point Unit (MPU). Notwithstanding this Alaska North Slope (ANS) specific categorization, we use the industry adopted, all-inclusive term “heavy oil” for all high viscosity oils. Resource characterization and additional details can be found in topical publications of Paskvan et al. (2016) and Targac et al. (2005). Despite the vast resource base, the development pace, vis-à-vis the production of heavy oils has been very slow and limited due to multiple factors such as cost, logistics, challenging arctic environment, poor waterflood sweep efficiency due to mobility contrasts, and significantly high minimum miscibility pressures (MMP). Most importantly, typical or standard thermal methods that are commonplace elsewhere (Canada, California) are inapplicable due to the continuous permafrost. As a consequence, cumulative production of heavy and viscous oils is a little over 1% of OOIP slope wide and currently, there is hardly any production from Ugnu. However, on a broader level, these unfavorable factors are outweighed by the fact that (1) these resources, within the established infrastructure, are too large to ignore because of their strategic importance to the Nation and the State of Alaska and (2) Prudhoe Bay type diluent crude oil is still available for heavy oil transport through the Trans Alaska Pipeline System (TAPS). Similarly, from a reservoir standpoint, the following factors also are important offsets: (1) favorable rock characteristics of Schrader Bluff; (2) the promise demonstrated by the initial scoping studies (Seright 2010, 2011) suggesting significant increase of heavy oil recovery using polymer flooding; (3) successful field implementation in Canada, China and elsewhere in the world, and (4) availability of the existing pairs of horizontal injector-producer in Schrader Bluff The foregoing was recognized as the best readily available opportunity for significant investment by the US Department of Energy and the field operator Hilcorp Alaska LLC to conduct the first ever field scale experiment to test the polymer flooding technology to unlock the vast heavy oil resources on ANS. With this primary goal in mind, the research team embarked on a ~4.5 years long project that focused on the field polymer pilot complemented by supporting laboratory and simulation studies. As documented in this final report, over the course of the project, many lessons have been learned and valuable field and supporting laboratory data has been collected, which also is complemented by numerical reservoir simulations. We have been able to establish the injectivity of polymer solution, evidence of significant reduction in the water cut of previously waterflooded pattern, effective propagation of a hydrolyzed polyacrylamide (HPAM), benefits of low salinity water, provide practical guidance on handling of produced fluids containing breakthrough polymer, fit-for-purpose forecast-worthy history matched simulation model, polymer EOR benefit of 700-1000 bopd over waterflood, and most importantly a low polymer utilization factor of ~1.7 lb/stb. In summary this project is deemed as a scientific, technical and economic success, having met all objectives, fulfilled deliverables and within budget, providing impetus to apply polymer EOR throughout the Milne Point Field paving the way for even heavier viscosity oils in the Ugnu area, eventually extending the economic life of TAPS.

02 PETROLEUM↗

A research program to measure the lifetime of spin polarized fuel

The use of spin polarized fuel could increase the deuterium-tritium (D-T) fusion cross section by a factor of 1.5 and, owing to alpha heating, increase the fusion power by an even larger factor. Issues associated with the use of polarized fuel in a reactor are identified. Theoretically, nuclei remain polarized in a hot fusion plasma. The similarity between the Lorentz force law and the Bloch equations suggests polarization can be preserved despite the rich electromagnetic spectrum present in a magnetic fusion device. The most important depolarization mechanisms can be tested in existing devices. The use of polarized deuterium and 3 He in an experiment avoids the complexities of handling tritium, while encompassing the same nuclear reaction spin-physics, making it a useful proxy to study issues associated with full D-T implementation. 3 He fuel with 65% polarization can be prepared by permeating optically-pumped 3 He into a shell pellet. Dynamically polarized 7 Li-D pellets can achieve 70% vector polarization for the deuterium. Cryogenically-frozen pellets can be injected into fusion facilities by special injectors that minimize depolarizing field gradients. Alternatively, polarized nuclei could be injected as a neutral beam. Once injected, the lifetime of the polarized fuel is monitored through measurements of escaping charged fusion products. Multiple experimental scenarios to measure the polarization lifetime in the DIII-D tokamak and other magnetic-confinement facilities are discussed, followed by outstanding issues that warrant further study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing Signal-Based Inertia and Frequency Response Estimation for Power Systems with High Penetration of Inverter-Based Resources: Preprint

Power system inertia is the inherent capability of a power system to resist changes in its frequency during disturbances. Real-time inertia estimation technology has become more important due to the low-inertia issues caused by the increasing integration levels of inverter-based resources (IBRs) from renewable energy; however, existing inertia estimation methods hardly consider multiple frequency response controls that act within the same time frame as conventional inertial response, thus making measured inertia values vary under different testing conditions. To resolve this issue, this paper proposes a novel real-time estimation method to simultaneously estimate a power system's inertia constant and frequency response droop constant using a well-designed probing signal. First, we formulate the inertia and frequency response model of a power system with IBRs. Second, through the integration and manipulation of the developed model, we propose a multivariate linear regression- based estimation method that is resilient to measurement noise. Third, we design a probing signal that can be injected by IBRs to incite the required transients for estimation. Finally, we validate the proposed estimation method through comprehensive power- hardware-in-the-loop experiments using inverter hardware and a realistic island power system model. The results demonstrate that the proposed method can accurately estimate the inertia and droop value of the power system with grid-following IBRs and grid-forming IBRs with virtual synchronous machine control.

frequency response↗

Accelerator Stewardship Test Facility Program - Elliptical Twin Cavity for Accelerator Applications

Funding is being requested pursuant to the proposals entitled Elliptical Twin Cavity for Accelerator Applications that was submitted and reviewed through the Portfolio Analysis and Management System (PAMS). The PAMS proposal identifier number is 0000219731. The proposed new type of superconducting cavity, the Elliptical Twin Cavity, is capable of accelerating or decelerating beams in two separate beam pipes. This configuration is particularly effective for high-current, low energy electron beams that will be used for bunched beam cooling of high-energy protons or ions. Having the accelerated beam physically separated from the decelerated beam, but interacting with the same RF mode, means that the low energy beam from the gun can be injected into to the superconducting cavity without bends enabling a small beam emittance to be maintained. A staff engineer who has been working with non-standard complicated cavity structures replaces the senior engineer (in the original budget) who is moving on to be a project leader. This is reflected in a slightly increased engineer time and in reduced costs. The Indirect costs for FY16 are lower than the previous projection. As a result, there is no scope reduction.

99 GENERAL AND MISCELLANEOUS↗

Direct comparison of gamma, electron beam and X-ray irradiation doses on characteristics of low-density polyethylene, polypropylene homopolymer, polyolefin elastomer and chlorobutyl rubber medical device polymers

There is a growing need for increased efficiency in the sterilization of single use medical devices and other products that contain polymer components. Gamma radiation is widely used for devices suited for radiation sterilization; however, safety, throughput and cobalt-60 source availability are challenging the prospect of relying on gamma radiation to meet the anticipated needs of the industry. Use of electron beam (e-beam) and X-rays as alternatives to gamma for radiation sterilization have been hampered in part by a concern that these modalities may adversely affect polymer integrity and performance relative to the gamma method, for which the industry has had much more experience. To address this concern, the effects of sterilization-relevant doses of ebeam, X-ray and gamma radiation were directly compared using common medical device polymers found in two prototypical commercial devices currently sterilized using cobalt-60 gamma irradiation. The Becton, Dickinson and Company (BD) Vacutainer™ Plus tube contains low-density polyethylene and chlorobutyl rubber components, while the BD Vacutainer™ Push Button Blood Collection Set contains polypropylene homopolymer and polyolefin elastomer components. Injection-molded samples prepared from the polymers used in these products were exposed to target doses of 15, 35, 50 and 80 kGy using gamma, e-beam and X-ray radiation. Changes in coloration, tensile properties and hardness were measured for each condition, and the effects of e-beam and X-ray irradiation compared with the effects of gamma irradiation on these properties. Both e-beam and X-ray appear as viable alternatives to gamma irradiation for sterilization of the polymers tested.

sterilization, gamma, X-ray, electron beam, radiat↗

Probing Signal-Based Inertia and Frequency Response Estimation for Power Systems With High Levels of Inverter-Based Resources

Power system inertia is the inherent capability of a power system to resist changes in its frequency during disturbances. Real-time inertia estimation technology has become more important due to the low-inertia issues caused by the increasing integration levels of inverter-based resources (IBRs) from renewable energy; however, existing inertia estimation methods hardly consider multiple frequency response controls that act within the same time frame as conventional inertial response, thus making measured inertia values vary under different testing conditions. To resolve this issue, this paper proposes a novel real-time estimation method to simultaneously estimate a power system's inertia constant and frequency response droop constant using a well-designed probing signal. First, we formulate the inertia and frequency response model of a power system with IBRs. Second, through the integration and manipulation of the developed model, we propose a multivariate linear regression-based estimation method that is resilient to measurement noise. Third, we design a probing signal that can be injected by IBRs to incite the required transients for estimation. Finally, we validate the proposed estimation method through comprehensive power-hardware-in-the-loop experiments using inverter hardware and a realistic island power system model. The results demonstrate that the proposed method can accurately estimate the inertia and droop value of the power system with grid-following IBRs and grid-forming IBRs with virtual synchronous machine control.

frequency response↗

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗

A Complete Machine-Learning-Based Workflow to Illuminate Earthquake Processes

Under this grant we developed, tested, and made available, machine learning models to improve the tasks in the earthquake monitoring workflow (Figure 1). We implemented these models as a part of an end-to-end workflow for seismic network processing and demonstrated, in a variety of settings, that these methods generalize and that they result in dramatically more comprehensive earthquake catalogs. These catalogs illuminate earthquake processes in detail and to an extent that had previously not been possible, and they do so for both tectonic seismicity and seismicity induced by fluid injection related to unconventional hydrocarbon development. This report summarizes the results from the 15 publications that resulted from this grant. Those contributes are divided into: (1) the development of specific tasks related to monitoring (7 publications), (2) the organization of those tasks into workflows for seismic monitoring (2 publications), and (3) applications of those workflows to data (6 publications).

58 GEOSCIENCES↗

Advanced Gasoline Turbocharged Direction Injection (GTDI) Engine Development

This program was undertaken in response to US Department of Energy Solicitation DE-FOA-0000079, resulting in a cooperative agreement with Ford and MTU to demonstrate improvement of fuel efficiency in a vehicle equipped with an advanced GTDI engine. Ford Motor Company has invested significantly in GTDI engine technology as a cost effective, high volume, fuel economy solution, marketed globally as EcoBoost technology. Ford envisions additional fuel economy improvement in the medium and long term by further advancing EcoBoost technology. The approach for the project was to engineer a comprehensive suite of gasoline engine systems technologies to achieve the project objectives, and to progressively demonstrate the objectives via concept analysis / computer modeling, single-cylinder and multi-cylinder engine testing on engine dynamometer, and vehicle level testing on chassis rolls.

33 ADVANCED PROPULSION SYSTEMS↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

Molybdenum isotope fractionation and speciation in a euxinic lake—Testing ways to discern isotope fractionation processes in a sulfidic setting

The molybdenum (Mo) isotope composition in euxinic shales has been used as a proxy for the global distribution of anoxic conditions in ancient oceans, and since more recently also as a proxy for sulfide concentrations in depositional environments. However, there is currently no way to distinguish isotope fractionation at low bottom water sulfide concentrations in ‘local’ basins from ‘global’ secular isotope variations associated with changing seawater composition. This uncertainty is challenging the use of Mo isotopes for paleoceanographic reconstructions. To explore this further, we present new data from sediments deposited over the past ~ 9800 years in one of the best studied euxinic localities in the world: Lake Cadagno in Switzerland. The sample set allows us to test ways to discern isotope fractionation processes at play in a highly restricted euxinic basin. Most of our drill core samples (n = 18) show high δ98Mo values similar to previously studied shallow sediments, indicative of quantitative Mo removal from the water column (Dahl et al. 2010a). However, a few samples (n = 3) deposited between about 1200 and 3400 years ago carry low δ98Mo values and have been isotopically fractionated in the lake. Sedimentological and geochemical characterizations show that these δ98Mo-fractionated sediments formed during times of frequent injection of O2- and sediment-rich river water into the deep sulfidic water column. A positive correlation between δ98Mo and sedimentary Mo contents suggests that isotope fractionation occurred during times of non-quantitative Mo removal, although Mn-oxide cycling at the chemocline might also contribute a subordinate proportion of (98Mo-depleted) molybdenum into the sulfidic zone. Sedimentary Mo/U enrichments relative to oxic lake water further supports the hypothesis that a particulate Mo shuttle was most efficient during times of quantitative Mo removal. Therefore, periods with inefficient Mo capture are ascribed to incomplete conversion of molybdate to particle reactive Mo species when bottom water H2S levels were low or less stable than today. Using XAFS spectroscopy, we found that the two distinct Mo compounds predominating in the sediments (MoIV-S and MoVI-OS) are not diagnostic for isotope fractionation that has occurred in Lake Cadagno. Instead, we infer that δ98Mo-fractionated products (forming via a low-sulfide Mo pathway) can be subsequently altered with little or no isotopic imprint during remobilization and re-precipitation (e.g., at higher sulfide levels in the sediments) as well as during post-depositional oxidation. Future work could investigate local δ98Mo-fractionation processes expressed in other euxinic settings and explore other sedimentary metrics to constrain the steps involved in the euxinic burial pathway(s). One tantalizing prospect of this is to distinguish between local bottomwater sulfide levels and variations in the fraction of global seafloor anoxia from the Mo isotope composition in ancient euxinic mudrocks.

58 GEOSCIENCES↗