Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID INJECTION”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

SPI Conformance Gel Applications in Geothermal Zonal Isolation

Zonal isolation in geothermal injection and producing wells is important while drilling the wells when highly fractured geothermal zones are encountered and there is a need to keep the fluids from interfering with the drilling operation. Department of Energy’s (DOE) Energy Efficiency and Renewable Energy (EERE) objectives are to advance technologies to make it more cost effective to develop, produce, and monitor geothermal reservoirs and produce geothermal energy. Thus, zonal isolation is critical to well cost, reservoir evaluation and operations. Traditional cementing off of the lost circulation or thief zones during drilling is often done to stem the drilling mud losses. This is an expensive and generally unsuccessful technique losing the potential of the remaining fracture system. Selective placement of strong SPI gels into only the offending fractures can maintain and even improve operational efficiency and resource life. The SPI gel system is a unique silicate based gel system that offers a promising solution to thief zones and conformance problems with water and CO2 floods and potentially geothermal operations. This gel system remains a low viscosity fluid until an initiator (either internal such as an additive or external such as CO2) triggers gelation. This is a clear improvement over current mechanical methods of using packers, plugs, liners and cementing technologies that often severely damage the highly fractured area that is isolated. In the SPI gels, the initiator sets up the fluid into a water-like (not a precipitate) gel and when the isolated zone needs to be reopened, the SPI gel may be removed with an alkaline solution without formation damage occurring. In addition, the SPI gel in commercial quantities is expected to be less expensive than competing mechanical systems and has unique deep placement possibilities. This project seeks to improve upon the SPI gel integrity by modifying the various components to impart temperature stability for use in geothermal.

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

The Big Group of People Looking at How to Control Putting the Parts of the Air That Are the Same as What You Breathe Out Into Small Spaces in Rocks

Representing the Nanoscale Control of Geologic CO2 (NCGC), this document is one of the entries in the Ten Hundred and One Word Challenge. As part of the challenge, the 46 Energy Frontier Research Centers were invited to represent their science in images, cartoons, photos, words and original paintings, but any descriptions or words could only use the 1000 most commonly used words in the English language, with the addition of one word important to each of the EFRCs and the mission of DOE energy. The mission of NCGC is to build a fundamental understanding of molecular-to-pore-scale processes in fluid-rock systems, and to demonstrate the ability to control critical aspects of flow, transport, and mineralization in porous rock media as applied to the injection and storage of carbon dioxide (CO2) in subsurface reservoirs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Effects of a CFD-improved dimple stepped-lip piston on thermal efficiency and emissions in a medium-duty diesel engine

Diesel piston-bowl shape is a key design parameter that affects spray-wall interactions and turbulent flow development, and in turn affects the engine’s thermal efficiency and emissions. It is hypothesized that thermal efficiency can be improved by enhancing squish-region vortices as they are hypothesized to promote fuel-air mixing, leading to faster heat-release rates. However, the strength and longevity of these vortices decrease with advanced injection timings for typical stepped-lip (SL) piston geometries. Dimple stepped-lip (DSL) pistons enhance vortex formation at early injection timings. Previous engine experiments with such a bowl show 1.4% thermal efficiency gains over an SL piston. However, soot was increased dramatically [SAE 2022-01-0400]. In a previous study, a new DSL bowl was designed using non-combusting computational fluid dynamic simulations. This improved DSL bowl is predicted to promote stronger, more rotationally energetic vortices than the baseline DSL piston: it employs shallower, narrower, and steeper-curved dimples that are placed further out into the squish region. In the current experimental study, this improved bowl is tested in a medium-duty diesel engine and compared against the SL piston over an injection timing sweep at low-load and part-load operating conditions. No substantial thermal efficiency gains are achieved at the early injection timing with the improved DSL design, but soot emissions are lowered by 45% relative to the production SL piston, likely due to improved air utilization and soot oxidation. However, these benefits are lost at late injection timings, where the DSL piston renders a lower thermal efficiency than that of the SL piston. Energy balance analyses show higher wall heat transfer with the DSL piston than with the SL piston despite a 1.3% reduction in the piston surface area. Vortex enhancement may not necessarily lead to improved efficiency as more energetic squish-region vortices can lead to higher convective heat transfer losses.

42 ENGINEERING↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

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

Investigating Liquid Water Transport in Different Pore Structure of Gas Diffusion Layers for PEMFC Using Lattice Boltzmann Method

