Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing thermodynamics”

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 37 records · Page 2

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗

Enhanced Energy Transfer to an Optomechanical Piston from Indistinguishable Photons

Thought experiments involving gases and pistons, such as Maxwell’s demon and Gibbs’ mixing, are central to our understanding of thermodynamics. In this paper, we present a quantum thermodynamic thought experiment in which the energy transfer from two photonic gases to a piston membrane grows quadratically with the number of photons for indistinguishable gases, while it grows linearly for distinguishable gases. This signature of bosonic bunching may be observed in optomechanical experiments, highlighting the potential of these systems for the realization of thermodynamic thought experiments in the quantum realm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermodynamic Band Gap Model for Photoinduced Phase Segregation in Mixed-Halide Perovskites

Provided is a comprehensive description of a band gap thermodynamic model, which predicts and explains key features of photosegregation in lead-based, mixed-halide perovskites. The model gives a prescription for illustrating halide migration driven by photocarrier energies. Where possible, model predictions are compared with experimental results. Free energy derivations are provided for three assumptions: (1) halide mixing in the dark, (2) a fixed number of photogenerated carriers funneling to and localizing in low band gap inclusions of the alloy, and (3) the statistical occupancy of said inclusions from a bath of thermalized photocarriers in the parent material. Model predictions include excitation intensity (I exc )-dependent terminal halide stoichiometries (x terminal ), temperature-independent x terminal , excitation intensity thresholds (I exc,threshold ) below which photosegregation is suppressed, and reduced segregation in nanocrystals as compared to thin films. Here, the model also offers insight into kinetically manipulating photosegregation rates via control of underlying mediators and rationalizes asymmetries in forward and reverse photosegregation rate constants/activation energies. What emerges is a cohesive framework for understanding ubiquitous photosegregation in mixed-halide perovskites and a rational basis by which to manage the phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implications of reduced-complexity aerosol thermodynamics on organic aerosol mass concentration and composition over North America

Atmospheric organic aerosol (OA) mass concentrations can be affected by water uptake through its impact on the gas–particle partitioning of semi-volatile compounds. Current chemical transport models (CTMs) neglect this process. We have implemented the Binary Activity Thermodynamics model coupled to a volatility basis set partitioning scheme in the GEOS-Chem CTM, providing an efficient reduced-complexity OA model that predicts relative-humidity-dependent mixing and partitioning thermodynamics, while limiting the impact on computational efficiency. We provide a quantitative assessment of this water-sensitive OA treatment, focusing on a subdomain over North America. The updated OA scheme predicts a spatiotemporal mean enhancement in surface-level OA mass concentration of 145 % for January 2019 and 76 % for July 2019 compared to GEOS-Chem's most advanced OA scheme. The temporal mean surface-level OA organic mass concentration can increase by up to ∼590 % for January 2019 and ∼280 % for July 2019, with the greatest enhancements occurring over the ocean. The updated OA scheme also quantifies the OA-associated water content. The simulations show how different OA precursors and related OA surrogates contribute and respond to water uptake, including that due to changes in temperature and relative humidity over the diurnal cycle in selected winter and summer months. These results are independent of future CTM improvements involving updates to chemical reaction schemes and emission inventories. Our water-sensitive OA scheme allows for a better representation of the seasonal and regional variations in OA mass concentration in CTMs.

54 ENVIRONMENTAL SCIENCES↗

Evaluating diffusion and the thermodynamic factor for binary ionic mixtures

Molecular dynamics (MD) simulations are a powerful tool for the calculation of transport properties in mixtures. Not only are MD simulations capable of treating multicomponent systems, they are also applicable over a wide range of temperatures and densities. In plasma physics, this is particularly important for applications such as inertial confinement fusion. While many studies have focused on the effect of plasma coupling on transport properties, here we focus on the effects of mixing. We compute the thermodynamic factor, a measure of ideal/non-ideal mixing, for three binary ionic mixtures. Here, we consider mixtures of hydrogen and carbon, hydrogen and argon, and argon and carbon, each at 500 randomly generated state points in the warm dense matter and plasma regimes. The calculated thermodynamic factors indicate different mixing behavior across phase space, which can significantly affect the corresponding mutual diffusion coefficients. As MD simulations are still computationally expensive, we apply modern data science tools to predict the thermodynamic factor over a large phase space. Further, we propose a more accurate approximation to the mutual diffusion coefficient than the commonly applied Darken relation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tethered Small-Molecule Acceptors Simultaneously Enhance the Efficiency and Stability of Polymer Solar Cells

