Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Impact 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ionization Yield in Silicon for eV-Scale Electron-Recoil Processes

The development of single charge resolving, macroscopic silicon detectors has opened a window into rare processes at the O(eV) scale. In order to reconstruct the energy of a given event, or model the charge signal obtained for a given amount of energy absorbed by the electrons in a detector, an accurate charge yield model is needed. In this paper we review existing measurements of charge yield in Silicon, focusing in particular on the region below 1 keV. We highlight a calibration gap between 12-50 eV (referred to as the "UV-gap") and employ a phenomenological model of impact ionization to explore the likely charge yield in this energy regime. Finally, we explore the impact of variations in this model on a test case, that of dark matter scattering off electrons, to illustrate the scientific impact of uncertainties in charge yield.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Visual Systems Mapping to Define and Compare Woody Biomass LCAs for Sustainable Systems

The challenge addressed in this research centres on the need to choose between several biomass sources and energy production processes, while supporting rural economies and resilience of forest systems. A key barrier to effective decision-making for strategies using biomass is the lack of standardized and transparent life cycle assessment (LCA) baselines. These baselines are critical for assessing the impacts of biomass strategies but often vary due to regional factors and chosen simplifying assumptions of the LCAs. However, omitting key variables can mean the LCA omits key feedback and balancing loops relevant to fully assessing impacts of the change or test scenario. To address these complexities, this project employs a systems engineering approach: visual systems mapping. This technique is used to define the boundaries and dynamic behaviours of LCA baselines, enhancing transparency. By examining five literature sources and their documented baseline scenarios, the systems mapping case-studies demonstrates an approach to documenting and archiving these baselines. Recommendations are that visual systems mapping should be used to document key assumptions, such as baselines, of LCAs. Further, where possible open data repositories should hold key information about LCA baselines and reproducible workflows (e.g., using open-source tools) should be used to improve transparency and comparability in LCAs. Given the consensus within the broader scientific community on the importance of replicable data practices, this research reinforces the need for standardized frameworks and systems engineering tools in LCAs. This research demonstrates a pathway to more transparent, standardized, and comparable LCAs, that may bolster decisions for biomass systems.

Davis, Maggie [ORNL] (ORCID:0000000181319328)↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2021

Production of desensitized, ultrafine PETN powder

Pentaerythritol tetranitrate (PETN) is a widely studied high explosive (HE), most commonly used in detonator applications. In this work, we use a spray drying technique to manufacture a new “nano-powder” morphology of PETN. The spray dried PETN is several orders of magnitude smaller in particle size than traditionally prepared PETN powders. The spray dried PETN has a mostly spherical and smooth morphology, while traditional crystalline PETN typically has sharp edges and crystal facets. Small-scale sensitivity tests including drop-weight impact, friction, and electrostatic discharge (ESD) indicate the spray dried PETN is less sensitive than traditionally used forms of PETN powders. Furthermore, we also observed no changes in chemical properties (melt temperature and onset of decomposition) after spray drying, and the material remains in the tetragonal, crystalline phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Seasonality and Albedo Dependence of Cloud Radiative Forcing in the Upper Colorado River Basin

Mountains create and enhance their own clouds, which both scatter and absorb shortwave radiation from the sun and absorb and re-emit land surface and atmospheric longwave radiation. However, the impacts of clouds on the surface radiation balance in high elevation snowy mountain terrain are poorly explored. In this study, we use data collected by the SAIL field campaign and partner organizations in the upper elevations (2,880 m.a.s.l) of the Upper Colorado River Basin (UCRB) over a 21-month period from September 2021 to June 2023 to estimate Cloud Radiative Forcing (CRF) in the shortwave, longwave, and the net effect. Longwave warming effects dominate during the winter when snow albedos are high (0.8–0.9) and the background atmospheric precipitable water vapor is low (<0.5 cm), yielding a maximum monthly average net CRF of +34.7 W·m -2 , meaning that clouds increase the net radiation relative to clear skies during this time period. The sign of net CRF switches in the warm season as snow recedes, sun-angles increase, and the North American monsoon arrives, yielding a minimum monthly average net CRF of -47.6 W·m -2 with hourly minima of -600 W·m -2 . The sign of net CRF is typically positive, even at solar noon, when the surface is snow covered, except for a brief period over melting, low-albedo snow (0.5–0.6) impacted by dust impurities. Sensitivity tests elucidate the role of the surface albedo on the net CRF. The results suggest that net CRF will increase in magnitude and lead to a more persistent cooling effect on the surface net radiation budget as the snow cover declines.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Threadsafe Dynamic Neighbor Lists for Monte Carlo Ray Tracing

