Search NASA⌕ Search

SEARCH · Search NASA

Results for “pressure diffusion”

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

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Resonant interactions involving local vibrational modes in crystals

When an impurity with a light mass is inserted into a crystal, it can undergo a high-frequency oscillation referred to as a local vibrational mode (LVM). A Fermi resonance may occur between the LVM and lower-frequency modes of the defect. The LVM may also interact with phonons or the electromagnetic field. Understanding these interactions can help model and control diffusion, defect reactions, and thermal conductivity. LVMs have been probed in semiconductors using pressure and alloying as experimental parameters, resulting in anticrossing between localized and extended vibrational modes. Furthermore, these types of vibrational interactions could play an important role in the stability and thermal properties of organic–inorganic hybrid semiconductors. The coupling between an LVM and electromagnetic wave yields an “LVM polariton,” an excitation that has significant vibrational and electric-field amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Dynamically Downscaled Hourly Future Weather Data with 12-km Resolution Covering Most of North America

This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.

Array↗
DOE OSTI2020

Facilitated Transport Membranes with Tunable Amine-CO 2 Chemistry for CO 2 /H 2 Separation

CO 2 -selective, amine-containing facilitated transport membranes are of great interest for syngas purification since high-pressure H 2 can be retained upon CO 2 removal. A variety of amine-containing polymers have shown decent chemical and thermal stability at aggressive conditions, but their CO 2 /H 2 separation properties are largely limited by the hydrostatic compaction and severe carrier saturation associated with the high syngas pressure. Herein, we report a new approach to enhance the CO 2 permeability by manipulating the steric hindrance of the amine carriers. A series of α-aminoacids with different alkyl or hydroxyethyl substituents are deprotonated by 2-(1-piperazinyl)ethylamine, leading to nonvolatile amine carriers with different degrees of steric hindrance. For hosting the low MW amine carriers, a highly crosslinked poly(vinyl acetal) is synthesized as a water-swellable polymer network. In order to avoid membrane compaction, perforated graphene oxide mono-sheets are dispersed as reinforcement fillers. In the presence of moisture, a bulkier alkyl substituent to the amino site (increasing steric hindrance) destabilizes the carbamate adduct and thus drastically increases the chemisorption of CO 2 , while the incorporation of ethylene oxide groups provides additional physisorption of CO 2 . The enhanced CO 2 solubility significantly mitigates the carrier saturation behavior, and an unprecedented CO 2 /H 2 selectivity greater than 100 is demonstrated at 107 °C and 12.5 atm of CO 2 partial pressure. As the CO 2 partial pressure reduces to 0.4 atm, a less hindered amine yields a higher reactive diffusivity of CO 2 , resulting in a CO 2 permeance of 435 GPU with a selectivity greater than 500. These reaction-mediated polymeric membranes are well above the theoretical upper bound, and they are of great interest for designing a highly-selective membrane process for syngas purification.

08 HYDROGEN↗
DOE OSTI2025

Microstructure-Based Degradation Modeling in Solid Oxide Cells

The presentation is intended for a graduate class in Georgia Southern University taught by Dr. Hayri Sezer on electrochemistry, SOFC and Li ion battery. It summaries NETL's recent works on microstructure-based modeling of degradation in solid oxide cells. It covers the degradation mechanism of Ni coarsening, Ni migration and the build-up of oxygen partial pressure inside the electrolyte. The phase-field models developed for Ni coarsening and migration are reviewed and the results are compared to experiments in literature to examine different diffusion mechanisms (Ni self-diffusion, Ni(OH)2 diffusion in pores and Ni(OH) surface diffusion) and driving forces (Ni(OH)x concentration and Ni-YSZ contact angle change). It also demonstrates the developed model for oxygen partial pressure and shows that the predicted oxygen partial pressure can explain the Ni oxidation in the hydrogen electrode and crack growth on the electrolyte-oxygen electrode interface under certain operating conditions.

nickel coarsening↗
DOE OSTI2021

