Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variability analysis”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restrains for non-protein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in Chlorophyll (Chl) proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C 2 symmetry of Photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and Cryogenic electron microscopy (CryoEM) PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restrains play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested, from cyanobacteria to algae. Furthermore, these results highlight both the opportunities and limitations inherent in the structure-based analysis of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

14 SOLAR ENERGY↗
DOE OSTI2024

Informing forest carbon inventories under the Paris Agreement using ground-based forest monitoring data

Human interactions with forests have shaped Earth's climate for millennia and will continue to do so as we target net-zero emission goals. Accurately characterizing these climate impacts requires making reliable forest carbon data available for forest monitoring and planning. Here, we develop a semi-automated process for submitting forest carbon measurements from the largest relevant scientific database to the International Panel on Climate Change's Emission Factor Database, which currently has sparse forest carbon data. Building this bridge from scientific research to international policy is an important step towards managing forests in a net-zero motivated future. Humans have been influencing Earth's climate via transformative impacts on forests for millennia, and forests are now recognized as critical to climate change mitigation under the Paris Agreement. The efficacy of climate change mitigation planning and reporting depends on quality data on forest carbon (C) stocks and changes. The Emission Factor Database (EFDB) of the International Panel on Climate Change (IPCC) is intended to be a definitive source for such data, but needs comprehensive and well-documented data to be so. To facilitate submission of forest C estimates from scientific studies to EFDB, we develop and document a process for semi-automated data submission from the Global Forest C database (ForC v4.0), which is the largest compilation of ground-based forest C estimates. We then assess the data currently available through ForC and provide recommendations for improving forest data collection, analysis, and reporting. As of September 2024, ForC contained ~19,286 records potentially relevant to EFDB, 1068 of which had been submitted and posted to EFDB. These represented 19% of the total EFDB records for forest land. Records were unevenly distributed across variables and geographic regions. ForC records (37%) reviewed could not be submitted because the original publication lacked required information. In the future, ground-based forest C estimates should target gaps in the record, and studies should ensure that they report all information necessary for inclusion in EFDB. Given that climate change is rapidly impacting the world's forests, timely reporting of recent estimates will be critical to accurate forest C inventories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using occupant interaction with advanced lighting systems to understand opportunities for energy optimization: Control data from a hospital NICU

For this study, lighting system control data from five patient rooms in a neonatal intensive care unit (NICU) over a 25-week monitoring period were analyzed to better understand how the NICU staff and families interacted with an automatic tunable lighting system (variable spectrum and intensity) with manual override options. The data were analyzed to determine the amount of time spent in available lighting modes, estimate energy consumption, and to observe any patterns in how the lighting system was being used by occupants. Further analysis revealed opportunities for optimizing the system for occupant needs and energy savings. The patient room occupants were clearly engaged in using the lighting system throughout the 24-hour day to meet needs and preferences. Giving occupants control over the lighting system did not result in a considerable increase in energy use. Implementing a time-out would minimize the extended periods of time when the lights were left on without need at the highest brightness level, further decreasing energy use. This analysis reveals the potential of lighting control data to contribute to energy savings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Meteorological environments associated with California wildfires and their potential roles in wildfire changes during 1984-2017

California has been seeing more wildfires in recent years, resulting in huge economic losses and threatening human health. Clarifying the meteorological environments of wildfires is foundational to improving the understanding and prediction of wildfires and their impacts. Here, 1535 California wildfires during 1984-2017 are systematically investigated. Based on two key meteorological factors - temperature and moisture anomalies - all wildfires are classified into four groups: hot-dry, hot-wet, cold-dry, cold-wet. Most (~60%) wildfires occurred on hot-dry days. Compositing the large-scale environments of the four groups shows that persistent high pressure and strong northeasterly wind descending from inland favor hot-dry conditions for wildfires. This analysis also reveals an important role of anomalies in southerly onshore flow that supports stronger convection, accompanied by more lightning flashes that provide a triggering mechanism for wildfires under hot-wet conditions. Self-organizing map analysis lends confidence in the large-scale meteorological pattern for dominating hot-dry wildfires in California. Besides wildfire occurrence, wildfire size is also influenced by meteorological anomalies through their magnitudes. Among them, moisture anomaly explains the largest fraction (~69%) of variability in wildfire sizes. Large-scale meteorological anomalies are found to play an important role in the devastating 2018 wildfire season in California. During 1984-2017, wildfire burned area has significantly increased by ~3.6% per year, indicating a doubling of burned area in 2017 relative to 1984, with the trend dominated by hot-dry wildfires in summer. Drying and warming in conjunction with strengthening of the high pressure in summer support more frequent and larger wildfires in California.

