Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗
DOE OSTI2024

The Potential of Stratospheric Aerosol Injection to Reduce the Climatic Risks of Explosive Volcanic Eruptions

Sulfur-rich volcanic eruptions happen sporadically. If Stratospheric Aerosol Injection (SAI) were to be deployed, it is likely that explosive volcanic eruptions would happen during such a deployment. Here we use an ensemble of Earth System Model simulations to show how changing the injection strategy post-eruption could be used to reduce the climate risks of a large volcanic eruption; the risks are also modified even without any change to the strategy. For a medium-size eruption (10 Tg-SO 2 ) comparable to the SAI injection rate, the volcanic-induced cooling would be reduced if it occurs under SAI, especially if artificial sulfur dioxide injections were immediately suspended. Alternatively, suspending injection only in the eruption hemisphere and continuing injection in the opposite would reduce shifts in precipitation in the tropical belt and thus mitigate eruption-induced drought. Finally, we show that for eruptions much larger than the SAI deployment, changes in SAI strategy would have minimal effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗
DOE OSTI2021

Characterization of injection and confinement improvement through impurity induced profile modifications on the Wendelstein 7-X stellarator

Pulsed injections of boron carbide granules into Wendelstein 7-X stellarator (W7-X) plasmas transiently increase the plasma stored energy and core ion temperatures above the reference W7-X experimental programs by up to 30%. In a series of 4 MW electron cyclotron resonance heating experiments, the PPPL Probe Mounted Powder Injector provided 50 ms bursts of 100 μm granules every 350 ms at estimated quantities ranging from approximately 1 mg/pulse to over 30 mg/pulse. For each injection, the stored energy was observed to initially drop and the radiated power transiently increased, while the radial electron density profile rose at the edge as material was assimilated. Once the injected boron carbide was fully absorbed, the density rise transitioned to the core while the stored energy increased above the previous baseline level by an amount linearly correlated with the injection quantity. During the injection, the ion temperature gradient steepened with peak core ion temperatures observed to increase from a nominal 1.7 keV to over 2.6 keV for the largest injection amounts. Enhanced performance is accompanied by a reversal of the radial electric field at ρ < 0.3, indicating that the core transport has switched to the ion root. These observations are suggestive of a change in transport and provide further evidence that externally induced profile modifications provide a possible path to enhanced W7-X performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The effect of pellet injection on turbulent transport in TJ-II

An improvement in confinement associated with the injection of pellets has been observed in TJ-II and other devices. By means of a simple model, we study here the modification of turbulent transport by the injection of pellets in TJ-II and how this modification affects the particle confinement time. The analysis of the results shows the relationship of the confinement improvement with the evolution of the shear flow due to turbulence, especially near low order rational surfaces. This work is based on pellet injections into the neutral beam injection heated phase of TJ-II plasmas. For this, single, double, or triple hydrogen pellet injection are considered and the post-injection evolution of selected plasma parameters are analyzed and compared.

Physics↗
DOE OSTI2025

Impact of threading dislocations on the V-defect assisted lateral carrier injection and recombination in InGaN quantum well LEDs

The nonuniform hole distribution between InGaN quantum wells (QWs) of light emitting diodes (LEDs) has a negative impact on LED efficiency. The uniformity can be increased by using lateral hole injection through sidewalls of V-defects, which form at threading dislocations. However, the inherent coupling between the V-defects and dislocations might affect efficiency of the hole injection and nonradiative recombination. In this work, we have tested the possible impact of the dislocations on the injection and recombination by means of scanning near-field electroluminescence and photoluminescence spectroscopy on single green-emitting InGaN QW LEDs containing large (∼0.5 μm) V-defects. The measurements have not provided any evidence of a lower hole injection efficiency or enhanced nonradiative recombination at the dislocations located at the V-defect facets or their apexes. This shows that large V-defects are excellent volumetric injectors for long wavelength InGaN LEDs. Furthermore, it was established that V-defects are preferential hole injectors even in single quantum well devices. Compared to vertical injection, the V-defect injection allows lowering the operating voltage, which should contribute to an enhanced wall plug efficiency.