Experimental study of sprays with isooctane, hexane, ethanol and their binary mixtures under different flash boiling intensities

Here, this work investigated the behaviors of single-component and multi-component sprays under different flash boiling intensities. Isooctane, hexane, ethanol, isooctane-hexane, and isooctane-ethanol mixtures were used as test fuels. The pressure ratio of the ambient pressure over saturation pressure (R p ) was used as the indicator of the superheated degree. Macroscopic characteristics of the sprays under 10 different R p (1.1 to 0.14) were investigated using Diffused Backlight Imaging (DBI). Droplet size, which revealed the spray microscopic characteristics, was also investigated using Phase Doppler Anemometry (PDA). The results showed that ambient pressure reduction led to higher flash boiling intensity, significant spray morphology changes, and smaller droplet size, regardless of the fuel type and blend ratios. Spray collapse was observed for all sprays once R p dropped below 0.3. Spray collapse dramatically increased the liquid penetration length, expanded the near-field spray angle, and reduced the far-field angle for sprays with all the fuels. Fuel properties, such as the latent heat of vaporization, had significant impacts on droplet size distributions, thus altering the liquid penetration length and spray angles. Under similar R p , sprays with high ethanol content generally had larger droplet sizes compared with other fuel sprays due to the high latent heat of vaporization. The larger droplets caused the high ethanol content sprays to be more resistant to the plume trajectory change, resulting in a smaller far-field spray angle variations with changing R p . At R p of 0.28, sprays with the similar latent heat of vaporization have similar droplet sizes, and droplet sizes at different radial locations were also similar to each other.

42 ENGINEERING↗
DOE OSTI2022

High-temperature and/or high pressure gas enclosure

A gas enclosure can include a refractory metal liner; a ceramic matrix composite cladding; and a diffusion barrier layer. The refractory metal liner is adapted to surround and enclose a gas to be contained within the gas enclosure. The diffusion barrier layer is disposed between the refractory metal liner and the ceramic matrix composite cladding.

Miao, Yinbin↗
DOE OSTI2021

Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors

Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.

Braun, James↗
DOE OSTI2022

Synthesis of wafer-scale SWCNT forests with remarkably invariant structural properties in a bulk-diffusion-controlled kinetic regime

Robust synthesis of vertically aligned carbon nanotubes (VACNT) at large scale is required to accelerate deployment of numerous cutting-edge devices to emerging commercial applications. To address this need, this work demonstrates that the structural characteristics of single-walled CNTs (SWCNTs) produced in a growth regime dominated by bulk diffusion of the gaseous carbon precursor are remarkably invariant over a broad range of process conditions (precursor concentration increased up to 30-fold, catalyst substrate area from 1 cm 2 to 180 cm 2 , growth pressure from 20 to 790 mbar, and gas flowrates up to 8-fold). A simple growth kinetics model is derived that accounts for both reactant bulk diffusion and competing byproduct formation in the presence of excess hydrogen, which quantitatively reproduces the experimental data in the entire isothermal parameter space. The model enables critical predictions for process scale-up optimization including a potential 6-fold increase in CNT production rate, a ~90% carbon conversion efficiency in select conditions, and elimination of reaction rate decay observed at high pressures by using a hydrogen-free growth environment. This model-guided opportunity to substantially improve synthesis throughput and efficiency, in conjunction with the remarkably invariant CNT properties, is appealing for reliable VACNT device fabrication at both small and industrial scales.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2022

Oxygen diffusion in garnet: Experimental calibration and implications for timescales of metamorphic processes and retention of primary O isotopic signatures

