Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction Kinetics”

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

Facile Synthesis of Rhodium‐Based Nanocrystals in a Metastable Phase and Evaluation of Their Thermal and Catalytic Properties

Abstract Controlling the polymorphism of metal nanocrystals is a promising strategy for enhancing properties and discovering new phenomena. However, previous studies on Rh nanocrystals have focused on their thermodynamically stable face‐centered‐cubic (fcc) phase. Herein, a facile synthesis of Rh‐based nanocrystals featuring the metastable hexagonal close‐packed (hcp) phase is reported by using Ru seeds in their native hcp phase to template the deposition of Rh atoms. The success of such phase‐controlled synthesis relies on the templating effect promoted by the small lattice mismatch between Ru and Rh and the slow dropwise titration of the precursor at an elevated temperature, ensuring the layer‐by‐layer growth mode and thus the formation of a conformal hcp‐Rh shell. Faster injection rate of Rh(III) precursor leads to the formation of a rough Rh shell in the conventional fcc phase due to accelerated reaction kinetics. Considering both thermodynamic and kinetic aspects of this system, the hcp‐Rh phase is favored when the low surface energy from smooth overlayers balances the high bulk energy of the metastable phase, achieved through tight control of reaction rates and deposition patterns. These Ru hcp @Rh hcp core–shell nanocrystals demonstrate thermal stability up to 400 °C, while exhibiting higher catalytic activity toward ethanol oxidation reaction compared to Ru hcp @Rh fcc counterparts.

Nguyen, Quynh N.↗
DOE OSTI2025

Coupling of Lipid Peroxidation and Criegee Intermediate Mediated Autoxidation in the Heterogeneous Oxidation of Linoleic Acid Aerosols

Autoxidation is an established mechanism for the degradation of organic molecules, which is relevant in the atmosphere, combustion processes, the environment, and the rancidification of lipids (commonly known as lipid peroxidation). Autoxidation proceeds via radical chain reactions involving hydroxyl (•OH), peroxy (RO2•), and alkoxy radicals, which are also prominent oxidants in the atmosphere. Recent reports have provided evidence for an alternative autoxidation mechanism driven instead by Criegee intermediates (CIs), which are produced from the reaction of β-hydroxy peroxy radicals (β-OH-RO2•). This work evaluates the contributions of these two mechanisms in the •OH initiated heterogeneous oxidation of linoleic acid (LA) aerosols. Reaction kinetics and product distributions are monitored using a vacuum ultraviolet photoionization aerosol time-of-flight mass spectrometer. To explain the observed kinetics, a kinetic model is developed that incorporates both the conventional peroxidation and alternative CI-mediated autoxidation mechanisms. We observe that the CI-mediated autoxidation pathways enhance the heterogeneous autoxidation rate, while the peroxidation reactions, although present, contributes less to the overall oxidation rate. α-Acyloxyalkyl hydroperoxides (AAHPs) are identified as key indicators for bimolecular reactions of CI with LA, highlighting the role of LA as a CI scavenger. Moreover, the measured functionalized LA products with hydroxyl or carbonyl group(s), serve as markers for the peroxidation reactions. In summary, this work presents a quantitative framework to understand the coupled reaction network of •OH, RO2•, β-OH-RO2• radicals, and CI in driving heterogeneous autoxidation, which is crucial for understanding degradation mechanisms of organic molecules in the environment and atmosphere.

Criegee intermediate↗
DOE OSTI2026

Transient Catalytic Reaction Analysis Through Signal Defragmentation

The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst↗
DOE OSTI2025

Photochemistry of Hypervalent Iodoazide Derivatives

The photochemical properties and reaction mechanisms of a series of hypervalent iodoazide compounds (R–IN3) were investigated, with substituents (−CH 3 , −H, −CF 3 ) tuning the electronic density of the phenyl ring. With the help of UV irradiation, ultrafast time-resolved spectroscopy, and density functional theory (DFT) calculations, we elucidated the mechanisms of azide radical (N 3 • ) release and its subsequent reactivity. UV–vis spectroscopy reveals that the photoconversion rates follow the trend CH3 > H > CF3, aligning with DFT-calculated ΔG values for ring-opening transitions. Homolytic cleavage of the I–N bond is identified as the dominant pathway for N3• generation upon UV irradiation, occurring within 400 fs. The released azide radicals react with the solvent molecules, while the iodo radical R–I C • fragments undergo a thermodynamically uphill lactone ring-opening (RO) reaction and subsequent hydrogen atom abstraction from the solvent to form carboxylic acids R–I–COOH, as validated by NMR and IR spectroscopy. The study also highlights the role of substituents in influencing reaction kinetics and intermediate stability, with electron-donating groups accelerating the release of the N 3 • species. This work bridges experimental observations with computational predictions, offering a foundation for future advancements in azide-based reactions and materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Role of Transport and Kinetics in Growth of Renal Stones

