Search NASA⌕ Search

SEARCH · Search NASA

Results for “zone of influence”

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 163 records · Page 9

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

Effect of the Quasi-Biennial Oscillation on the Madden Julian Oscillation Teleconnections in the Southern Hemisphere

The Madden Julian Oscillation (MJO) is the main source of intraseasonal variability in the tropics. MJO modulates rainfall in remote areas in the Southern Hemisphere by exciting tropical and extratropical wave trains. We use newly released reanalysis data to analyze how the Quasi-Biennial Oscillation (QBO) can influence the MJO's effects over the Southern Hemisphere, focusing on precipitation anomalies over South America. The anomalies in the intensity of the South Atlantic Convergence Zone (SACZ) and precipitation over Southeastern South America during MJO phases 1 and 4 are intensified during the easterly QBO. The extratropical wave train excited by the anomalous convection over the maritime continent during MJO phase 4 is affected by the QBO, with a stronger, better-defined pattern during westerly QBO. The conclusions are supported by a perturbation experiment of a case study, conducted using a high-top atmospheric global climate model where the QBO and MJO are controlled.

58 GEOSCIENCES↗

PARTICLE RESOLVED DNS STUDY OF TURBULENCE EFFECTS ON HYPORHEIC MIXING IN RANDOMLY PACKED SEDIMENT BEDS

Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.

Karra, Shashank↗

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗

Dynamic coevolution of baseflow and multiscale groundwater flow system during prolonged droughts

Field and numerical studies suggest that baseflow is composed of waters from a spectrum of groundwater flow paths termed the Groundwater Flow System (GWFS) – from shallow hillslope contributions to watershed-scale deep circulation originating in headwaters and discharging into lowland rivers. Here, we explore the evolution of the GWFS under prolonged droughts to understand its dynamics and multiscale nature, and to elucidate its role in baseflow generation and recession at the watershed scale. In this work, we consider three drought scenarios of varying severity and simulate groundwater flow in a 2-D cross-section of an idealized watershed with deep permeable bedrock, tracking the evolution of flow paths, baseflow, and residence times during the recession process. We find that baseflow generation at different drainage stages, and within different subwatersheds, is influenced distinctly by flow paths of different scales, depending on the relative strength of the flow paths and the position of the subwatersheds relative to the recharge/discharge zones of the deeper watershed-scale groundwater circulation. Despite having the same local relief, geology, and climate, baseflow from each subwatershed has a distinct recession behavior and time-dependent residence time distribution. Also, the hydraulic and transport characteristics of baseflow generation co-evolve and are strongly affected by the connection state of the water table to subwatersheds. These findings suggest that asynchrony and dissimilarity of baseflow generation from hillslopes under the impact of the watershed-scale groundwater flow, and interactions with local-scale and intermediate-scale groundwater flow, must be taken into account when interpreting baseflow recession data and building conceptual baseflow models at the watershed scale.

54 ENVIRONMENTAL SCIENCES↗

Low–severity fire as a mechanism of organic matter protection in global peatlands: Thermal alteration slows decomposition

Worldwide, regularly-recurring wildfires shape many peatland ecosystems to the extent that fire-adapted species often dominate plant communities, suggesting that wildfire is an integral part of peatland ecology rather than an anomaly. The most destructive blazes are smoldering fires that are usually initiated in periods of drought and can combust entire peatland carbon stores. However, peatland wildfires more typically occur as low-severity surface-burns that arise in the dormant season when vegetation is desiccated, and soil moisture is high. In such low-severity fires, surface layers experience flash heating, but there is little loss of underlying peat to combustion. This study examines the potential importance of such processes in several peatlands that span a gradient from hemiboreal to tropical ecozones and experience a wide range of fire return intervals. We show that low-severity fires can increase the pool of stable soil carbon by thermally altering the chemistry of soil organic matter (SOM), thereby reducing rates of microbial respiration. Using XPS and FTIR, we demonstrate that low-severity fires significantly increase the degree of carbon-condensation and aromatization of SOM functional groups, particularly on the surface of peat aggregates. Laboratory incubations show lower CO 2 emissions from peat subjected to low-severity fire, and predict lower cumulative CO 2 emissions from burned peat after one to three years. Also, low-severity fires reduce the temperature sensitivity (Q 10 ) of peat, indicating that these fires can inhibit microbial access to SOM. The increased stability of thermally-altered SOM may allow a greater proportion of organic matter to survive vertical migration into saturated and anaerobic zones of peatlands where environmental conditions physiochemically protect carbon stores from decomposition for thousands of years. Furthermore, across latitudes, low-severity fire is an overlooked factor influencing carbon cycling in peatlands, which is relevant to global carbon budgets as climate-change alters fire regimes worldwide.