Knowledge of oxygen diffusion in garnet is crucial for a correct interpretation of oxygen isotope signatures in natural samples. A series of experiments was undertaken to determine the diffusivity of oxygen in garnet, which remains poorly constrained. The first suite included high-pressure (HP), nominally dry experiments performed in piston-cylinder apparatus at: (1) T = 1050–1600 °C and P = 1.5 GPa and (2) T = 1500 °C and P = 2.5 GPa using yttrium aluminum garnet (YAG; Y3Al5O12) cubes. Second, HP H2O-saturated experiments were conducted at T = 900 °C and P = 1.0–1.5 GPa, wherein YAG crystals were packed into a YAG + Corundum powder, along with 18O-enriched H2O. Third, 1 atm experiments with YAG cubes were performed in a gas-mixing furnace at T = 1500–1600 °C under Ar flux. Finally, an experiment at T = 900 °C and P = 1.0 GPa was done using a pyrope cube embedded into pyrope powder and 18O-enriched H2O. Experiments using grossular were not successful. Profiles of 18O/(18O+16O) in the experimental charges were analyzed with three different secondary ion mass spectrometers (SIMS): sensitive high-resolution ion microprobe (SHRIMP II and SI), CAMECA IMS-1280, and NanoSIMS. Considering only the measured length of 18O diffusion profiles, similar results were obtained for YAG and pyrope annealed at 900 °C, suggesting limited effects of chemical composition on oxygen diffusivity. However, in both garnet types, several profiles deviate from the error function geometry, suggesting that the behavior of O in garnet cannot be fully described as simple concentration-independent diffusion, certainly in YAG and likely in natural pyrope as well. The experimental results are better described by invoking O diffusion via two distinct pathways with an inter-site reaction allowing O to move between these pathways. Modeling this process yields two diffusion coefficients (D values) for O, one of which is approximately two orders of magnitude higher than the other. Taken together, Arrhenius relationships are:logDm2s-1=-7.2(±1.3)+(-321(±32)kJmol-12.303RT) for the slow pathway, andlogDm2s-1=-5.4(±0.7)+(-321(±20)kJmol-12.303RT) for the fast pathway. We interpret the two pathways as representing diffusion following vacancy and inter-stitial mechanisms, respectively. Regardless, our new data suggest that the slow mechanism is prevalent in garnet with natural compositions, and thus is likely to control the retentivity of oxygen isotopic signatures in natural samples. The diffusivity of oxygen is similar to Fe-Mn diffusivity in garnet at 1000–1100 °C and Ca diffusivity at 850 °C. However, the activation energy for O diffusion is larger, leading to lower diffusivities at P-T conditions characterizing crustal metamorphism. Therefore, original O isotopic signatures can be retained in garnets showing major element zoning partially re-equilibrated by diffusion, with the uncertainty caveat of extrapolating the experimental data to lower temperature conditions.

Geochemistry & Geophysics↗
DOE OSTI2023

Non-equilibrium molecular dynamics studies of thermal diffusion of hydrogen isotopes in low concentration zirconium hydrides

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Keyhole fluctuation and pore formation mechanisms during laser powder bed fusion additive manufacturing

Keyhole porosity is a key concern in laser powder-bed fusion (LPBF), potentially impacting component fatigue life. However, some keyhole porosity formation mechanisms, e.g., keyhole fluctuation, collapse and bubble growth and shrinkage, remain unclear. Using synchrotron X-ray imaging we reveal keyhole and bubble behaviour, quantifying their formation dynamics. The findings support the hypotheses that: (i) keyhole porosity can initiate not only in unstable, but also in the transition keyhole regimes created by high laser power-velocity conditions, causing fast radial keyhole fluctuations (2.5-10 kHz); (ii) transition regime collapse tends to occur part way up the rear-wall; and (iii) immediately after keyhole collapse, bubbles undergo rapid growth due to pressure equilibration, then shrink due to metal-vapour condensation. Concurrent with condensation, hydrogen diffusion into the bubble slows the shrinkage and stabilises the bubble size. The keyhole fluctuation and bubble evolution mechanisms revealed here may guide the development of control systems for minimising porosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

High frequency impedance characteristics of a tunable microplasma device