42 ENGINEERING↗
DOE OSTI2020

Numerical simulation of Li-pellet injection experiments for ELM-pacing in EAST

The injection of light impurity granules into fusion research discharges has shown a great potential for ELM pacing in future nuclear devices. The present work employed a model considering both pellet ablation and transport to systematically study for the first time the evolution of both impurity Li species and target plasma in the edge plasma during single and repeated Li pellet injections, with all initial information extracted from the EAST experiments. The evolutionary distributions of atomic and ionized Li species were presented and show good consistency with experimental CCD camera images. Simulation studies were carried out for single and repeated pellet injections at the outer mid-plane, and statistical uncertainty of ELM pacing rates triggering ELMs is discussed with regard to related injection conditions. It is found that the pressure locally caused by an ELM pacing pellet (with a diameter of 0.6 ~ 0.7 mm) was 1.8–2.4 times its unperturbed magnitude in the EAST shots of interest. Accumulation from repeated pellet injections results in the rise of pedestal height at a high injection frequency and this effect becomes stronger with increasing pellet size. Pressure changes are also evident in a scenario of multiple small pellets per stroke, which could lead to a nonnegligible, probabilistic ELM triggering rate observed experimentally. The effects of other possible factors are also discussed. With the assumption that the ELM trigger rate of a Li pellet of 0.7 mm diameter is 100%, the predicted trigger rates are in good agreement with those observed in experiments. This work can help understand the statistically probabilistic behavior of ELM triggering efficiency observed in previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Toroidal injection angle dependence of EC assisted plasma initiation at DIII-D

An experimental scan of the electron cyclotron waves (EC) toroidal injection angle in plasma breakdown is performed at the DIII-D tokamak. The second harmonic, extraordinary mode EC is used for the study. The dependence of n e and T e on the EC injection angle cannot be conclusively drawn from this study due to the large error bars in the n e and T e measurements. On the other hand, high T e data points are observed in some discharges which can be explained by nonlinear heating. The D α emission measurement shows a clear relation between the breakdown time and the injection angle. An experimental investigation of the cause of the dependence of breakdown delay on the EC injection angle suggests that when the injection angle is oblique, the EC heating after the reflection at the inboard wall may become ineffective and cause the breakdown delay even when the EC heating before and directly upon injection remains effective. A preliminary run of the heat and transport balance code DYON indicates that the obtained dataset is suitable for a quantitative validation of EC absorption models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Imaging of a fluid injection process using geophysical data — A didactic example

In many subsurface industrial applications, fluids are injected into or withdrawn from a geologic formation. It is of practical interest to quantify precisely where, when, and by how much the injected fluid alters the state of the subsurface. Routine geophysical monitoring of such processes attempts to image the way that geophysical properties, such as seismic velocities or electrical conductivity, change through time and space and to then make qualitative inferences as to where the injected fluid has migrated. The more rigorous formulation of the time-lapse geophysical inverse problem forecasts how the subsurface evolves during the course of a fluid-injection application. Using time-lapse geophysical signals as the data to be matched, the model unknowns to be estimated are the multiphysics forward-modeling parameters controlling the fluid-injection process. Properly reproducing the geophysical signature of the flow process, subsequent simulations can predict the fluid migration and alteration in the subsurface. The dynamic nature of fluid-injection processes renders imaging problems more complex than conventional geophysical imaging for static targets. This work intents to clarify the related hydrogeophysical parameter estimation concepts.

58 GEOSCIENCES↗
DOE OSTI2023

Process Variability Effects on Tensile Response in Injection Molded, Fluorinated Thermoplastics

