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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Estimating the Impact of Seep Methane Oxidation on Ocean pH and Dissolved Inorganic Radiocarbon Along the U.S. Mid-Atlantic Bight

Ongoing ocean warming can release methane (CH 4 ) currently stored in ocean sediments asfree gas and gas hydrates. Once dissolved in ocean waters, this CH 4 can be oxidized to carbon dioxide(CO 2 ). While it has been hypothesized that the CO 2 produced from aerobic CH4oxidation could enhanceocean acidification, a previous study conducted in Hudson Canyon shows that CH 4 oxidation has a smallshort-term influence on ocean pH and dissolved inorganic radiocarbon. Here we expand upon thatinvestigation to assess the impact of widespread CH 4 seepage on CO 2 chemistry and possible accumulationof this carbon injection along 234 km of the U.S. Mid-Atlantic Bight. Consistent with the estimates fromHudson Canyon, we demonstrate that a small fraction of ancient CH 4 -derived carbon is being assimilatedinto the dissolved inorganic radiocarbon (mean fraction of 0.5 ± 0.4%). The areas with the highest fractionsof ancient carbon coincide with elevated CH 4 concentration and active gas seepage. This suggests thataerobic CH 4 oxidation has a greater influence on the dissolved inorganic pool in areas where CH 4 concentrations are locally elevated, instead of displaying a cumulative effect downcurrent from widespreadgroupings of CH 4 seeps. Afirst-order approximation of the input rate of ancient-derived dissolved inorganiccarbon (DIC) into the waters overlying the northern U.S. Mid-Atlantic Bight further suggests thatoxidation of ancient CH 4 -derived carbon is not negligible on the global scale and could contribute todeepwater acidification over longer time scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Detectability of subsolar mass neutron stars through a template bank search

Here, we study the detectability of gravitational-wave signals from subsolar-mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as 78.4% for an equal mass binary system of chirp mass 0.17M ⊙ , in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for subsolar-mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of subsolar-mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for subsolar-mass binary neutron star signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A fusion relevant data-driven engineering void swelling model for 9Cr tempered martensitic steels

The UCSB database on cavity evolution in 9-12Cr tempered martensitic steels (TMS), includes the results for both dual heavy and helium ion (DII), and High Flux Isotope Reactor (HFIR) in situ helium injection (ISHI) neutron irradiations at 500°C. These results were combined with literature single ion and fission neutron irradiation data to derive a model for the void volume fraction, f v , as a function of displacements per atom (dpa) and transmutant helium concentrations in atomic parts per million (appm). The scientific foundation for the paper is described in a companion paper entitled “Cavity Evolution and Void Swelling in Dual Ion Irradiated Tempered Martensitic Steels”. Here, in this study, we show that f v (dpa, He/dpa) is described by the incubation dose, dpa i , for the onset of void growth, and the post-incubation growth rate, f v ’(%/dpa). Both dpa i and f v ’ decrease with increasing He/dpa at > ~ 5. The dpa i is also lower for the ISHI neutron irradiations at the same He/dpa. Single heavy ion and fission reactor neutron irradiations, with low He/dpa ratios, have a much larger dpa i . Based on a combined analysis of DII, single ion, ISHI and fission neutron data, we further show that the post-incubation f v data analyzed here have a common empirical curve shape, with f v ’ reaching up to ~ 0.2%/dpa at very high dpa. We also show that f v ’ can be predicted based on a physical model of defect partitioning between evolving sinks. At 500°C and fusion relevant He/dpa ≈ 10, the best-fit model predicts nominal swelling, S = f v /(1-f v ), of ~ 1.1, 4.9 and 16% at 50, 100 and 200 dpa, respectively. The physically motivated, data-driven model includes estimated uncertainties for both dpa i and f v ’.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A single-zone zero-dimensional study of HCCI combustion of methanol dehydration products to enable ignition of direct-injected methanol