Monte Carlo (MC) transport codes offer high-fidelity modeling of particle transport physics, but their high computational cost makes them impractical for many applications. For some applications such as multiphysics and depletion that use finely discretized geometries, a large portion of this computational cost is attributable to ray tracing. Neighbor lists are a well-known method for accelerating ray-tracing calculations in a MC code, but despite their prevalence, little work has been published on the details of their implementation. The fine details can have a significant impact on performance, particularly when using shared-memory parallelism. This paper addresses these details of implementation with a discussion of different neighbor list schemes and their impact on software runtime. Performance tests were run by using OpenMC on a pin-cell problem discretized with up to 200 axial regions. The results demonstrate that switching from surface-based to cell-based neighbor lists leads to a 10 faster calculation rate for the most fine discretization. Finally, using a threadsafe shared-memory data structure results in a 20% faster calculation rate versus simple threadprivate neighbor lists. Results here show that a data structure that is contiguous in memory improves performance by only 1% to 2% over noncontiguous linked lists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Power-to-Gas Systems for Active Load Management at EV Charging Sites with High DER Penetration

Intermittent renewable systems and increasing electric vehicle (EV) penetrations increase load intermittency—and will require expensive distribution system infrastructure upgrades. Power-to-gas (P2G) systems are highly controllable loads that can be modulated to regulate the net load behind the meter. In this work, we present an active load management approach using P2G systems to regulate the net load for an EV charging station with high penetrations of renewable generation. The system also validates the feasibility of the proposed approach using a 750-kW electrolyzer operated with the distribution system in a real-time power-hardware-in-the-loop test. Finally, we evaluate the impact of the size of the P2G system on its ability to regulate the net load.

ADVANCED PROPULSION SYSTEMS,POWER TRANSMISSION AND↗
DOE OSTI2021

Investigation of Surface Skidding Response of PBX 9501 and PBX 9502

Introduction: Assess the hazards of handling polymer bonder explosives(PBX) for the case of skidding with grit particles in the frictional interface. High explosive (HE) materials tested: PBX 9501 and 9502. High melting point materials at the skidding interface between HE and a substrate is known to increase the probability of ignition by creating hot spots not subject to reduced heating from the melting of HE material. Skidding occurred at peak forces of 6-26 kN normal force, and 0.6-3.4 kN of frictional force. Velocities for skidding range from 5-30 m/s, see figure 7. Sparse amounts of 300 μm silica grit is applied loosely to the HE surface before testing. Previous works have tested ignition by grit. Testing parameters include angle of impact, drop height, grit size, grit density, substrate material. This study focuses on pre ignition conditions and detailed thermal imaging of the frictional interface. Conclusions: Temperature hot spots were observed at high points in HE surface topology. These high points may be created by grit particles. Since these points impact first high frictional work is concentrated and quickly increases temperature. Small amounts of grit do not affect the overall frictional force but do concentrate the frictional work to create hot spots.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical durability assessment of enhanced low-activity waste glasses through EPA method 1313

In this work, we report the progress of the Glass Leaching Assessment for Durability (GLAD) program on the implementation of the United States Environmental Protection Agency (EPA) Leaching Environmental Assessment Framework pH-dependent leach test (EPA Method 1313) to low-activity nuclear waste (LAW) glasses. The GLAD program seeks to develop new strategies to understand the chemical durability of nuclear waste glasses for the disposal in near-surface conditions. A series of 16 high-waste loading LAW glasses, currently under development, were selected using machine learning methods to study the corrosion behavior using EPA Method 1313. Reacted glass powders were examined using scanning electron microscopy and the eluate compositions were examined using inductively coupled plasma-optical emission spectroscopy. Compositional modeling was used to fit the measured elemental releases from EPA Method 1313. The compositional models demonstrated that elements such as Si reduce elemental release while B can increase elemental release (consistent with elemental modeling of the Product Consistency Test and Vapor Hydration Test) while other elements, such as Fe, exhibit pH-dependent behavior. The amount of acid added during the EPA testing was found to significantly impact the observed result, which was only apparent after preforming the present matrix study. The overall titration curves were able to be compositionally modeled for future process optimization.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Wake steering optimization under uncertainty

Turbines in wind power plants experience significant power losses when wakes from upstream turbines affect the energy production of downstream turbines. A promising plant-level control strategy to reduce these losses is wake steering, where upstream turbines are yawed to direct wakes away from downstream turbines. However, there are significant uncertainties in many aspects of the wake steering problem. For example, infield sensors do not give perfect information, and inflow to the plant is complex and difficult to forecast with available information, even over short time periods. Here, we formulate and solve an optimization under uncertainty (OUU) problem for determining optimal plant-level wake steering strategies in the presence of independent uncertainties in the direction, speed, turbulence intensity, and shear of the incoming wind, as well as in turbine yaw positions. The OUU wake steering strategy is first examined for a two-turbine test case to explore the impacts of different types of inflow uncertainties, and it is then demonstrated for a more realistic 11-turbine wind power plant. Of the sources of uncertainty considered, we find that wake steering strategies are most sensitive to uncertainties in the wind speed and direction. When maximizing expected power production, the OUU strategy also tends to favor smaller yaw angles, which have been shown in previous work to reduce turbine loading. Ultimately, the plant-level wake steering strategy formulated using an OUU approach yields 0.48?% more expected annual energy production for the 11-turbine wind plant than a strategy that neglects uncertainty when considering stochastic inputs. Thus, not only does the present OUU strategy produce more power in realistic conditions, but it also reduces risk by prescribing strategies that call for less extreme yaw angles.

