Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUIDS”

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 OSTI2020

Alternative, Less Water-Intensive, Fracturing Fluids for Enhanced Geothermal Systems

Four fracturing fluids; water, CO2, CO2 with water, and CO2 with an aqueous solution of a CO2-reactive polymer, poly(allylamine) (PAA, 1wt%) were evaluated using a high-temperature true-triaxial fracturing apparatus and ½ foot side granite cubic samples. All three CO2-based fracturing fluids, CO2, CO2 with water, and CO2 with aqueous PAA fractured granite at higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When partially saturating the rock sample with PAA solution followed by fracturing with CO2, the volume expansion caused by CO2-induced cross-linking of PAA leads to a faster pressure increase due to the associated stress generated and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. It is concluded that CO2 (when injected in HDR) and CO2/PAA are the fluids of choice for stimulation of granitic rock samples under the studied geothermal P/T conditions. CO2/PAA, however, offer the following three additional advantages; 1) it requires significantly lower volumes of CO2, 2) fracture permeability is independent of injection flow rate, and 3) the reversible (previously reported) viscosity increase is beneficial to transport proppants if they are ever developed enhanced geothermal systems.

Jian, Guoqing↗
DOE OSTI2024

Open aperture flow cells for on-line optical analysis of process fluids

A flow cell system includes a vessel and a fluid located in the vessel. A fluid surface of the fluid can be vented to a first gas pressure. The fluid surface can have a first cross-sectional area. The flow cell system includes a conduit in fluid communication with the vessel and positioned downstream of the vessel. The conduit can have a region that includes one or more orifices and has a second cross-sectional area. The second cross-sectional area can be less than the first cross-sectional area. The one or more orifices can be vented to a second gas pressure. The second gas pressure can be equal to or greater than the first gas pressure. Methods for analyzing a process fluid can include characterizing the fluid in the conduit.

Launiere, Cari A.↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗
DOE OSTI2023

Systematic study on the phase transition of confined fluid mixture up to the critical region

Systematic experimental study on the phase transition of confined fluids up to the critical region has been performed, in which many new insights have been obtained. For confined pure fluids, compared to the bulk phase, the critical point shifts to lower pressure and temperature, but the relative shift of critical pressure is much larger than that of critical temperature. The behavior of the shift, i.e., the shift as a function of pore size, is dependent on the pore type. For confined fluid mixtures, compared to the bulk phase, the critical point also shifts to lower pressure and temperature, and the relative shift of critical pressure is also much larger than that of critical temperature. Confined fluid mixture behaves similarly to confined pure fluid, i.e., it exhibits hysteresis at lower temperature range and does not exhibit phase coexistence or phase envelope. An equation of state for confined pure fluids, based on the Generalized van der Waals partition function, has also been developed and demonstrated to accurately predict the capillary condensation of simple molecules in MCM-41. The parameters of the equation of state are derived from experimental critical point.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Tabulated Fluid Properties Research Report

The Multiphysics Object-Oriented Simulation Environment (MOOSE) enables a wide range of advanced nuclear reactor simulations.[6] Under the guidance of MOOSE’s Thermal Hydraulics Team,I worked to expand the capabilities of Tabulated Fluid Properties (TFP) in the fluid properties module. The fluid properties module allows the user to determine a variety of fluid properties by interpolating points between tabulated data. I implemented the ability to use bilinear interpolation instead of bicubic interpolation for interpolating tabulated data. I also changed the method of variable set inversions to use a 2-dimensional Newton’s Method utility that I created. Variable set inversions are often done from (v,e) to (p,T), where v is specific volume, e is specific internal energy, p is pressure and T is temperature. New routines have also been added into TFP such that it can be used with more applications, such as the Navier Stokes and Thermal Hydraulics modules in MOOSE for Pronghorn[5] and RELAP-7[1] respectively. This work was spurred by interest from NASA in testing a Nuclear Thermal Propulsion (NTP) engine system. NTP engines have drastically different fluid properties throughout the engine and Tabulated Fluid Properties provides the flexibility needed to properly simulate and test these engines.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Assessment of the Utility of the Oral Fluid and Plasma Proteome for Hydrocodone Usage

Non-medical use and abuse of prescription opioids is a growing problem in both the civilian and military communities, with minimal technologies for detecting hydrocodone use. This study explored the proteomic changes that occur in the oral fluid and blood plasma following controlled hydrocodone administration in twenty subjects. The global proteomic profile was determined for samples taken at four time points per subject of pre-exposure and 4, 6, or 168 hours post-exposure. The oral fluid samples analyzed herein provided greater differentiation between “baseline” and “response” time points than was observed with blood plasma, at least partially due to significant person-to-person variability in the plasma proteome. A total of 399 proteins were identified from oral fluid samples, and the abundance of 102 of those proteins was determined to be significantly different upon metabolism of hydrocodone (4 and 6 hour time points) as compared to baseline levels in the oral fluid (pre-dose and 168 hour). This is the first, to our knowledge, assessment of the oral fluid and blood plasma proteome following hydrocodone administration, which demonstrates the potential of oral fluid as a non-invasive sample that may reveal features of hydrocodone or other opioid use and/or misuse.