Methanol is an alternative fuel gaining traction in the maritime sector. Its direct adoption, however, is accompanied by a unique set of technical challenges, such as low cetane number and high latent heat of vaporization. An approach to overcome these challenges is being developed at the US Department of Energy’s Oak Ridge National Laboratory, where onboard generation of dimethyl ether (DME) via catalytic dehydration of methanol can be used to assist in the mixing controlled combustion of direct-injected (DI) methanol. The generated mixture from this dehydration process can be premixed with intake air to condition the cylinder via. homogeneous charge compression ignition (HCCI) for subsequent DI methanol. In this preliminary work, various catalyst or reactor conversion efficiencies were simulated (using bottles) at constant DME and water flow at low load on a single-cylinder marine-variant of a CAT® C18 18 L engine with a 145 mm bore. To substantiate the experimental findings, a zero-dimensional engine model was developed in Cantera using a DME mechanism with 79 species and 658 reactions. Results presented include experimental and simulation heat release rate comparisons, species evolution information, and constant volume ignition delay (ID) for DI methanol with and without background species from HCCI of the premixed products from different reactor efficiencies. The results suggest that thermal effects dominate the DI methanol ignition process, and this work provides a chemical kinetic foundation or guideline for developing future control schemes.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗
DOE OSTI2021

Effects of Inlet Area Ratio on Operability of an Axial Air Inlet Rotating Detonation Combustor

This research presents the first results from a new Rotating Detonation Combustor (RDC) design with an axial air inlet and rear-facing angled fuel injection. This design features a variable air inlet that allows for the study of the effect of inlet area ratio on operation and to directly study the coupling between the air plenum and channel of the RDC while minimizing inlet pressure losses. This device maintains the same downstream channel configuration as those previously tested in our laboratory, which includes a 154 mm outer channel diameter and a constant area channel with a width of 7.6 mm. All tests were completed with a downstream hydrogen-air afterburner capable of flashing back in the channel to initiate detonation operation. Initial tests demonstrated successful operation at air mass flow rates from 160 to 770 g/s and equivalence ratios from 0.53 to 1.05 at multiple air inlet to channel area ratios, indicated that sensors installed within the air feed plenum are able to characterize the dynamic characteristics of waves within the detonation channel, and verified that the design goals for this configuration were achieved.

Shepard, Joshua↗
DOE OSTI2020

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗
DOE OSTI2022

Modelling of sawtooth-induced fast ion transport in positive and negative triangularity in TCV

Internal kinks are a common magneto hydro-dynamic (MHD) instability observed in tokamak operation when the q profile in the plasma core is close to unity. This MHD instability impacts both the transport of the bulk plasma (current, particle and energy transport) and minority species, such as fast ions. In TCV (R 0 /a = 0.88 m/0.25 m) the fast ion population is generated in the plasma by neutral beam tangential injection of energies up to 28 keV. TCV features 16 active shaping coils permitting a great flexibility in plasma shape, including negative triangularity (δ) configurations that show surprisingly high confinement. This study focuses on the transport of fast ions induced by sawteeth, by comparing two triangularity cases and simulation results with experimental data. Comparison of two equilibria with opposite δ shows that the fast ion drifts are larger for δ < 0. Furthermore, the sawtooth-induced transport in this case is larger than δ > 0 in similar conditions. Additionally, comparison with experimental data confirms the dominance of the modification of thermal kinetic profiles following the sawtooth crash in explaining drops in the neutron rates and FIDA signals. Additional fast ion diffusion, however, improves the interpretation of the experimental data. For δ < 0, the amplitude of the perturbation better representing the experimental data is larger. Finaly, an exploratory study for 50 keV particles (soon available in TCV) shows that the situation does not worsen for such particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Exploring Whether Subsurface Fluid Production Can Minimize Triggered Seismicity in Geothermal Fields

