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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fuel injection assembly for a turbomachine

Fuel injection assemblies and combustors are provided. A fuel injection assembly includes a fuel injector having a first annular wall that about an axial centerline and extends from a primary inlet disposed at a first end to a primary outlet disposed at a second end. A second annular wall surrounds the first annular wall. A fuel plenum is defined between the first annular wall and the second annular wall. A fuel duct extends from a fuel outlet defined in the second annular wall to a fuel inlet. wherein the fuel duct is in fluid communication with the fuel plenum. The fuel duct includes a polygonal segment and a cylindrical inlet segment. The polygonal segment extends from the fuel outlet to the cylindrical inlet segment.

Griffis, Clay Thomas↗
DOE OSTI2023

Impact of organic solutes on capillary phenomena in water-CO2-quartz systems

Hypothesis: The migration of supercritical CO2 (scCO2) injected into underground reservoirs as part of carbon capture and storage is influenced by organic contamination affecting mineral wettability. Molecular dynamics (MD) simulations of relevant systems that incorporate representative organic solutes allow detailed investigation of changes in fundamental interfacial and capillary properties. Experiments: We use MD simulations to explore the effects of four organic solutes (quinoline, decanoic acid, coronene, sorgoleone) on the wettability of quartz by water in the presence of scCO2. We examine the impacts of polar, alkyl, and aromatic moieties as well as fluid flow velocity at elevated temperatures and pressures. Findings: Organic molecules accumulate at the water-CO2 interface, where they distribute according to their size and functional groups. Certain organics penetrate the adsorbed water film at the quartz-CO2 interface, revealing two modes of hydrogen bonding between polar organic functional group, water, and quartz surface –OH groups. Interfacial energies and contact angles are minimally impacted by organic adsorption at the water-CO2 interface, possibly due to simultaneous CO2 desorption. Non- equilibrium MD simulations reveal that flow-induced redistribution of organic compounds modulates the radii of curvature of the advancing and receding water-CO2 interfaces. Our results indicate that organic adsorption on water surfaces is likely ubiquitous during multi-phase flow in soils and sedimen- tary rocks, with implications for the mobilization and transport of organic compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Material Plasma Exposure eXperiment High Heat Flux Microwave Absorber Design, Manufacture, and Articles Test

The Material Plasma Exposure eXperiment (MPEX) at Oak Ridge National Laboratory is in the final design phase. MPEX will be capable of exposing neutron-irradiated materials to plasmas for the study of plasma-material interaction. This facility will provide information about the complex effects of plasmas on materials and contribute to examining new materials that can withstand high heat fluxes and high ion fluences for future fusion devices. MPEX plasma is heated by 70-GHz or 105-GHz electron Bernstein wave/electron cyclotron heating (ECH), and the high-frequency microwaves are prone to scattering microwave power, which can have detrimental effects, especially on diagnostic components. A large portion of the injected ECH power is expected to be absorbed by plasma, but the remainder requires that microwave absorbers be placed immediately upstream and downstream of the ECH launcher to minimize stray microwaves leaving the ECH region. These microwaves can inadvertently heat components that cannot be shielded or otherwise protected. The microwave absorber design is based on an array of pyramid-shaped ceramic tiles brazed to a water-cooled explosion-bonded heat sink and a stainless steel plate to produce one tile module. Computational fluid dynamics and structural analyses were performed to optimize and validate the design. Multiple test coupons were produced to validate the process for brazing the two different tile materials to the Glidcop AL-15 baseplate. The articles were tested to evaluate the reliability and thermal performance through exposure to an electron beam with a heat flux of up to 1.5 MW/m2. Nondestructive testing was performed before and after testing to identify voids or separations that may have been introduced by the high heat flux. This paper discusses the details of high heat flux microwave absorber design, manufacturing details and associated challenges, and test results, demonstrating the effectiveness of the proposed design.

