Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Comparative Experimental Analysis of Natural Gas Dual Fuel Combustion Ignited by Diesel and Poly OxyMethylene Dimethyl Ether

Dual-fuel low-temperature combustion is a possible solution for alleviating the tradeoff between oxides of nitrogen and soot emissions in conventional diesel combustion, albeit with poor combustion stability, high carbon monoxide, and unburned hydrocarbon emissions at low engine loads. The present work compares emissions and combustion (heat release and other metrics) of both diesel and poly-oxy methylene dimethyl ether as high-reactivity fuels to ignite natural gas while leveraging spray-targeted reactivity stratification, which involved multiple injections of the high-reactivity fuels. The experiments included six parametric sweeps of: (1) start of first injection, (2) start of second injection, (3) percentage of energy substitution of natural gas, (4) commanded injection duration ratio, (5) rail pressure, and (6) intake pressure. The experiments were performed on a 1.8 L heavy-duty single-cylinder research engine operating at a medium speed of 1339 rev/min. Not-to-exceed limits for the indicated oxides of nitrogen emissions, maximum pressure rise rate, and the coefficient of variation of the indicated mean effective pressure were set to 1 g/kWh, 10 bar/CAD, and 10%, respectively. The indicated emissions decreased and combustion improved significantly for both fueling combinations when the experimental procedure was applied.

Partridge, Kendyl Ryan↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2021

Analysis and design of DC-link capacitor in cascaded H-bridge multilevel active front-end converter

Medium-voltage grid-tied systems often use a cascaded H-bridge multi-level active front-end. In this converter, dc link bus capacitors play an important role in stabilizing the converter and enabling both active and reactive power injections. The present work provides analytical expressions for the capacitor current, which are essential for optimizing system design (especially capacitor size vs. lifetime). Then, the expression is incorporated into the grid connected bidirectional power system model. Consequently, this work contributes to the guiding principles to choose accurate dc link capacitor ratings against grid-side power delivery requirements. DC-link capacitors come with an expiration date which is dependent on the operating power and voltage values. Lack of information about the lifetime of the dc-link capacitor creates uncertainty in the duration of online operation of the converter, thus increasing the probability of a contingency outage. The lifetime of the capacitor can range from only a few seconds of online operation of the station due to under-sizing, to a compromised quality of power flowing through other solid state devices which may not have been designed for the transients that are caused by over-sizing of the dc-link capacitor. A better estimate of the lifetime of these capacitors can help maximize usage and maintain scheduled outages more efficiently, without disruption in service.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Epitaxial heterointerfacial electron bridge synchronizes oxygen evolution activity and stability on a layered double hydroxide surface

Scalable green hydrogen production via electrocatalytic water splitting is largely restricted by the insufficient activity and stability of oxygen evolution reaction (OER) catalysts at the anode. As a class of the most active OER catalysts in alkaline electrolyzers, the application of layered double hydroxides (LDHs) remains a main challenge owing to the unstable lattice oxygen dissolution due to the dominant lattice oxygen-involving OER mechanism during long-term operation. Herein, we found that using an epitaxial hetero-interfacing nickel hydroxide (namely Ni(OH) 2 ) as an electron bridge between an active FeCo LDH and Ni foam support to form an LDH*/NFO catalyst, the electronic storage capacity around the Fermi level (-0.5 to +0.5 eV, e-D FE ) sharply increases from 0.93 per cell to 1.51 per cell. Subsequently, we demonstrate that this high e-D FE enables ceaseless and fast power injection into the kinetic process of intermediate species conversion and inhibits lattice oxygen dissolution in the active FeCo LDH. Consequently, it demonstrated a low OER overpotential of 246 mV at a current density of 100 mA cm -2 and ultrahigh stability for up to 3500 hours with an ultraslow overpotential increase rate of 9.4 × 10 -3 mV h -1 . Therefore, we developed an epitaxial hetero-interfacial electron bridging strategy to synchronize the activity and stability of available catalysts for scalable green hydrogen production via electrocatalytic water splitting.