Kaiser, Brooke LD↗
DOE OSTI2024

Mineral and fluid transformation of hydraulically fractured shale: case study of Caney Shale in Southern Oklahoma

This study explores the geochemical reactions that can cause permeability loss in hydraulically fractured reservoirs. The experiments involved the reaction of powdered-rock samples with produced brines in batch reactor system at temperature of 95 °C and atmospheric pressure for 7-days and 30-days respectively. Results show changes in mineralogy and chemistry of rock and fluid samples respectively, therefore confirming chemical reactions between the two during the experiments. The mineralogical changes of the rock included decreases of pyrite and feldspar content, whilst carbonate and illite content showed an initial stability and increase respectively before decreasing. Results from analyses of post-reaction fluids generally corroborate the results obtained from mineralogical analyses. Integrating the results obtained from both rocks and fluids reveal a complex trend of reactions between rock and fluid samples which is summarized as follows. Dissolution of pyrite by oxygenated fluid causes transient and localized acidity which triggers the dissolution of feldspar, carbonates, and other minerals susceptible to dissolution under acidic conditions. The dissolution of minerals releases high concentrations of ions, some of which subsequently precipitate secondary minerals. On the field scale, the formation of secondary minerals in the pores and flow paths of hydrocarbons can cause significant reduction in the permeability of the reservoir, which will culminate in rapid productivity decline. This study provides an understanding of the geochemical rock–fluid reactions that impact long term permeability of shale reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Phase diagram to demarcate supercritical, transcritical, and continuous phase regimes for binary fluid equilibrium mixing relevant to combustion applications

Here, a robust methodology to develop phase diagrams of binary fluid mixtures at fixed thermo dynamic conditions (pressure, temperature, and mole fraction) as well as of two initially separated fluids undergoing mixing near critical conditions are presented for fluids and con ditions relevant to rockets, gas turbines, and diesel engine applications. Phase equilibria of mixtures is first examined to provide insight into the continuous-phase mixing behavior (including but not limited to supercritical behavior), and to develop a broadly applicable phase-diagram for binary fluid mixtures at fixed conditions. Next, adiabatic mixing theory and reduced Helmholtz equations of state are used to predict the thermodynamic conditions required to attain continuous-phase binary fluid mixing near critical conditions. Then, a 3D surface diagram (P,T fuel ,T amb ) separating single and two-phase regions is constructed by varying the ambient pressure and the initial temperatures of the two fluids. Polynomial fits of the 3D surfaces for 10 different binary mixtures are tabulated for nitrogen-alkane and methane-oxygen blends relevant to air-breathing and propellant based engines, respectively.

42 ENGINEERING↗
DOE OSTI2021

Anomalous normal fluid response in a chiral superconductor UTe2

Abstract Chiral superconductors have been proposed as one pathway to realize Majorana normal fluid at its boundary. However, the long-sought 2D and 3D chiral superconductors with edge and surface Majorana normal fluid are yet to be conclusively found. Here, we report evidence for a chiral spin-triplet pairing state of UTe 2 with surface normal fluid response. The microwave surface impedance of the UTe 2 crystal was measured and converted to complex conductivity, which is sensitive to both normal and superfluid responses. The anomalous residual normal fluid conductivity supports the presence of a significant normal fluid response. The superfluid conductivity follows the temperature behavior predicted for an axial spin-triplet state, which is further narrowed down to a chiral spin-triplet state with evidence of broken time-reversal symmetry. Further analysis excludes trivial origins for the observed normal fluid response. Our findings suggest that UTe 2 can be a new platform to study exotic topological excitations in higher dimension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Development and validation of a mathematical model of heart rate response to fluid perturbation