Hussain, Aftab↗
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Long Hole Film Cooling Dataset for CFD Development : Infrared Thermography and Thermocouple Surveys - Part 1

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (0.01905 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time-Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗
NASA NTRS2014

Long Hole Film Cooling Dataset for CFD Development - Flow and Film Effectiveness

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30 deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (approx. 0.02 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
DOE OSTI2011

Install or Extend Additional Infrastructure (SECARB Work Product 6.1.b)

Infrastructure needs at the SECARB “early” Cranfield test are different from many other CCS sites, because this test is strongly supported by the commercial activities of site host Denbury Onshore LLC. Commercial development by Denbury provided the substantive cost-share required to make the project viable. Denbury contributions fundamental for the project included the retrofit of the former Sonat natural gas pipeline from the Jackson Dome CO 2 source to the center of Cranfield Field, construction of distributory pipelines to take CO 2 to injection wells, and construction of roads and well pads that allowed access for drilling and maintenance of injection wells. In addition, as part of the EOR operations, Denbury constructed a CO 2 separation plant to increase pressure from pipeline pressure to field pressure and support CO 2 recycle, construction of production lines to bring oil, water, CO 2 and gas to the separation plant, and worked over plugged and abandoned (P&A) wells to retrofit them as producers. Producers were used by the SECARB project as monitoring points to observe flood. Denbury’s well test facility allowed monthly quantification of the volumes of fluids produced and reinjected.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗
DOE OSTI2015

Enhanced enstrophy generation for turbulent convection in low-Prandtl-number fluids

Turbulent convection is often present in liquids with a kinematic viscosity much smaller than the diffusivity of the temperature. Here we reveal why these convection flows obey a much stronger level of fluid turbulence than those in which kinematic viscosity and thermal diffusivity are the same; i.e., the Prandtl number Pr is unity. We compare turbulent convection in air at Pr = 0.7 and in liquid mercury at Pr = 0.021. In this comparison the Prandtl number at constant Grashof number Gr is varied, rather than at constant Rayleigh number Ra as usually done. Our simulations demonstrate that the turbulent Kolmogorov-like cascade is extended both at the large- and small-scale ends with decreasing Pr. The kinetic energy injection into the flow takes place over the whole cascade range. In contrast to convection in air, the kinetic energy injection rate is particularly enhanced for liquid mercury for all scales larger than the characteristic width of thermal plumes. As a consequence, mean values and fluctuations of the local strain rates are increased, which in turn results in significantly enhanced enstrophy production by vortex stretching. The normalized distributions of enstrophy production in the bulk and the ratio of the principal strain rates are found to agree for both Prs. Finally, despite the different energy injection mechanisms, the principal strain rates also agree with those in homogeneous isotropic turbulence conducted at the same Reynolds numbers as for the convection flows. Thus, our results have interesting implications for small-scale turbulence modeling of liquid metal convection in astrophysical and technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Simulation of Effects of the Saffman Force and the Magnus Force on Sand Saltation in Turbulent Flow

The effects of both the Saffman force and Magnus force on sand saltation are investigated. Turbulent flows in a channel and over a barchans dune are considered with sand particles injected into the flow. The results show that both of the forces increase the height and skipping distance of sand saltation, with the Magnus force giving more significant effect on the height. These forces can also increase the sand settling at the lee side of the barchans dune.

Computational fluid dynamics↗
DOE OSTI2024

Application of quantitative risk assessment to address stakeholder questions in geologic carbon storage

Ambitious international greenhouse gas emissions reduction targets demand a rapid transformation to a low-carbon economy. This transformation includes the accelerated adoption of carbon dioxide (CO2) capture and storage (CCS) technology. However, as with any large-scale engineering enterprise, the widespread commercial-scale deployment of geologic carbon storage (GCS) raises important questions about technology and cost-effectiveness, safety, environmental risk, and long-term liability. Effectively assessing and managing risks and liability associated with GCS projects is a key technical need throughout the project life cycle-from site selection and permitting to monitoring design, operational risk management, and post-operational site closure. This presentation highlights recent advancements in tools for quantitative risk assessment, being developed by the National Risk Assessment Partnership (NRAP). NRAP is a multi-year, multinational laboratory research collaboration sponsored by the U.S. Department of Energy's Office of Fossil Energy and Carbon Management. Our focus will be on these tools' applications in addressing critical stakeholder questions related to supporting permitting to ensure secure and environmentally protective storage; designing effective and efficient monitoring plans; evaluating the effectiveness of remedial actions and risk management alternatives; and informing liability assessment and investment decisions. This paper will detail the key functionality of NRAP’s Open-Source Integrated Assessment Model (NRAP-Open-IAM), a computational framework for assessing leakage risk and containment assurance. This model features streamlined workflows for calculating leakage risk profiles, delineating risk-based area of review, and assessing contingency plans and post-injection site care requirements. ORION is an open-source, observation-based ensemble forecasting toolkit to help operators assess the seismic hazard at a carbon storage site. The State of Stress Analysis Tool (SOSAT), designed to assess subsurface stress conditions and evaluate geomechanical risk resulting from CO2 injection in an area of interest will also be presented. We will also introduce a prototype model to evaluate storage project costs and liability associated with risk management. The Technoeconomic and Liability Evaluation for Storage (TALES) model uses results from forecasts of leakage and induced seismicity risk to estimate the lifecycle cost of managing risk. Finally, a preliminary example of how the NRAP Risk-based Adaptive Monitoring Plan (RAMP) tool can be used to design efficient and effective site monitoring plans and estimate the detectability of fluid leakage will be provided. The relevance of these tools for addressing key stakeholder questions amidst uncertainty will be emphasized.

decision support↗
NASA NTRS1993

Experimental investigation of turbine disk cavity aerodynamics and heat transfer

An experimental investigation of turbine disk cavity aerodynamics and heat transfer was conducted to provide an experimental data base that can guide the aerodynamic and thermal design of turbine disks and blade attachments for flow conditions and geometries simulating those of the space shuttle main engine (SSME) turbopump drive turbines. Experiments were conducted to define the nature of the aerodynamics and heat transfer of the flow within the disk cavities and blade attachments of a large scale model simulating the SSME turbopump drive turbines. These experiments include flow between the main gas path and the disk cavities, flow within the disk cavities, and leakage flows through the blade attachments and labyrinth seals. Air was used to simulate the combustion products in the gas path. Air and carbon dioxide were used to simulate the coolants injected at three locations in the disk cavities. Trace amounts of carbon dioxide were used to determine the source of the gas at selected locations on the rotors, the cavity walls, and the interstage seal. The measurements on the rotor and stationary walls in the forward and aft cavities showed that the coolant effectiveness was 90 percent or greater when the coolant flow rate was greater than the local free disk entrainment flow rate and when room temperature air was used as both coolant and gas path fluid. When a coolant-to-gas-path density ratio of 1.51 was used in the aft cavity, the coolant effectiveness on the rotor was also 90 percent or greater at the aforementioned condition. However, the coolant concentration on the stationary wall was 60 to 80 percent at the aforementioned condition indicating a more rapid mixing of the coolant and flow through the rotor shank passages. This increased mixing rate was attributed to the destabilizing effects of the adverse density gradients.

Daniels, W. A.↗
NASA NTRS1993

Development of an Atomization Methodology for Spray Combustion

In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.

Seung, S. P.↗
NASA NTRS2002

Propulsion Simulations with the Unstructured-Grid CFD Tool TetrUSS

A computational investigation has been completed to assess the capability of the NASA Tetrahedral Unstructured Software System (TetrUSS) for simulation of exhaust nozzle flows. Three configurations were chosen for this study: (1) a fluidic jet effects model, (2) an isolated nacelle with a supersonic cruise nozzle, and (3) a fluidic pitchthrust- vectoring nozzle. These configurations were chosen because existing data provided a means for measuring the ability of the TetrUSS flow solver USM3D for simulating complex nozzle flows. Fluidic jet effects model simulations were compared with structured-grid CFD (computational fluid dynamics) data at Mach numbers from 0.3 to 1.2 at nozzle pressure ratios up to 7.2. Simulations of an isolated nacelle with a supersonic cruise nozzle were compared with wind tunnel experimental data and structured-grid CFD data at Mach numbers of 0.9 and 1.2, with a nozzle pressure ratio of 5. Fluidic pitch-thrust-vectoring nozzle simulations were compared with static experimental data and structured-grid CFD data at static freestream conditions and nozzle pressure ratios from 3 to 10. A fluidic injection case was computed with the third configuration at a nozzle pressure ratio of 4.6 and a secondary pressure ratio of 0.7. Results indicate that USM3D with the S-A turbulence model provides accurate exhaust nozzle simulations at on-design conditions, but does not predict internal shock location at overexpanded conditions or pressure recovery along a boattail at transonic conditions.

Deere, Karen A.↗
NASA NTRS1979

Autoignition of hydrogen injected transverse to a supersonic airstream

A parametric study of autoignition characteristics of transverse fuel injector patterns for a hydrogen-fueled, supersonic-combustion ramjet was conducted in direct-connect tests over a range of simulated hypersonic flight conditions. Autoignition limits, defined directly by visual observations and indirectly by measured wall pressure and temperature, are compared with a semi-empirical ignition limit model which accounts for first order effects on autoignition. It is shown that to initiate combustion in the immediate region of the fuel injector, ignition must occur in the separation region ahead of the fuel jet. Autoignition within this separation zone exhibits strong dependence on flow pressure and temperature, injector diameter, local boundary layer energy thickness, wall temperature, and on fluid dynamic interactions with adjacent downstream facing steps (such as inlet isolation steps). Secondary effects are also discussed. The empirical model is used to define limits of test temperature, pressure, and scale required for autoignition in a current subscale engine test program.

Mcclinton, C. R.↗
DOE OSTI2018

KLF2 attenuates bleomycin-induced pulmonary fibrosis and inflammation with regulation of AP-1

Highlights: • KLF2 is downregulated in BLM-treated rats. • KLF2 alleviates BLM-induced lung fibrosis. • KLF2 suppresses the production of pro-inflammatory cytokines. • KLF2 prevents TGF-β1 induced lung fibrosis in vitro. • Effects of KLF2 on AP-1. Pulmonary fibrosis (PF) is a chronic, fibrosing interstitial pneumonia and devastating disease. Here we investigated the potential roles of Kruppel-like factor 2 (KLF2) on pulmonary fibrosis and inflammation response. A mouse model of pulmonary fibrosis was established by intratracheal injection of bleomycin (BLM). The mRNA and protein levels of KLF2 were assayed by RT-PCR and Western blotting respectively. The extent of lung fibrosis was determined using hematoxylin and eosin (HE) staining and Masson's trichrome staining, and the hydroxyproline content was quantified. RT-PCR was used to evaluate the mRNA expression of collagen type 1a1 (col1a1), col3a1, α-SMA, TNF-α, IL-1β and IL-6. The concentrations of TNF-α, IL-1β, and IL-6 in bronchoalveolar lavage fluid (BALF) and lung tissue were examined by ELISA. Also, the effects of KLF2 on activator protein-1 (AP-1) were evaluated by measuring the c-Jun and c-Fos protein levels. We found that KLF2 was remarkably downregulated in BLM-treated rats, both in mRNA and protein levels. Additionally, overexpression of KLF2 attenuated the destruction of the alveolar space and pulmonary interstitial collagen hyperplasia, and deposition reduced the expression of col1a1, col3a1, and α-SMA, and blocked the production of TNF-α, IL-1β, and IL-6 in BALF and lung tissue in vivo. Moreover, adenoviral transduction of KLF2 inhibited TGF-β1-induced expression of col1a1, col3a1, and α-SMA in vitro. Mechanically, BLM up-regulated c-Jun and c-Fos expression, which was impeded by KLF2 overexpression. Taken together, our data indicate that KLF2 attenuates pulmonary fibrosis and inflammation, possibly through the regulation of AP-1.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Geomechanical effects on CO 2 leakage through fault zones during large-scale underground injection

The importance of geomechanics—including the potential for faults to reactivate during large-scale geologic carbon sequestration operations—has recently become more widely recognized. However, notwithstanding the potential for triggering notable (felt) seismic events, the potential for buoyancy-driven CO 2 to reach potable groundwater and the ground surface is actually more important from public safety and storage-efficiency perspectives. In this context, this paper extends the previous studies on the geomechanical modeling of fault responses during underground carbon dioxide injection, focusing on the short-term integrity of the sealing caprock, and hence on the potential for leakage of either brine or CO 2 to reach the shallow groundwater aquifers during active injection. We consider stress/strain-dependent permeability and study the leakage through the fault zone as its permeability changes during a reactivation, also causing seismicity. We analyze several scenarios related to the volume of CO 2 injected (and hence as a function of the overpressure), involving both minor and major faults, and analyze the profile risks of leakage for different stress/strain-permeability coupling functions. We conclude that whereas it is very difficult to predict how much fault permeability could change upon reactivation, this process can have a significant impact on the leakage rate. Moreover, our analysis shows that induced seismicity associated with fault reactivation may not necessarily open up a new flow path for leakage. Results show a poor correlation between magnitude and amount of fluid leakage, meaning that a single event is generally not enough to substantially change the permeability along the entire fault length. Finally, and consequently, even if some changes in permeability occur, this does not mean that the CO 2 will migrate up along the entire fault, breaking through the caprock to enter the overlying aquifer.

54 ENVIRONMENTAL SCIENCES↗