For polymer solar cells (PSCs), the mixture of polymer donors and small-molecule acceptors (SMAs) is fine-tuned to realize a favorable kinetically trapped morphology and thus a commercially viable device efficiency. However, the thermodynamic relaxation of the mixed domains within the blend raises concerns related to the long-term operational stability of the devices, especially in the record-holding Y-series SMAs. Here, a new class of dimeric Y6-based SMAs tethered with differential flexible spacers is reported to regulate their aggregation and relaxation behavior. In their polymer blends with PM6, it is found that they favor an improved structural order relative to that of Y6 counterpart. Most importantly, the tethered SMAs show large glass transition temperatures to suppress the thermodynamic relaxation in mixed domains. For the high-performing dimeric blend, an unprecedented open circuit voltage of 0.87 V is realized with a conversion efficiency of 17.85%, while those of regular Y6-base devices only reach 0.84 V and 16.93%, respectively. Most importantly, the dimer-based device possesses substantially reduced burn-in efficiency loss, retaining more than 80% of the initial efficiency after operating at the maximum power point under continuous illumination for 700 h. In conclusion, the tethering approach provides a new direction to develop PSCs with high efficiency and excellent operating stability.

14 SOLAR ENERGY↗

Accurate prediction of oxygen vacancy concentration with disordered A-site cations in high-entropy perovskite oxides

Abstract Entropic stabilized ABO 3 perovskite oxides promise many applications, including the two-step solar thermochemical hydrogen (STCH) production. Using binary and quaternary A-site mixed {A}FeO 3 as a model system, we reveal that as more cation types, especially above four, are mixed on the A-site, the cell lattice becomes more cubic-like but the local Fe–O octahedrons are more distorted. By comparing four different Density Functional Theory-informed statistical models with experiments, we show that the oxygen vacancy formation energies ( $${E}_{V}^{f}$$ E V f ) distribution and the vacancy interactions must be considered to predict the oxygen non-stoichiometry ( δ ) accurately. For STCH applications, the $${E}_{V}^{f}$$ E V f distribution, including both the average and the spread, can be optimized jointly to improve Δ δ (difference of δ between the two-step conditions) in some hydrogen production levels. This model can be used to predict the range of water splitting that can be thermodynamically improved by mixing cations in {A}FeO 3 perovskites.

08 HYDROGEN↗

Synergistic effects of mixing and strain in high entropy spinel oxides for oxygen evolution reaction

Developing stable and efficient electrocatalysts is vital for boosting oxygen evolution reaction (OER) rates in sustainable hydrogen production. High-entropy oxides (HEOs) consist of five or more metal cations, providing opportunities to tune their catalytic properties toward high OER efficiency. This work combines theoretical and experimental studies to scrutinize the OER activity and stability for spinel-type HEOs. Density functional theory confirms that randomly mixed metal sites show thermodynamic stability, with intermediate adsorption energies displaying wider distributions due to mixing-induced equatorial strain in active metal-oxygen bonds. The rapid sol-flame method is employed to synthesize HEO, comprising five 3d-transition metal cations, which exhibits superior OER activity and durability under alkaline conditions, outperforming lower-entropy oxides, even with partial surface oxidations. The study highlights that the enhanced activity of HEO is primarily attributed to the mixing of multiple elements, leading to strain effects near the active site, as well as surface composition and coverage.

08 HYDROGEN↗

Investigation of Isobaric Mixing as a Mechanism for Boundary‐Layer Cloud Formation

This study investigates the potential role of isobaric mixing in the formation of marine boundary layer clouds. Cloud formation theory emphasizes uplift and adiabatic cooling, but recent observations support the existence of small clouds forming at various altitudes, even below the lifting condensation level. Isobaric mixing of air with different thermodynamic properties can generate localized supersaturation. A Gaussian mixing model is employed to simulate this process, considering the correlation between temperature and water vapor. Cloud droplet size distributions from aircraft measurements show a persistent and prominent mode of small droplets at 9 m, and the size of this mode compares favorably with predictions from the model. The results suggests that isobaric mixing plausibly contributes to the formation of clouds, particularly those observed at multiple altitudes with narrow droplet size distributions. This finding highlights the importance of considering isobaric mixing processes in understanding and modeling cloud formation.

54 ENVIRONMENTAL SCIENCES↗

Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria

Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.

09 BIOMASS FUELS↗

Mixed gauge-global symmetries, elliptic modes, and black hole thermodynamics in Hořava-Lifshitz gravity