The mechanical properties of fluorinated thermoplastics (i.e., tensile strength and elongation) can vary with changes in injection molding processing parameters. Four fluoropolymers are examined: poly(vinylidene fluoride) (PVDF) and random poly(vinylidene fluoride-co-chlorotrifluoroethylene) (PVDF-CTFE) with three CTFE concentrations. Dog bones were manufactured with various cylinder dwell times and mold cooling times to assess the manufacturing sensitivity to the tensile response. Dwell and cooling times increasingly impact mechanical performance as CTFE concentration increases. Specimens exhibit higher tensile strength as a function of injection order. The first injected specimen exhibits the lowest tensile strength and highest elongation in all copolymers. This trend becomes more pronounced among fluoropolymers with higher CTFE concentration and lower weight-averaged molecular weight. Parallel plate rheology was used to obtain the zero-shear viscosity as a function of material type, process, and injection order. We found that in the copolymers, the first injected sample exhibited a lower zero-shear viscosity than the next, which indicates a lower molecular weight in the first injected specimen. This phenomenon was not present for the PVDF homopolymer. Copolymer mechanical uncertainties are hypothesized to result from the shorter molecular weight chains extruding out of the specimens’ sides as a flash due to higher mobility with CTFE segments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Progress of Gas Injection EOR Surveillance in the Bakken Unconventional Play—Technical Review and Machine Learning Study

Although considerable laboratory and modeling activities were performed to investigate the enhanced oil recovery (EOR) mechanisms and potential in unconventional reservoirs, only limited research has been reported to investigate actual EOR implementations and their surveillance in fields. Eleven EOR pilot tests that used CO2, rich gas, surfactant, water, etc., have been conducted in the Bakken unconventional play since 2008. Gas injection was involved in eight of these pilots with huff ‘n’ puff, flooding, and injectivity operations. Surveillance data, including daily production/injection rates, bottomhole injection pressure, gas composition, well logs, and tracer testing, were collected from these tests to generate time-series plots or analytics that can inform operators of downhole conditions. A technical review showed that pressure buildup, conformance issues, and timely gas breakthrough detection were some of the main challenges because of the interconnected fractures between injection and offset wells. The latest operation of co-injecting gas, water, and surfactant through the same injection well showed that these challenges could be mitigated by careful EOR design and continuous reservoir monitoring. Reservoir simulation and machine learning were then conducted for operators to rapidly predict EOR performance and take control actions to improve EOR outcomes in unconventional reservoirs.

Energy & Fuels↗
DOE OSTI2023

Performance, Combustion and Emissions Evaluation of Liquid Phase Port-Injected LPG on a Single Cylinder Heavy-Duty Spark Ignited Engine

Liquefied petroleum gas (LPG), like many other alternative fuels, has witnessed increased adoption in the last decade, and its use is projected to rise as stricter emissions regulations continue to be applied. However, much of its use is limited to dual fuel applications, gaseous phase injection, light-duty passenger vehicle applications, or scenarios that require conversion from gasoline engines. Therefore, to address these limitations and discover the most efficient means of harnessing its full potential, more research is required in the development of optimized fuel injection equipment for liquid port and direct injection, along with the implementation of advanced combustion strategies that will improve its thermal efficiency to the levels of conventional fuels. This paper focuses on the development of a liquid phase port-injection system for LPG, the design of a reference piston, and the baseline evaluation of the performance, combustion, and emissions characteristics of a single cylinder research engine to establish a benchmark comparable to existing LPG engines. A sweep of start of injection (SOI) timing is performed by injecting liquid LPG at several closed and open intake valve timings, which demonstrates no significant variation in engine performance, but accounts for a 10% reduction in bsCO with the optimal SOI timing. Spark timing sweep demonstrates the 50% burn crank angle location related to maximum brake torque (MBT) point with a brake thermal efficiency (BTE) of ~34% for the tested load case. The effect of equivalence ratio is also presented with optimal SOI timing at MBT condition. The engine starts exhibiting knocking combustion at 140kPa intake manifold air pressure (IMAP) with a peak torque of 253Nm and a 5% reduction in brake specific fuel consumption compared to the naturally aspirated scenario.