58 GEOSCIENCES↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Variational Autoencoder Model Toward Molecular Structure Representation Learning of Fuels

Here, in this work, a Variational Autoencoder (VAE)-based data-driven modeling framework is developed with the overarching goal of enabling fuel design. The VAE model is trained on a large dataset with several chemical species to learn a compressed latent space molecular representation. Chemical structure in the form of Simplified Molecular Input Line Entry System (SMILES) string is fed as input, encoded into the VAE latent space, and decoded back to the SMILES string using Long Short-Term Memory (LSTM) networks. Complexities of the VAE training loss function are thoroughly examined by varying the weightage (beta (𝜷) parameter) of the latent space regularization term, thereby assessing the balance between reconstruction accuracy and validity, and focusing on both accurate molecular structure reconstruction and latent space consistency. Two different strategies for 𝜷 variation are evaluated: linear annealing and cyclic annealing. In addition, the impact of total correlation adjustment and hierarchical priors is also studied with regard to the balance between reconstruction fidelity and latent space regularization, and potential issues such as posterior collapse, over-regularization, and poor disentanglement of latent variables. Overall, the best performance of the model is achieved with hierarchical priors and incrementally increasing 𝜷 from 0 to a threshold value of 0.25 over 75 epochs. The generative VAE model can be readily coupled with Quantitative Structure–Property Relationship (QSPR) analysis to develop an integrated end-to-end framework for fuel-property prediction and molecular design of novel promising fuels.

fuel design↗
DOE OSTI2024

Understanding GPU Memory Corruption at Extreme Scale: The Summit Case Study

GPU memory corruption and in particular double-bit errors (DBEs) remain one of the least understood aspects of HPC system reliability. Albeit rare, their occurrences always lead to job termination and can potentially cost thousands of node-hours, either from wasted computations or as the overhead from regular checkpointing needed to minimize the losses. As supercomputers and their components simultaneously grow in scale, density, failure rates, and environmental footprint, the efficiency of HPC operations becomes both an imperative and a challenge. We examine DBEs using system telemetry data and logs collected from the Summit supercomputer, equipped with 27,648 Tesla V100 GPUs with 2nd-generation high-bandwidth memory (HBM2). Using exploratory data analysis and statistical learning, we extract several insights about memory reliability in such GPUs. We find that GPUs with prior DBE occurrences are prone to experience them again due to otherwise harmless factors, correlate this phenomenon with GPU placement, and suggest manufacturing variability as a factor. On the general population of GPUs, we link DBEs to short- and long-term high power consumption modes while finding no significant correlation with higher temperatures. We also show that the workload type can be a factor in memory’s propensity to corruption.

Oles, Vlad↗
DOE OSTI2025

Misclassification in Workers’ Telecommuting Frequency Choices Using a Generalized Extreme Value Model

Telecommuting frequency is a response variable collected in travel surveys and is, therefore, prone to errors leading to mismeasurements or misclassification. Misclassification of explanatory variables is a common risk when using statistical modeling techniques. We define “misclassification” as a response reported or recorded in the wrong category; for example, a variable is recorded as a 1 when it should be 0. Here, in this context, this study aims to develop a statistical model to analyze telecommuting data which accounts for potential misclassification errors by building on existing literature in econometrics. The empirical analysis was undertaken using the 2017 National Household Travel Survey (NHTS) and the general extreme value (GEV) models available in the literature. Specifically, the frequency of telecommuting days was analyzed using the negative binomial (NB) model recast as the multinomial logit (MNL) model. By nature—and consistent with other studies—NHTS data are prone to errors that can be classified as intentional or unintentional misinformation provided by the person being interviewed. Ignoring these errors while modeling telecommuting frequencies using standard discrete count models can result in biased parameter estimates. The misclassification parameter was calculated for both over-reporting and under-reporting scenarios. The misclassification errors can be as high as 14% over-reported and 10% under-reported, particularly for the neighboring values. Statistical fit comparison between the models shows that models that ignore misclassification have worse data fit and biased parameter estimates with significant policy implications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