Renal stone disease is not only a concern on earth but could conceivably pose as a serious risk to the astronauts health and safety in Space. In this paper, a combined transport-kinetics model for growth of calcium oxalate crystals is presented. The model is used to parametrically investigate the growth of renal calculi in urine with a focus on the coupled effects of transport and surface reaction on the ionic concentrations at the surface of the crystal and their impact on the resulting growth rates. It is shown that under nominal conditions of low solution supersaturation and low Damkohler number that typically exist on Earth, the surface concentrations of calcium and oxalate approach their bulk solution values in the urine and the growth rate is most likely limited by the surface reaction kinetics. But for higher solution supersaturations and larger Damkohler numbers that may be prevalent in the microgravity environment of Space, the calcium and oxalate surface concentrations tend to shift more towards their equilibrium or saturation values and thus the growth process may be limited by the transport through the medium. Furthermore, parametric numerical studies suggest that changes to the renal biochemistry of astronauts due in space may promote development of renal calculi during long duration space expeditions.

Kassemi, Mohammad↗
NASA NTRS1979

Development of a Process for a High Capacity Arc Heater Production of Silicon for Solar Arrays

A program was established to develop a high temperature silicon production process using existing electric arc heater technology. Silicon tetrachloride and a reductant (sodium) are injected into an arc heated mixture of hydrogen and argon. Under these high temperature conditions, a very rapid reaction is expected to occur and proceed essentially to completion, yielding silicon and gaseous sodium chloride. Techniques for high temperature separation and collection were developed. Included in this report are: test system preparation; testing; injection techniques; kinetics; reaction demonstration; conclusions; and the project status.

Reed, W. H.↗
NASA NTRS2011

Evaluation of Bosch-Based Systems Using Non-Traditional Catalysts at Reduced Temperatures

Oxygen and water resupply make open loop atmosphere revitalization (AR) systems unfavorable for long-term missions beyond low Earth orbit. Crucial to closing the AR loop are carbon dioxide reduction systems with low mass and volume, minimal power requirements, and minimal consumables. For this purpose, NASA is exploring using Bosch-based systems. The Bosch process is favorable over state-of-the-art Sabatier-based processes due to complete loop closure. However, traditional operation of the Bosch required high reaction temperatures, high recycle rates, and significant consumables in the form of catalyst resupply due to carbon fouling. A number of configurations have been proposed for next-generation Bosch systems. First, alternative catalysts (catalysts other than steel wool) can be used in a traditional single-stage Bosch reactor to improve reaction kinetics and increase carbon packing density. Second, the Bosch reactor may be split into separate stages wherein the first reactor stage is dedicated to carbon monoxide and water formation via the reverse water-gas shift reaction and the second reactor stage is dedicated to carbon formation. A series system will enable maximum efficiency of both steps of the Bosch reaction, resulting in optimized operation and maximum carbon formation rate. This paper details the results of testing of both single-stage and two-stage Bosch systems with alternative catalysts at reduced temperatures. These results are compared to a traditional Bosch system operated with a steel wool catalyst.

Abney, Morgan B.↗
DOE OSTI2023

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗
DOE OSTI2023

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Arcjet thruster research and technology, phase 1

The objectives of Phase 1 were to evaluate analytically and experimentally the operation, performance, and lifetime of arcjet thrusters operating between 0.5 and 3.0 kW with catalytically decomposed hydrazine (N2H4) and to begin development of the requisite power control unit (PCU) technology. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. The VNAP2 flow code was used to analyze arcjet nozzle performance with non-uniform entrance profiles. Viscous losses become dominant beyond expansion ratios of 50:1 because of the low Reynolds numbers. A survey of vortex phenomena and analysis techniques identified viscous dissipation and vortex breakdown as two flow instabilities that could affect arcjet operation. The gas kinetics code CREK1D was used to study the gas kinetics of high temperature N2H4 decomposition products. The arc/gas energy transfer is a non-equilibrium process because of the reaction rate constants and the short gas residence times. A thermal analysis code was used to guide design work and to provide a means to back out power losses at the anode fall based on test thermocouple data. The low flow rate and large thermal masses made optimization of a regenerative heating scheme unnecessary.

Knowles, Steven C.↗
DOE OSTI2024

Dehydration kinetics of nanoconfined water in beryl probed by high temperature single crystal synchrotron X-ray diffraction

Understanding changes in material properties through external stimuli plays a key role in validating the expected performance of materials and engineering material properties in a controlled manner. Here, we introduce a fundamental protocol to deduce dehydration reactions kinetics of water confined in nanopore channels, with the cyclosilicate beryl as the scaffold of interest, using time-resolved synchrotron X-ray diffraction (SXRD), in the temperature interval of 298–1038 K. The temperature-dependent intensity (I) of the strongest reflection (112) was used as the crystallite variable. An estimation of an isobaric thermal crystallite coefficient, k, analogous with the isobaric thermal expansion coefficient, established the rate of relative crystallization as a function of temperature, $\frac{∂I}{∂T}$. A plot of lnk and $\frac{1}{T}$ gives rise to two kinetic steps, indicating a slow dehydration stage up to ~ 700 K and a fast dehydration stage up to the investigated temperature 1038 K. The crystal structure of beryl determined up to 1038 K, in temperature increment as small as 10 K, indicates the presence of channel ions Na and Fe and a gradual decrease of water upon heating.

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

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

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