17 WIND ENERGY↗
DOE OSTI2023

Heliostat Consortium: Gap Analysis on Wind Load for Achieving a Fully Competitive Heliostat Industry

The Heliostat Consortium (HelioCon) Wind Load Subtask was initiated with the aim of bringing research work pertaining to wind load measurement, characterization, and prediction taking place across several tasks, including Advanced Manufacturing, Components and Controls, and Field Deployment. The cross-cutting wind load subtopic in the HelioCon roadmap report [1] highlighted standardized methods and tools that are needed for a more detailed understanding of the static and dynamic loads on a heliostat. This will enable cost reduction of wind-dependent heliostat components to avoid unnecessarily conservative, overly constrained designs and increase field efficiency/reliability, to reduce the risk of component failures due to high-wind events (>15 m/s). Gaps related to heliostat wind load include site characterization for wind measurements, critical load cases for heliostat design, turbulence impacts on heliostat tracking error, testing of heliostat array configurations, understanding spatial variation of maximum loads across the solar field, and heliostat field layout and operating strategies. The recommended highest-priority pathway as first steps taken by HelioCon to address these gaps are to develop site characterization guidelines for heliostat design and field load measurements.

Emes, Matthew (ORCID:0000000341474387)↗
DOE OSTI2022

Slides for MS&T2022: Carbon-Carbon Composites from 3D Printed Preforms Graphitized by Electric Field Assisted Sintering

Formation of Carbon-Carbon (C-C) composite materials requires putting a carbon fiber preform through a cycle of densification, carbonization, and graphitization. This cycle is then repeated until the desired properties are achieved. The graphitization phase is particularly energy demanding as it involves heating the components to in excess of 2000°C. In traditional furnaces this is also quite time consuming as typical heating rates are low, and the actual heat treatment may be limited by component dimensions, as heat is conducted throughout the part primarily by thermal diffusion. Joule heating by means of Electric Field Assisted Sintering (EFAS) systems enables high heating rates and causes the part to be heated from the inside, which eliminates the thermal diffusion time constraint. In our process, C-C composites are formed by pyrolyzation of 3D-printed carbon fiber preforms, vacuum infiltration of naphthalene based mesophase pitch, and carbonization. The part is then graphitized in the EFAS system. By this method graphitization of the composite matrix was achieved in less than an hour. A part density of greater than 1.6g/cc (up to 1.9g/cc demonstrated) was achieved after only two rounds of infiltration and a single graphitization. Overall the process took less than a week which is significantly faster than traditional methods. Optical microscopy and XRD results are discussed to elucidate the internal porosity and characterize the degree of graphitization. Testing was performed showing the impact on fiber direction to electrical and thermal properties of the final part.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Updates and Validation for the n+ 63,65 Cu Cross Sections [Abstract]

The neutron induced total, elastic, and capture cross sections of 63,65 Cu isotopes were selected for evaluation in the resolved and unresolved resonance energy ranges by the National Criticality Safety Program to resolve discrepancies related to benchmark performance. This is especially evident for the series of ZEUS benchmarks in which copper is used as a reflector. Because copper is also used as structural material in both fission and fusion reactors, the need to address benchmark discrepancies linked to nuclear data deficiencies is a task of primary importance. The aim of this work is to describe the steps of evaluation work towards a consistent improvement of the benchmark performance. The R-matrix analysis with the SAMMY code focused on the 63 Cu(n,γ) reaction channel between 100-300 keV coupled to unresolved resonance region parameters up to 650 keV to fit average cross section data from a recent experiment. Due to the high sensitivity of many benchmarks to elastic scattering angular distribution data, especially for the 65 Cu isotope, the impact of these data was tested by generating Legendre coefficients from both resonance parameters and the Hauser-Feshbach model. Guided by the findings of Shaw et al., the performance of the current evaluation for 65 Cu was compared to that of ENDF/B-VII.1 and ENDF/B-VIII.0 by testing the reactivity coefficients corresponding to the validation suite of experimental criticality benchmarks for thermal, intermediate, and fast systems taken from the International Criticality Safety Benchmark Experiments Project Handbook. The benchmark performance is especially sensitive to 63 Cu(n,γ) and 65 Cu elastic scattering for neutron energies in the 100–500 keV region, whereas 100 keV is the upper limit of the resolved resonance region in the ENDF/B-VIII.0 evaluations for 63,65 Cu. The results highlight the need to handle the transition from the resolved resonance region to the high energy region carefully.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Exponential Backoff and Its Security Implications for Safety-Critical OT Protocols over TCP/IP Networks