Fluid injection and production related to energy recovery and other industrial operations alter the pressure and stress state of subsurface reservoirs which can lead to induced seismicity. Here, the primary goal is to investigate the hypothesis that the modulation of stress state in subsurface reservoirs through fluid injection/production operations can reduce the likelihood of inducing seismicity. Validation of this mitigation strategy will provide an active operational method to control induced seismicity in subsurface reservoir exploitations such as geothermal energy recovery and carbon storage. In this project, we analyze the relationship between extensive fluid production and the paucity of aftershock activity at the Coso Geothermal plant (CGP) following the 2019 Ridgecrest earthquake sequence, where high rates of aftershock triggering were expected. We developed a high-fidelity multiphase coupled thermo-hydro-mechanical (THM) model to simulate fluid injection/production activities at the CGP between 1986 and 2020. THM results of surface subsidence due to high rates of fluid production agree well with field observations from global positioning system (GPS) and interferometric synthetic aperture radar (InSAR). Subsequently, the variations of pore pressure, temperature and stress state were used to drive numerical earthquake simulations of the aftershock response to the 2019 Ridgecrest earthquakes. The earthquake simulations show that seismic quiescence may occur following the Ridgecrest event, depending on the initial stress state at the time of the mainshock and at the initiation of geothermal energy production. Seismic quiescence occurs in 20% of the cases we explored, where rates are decreased by 50% of the background rate in the two years following the mainshock. In circumstances where aftershock rates increase near the CGP following Ridgecrest, the average rate change is a factor of two larger than background rates from simulations with no operations. These findings indicate that unlike many other operations that bring faults closer to failure, operations at the CGP are acting to stabilize faults such that triggering in the current stress state is minimal.

58 GEOSCIENCES↗
DOE OSTI2020

An expanded kinetics model describing injector region desorption

During the initial stages of electron beam injection, absorbed or entrained material can be released into the nascent beam pulse by thermal desorption or electron beam impact. These processes, and subsequent ionization and dissociation, are examined assuming a prototypical temporal beam pulse. An expanded set of kinetic equations that includes thermal desorption, electron-stimulated ion emission, electron impact dissociation, electron impact ionization, and thermal decomposition of adsorbed water molecules is developed. Estimates of the three most important ionized species H 2 O + , H + and OH + are provided for typical injector conditions. Electron impact dissociative ionization of neutral H 2 O is the dominant source of H + under these conditions. Since the abundance of thermally desorbed neutral H 2 O has an exponential rate, the production of H+ is likewise controlled by this exponential desorption rate. Thus thermal management of the confining walls could suppress proton production in the injector region. Thermal decomposition of H 2 O will not be a contributing physical process in the injector region on time scales less than a μs. The formalism is sufficiently flexible to include entrained hydrogen.

42 ENGINEERING↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Resource Characterization to Estimate Potential for Electricity Co-Production at Blackburn Oil Field, Nevada

The U.S. Department of Energy estimates that an annual average of 25 billion barrels of hot water are produced from oil and gas wells within the United States. The thermal energy available in the co-produced water stream is usually discarded, as the produced waters are considered an inconvenience by the operators and are disposed of using injection wells. However, utilizing organic Rankine cycle (ORC) generators, a vast amount of thermal energy can be captured and converted into electricity (albeit at relatively low efficiency due to the low temperatures). The National Renewable Energy Laboratory (NREL), in collaboration with Transitional Energy and Grant Canyon Oil & Gas, evaluated the feasibility of geothermal co-production of electricity by utilizing existing oil wells in Blackburn oil field in Nevada. The once prolific Blackburn oil field is located in Pine Valley, approximately 45 miles east-southeast of Elko, Nevada. Currently, the wells targeting the highly fractured Devonian Nevada dolomite reservoir are operating at a water cut ratio of more than 99%, with individual fluid (oil and water) production rates reaching 7.4 L/s (4,021 BBL/day). Analysis of publicly available data showed that the combination of the suitable wells' maximum historical production rates reached 22.90 L/s. The production from these wells occurs naturally and the wells are choked (and even shut down) by the operator to mitigate excessive water production, indicating a strong reservoir recharge and future opportunity to increase the water production for geothermal electricity generation. The main goal of this study was to evaluate the productivity of the existing wells, the performance of the reservoir, the surface network, and the operational constraints in order to achieve 1 MWe of electricity production from the field's water production. Utilizing the GEOPHIRES tool, we have determined that a twofold to threefold increase in the total fluid production, compared to the historical production under artificial restraint (choke), is required to reach a 1-MWe net target output for a low-temperature ORC system with air-cooled condensers. Lower flow rates would be required when utilizing water-based condensers instead of air-cooled condensers. However, that would require a constant supply of cold water, which may be challenging given the arid environment of the project site.