The MOOSE Thermal Hydraulics Module

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source object-oriented finite element framework written in C++ (Lindsay et al., 2022). The Thermal Hydraulics Module (THM) is an optional MOOSE physics module that provides capabilities for studying thermal hydraulic systems. Its core capability lies in assembling a network of coupled components, for instance, pipes, junctions, and valves. THM provides several new systems to MOOSE to enable and facilitate thermal hydraulic simulations, most notably the Components system, which provides a higher-level syntax to MOOSE’s lower-level objects. This system is extensible by the user, but the current library primarily includes components based on a one-dimensional, single-phase, variable-area, compressible flow model, as well as heat conduction.

42 ENGINEERING↗
DOE OSTI2025

Using multimodal X-ray computed tomography to advance 3D petrography: A non-destructive investigation of olivine inside a carbonaceous chondrite

Rocks form in three dimensions through time and studying them provides information from inside dynamic systems we cannot otherwise observe. Yet how we typically access the interior of the rocks themselves to gain that information may limit our understanding and influence how we reconstruct the processes that formed them. Here, we demonstrate combined non-destructive 3D X-ray imaging techniques that produce quantitative densitometric and crystallographic maps of entire individual grains inside a rock. Olivine grains throughout a sample of the carbonaceous chondrite Northwest Africa (NWA) 11346 were each characterized by size, shape, composition, zoning intensity, and crystallographic orientation. The addition of 3D crystallographic mapping to calibrated 3D densitometric analysis—used to calculate chemical composition—demonstrates a fully non-destructive petrographic method and provides unique insight. For instance, in our case, using crystallographic data to delineate individual grains and then measuring the 3D size, shape, and composition of each distinguishes variably reset relict grains from those later crystallized after a melting event. Intersection in a 2D slice could not have led to this interpretation because the integration of three-dimensional size, rounding, composition, location, and crystallographic orientation measured from each grain forms the key patterns. Multimodal laboratory X-ray imaging has strong potential to advance 3D petrography.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vitrification of Hanford Tank 241-AP-107 with Recycled Condensate

During the vitrification of nuclear waste at the Hanford Waste Treatment and Immobilization Plant (WTP) – the primary mission of the U.S. Department of Energy Office of River Protection – the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, a test apparatus was designed to mimic the EMF evaporator and used to concentrate a volume of condensate that had been previously produced during the vitrification of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The test apparatus successfully concentrated the AP-107 condensate by a factor of 10 while retaining over 90 % of the technetium-99 ( 99 Tc), Cs, and I inventory. A second portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC, given to Pacific Northwest National Laboratory, and combined with the AP-107 condensate concentrate after undergoing solids filtration and cesium removal by ion exchange. This combination served to approximate the recycling action to be performed at the WTP. After the addition of glass-forming chemicals (GFCs), the combined AP-107 waste and AP-107 condensate concentrate were processed in the CLSM to produce a glass, called AP-107-1R, that was designed to satisfy the WTP baseline requirements (Kim et al. 2012). During the 8.87 hours of processing, 7.27 kg of AP-107-1R glass were produced for an average glass production rate of 1739 kg m 2 d -1 . Compared to the previous run in the CLSM without recycled condensate, the run with the recycle had a greater average glass production rate, but the rate was within the potential range of variability when processing melter feeds with similar composition in the CLSM. The glass produced from the AP-107 recycle run in the CLSM was within 10 % of the target AP-107-1R glass composition with respect to the primary glass components. Analysis of the minor component impurities revealed that their content in the glass product had approached their nominal target after 2 turnovers of the glass inventory in the CLSM while the activity of the minor radionuclides was retained in the glass product. The 99 Tc and total cesium content in the combined AP-107 waste and recycled condensate were maintained at concentrations expected to be experienced at the WTP. During processing in the CLSM, at discrete sampling time periods, the target 99 Tc/Cs mass ratio in the glass formulation varied from 0.9 to 62.9. Across this range, the Cs retention in the glass ranged from 53 to 60 %, while the retention from the entire runtime totaled 68 %, values which align with Cs retention in other scaled melter systems while processing LAW melter feeds at 99 Tc/Cs mass ratios ranging from 1 to 100. The 99 Tc retention in the glass ranged from 22 to 32 %, primarily due to the cold-cap coverage on the glass melt surface, the area covered by reacting melter feed, varying from ~80 % to ~95 % during processing, demonstrating greater volatility of 99 Tc from the glass while more surface was exposed, as expected based on previous 99 Tc volatility studies.

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