Hybrid Density Functional Theory Comparison of Oxygen Release and Solvent Decomposition Kinetics on Li x NiO 2 Surfaces

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

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

Radiation-Induced Effects on the Extraction Properties of Hexa-n-octylnitrilo-triacetamide (HONTA) Complexes of Americium and Europium

The candidate An(III)/Ln(III) separation ligand hexa-n-octylnitrilo-triacetamide (HONTA) was irradiated under envisioned SELECT (Solvent Extraction from Liquid waste using Extractants of CHON-type for Transmutation) process conditions (n-dodecane/0.1 M HNO 3 ) using a solvent test loop in conjunction with cobalt-60 gamma irradiation. The extent of HONTA radiolysis and complimentary degradation product formation was quantified by HPLC-ESI-MS/MS. Further, the impact of HONTA radiolysis on process performance was evaluated by measuring the change in 243 Am and 154 Eu distribution ratios as a function of absorbed gamma dose. HONTA was found to decay exponentially with increasing dose, affording a dose coefficient of d = (4.48 ± 0.19) × 10 –3 kGy –1 . Multiple degradation products were detected by HPLC-ESI-MS/MS with dioctylamine being the dominant quantifiable species. Both 243 Am and 154 Eu distribution ratios exhibited an induction period of ~70 kGy for extraction (0.1 M HNO 3 ) and back-extraction (4.0 M HNO 3 ) conditions, after which both values decreased with absorbed dose. The decrease in distribution ratios was attributed to a combination of the destruction of HONTA and ingrowth of dioctylamine, which is capable of interfering in metal ion complexation. Furthermore, the loss of HONTA with absorbed gamma dose was predominantly attributed to its reaction with the n-dodecane radical cation (R˙ + ). These R˙ + reaction kinetics were measured for HONTA and its 241 Am and 154 Eu complexes using picosecond pulsed electron radiolysis techniques. All three second-order rate coefficients (k) were essentially diffusion limited in n-dodecane indicating a significant reaction pathway: k(HONTA + R˙ + ) = (7.6 ± 0.8) × 10 9 M –1 s –1 , k(Am(HONTA) 2 + R˙ + ) = (7.1 ± 0.7) × 10 10 M –1 s –1 , and k(Eu(HONTA) 2 + R˙ + ) = (9.5 ± 0.5) × 10 10 M –1 s –1 . HONTA-metal ion complexation afforded an order-of-magnitude increase in rate coefficient. Additional nanosecond time-resolved measurements showed that both direct and indirect HONTA radiolysis yielded the short-lived (< 100 ns) HONTA radical cation and a second long-lived (µs) species identified as the HONTA triplet excited state. The latter was confirmed by a series of oxygen quenching picosecond pulsed electron measurements, affording a quenching rate coefficient of k( 3 [HONTA]* + O 2 ) = 2.2 × 10 8 M –1 s –1 . Overall, both the HONTA radical cation and triplet excited state are important precursors to the suite of measured HONTA degradation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Direct propylene epoxidation via water activation over Pd-Pt electrocatalysts

Direct electrochemical propylene epoxidation by means of water-oxidation intermediates presents a sustainable alternative to existing routes that involve hazardous chlorine or peroxide reagents. Here, we report an oxidized palladium-platinum alloy catalyst (PdPtO x /C), which reaches a Faradaic efficiency of 66 ± 5% toward propylene epoxidation at 50 milliamperes per square centimeter at ambient temperature and pressure. Embedding platinum into the palladium oxide crystal structure stabilized oxidized platinum species, resulting in improved catalyst performance. The reaction kinetics suggest that epoxidation on PdPtO x /C proceeds through electrophilic attack by metal-bound peroxo intermediates. This work demonstrates an effective strategy for selective electrochemical oxygen-atom transfer from water, without mediators, for diverse oxygenation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Novel laser gain and time-resolved FTIR studies of photochemistry

Several techniques are discussed which can be used to explore laboratory photochemical processes and kinetics relevant to planetary atmospheres; these include time-resolved laser gain-versus-absorption spectroscopy and time-resolved Fourier transform infrared (FTIR) emission studies. The laser gain-versus-absorption method employed tunable diode and F-center lasers to determine the yields of excited photofragments and their kinetics. The time-resolved FTIR technique synchronizes the sweep of a commercial FTIR with a pulsed source of light to obtain emission spectra of novel transient species in the infrared. These methods are presently being employed to investigate molecular photodissociation, the yields of excited states of fragments, their subsequent reaction kinetics, Doppler velocity distributions, and velocity-changing collisions of translationally fast atoms. Such techniques may be employed in future investigations of planetary atmospheres, for example to study polycyclic aromatic hydrocarbons related to cometary emissions, to analyze acetylene decomposition products and reactions, and to determine spectral features in the near infrared and infrared wavelength regions for planetary molecules and clusters.

Leone, Stephen R.↗
DOE OSTI2026

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS↗