The convergence of Operational Technology (OT) and Information Technology (IT) networks has become increasingly prevalent with the growth of Industrial Internet of Things (IIoT) applications. This shift, while enabling enhanced automation, remote monitoring, and data sharing, also introduces new challenges related to communication latency and cybersecurity. Oftentimes, legacy OT protocols were adapted to the TCP/IP stack without an extensive review of the ramifications to their robustness, performance, or safety objectives. To further accommodate the IT/OT convergence, protocol gateways were introduced to facilitate the migration from serial protocols to TCP/IP protocol stacks within modern IT/OT infrastructure. However, they often introduce additional vulnerabilities by exposing traditionally isolated protocols to external threats. This study investigates the security and reliability implications of migrating serial protocols to TCP/IP stacks and the impact of protocol gateways, utilizing two widely used OT protocols: Modbus TCP and DNP3. Our protocol analysis finds a significant safety-critical vulnerability resulting from this migration, and our subsequent tests clearly demonstrate its presence and impact. A multi-tiered testbed, consisting of both physical and emulated components, is used to evaluate protocol performance and the effects of device-specific implementation flaws. Through this analysis of specifications and behaviors during communication interruptions, we identify critical differences in fault handling and the impact on time-sensitive data delivery. The findings highlight how reliance on lower-level IT protocols can undermine OT system resilience, and they inform the development of mitigation strategies to enhance the robustness of industrial communication networks.

DNP3↗
DOE OSTI2023

Local Power Impact Experiment Design for a New Fuel Type for use in the Advanced Test Reactor

The Advanced Test Reactor (ATR), and complimentary zero-power ATR Critical (ATRC) reactor, located at Idaho National Labs (INL), are undergoing conversion from Highly Enriched Uranium (HEU) to Low Enriched Uranium (LEU). Both have a variety of testing locations that can receive large variations in flux due to its unique serpentine design, consisting of five lobes (see Figure 1). Initial criticality and power distribution throughout the core are controlled by core-external outer shim control cylinders (OSCCs). Distinct test loops allow for testing at specific temperatures, pressures, and irradiation conditions. The ATR is one of the key nuclear engineering research and testing facilities within the DOE National Laboratory Complex, and the ATRC supports its operation [1]. Currently, the Office of Material Management and Minimization (M3) within the National Nuclear Security Administration of the DOE is working to convert the remaining research reactors, including the ATR, from 93% HEU fuel to 19.75% LEU fuel (LEU) to support non-proliferation [2]. Extensive materials testing at INL and internationally has demonstrated that a high-density uranium molybdenum (U 10Mo) alloy can meet the performance requirements of the remaining high powered research reactors. However, there are many technical challenges to address before the conversion to LEU can be successful, including the accurate characterization of the reactor core physics with LEU fuel. Reactor physics safety evaluations currently use Monte Carlo for the 21st Century (MC21), a continuous-energy Monte Carlo radiation transport code [3]. Existing MC21 models of the ATR and ATRC cores have a validation basis for use in neutronics analyses with HEU fuel. The models are used to support safety analyses that include comparisons to the safety requirements for the reactors. However, the use of the LOWE element in the ATR and ATRC is not currently covered by the current model validation basis. To deploy the new fuel type, extensive computational reactor physics support is necessary to support the use of LOWE in the ATR and ATRC. Therefore, LOWE requires a rigorous validation basis, aligned with that of HEU fuel, that takes advantage of the existing software tools and processes currently used for the ATR and ATRC. The experiment to validate of the MC21 models for determining power, the Power Impact Validation Experiment, will consist of two flux runs in the ATRC, one with fully HEU loading and one with a single LOWE element. Both flux runs will be instrumented with 20 sets of azimuthal fission wires and 3 sets of axial fission wires, as shown in Figure 4. Standard flux run methodology will be used [4]. Power Impact Validation Experiment data will be compared against MC21 calculated data, both for absolute fission rate accuracy and to determine the relative change in fission rates between the two runs. The results of the Power Impact Validation Experiment and subsequent evaluations will provide the validation basis for MC21 for use with LOWE elements. Key features of the Power Impact Validation Experiment include: (1) Two flux runs to allow for LOWE perturbed measurements to be compared to already validated measurements taken from a full core of HEU fuel, (2) Optimization of instrumentation to balance analytical needs with practical considerations (e.g., limited time window to count beta particles from fission products), and (3) Standard ATRC core loading, including both driver positions and flux traps, to minimize cost while remaining representative of typical ATR core loading.

42 ENGINEERING↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