Computational studies on high frequency impedance characteristics of a microplasma device are reported. While microplasma is ignited using a primary excitation signal, frequency response of plasma impedance is determined by a secondary high frequency probe signal with significantly lower voltage amplitude such that it does not influence the plasma parameters. Here, the computational model utilizing the drift–diffusion approximation is first validated by comparing with experimental data for microplasmas ignited at pressures ranging from 1 to 5 Torr. In spite of quantitative discrepancies, good overall agreement is obtained between the measured frequency response of impedance of the discharge. Comparisons are also presented for various plasma parameters including mean electron number density, sheath thickness, mean electron temperature, and collision frequency that were inferred from the impedance measurements. The computational model is then used to perform simulations of near-atmospheric pressure microplasmas with the probe signal frequency ranging from 3 to 20 GHz. The simulations demonstrate the presence of a resonance frequency at which the impedance vanishes. More importantly, it is shown that this resonant frequency can be tuned effectively by suitably modifying the operating parameters (gap size, pressure, and excitation voltage). The simulated impedance characteristics are used to determine the effective plasma inductance and capacitance using a non-linear fitting approach, thereby showing the dependence of these electrical parameters on the plasma operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phase stability of a eutectic high entropy alloy under extremes of pressures and temperatures

Additively manufactured high-entropy alloys are of interest because of their unique combination of high yield strength and large ductility achieved with far-from-equilibrium crystalline phases and micro/nanostructure morphology. We report on the phase transformation and thermal equation of state of the eutectic high-entropy alloy (EHEA) Al 18 Co 20 Cr 10 Fe 10 Ni 40 W 2 , produced by laser powder-bed fusion (L-PBF). The EHEA was studied in a large-volume Paris–Edinburgh cell using energy-dispersive x-ray diffraction to a pressure of 5.5 GPa and a temperature of 1723 K. Static compression studies in diamond anvil cells using angle-dispersive x-ray diffraction extended the high-pressure structural data to 317 GPa at ambient temperature. The initial dual-phase nanolamellar face-centered cubic (FCC) and body-centered cubic (BCC) structure of Al 18 Co 20 Cr 10 Fe 10 Ni 40 W 2 transforms into a single FCC phase under high pressure, with the BCC-to-FCC phase transformation completing at 9 ± 2 GPa. The FCC phase remained stable up to the highest pressure of 317 GPa. The measured thermal equation of state for the FCC phase of Al 18 Co 20 Cr 10 Fe 10 Ni 40 W 2 is presented up to 5.5 GPa and 1473 K. We observed melting of the EHEA at 1698 ± 25 K at a pressure of 5.5 GPa, and the recrystallized sample shows an increased fraction of the CsCl-type (B2) phase at ambient conditions following release from the high-pressure high-temperature state. The BCC-to-FCC phase transition completion pressure is correlated with the nanolamellae thickness of the BCC layer in this diffusion-less transformation at ambient temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

NMR relaxometry characterization of water adsorption in corn stover anatomical fractions

Nuclear magnetic resonance (NMR) relaxometry is applied to provide direct measurement of water adsorption in anatomical fractions of corn stover. NMR transverse T 2 relaxation time distribution measurements indicate multiple water populations, which vary with anatomical fraction and water adsorption. Measured T 2 data are used to calculate thermodynamic properties of Brunauer-Emmet-Teller adsorption theory using a model to estimate mono and bilayer relaxation. T 2 data are used directly to determine rotational diffusion correlation times indicating adsorption interaction strength. T 1 -T 2 longitudinal-transverse relaxation time correlation measurements quantify differences in the molecular level structural order of the adsorbate surface water as a function of water activity, i.e. relative humidity or water vapor partial pressure. The T 1 /T 2 ratio provides a measure of the surface energy related to the adsorption strength and surface diffusive mobility of the water adsorbate and differentiates the anatomical fractions. Furthermore, the results indicate that direct measurement of NMR relaxation times can be used to characterize corn stover biomass water adsorption, data relevant to processing and handling considerations.

09 BIOMASS FUELS↗