In Hořava-Lifshitz gravity, a putative consistent theory of quantum gravity for which there is evidence for both black hole thermodynamics and a holographic construction, spacetime is endowed with a preferred dynamical spacelike foliation. The theory has a leaf reparameterization symmetry that is neither global nor local gauge, hyperbolic and elliptic equations of motion, a lack of splittability, and universal horizon black hole solutions. The reparameterization symmetry is “mixed”: it is a local symmetry in one coordinate yet global on each leaf. More broadly it is an example of both unfree and projectable gauge symmetries. The mixed symmetry and associated charge has not yet been accounted for in calculations of universal horizon thermodynamics in Hořava-Lifshitz gravity. This has led to problems, in particular the failure of the first law in a class of asymptotically AdS solutions where the normal to the leaves of the foliation is not aligned with the time translation Killing vector at infinity. We show how the dynamics of the charge corresponding to this symmetry coupled with the other features above resolves this issue. We then briefly comment how this mixed symmetry, the corresponding charge, and the elliptic equations of motion also conspire to evade recent holographic arguments for only local gauge fields in consistent theories of quantum gravity due to the lack of splittability of the elliptic equation and associated mode.

Global Symmetries↗

Excitation Intensity- and Size-Dependent Halide Photosegregation in CsPb(I 0.5 Br 0.5 ) 3 Perovskite Nanocrystals

Although broad consensus exists that photoirradiation of mixed-halide lead perovskites leads to anion segregation, no model today fully rationalizes all aspects of this near ubiquitous phenomenon. In this work, we quantitatively compare experimental, CsPb(I 0.5 Br 0.5 ) 3 nanocrystal (NC) terminal anion photosegregation stoichiometries and excitation intensity thresholds to a band gap-based, thermodynamic model of mixed-halide perovskite photosegregation. Mixed-halide NCs offer strict tests of theory given physical sizes, which dictate local photogenerated carrier densities. We observe that mixed-anion perovskite NCs exhibit significant robustness to photosegregation, with photosegregation propensity decreasing with decreasing NC size. Observed size- and excitation intensity-dependent photosegregation data agree with model predicted size- and excitation intensity-dependent terminal halide stoichiometries. Established correspondence between experiment and theory, in turn, suggests that mixed-halide perovskite photostabilities can be predicted a priori using local gradients of (empirical) Vegard’s law expressions of composition-dependent band gaps.

36 MATERIALS SCIENCE↗

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of New Additions to OLI Software in Predicting Mercuric and Mercurous Species in Liquid Waste Operations

Speciation of mercury during the pretreatment steps of tank waste processing is critical to successful mercury removal prior to vitrification during Liquid Waste Operations (LWO) at SRS. OLI software has been used to predict mercury speciation and activity throughout LWO. The OLI software operates based on a thermodynamic framework called the Mixed Solvent Electrolyte (MSE) framework. The MSE framework allows prediction in theoretically infinitely dilute to concentrated mixtures (e.g., purely solute solutions). Before modification to the MSE framework databanks, certain critical mercury species were missing in the MSE databank, and some thermodynamic data needed to be updated for the OLI software to accurately predict mercury chemical species in SRS waste tanks. To better reflect streams across LWO, new mercury species were integrated into the MSE database. To evaluate the changes to the OLI MSE framework per the Technical Task Request (TTR) and the Task Technical and Quality Assurance Plan (TTQAP), waste stream compositions from Tanks 38, 43, and Tank 50 decontaminated salt solution (DSS) were used as model inputs. Models were developed and executed using both the old and new databases. Compositional analyses from caustic Tank 50 DSS and caustic Tanks 38 and 43 were used as the input streams. These streams represent the most comprehensive chemical data sets where both mercury and tank constituents were measured together. Results for Tank 50 DSS predict HgO as the predominant species in both databases. Both methyl and dimethyl Hg species are present when the new database is ‘on’ and are not predicted with the new database turned ‘off’. The new database predicts a greater amount of HgO and a greater fraction of it in the solid phase. Pourbaix diagrams (potential vs. pH) generated for each Tank 50 DSS were identical regardless of which database was used. Elemental Hg and HgO were predicted in the water stable region under basic conditions. Tanks 38 and 43 follow similar trends as the Tank 50 DSS models. Unlike Tanks 38 and 50 DSS, the Tank 43 Pourbaix plot shows a region of stability for an aqueous HgOHCO3 - species between approximately pH 7-11. In all streams, when MeHg+ is included in the inputs, the new database predicts aqueous MeHgOH as the dominant species. If elemental or dimethyl mercury is in the waste stream, the new database model predicts they are unchanged and remain in those states and quantities. Additionally, the total mercury values are reported for both the measured input data and the OLI output data for all considered tanks. The summary indicates that the percentage error between the measured and calculated values is less than 1% in all cases The reconciliations and generation of the Pourbaix diagrams for Tank 50 DSS took approximately ten times longer with the new database ‘on’. In addition, over the course of that time, models with the new database ‘on’ were more likely to crash or display an error. Some modest performance improvements were noted when modeling with an i7 processor versus an i5. An example error is found in Appendix A. Furthermore, Appendix B provides V&V for two chemical systems analyzed with the OLI software, results were satisfactory. It is recommended to utilize the new databases (i.e., HCO.ddb and SR-Hg.ddb) in future Savannah River Mission Completion applications of OLI to represent pseudo steady-state. Furthermore, the integration and utilization of the new databases (i.e., HCO.ddb and SR-Hg.ddb) in modeling applications (e.g., Aspen) is also recommended.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Studies on the Encapsulation of Brine Generated from a Process Using Selective Electrodialysis Reversal