Proton exchange membrane fuel cells (PEMFC) require a gas diffusion layer (GDL) to aid in the transport of liquid fuel to the catalyst layer. In this work, direct modeling using the Lattice Boltzmann Method (LBM) was applied to X-ray CT scans of four different carbon gas diffusion layers to understand the mass transport properties through the samples. Three injection orientations were used to study local saturation levels, water evolution through the sample, and mass transport behavior at breakthrough conditions. The LBM, combined with computational fluid dynamic modeling techniques, can accurately predict liquid saturation at the macro and micro scale, which provides more insight into the mass transport phenomena through the GDL. The change of pore structure and orientation in both the in-plane and through-plane determines the path that liquid water must take, which could aid or impact PEMFC performance. The outcomes from this work will also benefit any research that needs knowledge of internal mass transport qualities of gas diffusion media.

Sepe, M.↗
DOE OSTI2023

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Utah FORGE: Well 16B(78)-32 Logs from Schlumberger Technologies

This dataset is a collection of well logs provided by Schlumberger Technologies from the Utah FORGE well 16B(78)-32 drilling project. Information here includes critical borehole information collected by an ultrasonic borehole imager (UBI) and a fullbore formation microimager (FMI). Well 16B(78)-32 serves as the production well for reservoir creation, fluid circulation, and demonstration of heat extraction for the FORGE project. It has been drilled as a doublet approximately 300 feet parallel to and above the injection well 16A(78)-32. The total depth measured 10,947 feet and the vertical depth measured 8,357 feet.

11.5 inch intermediate casing logs↗
DOE OSTI2014

Dual-Mode Propulsion System Enabling CubeSat Exploration of the Solar System: NASA Innovative Advanced Concepts (Phase I Final Report)

It is apparent the cost of planetary exploration is rising as mission budgets declining. Currently small scientific beds geared to performing limited tasks are being developed and launched into low earth orbit (LEO) in the form of small-scale satellite units, i.e., CubeSats. These micro- and nano-satellites are gaining popularity among the university and science communities due to their relatively low cost and design flexibility. To date these small units have been limited to performing tasks in LEO utilizing solar-based power. If a reasonable propulsion system could be developed, these CubeSat platforms could perform exploration of various extra-terrestrial bodies within the solar system engaging a broader range of researchers. Additionally, being mindful of mass, smaller cheaper launch vehicles (~1,000 kgs to LEO) can be targeted. This, in effect, allows for beneficial exploration to be conducted within limited budgets. Researchers at the Center for Space Nuclear Re-search (CSNR) are proposing a low mass, radioisotope-based, dual-mode propulsion system capable of extending the exploration realm of these CubeSats out of LEO. The proposed radioisotope-based system would leverage the high specific energies [J/kg] associated with radioisotope materials and enhance their inherent low specific powers [W/g]. This is accomplished by accumulating thermal energy from nuclear decay within a central core over time. This allows for significant amounts of power to be transferred to a flowing gas over short periods of time. In the proposed configuration the stored energy can be utilized in two ways: (1) with direct propellant injection to the core, the energy can be converted into thrust through the use of a converging-diverging nozzle and (2) by flowing a working fluid through the core and subsequent Brayton engine, energy within the core can be converted to electrical energy. The first scenario achieves moderate ranges of thrust, but at a higher Isp than traditional chemical-based systems. The second scenario allows for the production of electrical power, which is then available for electric-based propulsion. Additionally, once at location the production of electrical power can be dedicated to the payload’s communication system for data transfer. Ultimately, the proposed dual-mode propulsion platform capitalizes on the benefits of two types of propulsion methods – the thrust of thermal propulsion ideal for quick orbital maneuvers and the specific impulse of electric propulsion ideal for efficient inter-planetary travel. Previous versions of this RTR-based concept have been studied for various applications [NETS 1-3]. The current version of this concept is being matured through a NASA Innovative Advanced Concepts (NIAC) Phase I grant, awarded for FY 2014. In this study the RTR concept is being developed to deliver a 6U CubeSat payload to the orbit of the Saturnian moon - Enceladus. Additionally, this study will develop an entire mission architecture for Enceladus targeting a total allowable launch mass of 1,000 kg.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

Probabilistic Geomechanical Analysis of Compartmentalization at the Snøhvit CO 2 Sequestration Project