Solar Power + Electric Vehicle Charging: Capturing Synergies in Minnesota

Despite the state policy support, significant market barriers to decarbonization of transportation remain. Analysis conducted by the National Renewable Energy Laboratory (NREL) and others confirm that without routine access to charging infrastructure, both at home and in non-residential locations, bullish market forecasts are unlikely to materialize. Similarly, the expansion of EV charging and variable renewable energy deployment each have potentially disruptive and costly effects on the distribution grid. Finally, electrification of passenger vehicles is not enough to decarbonize transportation without pathways to carbon-free electrification, linking EV charging and clean energy generation. This document examines the opportunity to jointly deploy solar energy and EV charging infrastructure to capture synergistic value, mitigate for risks associated with integrating and scaling both solar and EV charging, and accelerate their mutual deployment.

14 SOLAR ENERGY↗
DOE OSTI2023

Techno-Economic Analysis and Optimization of a Compressed-Air Energy Storage System Integrated with a Natural Gas Combined-Cycle Plant

To address the rising electricity demand and greenhouse gas concentration in the environment, considerable effort is being carried out across the globe on installing and operating renewable energy sources. However, the renewable energy production is affected by diurnal and seasonal variability. To ensure that the electric grid remains reliable and resilient even for the high penetration of renewables into the grid, various types of energy storage systems are being investigated. In this paper, a compressed-air energy storage (CAES) system integrated with a natural gas combined-cycle (NGCC) power plant is investigated where air is extracted from the gas turbine compressor or injected back into the gas turbine combustor when it is optimal to do so. First-principles dynamic models of the NGCC plant and CAES are developed along with the development of an economic model. The dynamic optimization of the integrated system is undertaken in the Python/Pyomo platform for maximizing the net present value (NPV). NPV optimization is undertaken for 14 regions/cases considering year-long locational marginal price (LMP) data with a 1 h interval. Design variables such as the storage capacity and storage pressure, as well as the operating variables such as the power plant load, air injection rate, and air extraction rate, are optimized. Results show that the integrated CAES system has a higher NPV than the NGCC-only system for all 14 regions, thus indicating the potential deployment of the integrated system under the assumption of the availability of caverns in close proximity to the NGCC plant. The levelized cost of storage is found to be in the range of 136–145 $/MWh. Roundtrip efficiency is found to be between 74.6–82.5%. A sensitivity study with respect to LMP shows that the LMP profile has a significant impact on the extent of air injection/extraction while capital expenditure reduction has a negligible effect.

25 ENERGY STORAGE↗
DOE OSTI2020

Validation of Reanalysis-Based Offshore Wind Resource Characterization Using Lidar Buoy Observations

The offshore wind industry in the U.S. is gaining strong momentum to achieve sustainable energy goals, and the need for observations to provide resource characterization and model validation is greater than ever. Pacific Northwest National Laboratory (PNNL) operates two lidar buoys for the U.S. Department of Energy (DOE) in order to collect hub height wind data and associated meteorological and oceanographic information near the surface in areas of interest for offshore wind development. This work evaluates the performance of commonly used reanalysis products and spatial approximation techniques using lidar buoy observations off the coast of New Jersey and Virginia, USA. Reanalysis products are essential tools in order to set performance expectations and quantify the wind resource variability at a given site. Long-term accurate observations at typical wind turbine hub-heights have been lacking at offshore locations. Using wind speed observations from both lidar buoy deployments, biases and degrees of correspondence for the Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), and the analysis system of the Rapid Refresh (RAP) are examined both at hub height and near surface. Results provide insights on the performance and uncertainty of using reanalysis products for long-term wind resource characterization.

lidar buoy, model validation, offshore wind energy↗
DOE OSTI2022

Heavy Duty Natural Gas Single Cylinder Research Engine Installation, Commissioning, and Baseline Testing