Blackburn↗
DOE OSTI2020

Hot Electron Driven Photocatalysis on Plasmon-Resonant Grating Nanostructures

Here, we demonstrate the hot electron injection of photoexcited carriers in an Ag-based plasmon resonant grating structure. By varying the incident angle of irradiation, sharp dips are observed in the reflectance with p-polarized light (electric field perpendicular to grating lines) when there is wavevector matching between the incident light and the plasmon resonant modes of the grating and no angle dependence is observed with s-polarized light. This configuration enables us to compare photoelectrochemical current produced by plasmon resonant excitation with that of bulk metal interband absorption simply by rotating the polarization of the incident light while keeping all other parameters of the measurement fixed. With 633 nm light, we observed a 12-fold enhancement in the photocurrent (i.e., reaction rate) between resonant and nonresonant polarizations at incident angles of ±7.6° from normal. At 785 nm irradiation, we observed similar resonant profiles to those obtained with 633 nm wavelength light but with a 44-fold enhancement factor. Using 532 nm light, we observed two resonant peaks (with approximately 10× enhancement) in the photocurrent at 19.4° and 28.0° incident angles, each corresponding to higher order modes in the grating with more nodes per period. The lower enhancement factors observed at shorter wavelengths are attributed to interband transitions, which provide a damping mechanism for the plasmon resonance. Finite difference time domain (FDTD) simulations of these grating structures confirm the resonant profiles observed in the angle-dependent spectra of these gratings and provide a detailed picture of the electric field profiles on and off resonance.

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

Pore-Scale Oil Connectivity and Displacement by Controlled-Ionic-Composition Waterflooding Using Synchrotron X-Ray Microtomography

Controlled-ionic-composition waterflooding is an economic and effective method to improve oil recovery in carbonate oil reservoirs. Recent studies show controlling the salinity and ionic composition of injection water can alter the wettability of carbonate mineral surfaces. The pore-scale oil connectivity and displacement by controlled-ionic-composition waterflooding in heterogeneous carbonate reservoirs, especially at the early stage, is still unclear. The goal of this study is to examine the role of ion concentrations and types in the oil displacement efficiency and investigate the impact of the waterflooding on the pore-scale oil displacement using the national synchrotron facility. A carbonate rock sample was flooded with synthetic high-salinity water and other water solutions with different sulfate concentrations. The waterflooding processes were visualized with synchrotron X-ray microtomography to follow the evolution of pore-scale oil/brine interactions at typical field flow rates. Experimental results show that the water with lower sulfate concentration and higher salinity did not change the wettability of the pore surfaces. Higher sulfate ion concentrations in the water, in contrast, altered the wettability of carbonate pore surfaces from oil-wet to neutral-wet within the first few minutes of waterflooding. Finally, novel insight was gained on the ability of water with high-sulfate concentration to displace oil in the small pores and through abundant oil channels, which could consequently lead to higher oil recovery from the carbonate rock.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Techno-Economic Performance of Eavor Loop 2.0