Physiological closed-loop controlled (PCLC) medical devices monitor and automatically adjust the patient’s condition by using physiological variables as feedback, ideally with minimal human intervention to achieve the target levels set by a clinician. PCLC devices present a challenge when it comes to evaluating their performance, where conducting large clinical trials can be expensive. Virtual physiological patients simulated by validated mathematical models can be utilized to obtain preclinical evidence of safety and assess the performance of the PCLC medical device during normal and worst-case conditions that are unlikely to happen in a limited clinical trial. A physiological variable that plays a major role during fluid resuscitation is heart rate (HR). For in silico assessment of PCLC medical devices regarding fluid perturbation, there is currently no mathematical model of HR validated in terms of its predictive capability performance. This paper develops and validates a mathematical model of HR response using data collected from sheep subjects undergoing hemorrhage and fluid infusion. The model proved to be accurate in estimating the HR response to fluid perturbation, where averaged between 21 calibration datasets, the fitting performance showed a normalized root mean square error (NRMSE) of 7.41 ± 2.8%. The model was also evaluated in terms of model predictive capability performance via a leave-one-out procedure (21 subjects) and an independent validation dataset (6 subjects). Two different virtual cohort generation tools were used in each validation analysis. The generated envelope of virtual subjects robustly met the defined acceptance criteria, in which 95% of the testing datasets presented simulated HR patterns that were within a deviation of 50% from the observed data. In addition, out of 16000 and 18522 simulated subjects for the leave one-out and independent datasets, the model was able to generate at least one virtual subject that was close to the real subject within an error margin of 9.56 ± 3.15% and 11.1 ± 1.22% NRMSE, respectively. In conclusion, the model can generate valid virtual HR physiological responses to fluid perturbation and be incorporated into future non-clinical simulated testing setups for assessing PCLC devices intended for fluid resuscitation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

American eel resilience to simulated fluid shear associated with passage through hydroelectric turbines

American eel (Anguilla rostrata) populations have declined within their native range along the eastern coast of North America due to factors such as commercial fishing, habitat alteration, and dams. American eel are catadromous fish species, and high mortality rates (>40%) have been observed for freshwater life-stage adult eel passing downstream through hydropower turbines. Lacerations and sectioning of fish have been observed downstream of turbines and these injuries are commonly associated with direct contact with the turbine runner, whether through blade strike or pinching and grinding. Exposure to fluid shear may also be a source of injury, however, little is known about American eel susceptibility to this physical stressor. Eels are considerably flexible when compared to other fish species and lack other morphological characteristics that would make them susceptible to fluid shear, such as protruding eyes, large scales, and large operculum. European eel, which have previously been tested for susceptibility to fluid shear, were found to be resilient. To determine if American eel are also resilient to fluid shear, forty American eel were exposed to a water jet, simulating severe fluid shear (strain rate > 800 s –1 ) that fish may experience when passing downstream through turbines. No immediate or delayed (48 h) signs of injury were observed after exposure to severe fluid shear. Based on this study, and a previous study conducted on American eel susceptibility to barotrauma, the source of injury and mortality of American eel passing through turbines is likely attributed to blade strike or pinching and grinding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An incompressible Eulerian method for fluid–structure interaction with mixed soft and rigid solids

Here, we present a general simulation approach for incompressible fluid–structure interactions in a fully Eulerian framework using the reference map technique. The approach is suitable for modeling one or more rigid or finitely deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid–soft solver [Rycroft et al., “Reference map technique for incompressible fluid–structure interaction,” J. Fluid Mech. 898, A9 (2020)] and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.

42 ENGINEERING↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Oxygen isotope evidence for input of magmatic fluids and precipitation of Au-Ag-tellurides in an otherwise ordinary adularia-sericite epithermal system in NE China

Abstract Tellurium-rich (Te) adularia-sericite epithermal Au-Ag deposits are an important current and future source of precious and critical metals. However, the source and evolution of ore-forming fluids in these deposits are masked by traditional bulk analysis of quartz oxygen isotope ratios that homogenize fine-scale textures and growth zones. To advance understanding of the source of Te and precious metals, herein, we use petrographic and cathodoluminescence (CL) images of such textures and growth zones to guide high spatial resolution secondary ion mass spectroscopy (SIMS) oxygen isotope analyses (10 µm spot) and spatially correlated fluid inclusion microthermometric measurements on successive quartz bands in contemporary Te-rich and Te-poor adularia-sericite (-quartz) epithermal Au-Ag vein deposits in northeastern China. The results show that large positive oxygen isotope shifts from –7.1 to +7.7‰ in quartz rims are followed by precipitation of Au-Ag telluride minerals in the Te-rich deposit, whereas small oxygen isotope shifts of only 4‰ (–2.2 to +1.6‰) were detected in quartz associated with Au-Ag minerals in the Te-poor deposits. Moreover, fluid-inclusion homogenization temperatures are higher in comb quartz rims (avg. 266.4 to 277.5 °C) followed by Au-Ag telluride minerals than in previous stages (~250 °C) in the Te-rich deposit. The Te-poor deposit has a consistent temperature (~245 °C) in quartz that pre- and postdates Au-Ag minerals. Together, the coupled increase in oxygen isotope ratios and homogenization temperatures followed by precipitation of Au-Ag tellurides strongly supports that inputs of magmatic fluid containing Au, Ag, and Te into barren meteoric water-dominated flow systems are critical to the formation of Te-rich adularia-sericite epithermal Au-Ag deposits. In contrast, Te-poor adularia-sericite epithermal Au-Ag deposits show little or no oxygen isotope or fluid-inclusion evidence for inputs of magmatic fluid.