25 ENERGY STORAGE↗
DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗
DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Examination of entrainment, combustibility, and heat transfer due to cell venting inside simplified battery energy storage enclosures

Battery energy storage systems (ESS) assembled of modules containing lithium-ion batteries can pose fire and explosion hazards during thermal runaway because vented gases may accumulate, combust, or escape enclosures. Here, we propose a prototypical configuration composed of the main components of an ESS rack, where modules are represented as rectangles within an enclosure and vent gas is injected steadily through the top of a module. This idealized geometry allows for exploration of how rack-level characteristics can influence smaller single-cell or module-level scales. Through this geometry, we conduct numerical simulations to understand how variations in the vent velocity, temperature, and vertical location affect entrainment rates, fraction of uncombusted fuel, rich-gas coverage, and heat flux to other modules. At lower vent velocities and vertical positions, buoyancy drives enough entrainment and mixing to react with most of the vent gases on the converging plenum. As jet momentum increases and/or buoyancy decreases via higher vent position in the rack, vent gases can propagate into and burn in the divergence plenum. Heat flux to the impinged module scales strongly with vent momentum and temperature, while heat transfer to upper modules peaks when flames in the diverging plenum reach the surface and declines when combustion is suppressed.

Computational fluid dynamics↗
DOE OSTI2026

Characterization of high-flux cadmium zinc telluride detectors coupled to MM-PAD readout ASICs for high-flux X-ray applications

High-flux cadmium zinc telluride (HF-CZT) is a promising material for next-generation X-ray imaging due to its high stopping power, wide bandgap, and improved resistance to radiation-induced polarization [1]. We bonded 2 mm thick HF-CZT sensors with Cornell's Mixed-Mode Pixel Array ASICs, which combine analog charge integration and in-pixel digital counting to achieve single X-ray sensitivity, large full well capacity, and frame rates up to 10 kHz [2,3]. Detector characterization showed stable electron collection above 150 V bias, while hole collection remained limited by trapping even at 700 V. Temporal studies revealed non-negligible afterglow, with an unexpected increase in afterglow at higher biases, which was reduced by infrared illumination, suggesting that injected carriers reduce trapping near contact interfaces. Polarization measurements confirmed improved stability compared to conventional CZT, though polarization effects persisted, particularly with hole collection, with up to 5% signal loss and spatial distortion at high doses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Soil Environmental Conditions in the East River Watershed, CO

This dataset contains depth-resolved soil environmental variables collected between June 2016 and October 2018 from the Bradley Meadow site in the East River watershed. Datalogger output from 7 Campbell CS655 12 cm Soil Moisture and Temperature Sensors including depth, temperature, volumetric water content, and electrical conductivity are included in ‘BM_Merged_T_VWC_0616_10_18.csv’. Sensors were installed in June 2016. Data have been processed to merge raw sensor data and to remove additional variables. No additional QA/QC was performed. Measured pCO2 and depth are contained in ‘BM_EGM_WELL_CO2.csv’. These data represent gas well samples measured via injection using a PP Systems EGM-5 infra-red gas analyzer as described in the methods section. Measured CO2 surface flux values are contained in 'BM_EGM_Surface_Flux.csv'. These were measured using the PP Systems EGM-5 with the Soil Respiration Chamber attachment. Multiple reps are reported for each measurement including CO2 flux rate calculations using built-in linear and quadratic function methods, and average and standard deviations are taken across all reps and calculation methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2022

Correlation between the relative blob fraction and plasma parameters in NSTX

Here, we present a simple analysis of the fraction of blobs observed in the scrape-off layer of NSTX and calculates the correlation between this blob fraction and various plasma parameters. The measurements were made using the gas puff imaging (GPI) diagnostic on a database of 103 shots from 2010 with neutral beam power P NBI = 0–6 MW. The blob fraction ranged from f b = 0.1% to 4.8%, where f b is defined as the fraction of time the GPI signal was larger than three times its average value within 0–7.5 cm outside the separatrix. The blob fraction generally decreased with increasing neutral beam injection power and was lower on average in H-mode plasmas than Ohmic and L-mode plasmas. The blob fraction had its highest correlation with the measured poloidal turbulence velocity and its radial gradient just inside the separatrix. This is in part consistent with a model in which the blob generation rate increases with the poloidal velocity shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Synergy of Coupled Kink and Tearing Modes in Fast Ion Transport