This project evaluated techno-economic performance for a sample Eavor-Loop 2.0 design for electricity production and direct-use heating. The Eavor-Loop 2.0 design investigated is a 7.5-km deep closed-loop geothermal system consisting of 12 laterals for a total of more than 90 km of downhole well and lateral length. Both a high geothermal gradient scenario of 60 degrees C/km and a low geothermal gradient scenario of 30 degrees C/km were considered. With pure water injected at 60 degrees C and 80 kg/s, reservoir simulations with the Slender-Body Theory simulator indicate average production temperatures over a 30-year lifetime of ~125 degrees C and ~210 degrees C for the low and high geothermal gradient scenario, respectively. These correspond to heat production of ~22 M Wth and ~51 M Wth, respectively. Using IPSEpro simulations, we find average power production of ~2.2 M We and ~8.6 M We, respectively, for a subcritical organic Rankine cycle power plant with air-cooled condensers. Cost estimates indicate the overall capital and levelized costs are dominated by the lateral drilling cost. Obtaining a levelized cost of electricity below $70/M Wh requires a geothermal gradient of 60 degrees C/km, a discount rate below 9%, and lateral drilling cost below $400/m. A well cost model indicates that ~$400/m for the Eavor-Loop 2.0 design investigated can be obtained for a drilling rate of penetration about 40 ft/hr (with bit life of 50 hours), and omitting casing and cement. Traditional (geothermal) well drilling has achieved these drilling rate conditions, including the Utah FORGE project where the rate of penetration has exceeded 50 ft/hr in granite. However, it is unclear if these conditions are still valid for drilling the Eavor-Loop 2.0 laterals (i.e., ~82 km of laterals at 4 to 7.5-km vertical depth with rock temperatures up to 460 degrees C), as such downhole completion has never been developed before. Competitive levelized cost of heat values ($1.2-$8.2/GJ) are calculated, even for the low geothermal gradient scenario (30 degrees C/km) and lateral drilling cost of $600/m.

advanced geothermal system↗
DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗
DOE OSTI2024

Produced water sharing: Improved economics and reduced community impact – A Pennsylvania case study

Here, hydraulic fracturing for oil and gas extraction from unconventional reservoirs is water intensive. Between water sourced for fracturing purposes and water present in rock formations, operators often produce a greater volume of water than oil or gas. Historically, this surplus of produced water has mainly been disposed of via deep injection wells. Rising disposal costs, seasonally limited water availability, and concerns over induced seismicity have incentivized produced water recycling practices, where an operator uses produced water for hydraulic fracturing operations. The logistical challenges associated with produced water recycling have also encouraged operators to adopt ad-hoc water exchange practices, in which competing operators will exchange produced water for mutual cost savings. In this paper, we investigate the potential benefit from systematic produced water exchange among operators in Northeastern Pennsylvania. We leverage PARETO, a free and open-source modeling framework for produced water management optimization, to quantify the benefits of water exchange practices. In an example drawn from FracFocus data, we find that the adoption of systematic water sharing could improve produced water recycling rates from 49.2% to 99%, decreasing operating and trucking costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Small Polarons and Surface Defects in Metal Oxide Photocatalysts Studied Using XUV Reflection–Absorption Spectroscopy

In photocatalytic transition metal oxides, the surface electronic structure and carrier kinetics are complicated by both charge trapping at defects and self-trapping by small polaron formation. These effects are specific to both the material and the type of defect or surface states involved. Here we review recent studies on ultrafast photoexcited charge carrier dynamics at the surfaces of alpha-Fe 2 O 3 (hematite) and NiO and at the interface of a alpha-Fe 2 O 3 /NiO heterojunction using extreme ultraviolet (XUV) reflection-absorption spectroscopy. We study dynamics of small polaron formation at hematite surfaces showing that the rate of carrier self-trapping is slower than in the bulk due to a greater lattice reorganization energy required for surface polaron formation. We also show that surface dynamics of hematite can be systematically tuned to control carrier transport properties through surface functionalization. To better understand the impact of defects on carrier events, we study charge carrier recombination of NiO. It is shown that O vacancies have no detrimental effect on carrier lifetime while grain boundaries cause fast recombination. Finally, we study a model heterojunction consisting of alpha-Fe 2 O 3 and NiO and find that interfacial charge transfer across these two materials occurs by a two-step mechanism where the interfacial electric field drives fast exciton dissociation followed by hole injection from alpha-Fe2O3 to NiO. Finally, these examples illustrate future opportunities to observe surface electron dynamics in photocatalytic materials with element and chemical state resolution using ultrafast XUV spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