Geochemistry & Geophysics↗
DOE OSTI2020

High operating temperature heat transfer fluids for solar thermal power generation. Final report

This is the final technical report for the DOE MURI project on developing high-operating temperature heat transfer fluids for solar thermal power generation. This project aimed to evaluate liquid metals as a heat-transport fluid and potential energy storage medium for concentrated solar powered (CSP) systems. The project consisted of five phases, starting with fluid selection and screen and ending with potential CSP integration. The initial fluid selection arrived at heavy liquid metals as candidate fluids such as lead bismuth eutectic (LBE) as the prime liquid-alloy candidate to be studied, given the conditions above. Fifty-seven liquid compositions were initially identified with liquidus temperatures lower than 200°C, composed of elements that provide no driving force to the reduction of Cr 2 O 3 or Al 2 O 3 passivation layers in pipe material which is required for corrosion prevention. LBE was identified as the primary candidate fluid after these compositions were compared against computational thermodynamic modeling to screen all of the selected systems that would be compatible with our design objectives. Twenty-eight of these compositions were also investigated with the rapid characterization methods; the accuracies of the liquidus temperature and the specific heats predicted by the computations were better than 3 percent and 9 percent, respectively. We conducted numerus corrosion tests of various Fe-Cr-Al based steels in LBE and its variants at temperate as high as 800°C. It was found that minor additions to the liquid did not provide sufficient corrosion protection to warrant their use in real applications. Furthermore, we investigated the detailed mechanism of corrosion between heavy liquid metal and the pipe material at these very high temperatures. A liquid metal loop was constructed to operate at these high temperatures. The loop construction and operation lead to valuable lessons learned that will be needed for industrial scale loops. Furthermore, we learned that electromagnetic pumps at these extreme temperatures will have to be reengineered and redesigned to be used for CSP applications.

14 SOLAR ENERGY↗
DOE OSTI2021

Development and Field Testing Novel Natural Gas Surface Process Equipment for Replacement of Water as Primary Hydraulic Fracturing Fluid (Final Report)

Southwest Research Institute® (SwRI ® ), Schlumberger Technology Corporation (SLB), and Chevron Corporation ® (Chevron) developed a novel, optimized, and lightweight process for natural gas (NG) to replace water as a low-cost fracturing medium with a low environmental impact. Hydraulic fracturing is used to increase oil and NG production by injecting high-pressure fluid, primarily water, into a rock formation, which fractures the rock and releases trapped oil and NG. This method was developed to increase yield and make feasible production areas that would not otherwise be viable for large-scale oil and NG extraction using traditional drilling technologies. Since the fracturing fluid is composed of approximately 90% water, one of the principal drawbacks to hydraulic fracturing is its excessive water use and associated large environmental footprint. According to recent data collected at fracture sites within the United States, fracturing applications in North America can consume as much as 9 million gallons of water per well. During the fracturing process, some of the fracturing fluid is permanently lost and the portion that is recovered is contaminated by both fracturing chemicals and dissolved solids from the formation. The recovered water or flow-back represents a significant environmental challenge, as it must be treated before it can be reintroduced into the natural water system. Although there is some recycling for future fracturing, the majority of the flow-back water is hauled from the well site to a treatment facility or to an injection well for permanent underground disposal. To mitigate these issues, an optimized, lightweight and modular surface process using NG to replace a majority of the water was developed as a cost-effective and environmentally clean fracturing fluid. Using NG will result in significantly less consumption since the gas that is injected as a fracturing fluid will be mixed with the formation gas and extracted as if it were from the formation itself. This process will minimize the collection, waste, and treatment of large amounts of water and will reduce the environmental impact of transporting and storing the fracturing fluid. There are two major steps involved in utilizing NG as the primary fracturing medium: (1) increasing the supply pressure of NG to wellhead pressures suitable for fracturing and (2) mixing the required chemicals and proppant needed for the fracturing process at these elevated pressures. The second step (NG-proppant mixing at elevated pressures) still requires demonstration; but very similar processes have been demonstrated in the field with other gases such as nitrogen (N 2 ) and carbon dioxide (CO 2 ). However, the first step (a compact, on-site unit for generating high-pressure NG at costs feasible for fracturing) has not been developed and is currently not commercially available. Due to the inherent compressibility of NG, more energy is required to compress the gas than what is required to pump water (or other incompressible liquids) to the very high pressure required for downhole injection. This project aimed to develop a novel, hybrid method to overcome this challenge. The project accomplishments and findings are discussed in this report. Ultimately, the research and development efforts described herein demonstrate that fracturing with NG foam is a feasible alternative to using water.

02 PETROLEUM↗