One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the ''kick" model has been developed to compute the evolution of the fast ion distribution from the neutral beam injection using instabilities as phase space resonance sources. The goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft X-ray diagnostics are used to identify the mode parameters in NSTX. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transpor

Energetic Particles↗
DOE OSTI2021

Resonant and random excitations on the proton beam in the Large Hadron Collider for active halo control with pulsed hollow electron lenses

We present the results of numerical simulations and experimental studies about the effects of resonant and random excitations on proton losses, emittances, and beam distributions in the Large Hadron Collider (LHC). In addition to shedding light on complex nonlinear effects, these studies are applied to the design of hollow electron lenses (HEL) for active beam halo control. In the High-Luminosity Large Hadron Collider (HL-LHC), a considerable amount of energy will be stored in the beam tails. To control and clean the beam halo, the installation of two hollow electron lenses, one per beam, is being considered. In standard electron-lens operation, a proton bunch sees the same electron current at every revolution. Pulsed electron beam operation (i.e., different currents for different turns) is also considered, because it can widen the range of achievable halo removal rates. For an axially symmetric electron beam, only protons in the halo are excited. If a residual field is present at the location of the beam core, these particles are exposed to time-dependent transverse kicks and to noise. We discuss the numerical simulations and the experiments conducted in 2016 and 2017 at injection energy in the LHC. The excitation patterns were generated by the transverse feedback and damping system, which acted as a flexible source of dipole kicks. Proton beam losses, emittances, and transverse distributions were recorded as a function of excitation patterns and strengths. The resonant excitations induced rich dynamical effects and nontrivial changes of the beam distributions, which, to our knowledge, have not previously been observed and studied in this detail. We conclude with a discussion of the tolerable and achievable residual fields and proposals for further studies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

UV-Induced Degradation and Associated Metastability in TOPCon Photovoltaic Modules: Understanding Kinetics and Cell Variance

Tunnel oxide passivated contact (TOPCon) silicon photovoltaic (PV) modules are dominating the PV market, but they may be susceptible to degradation under ultraviolet (UV)-containing light. Quantifying the impacts of UV-induced degradation (UVID) is complicated by an associated metastability causing further degradation under dark storage and rapid recovery under sunlight. Here, we study modules that have -2.3% to -3.2% nonrecoverable UVID loss after 60 kWh/m2 dose of 340 nm light and additional recoverable loss under dark storage. We use in situ electroluminescence (EL) imaging to characterize the post-UVID metastability at the module level. The cell-by-cell dark degradation and recovery kinetics span a wide range from +6% to -70% changes in EL intensity after 520 h of dark storage, which returns to +- 4% of the initial post-UVID EL intensity after illumination. The kinetics follow double exponential rates with dark storage degradation time constants of 345 and 45 h, and UV light recovery time constants of 5 min and 36 s. We propose that this is consistent with prior reports of kinetics for light-soaking metastability in Al2O3 passivation. Finally, we further show that cells having high UVID also have injection-dependent effective carrier lifetimes and significant intra-cell variance, suggesting possible origins of processing inconsistency.

14 SOLAR ENERGY↗
DOE OSTI2022

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗
DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗
DOE OSTI2020

Impacts of Mixed-Wettability on Brine Drainage and Supercritical CO 2 Storage Efficiency in a 2.5-D Heterogeneous Micromodel

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO 2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO 2 (scCO 2 ) is used in enhanced oil recovery. Here, we performed scCO 2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from –8.1 to –4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO 2 saturation increase at both capillary fingering and crossover flow regimes (–8.1 ≤ logCa ≤ –6.1), the one heterogeneous wetting to scCO 2 results in ~0.09 saturation increase only at the crossover flow regime (–7.1 ≤ logCa ≤ –6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO 2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > –6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO 2 saturation and distribution suggest the importance of wettability on CO 2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

58 GEOSCIENCES↗