This topical report summarizes the laboratory testing, geochemical modeling, and mineralogical testing conducted to examine the efficacy of wastewater encapsulation for the brine generated from the Saltworks Flex electrodialysis reversal technology. Within the laboratory testing, the short- and long-term physical properties of seven separate mix designs were examined along with the leaching of each mix design. Additionally, geochemical thermodynamic modeling was used along with mineralogical analyses to help gain a fundamental understanding of the chemistry involved in the encapsulation process. Encapsulation is a process that results in the immobilization of constituents of interest by creating an engineered material that contains the waste within its matrix, both physically and chemically, in such a manner that the receiving environment is protected from adverse effects. Within the context of this study, encapsulation mix designs involved mixing the Saltworks brine with fly ash and various additives, principally lime (calcium oxide) and cement, to form a paste that eventually sets and hardens into a monolith upon placement in a landfill. Moisture conditioning fly ash with wastewater sources for on-site landfill disposal is a process that is in use at some sites and has been presented as a viable option to dispose of liquid wastes. However, because chemical reactions that can lead to sequestration of constituents of interest are not optimized, leaching can occur. Therefore, in this study, mixes were created to mimic moisture conditioning, as well as engineered pastes of varying mix designs. The study found that engineered encapsulation mixes provide significant resistance to leaching compared to moisture conditioning alone. The best performing encapsulation mix used lime as the additive, which proved superior to the use of cement. After 90 days of curing, this encapsulated mix was 5.5 times stronger than the moisture conditioned mix. Additionally, this mix exhibited an approximate 2.5 order of magnitude drop in hydraulic conductivity compared to the moisture conditioned mix. When comparing leaching, nearly all the chloride was found to evolve from the moisture conditioned mix, whereas only 6.7% of the chlorine contained in the best performing encapsulated mix leached out. Mineralogical analysis via X-ray diffraction (XRD) showed a significant concentration of the chloride sequestering mineral, Friedel’s salt, in the best performing encapsulation mix, whereas none was observed in the moisture conditioned mix, pointing to it being the primary chemical mechanism for chloride retention, along with the physical impermeability of the material. With respect to the modeling study, the model was successful in predicting the presence of most mineral phases observed in the XRD analysis.

01 COAL, LIGNITE, AND PEAT↗

Ship-based Observations and Climate Model Simulations of Cloud Phase over the Southern Ocean

The Southern Ocean (SO) clouds exert a significant influence on the Earth’s radiation budget. Here we analyzed ship-based remote sensing observations of SO clouds over a five-month long DOE ARM Measurements of Aerosols, Radiation and Clouds over the Southern Oceans (MARCUS) field campaign to better understand cloud phase variability. We developed a method to classify eight categories of hydrometeors (ice, liquid, mixed phase, rain, drizzle, snow, aerosols, and clear sky) based on measurements of lidar, radar and radiosondes. Cloud thermodynamic phases (liquid, ice and mixed phase) at coarser scales were further derived to compare with the DOE Energy Exascale Earth System Model version 1 (E3SMv1) simulation and the Earth Model Column Collaboratory (EMC 2 ) instrument simulator. For a scale-aware comparison with climate model simulations, we found that spatially averaging the raw remote sensing data (e.g., backscatter, reflectivity) results in increased cloud cover and cloud liquid, whereas, directly averaging cloud phase from higher to lower resolution maintains clear air regions and is thus recommended for in-cloud frequency comparisons. For cloud thermodynamic phases in stratiform clouds, the E3SM underestimates cloud ice and overestimates cloud liquid at temperatures between -40 and 0°C. When latitudes increase, both observations and simulations show a transition of dominant phase from liquid to ice for cloud tops as well as for the entire cloud columns, but the model underestimation of ice phase is more severe at higher latitudes. Such model bias is unlikely caused by spatial scale differences or lack of heterogeneity in cloud vertical structure in the simulation.

54 ENVIRONMENTAL SCIENCES↗

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