Natural Gas (NG) Internal Combustion Engines (ICE) are a promising alternative to diesel engines for on-road heavy-duty applications to reduce greenhouse gas and harmful pollutant emissions. NG engines have not been widely adopted due to the lower thermal efficiency compared with diesel engine counterparts. To develop the base knowledge required to reach the desired efficiency, a Single Cylinder Engine (SCE) is the most effective platform to acquire reliable and repeatable data. A SCE test cell was developed using a Cummins 15-liter six-cylinder heavy-duty engine block modified to fire one cylinder (2.5-liter displacement). A Woodward Large Engine Control Module (LECM) is integrated to permit implementation of real-time advanced combustion control. Intake and exhaust characteristics, fuel composition, and exhaust gas recirculated substitution rate (EGR) are fully adjustable. A high-speed data acquisition system acquires in-cylinder, intake, and exhaust pressure for combustion analysis. The baseline testing shows reliable and consistent results for engine thermal efficiency, indicated mean effective pressure (IMEP), and coefficient of variance of the IMEP over a wide range of operating conditions while achieving effective control of all engine control and operation variables. This test cell will be used to conduct a research program to develop new and innovative control algorithms and CFD optimized combustion chamber designs, allowing ultra-high efficiency and low emissions for NG ICE heavy-duty on-road applications.

03 NATURAL GAS↗
DOE OSTI2023

Decision Support Tool for Solar Energy Cybersecurity Policy and Regulation

The Decision Support Tool helps users address four discrete challenges: (1) complex requirements of relevant codes and standards, (2) technical complexity of solar assets, (3) undefined cyber risk severity, and (4) unclear roles and responsibilities. Content includes an analysis of cyber vulnerability risks, a decision support resource for policymakers to mitigate cyber risks to solar photovoltaic systems, and background resources for informing policy development. A key component of this tool is the Probable Risk Assessment (PRA), based on established, formal variables, models, and consequences, which helps users understand the determined risk and assigned ownership of the physical components. It helps states draw logical lines between vulnerabilities and mitigative solutions, of which some have been further detailed as Cybersecurity Advisory Team for State Solar (CATSS) tools.

14 SOLAR ENERGY↗
DOE OSTI2024

Innovative Joining Method for Hybrid Composites with Tailored Performance

In various industries, such as renewables and aerospace, the demand for carbon fiber is rapidly growing to meet the requirements of high-performance applications with optimized designs. Traditionally, these optimized designs have been limited to single-material composites due to process constraints. However, this study aims to explore the feasibility of transitioning between two different fiber types for pultrusion and filament winding processes, thereby enhancing design flexibility, and optimizing performance. This research introduces a novel method of hybrid reinforcement through splicing techniques. By employing a splicing method, fibers are merged prior to infusion, resulting in a transition of higher strength compared to traditional joints. Moreover, this technique enables the creation of composites with variable compositions, allowing for the selective placement of properties in a process previously characterized by static properties. The objective of this paper is to evaluate the mechanical properties of dry fiber splicing and its impact on the resulting composite material. Through comprehensive analysis, we aim to provide insights into the feasibility and advantages of employing this innovative splicing technique in composite manufacturing processes.

Guzorek, Steven↗
DOE OSTI2025

Influence of particle size on NIR spectroscopic characterization of sorghum biomass for the biofuel industry

NIR spectroscopy is a rapid and accurate green technology for high-throughput biomass characterization, including sorghum (Sorghum bicolor), a promising energy crop for the biofuel industry. This study assessed the influence of particle size on NIR spectroscopic analysis (wavelength range: 867–2535 nm) of sorghum biomass composition. Grown under field conditions, a total of 113 types of genetically diverse sorghum accessions were dried, ground, and sieved (<250, 250–600, 600–850, and > 850 µm particle size) for developing partial least square regression (PLSR) prediction models for moisture, ash, extractive, glucan, xylan, acid-soluble lignin (ASL), acid-insoluble lignin (AIL), and total lignin (ASL + AIL). Overall, smaller particle sizes provided better model performance, while no single particle size provided the best performance for all the selected components. With only 9 selected bands and 4 latent variables (LVs), the best PLSR model was obtained for moisture with particle size of 600–850 µm with the square root of the coefficient of determination (R) of 0.85, the ratio of prediction to deviation (RPD) of 2.2, and the root mean square error (RMSE) of 0.46 % in external validation. Similar model performances were also obtained for ash, extractive, glucan, and xylan. This study showed that size reduction could effectively improve NIR spectroscopic analysis for lipid-producing sorghum biomass for the biofuel industry.

09 BIOMASS FUELS↗