Fosudo, Toluwalase↗
DOE OSTI2020

Results of Multiple Tracer Injections into Fractures in the EGS Collab Testbed-1

The EGS Collab project constructed an intermediate scale (~10-20 m) testbed at the 4850 level of the Stanford Underground Research Facility (SURF) in South Dakota for testing and validating fracture stimulation and flow/transport models. This testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the crystalline rocks of the Poorman Formation from the West Access Drift tunnel. Of the eight boreholes, one borehole is used as an injection/stimulation well, while another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as geophysical/fluid sampling monitoring wells. Hydraulic stimulation activities were conducted at three locations along the injection hole in an attempt to create direct fracture connections to the production hole. A flow system has been established between injection and production boreholes through a set of hydraulically stimulated fractures propagated from a notch located at 164 ft in the injection hole. Although we planned for a single production well, a set of natural fractures in the testbed is believed to have intersected the stimulated hydraulic fractures and provided additional flow paths for water to be transported to the drift through multiple monitoring boreholes and weep zones. As the flow tests continued after stimulation activities, single or a combination of two or three producers become dominant producers at different times, mostly as a response to intersection of hydraulic fracture and testbed wells, activation of natural fractures and making new leak points to the monitoring wells, resealing of leaky wells, and so on. Since late October 2018, multiple tracers were injected into the fracture system at the 164 ft location that involves both stimulated and natural fractures, and tracers were recovered from multiple locations in nearby wells and weep. The cumulative water recoveries over the time have ranged from 50 to 90%; however, the injected tracer recoveries from various production sources are much less (ranging from a few percentages to 38%). Changes in water/tracer recoveries, shifting of major producing wells from one well to the others, and other observations (e.g., microseismic, electric resistivity tomography, etc.) indicate a testbed that has undergone several changes since October 2018. In this paper, we present tracer recovery data accumulated during several tracer campaigns, provide conceptual testbed flow pathways, and provide observations that suggest a evolutionary nature of fracture volume and fracture geometry in the testbed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Fuel injection module for segmented annular combustion system

The present disclosure is directed to a fuel injection module for a segmented annular combustion system. The fuel injection module includes a housing body, a fuel nozzle portion, and at least one fuel injection lance. The fuel nozzle portion is fluidly coupled to a fuel nozzle plenum within the housing body, and the at least one fuel injection lance is fluidly coupled to an injector fuel plenum within the housing body. In some cases, the fuel nozzle portion is a bundled tube fuel nozzle having one or more subsets of tubes. The fuel injection lances are positioned along a radial side of the housing body or circumferentially between two subsets of tubes. Liquid fuel cartridges extend through the fuel nozzle portion, the fuel injection lances, or both.

42 ENGINEERING↗
DOE OSTI2025

Injection Locking of Gigahertz‐Frequency Surface Acoustic Wave Phononic Crystal Oscillator

Low-noise gigahertz (GHz) frequency sources are essential for applications in signal processing, sensing, and telecommunications. Surface acoustic wave (SAW) resonator-based oscillators offer compact form factors and low-phase noise due to their short mechanical wavelengths and high-quality (Q) factors. However, their small footprint makes them vulnerable to environmental variation, resulting in their poor long-term frequency stability. Injection locking is widely used to suppress frequency drift of lasers and oscillators by synchronizing to an ultra-stable reference. Here, injection locking of a 1-GHz SAW phononic-crystal oscillator is demonstrated, achieving 40-dB phase noise reduction at low offset frequencies and unperturbed low noise at large offset frequencies. Compared to a free-running SAW oscillator, which typically exhibits frequency drifts of several hundred hertz over minutes, the injection-locked oscillator reduces the frequency deviation to below 0.35 Hz. The locking range and oscillator dynamics is also investigated in the injection pulling region. The demonstrated injection-locked SAW oscillator could find applications in high-performance portable telecommunications and sensing systems.

injection locking↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