54 ENVIRONMENTAL SCIENCES↗

Interpretation of HFTS 2 Microseismic Data Using Bedding-Plane-Slip Mechanism

Summary The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (low-frequency distributed acoustic sensing, or LF-DAS, strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data have three unique observations that have not been commonly observed in other datasets: (1) Due to the influence of offset pad depletion, microseismic data show that hydraulic fractures from the child well can propagate over 3,000 ft into the depleted low-stress zone. (2) By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that the microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages is much shorter than what is shown with LF-DAS strain data. (3) Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. This model successfully explains the above field observations from the HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

Engineering↗

Process Intensification by One-Step, Plasma-Assisted Catalytic Synthesis of Liquid Chemicals from Light Hydrocarbons

In this project, we present a plasma-assisted, catalytic process that can handle variations in production rates and gas compositions of producing wells and reliably converts the hydrocarbons to gas phase olefins, high molecular weight alkanes, and liquid chemicals. The low-temperature plasma serves as a reactive chemical environment that activates and converts the light hydrocarbons to these valuable products. We present results from the integration of a catalyst into the low-temperature plasma zone to improve reaction efficiency and product selectivity. We show that the gas composition, bulk gas temperature, and input power of the plasma, with and without catalysts, influence the production rates of liquids from natural gas feeds.

03 NATURAL GAS↗

Characterizing the Non-linear Interactions Between Tide, Storm Surge, and River Flow in the Delaware Bay Estuary, United States

Low-lying coastal areas in the mid-Atlantic region are prone to compound flooding resulting from the co-occurrence of river floods and coastal storm surges. To better understand the contribution of non-linear tide-surge-river interactions to compound flooding, the unstructured-grid Finite Volume Community Ocean Model was applied to simulate coastal storm surge and flooding in the Delaware Bay Estuary in the United States. The model was validated with tide gauge data in the estuary for selected hurricane events. Non-linear interactions between tide-surge-river were investigated using a non-stationary tidal analysis method, which decomposes the interactions’ components at the frequency domain. Model results indicated that tide-river interactions damped semidiurnal tides, while the tide-surge interactions mainly influenced diurnal tides. Tide-river interactions suppressed the water level upstream while tide-surge interaction increased the water level downstream, which resulted in a transition zone of damping and enhancing effects where the tide-surge-river interaction was prominent. Evident compound flooding was observed as a result of non-linear tide-surge-river interactions. Furthermore, sensitivity analysis was carried out to evaluate the effect of river flooding on the non-linear interactions. The transition zone of damping and enhancing effects shifted downstream as the river flow rate increased.

54 ENVIRONMENTAL SCIENCES↗

Rock weathering controls the potential for soil carbon storage at a continental scale

Abstract As rock-derived primary minerals weather to form soil, they create reactive, poorly crystalline minerals that bind and store organic carbon. By implication, the abundance of primary minerals in soil might influence the abundance of poorly crystalline minerals, and hence soil organic carbon storage. However, the link between primary mineral weathering, poorly crystalline minerals, and soil carbon has not been fully tested, particularly at large spatial scales. To close this knowledge gap, we designed a model that links primary mineral weathering rates to the geographic distribution of poorly crystalline minerals across the USA, and then used this model to evaluate the effect of rock weathering on soil organic carbon. We found that poorly crystalline minerals are most abundant and most strongly correlated with organic carbon in geographically limited zones that sustain enhanced weathering rates, where humid climate and abundant primary minerals co-occur. This finding confirms that rock weathering alters soil mineralogy to enhance soil organic carbon storage at continental scales, but also indicates that the influence of active weathering on soil carbon storage is limited by low weathering rates across vast areas.

54 ENVIRONMENTAL SCIENCES↗

Economic analysis of CCUS: Accelerated development for CO 2 EOR and storage in residual oil zones under the context of 45Q tax credit

Residual oil zones (ROZ) undergoing CO 2 Enhanced Oil Recovery (CO 2 -EOR) may benefit from specific strategies to maximize their value. We evaluated several strategies for producing from a Permian Basin, West Texas, USA field’s ROZ. This ROZ lies below the main pay zone (MPZ) of the field. Such brownfield ROZs occur in the Permian Basin and elsewhere. Since brownfield ROZs are hydraulically connected to the MPZs, development sequences and schemes influence oil production, CO 2 storage, and net present value (NPV). We conducted economic assessments of various CO 2 injection/production schemes in the stacked ROZ-MPZ reservoir based on flow simulations of a high-resolution geocellular model built from wireline logs and core data and calibrated through production history matching. Flow simulations of water alternating gas (WAG) injection, such as switching injection from the MPZ to the ROZ and commingled production, were studied. Simulation results showed that simultaneous CO 2 injection into the MPZ and ROZ lead to the largest oil production and, generally, the largest NPV. If instead, CO 2 was simultaneously injected into the MPZ and ROZ, then into the ROZ alone, this maximized CO 2 storage. CO 2 storage can be used as a tax credit under the Internal Revenue Code, Section 45Q. Storage performance depends on the development approach and WAG ratio. Developing the ROZ increased storage compared to only producing from the MPZ. The WAG ratio to maximize oil production did not always yield the largest NPV. These findings are potentially applied to other Brownfield ROZs, which are common below San Andres reservoirs in the Permian Basin and other basins. ROZ development can increase oilfields’ NPV and carbon storage potential. Our study can serve as an analog for similar reservoirs. Here this work provides valuable insights into the further optimization of brownfield ROZ development and information for operators to plan to develop stacked ROZ-MPZ reservoirs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

CO 2 rock physics modeling for reliable monitoring of geologic carbon storage

Monitoring, verification, and accounting (MVA) are crucial to ensure safe and long-term geologic carbon storage. Seismic monitoring is a key MVA technique that utilizes seismic data to infer elastic properties of CO 2 -saturated rocks. Reliable accounting of CO 2 in subsurface storage reservoirs and potential leakage zones requires an accurate rock physics model. However, the widely used CO 2 rock physics model based on the conventional Biot-Gassmann equation can substantially underestimate the influence of CO 2 saturation on seismic waves, leading to inaccurate accounting. We develop an accurate CO 2 rock physics model by accounting for both effects of the stress dependence of seismic velocities in porous rocks and CO 2 weakening on the rock framework. We validate our CO 2 rock physics model using the Kimberlina-1.2 model (a previously proposed geologic carbon storage site in California) and create time-lapse elastic property models with our new rock physics method. We compare the results with those obtained using the conventional Biot-Gassmann equation. Our innovative approach produces larger changes in elastic properties than the Biot-Gassmann results. Using our CO 2 rock physics model can replicate shear-wave speed reductions observed in the laboratory. Our rock physics model enhances the accuracy of time-lapse elastic-wave modeling and enables reliable CO 2 accounting using seismic monitoring.

58 GEOSCIENCES↗

Particulate organic matter (POM) transport and transformation at the terrestrial-aquatic interface (Final Report)

This project investigates the input, transport, and degradation of particulate organic matter (POM) in near-surface riverbed sediments at the Hanford 300 Area of the Columbia River, a dynamic, regulated river system influenced by upstream dam operations. Riverbed sediments are biogeochemical hot spots where organic carbon inputs stimulate intense microbial activity, affecting nutrient cycling and redox transformations in the hyporheic zone (HZ). While dissolved organic matter (DOM) cycling has been studied extensively, little is known about the infiltration and transformation of POM—particularly under variable flow regimes common in large, regulated rivers.

54 ENVIRONMENTAL SCIENCES↗

Role of SaPCR2 in Zn Uptake in the Root Elongation Zone of the Zn/Cd Hyperaccumulator Sedum alfredii

Zn pollution is a potential toxicant for agriculture and the environment. Sedum alfredii is a Zn/Cd hyperaccumulator found in China and has been proven as a useful resource for the phytoremediation of Zn-contaminated sites. However, the molecular mechanism of Zn uptake in S. alfredii is limited. In this study, the function of SaPCR2 on Zn uptake in S. alfredii was identified by gene expression analysis, yeast function assays, Zn accumulation and root morphology analysis in transgenic lines to further elucidate the mechanisms of uptake and translocation of Zn in S. alfredii. The results showed that SaPCR2 was highly expressed in the root elongation zone of the hyperaccumulating ecotype (HE) S. alfredii, and high Zn exposure downregulated the expression of SaPCR2 in the HE S. alfredii root. The heterologous expression of SaPCR2 in yeast suggested that SaPCR2 was responsible for Zn influx. The overexpression of SaPCR2 in the non-hyperaccumulating ecotype (NHE) S. alfredii significantly increased the root uptake of Zn, but did not influence Mn, Cu or Fe. SR-μ-XRF technology showed that more Zn was distributed in the vascular buddle tissues, as well as in the cortex and epidermis in the transgenic lines. Root morphology was also altered after SaPCR2 overexpression, and a severe inhibition was observed. In the transgenic lines, the meristematic and elongation zones of the root were lower compared to the WT, and Zn accumulation in meristem cells was also reduced. These results indicate that SaPCR2 is responsible for Zn uptake, and mainly functions in the root elongation zone. This research on SaPCR2 could provide a theoretical basis for the use of genetic engineering technology in the modification of crops for their safe production and biological enhancement.

59 BASIC BIOLOGICAL SCIENCES↗

Interfacial microstructure and mechanical properties of rotary inertia friction welded dissimilar 422 martensitic stainless steel to 4140 low alloy steel joints

In this work, dissimilar rotary inertia friction welds between 422 martensitic stainless steel and 4140 martensitic low-alloy steel were made to fabricate prototype heavy-duty diesel engine pistons. The influence of the inertia friction welding process and post weld heat treatment (PWHT) temperature on the interfacial microstructure evolutions and corresponding effects on mechanical properties of the 422/4140 welds were evaluated in detail. Carbon diffused from the 4140 side to the 422 side during PWHT at 650 °C for 1.5 h, causing the formation of a hard carbide-rich layer on the 422 side, and a softer but discontinuous C-depleted layer the 4140 side. PWHT at 700 °C for 1.5 h greatly accelerated C diffusion across the interface relative to 650 °C, resulting in a thicker hard carbide-rich layer and a relatively thick and continuous layer of coarse C-depleted grains (ferrite) on the 4140 side. In addition, the PWHT temperature greatly influenced the tensile properties and fracture behavior of the welds, with the 650 °C PWHT-ed samples failing predominately in a ductile manner in the 4140 heat affected zone during tensile testing. Conversely, the 700 °C PWHT specimens exhibited a strength reduction compared with the 650 °C PWHT specimens because of additional coarsening of the interfacial ferrite layer and softening of the base materials during PWHT, with brittle fracture between the hard and soft layers the predominate failure mechanism. Based on the findings, a reduced PWHT temperature and/or time, minimizing the hardness differential of the base metals, and pre-heating the 422 steel prior to welding are the potential pathways to achieve a more optimal balance between desirable tempering and stress relief of the weld microstructure and undesirable C migration across the weld interface, and to reduce the strength mismatch across the weld.

36 MATERIALS SCIENCE↗

A comprehensive analysis of the crop effect on the urban-rural differences in land surface phenology

The response of land surface phenology (LSP) to the urban heat island effect (UHI) is a useful biological indicator for understanding how vegetated ecosystems will be affected by future climate warming. However, vegetation cover in rural areas is often dominated by cultivated land, whose phenological timing is considerably influenced by agricultural managements (e.g., timing of sowing and harvesting), leading to biased conclusions derived from the urban-rural LSP differences. To demonstrate this problem, we investigated the crop influence on the phenological response to a warmer environment resulting from the UHI effect. We partitioned cities in the United States into cultivated and non-cultivated categories according to the proportion of crops in rural areas. We then built continuous buffer zones starting from the urban boundary to explore the urban-rural LSP differences considering the UHI effect on them. In this work, the results suggest crop inclusion is likely to lead to >14 days of urban-rural differences at both the start of the season (SOS) and the end of the season (EOS) between cultivated and non-cultivated cities. The temperature sensitivity (ST) of SOS is overestimated by approximately 2.7 days/°C, whereas the EOS is underestimated by 3.6 days/°C. Removing crop-dominated pixels (i.e., above 50 %) can minimize the influence of crop planting/harvesting on LSP and derive reliable results. We, therefore, suggest explicit consideration of crop impacts in future studies of phenological differences between urban and rural areas and the UHI effect on LSP in urban domains, as presented by this comprehensive study.

54 ENVIRONMENTAL SCIENCES↗

Solution-Phase Growth of Cu Nanowires with Aspect Ratios Greater Than 1000: Multiscale Theory

Penta-twinned metal nanowires are finding widespread application in existing and emerging technologies. However, little is known about their growth mechanisms. Here, we probe the origins of chloride- and alkylamine-mediated, solution-phase growth of penta-twinned Cu nanowires from first principles using multi-scale theory. Using quantum density functional theory (DFT) calculations, we characterize the binding and surface diffusion of Cu atoms on chlorine-covered Cu(100) and Cu(111) surfaces. We find stronger binding and slower diffusion of Cu atoms on chlorinated Cu(111) than on chlorinated Cu(100) – a reversal of the trend for bare Cu surfaces. We also probe inter-facet diffusion and find this proceeds faster from Cu(100) to Cu(111) than the reverse. Using the DFT rates for hopping between individual sites at Ångstrom scales, we calculate coarse-grained, inter-facet rates for nanowires of various lengths – up to hundreds of μm – and diameters in the 10-nm range. We predict nanowires with aspect ratios around 100 based on surface diffusion, alone. We also account for the influence of a self-assembled alkylamine layer that covers most of the {100} facets, but is absent or thin and disordered on the {111} facets and in an “end zone” near the {100}-{111} boundary. With an end zone, we predict a wide range of nanowire aspect ratios in the experimental ranges. Our work reveals the mechanisms by which a halide – chloride – promotes the growth of high-aspect-ratio nanowires.

36 MATERIALS SCIENCE↗

Investigating the Role of Accident Tolerant Cladding on Source Term Reduction for High-Burnup PWRs Using MELCOR

The use of accident tolerant fuel (ATF) cladding can increase coping times during and beyond design basis accidents. While such gains may be incremental, they provide a margin that can potentially be recovered to enable high-burnup (HBU) operation. Realizing such a margin requires demonstrating that the combination of HBU and ATF has not led to an overall increase in source term. This study investigates the influence of cladding technology (Zr-based, Cr-coated Zr, and FeCrAl) and fuel cycle length (18 and 24 months) on radiological dose at the boundary of the exclusion zone for a four-loop pressurized water reactor to investigate whether ATF claddings can provide such benefits. We analyze a recovered large break loss-of-coolant accident scenario to investigate the impact of transient timescale on the benefits of such coping time increases. The simulations have been performed using the MELCOR and MELCOR Accident Consequence Code System codes. For the cases analyzed, increased fuel cycle length did not necessarily increase radionuclide release and hydrogen generation, as these were found to be sensitive to the core power distribution. Similarly, off-site dose consequence is dominated by short-lived radionuclides that tend to saturate earlier in the burnup, so higher burnup operation did not necessarily increase the source term for the phenomena and transients analyzed here. Delays in recovery of the lowpressure safety injection system increase hydrogen production and radionuclide release, especially between 780 s and 1620 s, due to the nonlinear oxidation and core degradation behavior. Results show that Cr-coated Zr enhances safety by delaying heatup and gap release. Here, when uncertainty propagation on oxidation properties is considered, FeCrAl exhibits the lowest overall radionuclide release and off-site dose throughout the spectrum. However, while the considered “base model” performance is superior under delayed injection scenarios, upper-bound cases display hydrogen generation risk comparable to the Zr-based cladding.

Accident Tolerant Fuel↗