Pressure build-up caused by large-scale CO 2 injection is one of the key concerns during a carbon sequestration project. Large overpressures can potentially compromise seal integrity, reactivate faults, and induce concerning seismicity, among other hazards. Furthermore, pressure build-up is directly related with storage capacity. In this work we study the geomechanical response to CO 2 injection in the Tubåen Fm at Snøhvit, with the objective of understanding the potential for fault reactivation, leakage, and contamination of the producing interval in an adjacent block through the bounding faults. Furthermore, we evaluate the potential contribution of a structural component to reservoir compartmentalization. We combine simplified analytical models, based on critically stressed fracture theory and a Mohr-Coulomb failure criterion, with a quantitative risk assessment methodology to perform a rigorous uncertainty quantification. Large stress uncertainties are present and are reflected in the modeling results. It was found that under the “most likely” stress state (Strike-Slip and North-South maximum horizontal stress (S Hmax ) azimuth) the faults are fairly stable, and tensile hydrofracturing of the caprock would be expected before fault reactivation. In most of the analyzed cases, the critical pressure perturbation needed for reactivation is above 13 MPa, which was the limiting pressure increase before reaching the fracture pressure. Faults were found to be ~ 20% less stable when considering plausible variations in the orientation of S Hmax . In those cases, fault reactivation could be expected before caprock failure if injection continued, posing risks of induced seismicity, CO 2 leakage and/or contamination of the gas producing reservoir. Finally, a potential structural component in the compartmentalization and fluid migration is difficult to assess due to stress orientation uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Remedial Amendment Delivery near the Water Table Using Shear Thinning Fluids: Experiments and Numerical Simulations

The use of shear thinning fluids (STFs) containing xanthan is a potential enhancement for emplacing a solute amendment near the water table and within the capillary fringe. Most research to date related to STF behavior has involved saturated and confined conditions. A series of flow cell experiments were conducted to investigate STF emplacement in variable saturated homogeneous and layered heterogeneous systems. Besides flow visualization using dyes, amendment concentrations and pressure data were obtained at several locations. The experiments showed that injection of STFs considerably improved the subsurface distribution near the water table by mitigating preferential flow through higher permeability zones compared to no-polymer injections. The phosphate amendment migrated with the xanthan SFT without retardation. Despite the high viscosity of the STF, no excessive mounding or preferential flow were observed in the unsaturated zone. The STOMP simulator was able to predict the experimentally observed fluid displacement and amendment concentrations reasonably well. Cross flow between layers could be interpreted as the main mechanism to transport STFs into lower permeability layers based on the observed pressure gradient and concentration data in layers of differing hydraulic conductivity.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2024

Risk Considerations of Transitioning CO2-EOR Field to CO2 storage Field: Case Study

In the United States (U.S.), carbon dioxide (CO2) injection wells at EOR sites are currently regulated as Class II wells under the U.S. Environmental Protection Agency’s (EPA) Underground Injection Control (UIC) program while dedicated geological CO2 storage (GCS) wells are considered Class VI wells. A CO2-EOR facility considering a transitioning from tertiary oil recovery to injecting CO2 for the primary purpose of long-term storage is required to obtain a Class VI permit where this transition poses an increased risk to underground sources of drinking water. This study considers how transitioning operations from CO2-EOR to storage can impact reservoir plume and pressure transient in the storage envelope, and how these changes could impact area of review and potential unwanted fluid migration. We developed a case study to assess subsurface response and leakage risks associated with a representative, hypothetical operation in a carbonate reservoir. This reservoir transitions from tertiary hydrocarbon recovery to dedicated GCS. Reservoir simulations were run for a set of credible CO2-EOR scenarios to estimate distributions of fluids phases and pressures throughout the model domain after CO2 flooding as well as forecasting the behavior of the reservoir after the transition to a dedicated storage phase. The evolutions from all simulated scenarios were used as the basis for leakage risk quantification using the National Risk Assessment Partnership’s Open-Source Integrated Assessment Model (NRAP-Open-IAM) with a novel reduced-order model to estimate time-dependent leakage of CO2, brine, and hydrocarbon fluids through potentially leaky wells. Results include a description of reservoir response, an estimate of the areal extent that could potentially be impacted by leakage to underground sources of drinking water, and estimates of the magnitude of potential leakage. Considerations for dedicated storage injection well selection, injectivity, and injection scheme performance and potential leakage risk are presented, with implications for risk assessment of well transition discussed. This study presents a risk-based workflow for the Class II to Class VI well transition. Integrating credible numerical simulation of viable CO2-EOR to dedicated CO2 storage with quantitative risk assessment tools, such as the NRAP-Open-IAM, will provide a valuable means to devise operational scenarios and inform decision-making related to storage benefit, leakage risk, and liability. Presented at the SPE/AAPG/SEG Carbon Capture Utilization and Storage Conference in Houston, TX, March 11-13, 2024.

Liu, Guoxiang↗
DOE OSTI2023

Numerical Optimization of Spray-Guided Spark Assistance for Cold Idle Operation in a Heavy-Duty Gasoline Compression Ignition Engine

This article describes the results of a response surface model (RSM)-based numerical optimization campaign for spray-guided spark assistance at cold operations in a heavy-duty gasoline compression ignition (GCI) engine. On the basis of an earlier work on spark-assisted GCI cold combustion, a space-filling design of experiments (DoE) method was first undertaken to investigate a multitude of hardware design variables and engine operating parameters. The main design variables included the number of injector nozzles, fuel split quantities and injection timings, and spark timing. The objective variables were engine combustion efficiency (ŋ c ), maximum pressure rise rate (MPRR), and engine-out nitrogen oxide (NO x ) emissions. A total of 150 design candidates were automatically generated using the Sobol sequence method provided by the commercial software package, CAESES. Then, closed-cycle computational fluid dynamic (CFD) spark-assisted GCI simulations under cold idling operations were performed. The outcomes from the CFD-DoE design campaign were utilized to construct high-fidelity RSMs that allowed for further design optimization of the spark plug- and fuel injector-related design variables, along with fuel injection strategy parameters. A merit function with respect to objective variables was formulated with an appropriate weight assignment on each objective variable. Finally, the best design candidate was identified from the RSM-based optimization process and further validated in the CFD analysis. The best design candidate showed the potential to significantly improve combustion efficiency (ŋ c > 90%) over the baseline at cold idle while satisfying MPRR and NOx emissions constraints (MPRR < 5 bar/CAD and NO x < 4.5 g/kWh).

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Finite Element Analysis (FEA) for Water-Foam Fracturing of Granite Rock

In addition to the foam data that were obtained from literature and that were collected from the current study, simulation data was also generated from finite element analysis (FEA) conducted in this study using COMSOL Multiphysics software. The FEA models were built to simulate the experiments conducted at Oak Ridge National Laboratory (ORNL) on cement and granite samples. In these FEA models, temperature was kept at ambient while the pressure profile resembled the loading conditions during the ORNL experiments, where pressure was either monotonically increased or applied cyclically. The cement material was used as a model material and was used to study Von Mises stress and tensile stress distribution for different bore hole length geometry using a parametric sweep with water as fracturing fluid using solid-fluid interaction module. For the granite material, FEA models were developed for stress analysis of cylindrical samples with water or foam fluids. The solid mechanics module in COMSOL was implemented to solve for Von Mises stress and tensile stress. The fluid-structure interaction module was implemented to solve for water-foam interaction on granite cylinder with addition of fluid-loading on structure, i.e., large deformation in solid mechanics with no impact on fluid deformation. Foam was considered as a pseudo single-phase compressible fluid for which material properties were calculated from water and gas (nitrogen) phases. The density of foam is calculated as a function of the densities of water and nitrogen, while viscosity is a function of temperature. Four types of FEA analyses were modelled: 1. Monotonic injection with water 2. Monotonic injection with foam 3. Cyclic injection with water 4. Cyclic injection with foam All the COMSOL files are converted to a zip file which is save in .mph.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Modification of the Alfvén wave spectrum by pellet injection

Alfvén eigenmodes driven by energetic particles are routinely observed in tokamak plasmas. These modes consist of poloidal harmonics of shear Alfvén waves coupled by inhomogeneity in the magnetic field. In this work, further coupling is introduced by 3D inhomogeneities in the ion density during the assimilation of injected pellets. This additional coupling modifies the Alfvén continuum and discrete eigenmode spectrum. We find that the frequencies of Alfvén eigenmodes drop dramatically when a pellet is injected in JET. From these observations, information about the changes in the ion density caused by a pellet can be inferred. To use Alfvén eigenmodes for MHD spectroscopy of pellet injected plasmas, the 3D MHD codes Stellgap and AE3D were generalised to incorporate 3D density profiles. Additionally, a model for the expansion of the ionised pellet plasmoid along a magnetic field line was derived from the fluid equations. Thereby, the time evolution of the Alfvén eigenfrequency is reproduced. By comparing the numerical frequency drop of a toroidal Alfvén eigenmode (TAE) to experimental observations, the initial ion density of a cigar-shaped ablation region of length 4 cm is estimated to be n * =6.8×10 22 m -3 at the TAE location (r/a≈0.75). The frequency sweeping of an Alfvén eigenmode ends when the ion density homogenises poloidally. Modelling suggests that the time for poloidal homogenisation of the ion density at the TAE position is τ h =18±4 ms for inboard pellet injection, and τ h =26±2 ms for outboard pellet injection. By reproducing the frequency evolution of the elliptical Alfvén eigenmode (EAE), the initial ion density at the EAE location (r/a≈0.9) can be estimated to be n * =4.8×10 22 m -3 . Poloidal homogenisation of the ion density takes 2.7 times longer at the EAE location than at the TAE location for both inboard and outboard pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