Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “combined indicator”

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 307 records · Page 17

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗

Characterizing the Spread of COVID-19 from Human Mobility Patterns and SocioDemographic Indicators

Mobility is an indicator of human movement through space and time. With the increasing availability of geolocated data (from GPS, accelerometers, etc.), it is now possible to examine individual as well as group human mobility patterns. Human mobility is influenced by both intrinsic (i.e. personal motivations) and extrinsic (i.e., events like natural hazards or a pandemic like the COVID-19) factors. However, the intricate relationships between human mobility patterns and sociodemographic characteristics in the context of a pandemic are yet to be fully explored. Our goal is to overcome this gap by using human mobility data at the census block group level from mobile phones and combining those with social vulnerability indicators to examine the overall spread of COVID-19 at local spatial scales. We used 585,878 weekly visits to 37,871 points of interests (POIs) from Safegraph to quantify mobility indices and social distancing metrics in 2,820 census block groups in the city of Los Angeles (LA) - before and during lockdown as well as during the phase1 and phase 2 reopening. Finally, using supervised machine learning algorithms, we classified the census block groups in LA into High, Medium and Low categories that represented the vulnerability of these block groups based on the cumulative number of occurrences of COVID-19 cases till July 24, 2020. Our results indicate that the tree-based classifiers performed well in comparison to the Support Vector Machines and Multinomial Logit models. Gradient Boosting had the highest classification accuracy of 97.4% COVID-19 with an AUC score of 0.987. The block groups with high COVID-19 cases also had a high concentration of socially vulnerable populations, high human mobility index and a low social distancing index.

Roy, Avipsa↗

GPM Ground Validation at NASA Wallops Precipitation Research Facility

The Wallops Precipitation Research Facility (WPRF) at NASA Goddard Space Flight Center, Wallops Island, VA has been established as a semi-permanent super-site for the Global Precipitation Measurement (GPM) Ground Validation (GV) program. WPRF is home to research quality precipitation instruments, including NASA’s S-band dual-polarimetric radar (NPOL), and a network of profiling radars, disdrometers, and rain gauges. This study investigates the statistical agreement of the GPM Core Observatory Dual Frequency Precipitation Radar (DPR), combined DPR-GPM Microwave Imager (GMI) and GMI Level II precipitation retrievals compared to WPRF ground observations from a six-year collection of satellite overpasses. Multi-sensor observations are integrated using the System for Integrating Multiplatform Data to Build the Atmospheric Column (SIMBA) software package. SIMBA ensures measurements recorded in a variety of formats are synthesized into a common reference frame for ease in comparison and analysis. Given that instantaneous satellite measurements are observed above ground level, this study investigates the possibility of a time lag between satellite and surface mass-weighted mean diameter (D(sub m)), reflectivity (Z), and precipitation rate (R) observations. Results indicate that time lags vary up to 30 minutes after overpass time but are not consistent between cases. In addition, GPM Core D(sub m) retrievals are within Level I mission science requirements as compared to WPRF ground observations. Results also indicate GPM algorithms overestimate light rain (< 1.0 mm hr(exp -1)). Two very different stratiform rain vertical profiles show differing results when compared to ground reference data. A key finding of this study indicates multi-sensor DPR/GMI combined algorithms outperform single sensor DPR algorithm.

Charanjit S Pabla↗

Ka-Band Wide-Bandgap Solid-State Power Amplifier: Hardware Validation

Motivated by recent advances in wide-bandgap (WBG) gallium nitride (GaN) semiconductor technology, there is considerable interest in developing efficient solid-state power amplifiers (SSPAs) as an alternative to the traveling-wave tube amplifier (TWTA) for space applications. This article documents proof-of-concept hardware used to validate power-combining technologies that may enable a 120-W, 40 percent power-added efficiency (PAE) SSPA. Results in previous articles [1-3] indicate that architectures based on at least three power combiner designs are likely to enable the target SSPA. Previous architecture performance analyses and estimates indicate that the proposed architectures can power combine 16 to 32 individual monolithic microwave integrated circuits (MMICs) with >80 percent combining efficiency. This combining efficiency would correspond to MMIC requirements of 5- to 10-W output power and >48 percent PAE. In order to validate the performance estimates of the three proposed architectures, measurements of proof-of-concept hardware are reported here.

Epp, L.↗

Classifying Agnostic Biosignatures using Raman, VNIR, and Elemental Data

How can we use our current wealth of terrestrial data, encompassing biogenic and abiogenic systems, to determine the distinguishing properties of life? SCOBI (Statistical Classification of Biosignature Information) uses machine learning techniques to algorithmically identify combinations of measurements that are “indicative of life”. A set of ~1000 observations, comprising elemental abundance, isotopic fractionation, VNIR reflectance, and (in progress) Raman spectra, have been assembled from existing literature and databases. The observations cover systems classified as “indicative alive” (e.g., cells, vegetation), “indicative non-alive” (e.g., fossils, teeth), “mixed indicative” (e.g., soil, pond water), or “non-indicative” (e.g., rocks, meteorites). VNIR data was preprocessed by linear interpolation from 400-2100 nm and smoothed with a Savitzky-Golay filter. To limit the amount of Earth-biochemistry-specific (non-agnostic) information included, the first five spectral features extracted were number of peaks, number of troughs, mean reflectance, mean peak width, and broadest peak width. To help further emphasize agnostic biosignatures, Earth-specific features such as chlorophylls have been manually flagged so that feature importance with and without them can be compared. Classifiers including k-nearest neighbors (KNN), Gaussian Naïve Bayes (GNB), logistic regression (LR), random forest (RF), and support vector machine (SVM) were implemented, as was a combination voting classifier. Performance metrics included false positive rates, false negative rates, and AUC with 50-50 test/train splits (Monte Carlo simulations). Key takeaways from this stage, prior to the inclusion of Raman spectra, are (1) the overall success rate of 0.933 AUC was most heavily influenced by the elemental abundance data; and (2) VNIR reflectance had the lowest classification performance with 0.52 AUC (58% of objects correctly classified). The next steps are to complete integration of Raman spectral data and to improve the approach to pre-processing and feature extraction for both types of spectral data, such as automated baseline removal, whole spectrum matching, and dimensionality reduction.

Biosignatures↗

Classifying Forest Type in the National Forest Inventory Context with Airborne Hyperspectral and Lidar Data

Forest structure and composition regulate a range of ecosystem services, including biodiversity, water and nutrient cycling, and wood volume for resource extraction. Forest type is an important metric measured in the US Forest Service Forest Inventory and Analysis (FIA) program, the national forest inventory of the USA. Forest type information can be used to quantify carbon and other forest resources within specific domains to support ecological analysis and forest management decisions, such as managing for disease and pests. In this study, we developed a methodology that uses a combination of airborne hyperspectral and lidar data to map FIA-defined forest type between sparsely sampled FIA plot data collected in interior Alaska. To determine the best classification algorithm and remote sensing data for this task, five classification algorithms were tested with six different combinations of raw hyperspectral data, hyperspectral vegetation indices, and lidar-derived canopy and topography metrics. Models were trained using forest type information from 632 FIA subplots collected in interior Alaska. Of the thirty model and input combinations tested, the random forest classification algorithm with hyperspectral vegetation indices and lidar-derived topography and canopy height metrics had the highest accuracy (78% overall accuracy). This study supports random forest as a powerful classifier for natural resource data. It also demonstrates the benefits from combining both structural (lidar) and spectral (imagery) data for forest type classification.

random forest↗

Patterns in Mobility and Modification of Middle and High-Latitude Southern Hemisphere Dunes on Mars

Change detection analyses of aeolian bedforms (dunes and ripples), using multitemporal imagesacquired by the Mars Reconnaissance Orbiter High Resolution Imaging Science Experiment (HiRISE), canreveal migration of bedforms on Mars. Here we investigated bedform mobility (evidence of wind-drivenmigration or activity), from analysis of HiRISE temporal image pairs, and dune field modification (i.e., apparentpresence/lack of changes or degradation due to nonaeolian processes) through use of a dune stability indexor SI (1–6; higher numbers indicating increasing evidence of stability/modification). Combining mobility dataand SI for 70 dune fields south of 40°S latitude, we observed a clear trend of decreasing bedform mobility withincreasing SI and latitude. Both dunes and ripples were more commonly active at lower latitudes, althoughsome high-latitude ripples are migrating. Most dune fields with lower SIs (3) were found to be active whilethose with higher SIs were primarily found to be inactive. A shift in prevalence of active to apparently inactivebedforms and to dune fields with SI 2 occurs at ~60°S latitude, coincident with the edge of highconcentrations of H2O-equivalent hydrogen observed by the Mars Odyssey Neutron Spectrometer. This result isconsistent with previous studies suggesting that stabilizing agents, such as ground ice, likely stabilize bedformsand limit sediment availability. Observations of active dune fields with morphologies indicative of stability(i.e., migrating ripples in SI = 3 dune fields) may have implications for episodic phases of reworking or dunebuilding, and possibly geologically recent activation or stabilization corresponding to shifts in climate

Banks, Maria E.↗

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗

Radiation Heat Transfer Procedures for Space-Related Applications

Over the last contract year, a numerical procedure for combined conduction-radiation heat transfer using unstructured grids has been developed. As a result of this research, one paper has been published in the Numerical Heat Transfer Journal. One paper has been accepted for presentation at the International Center for Heat and Mass Transfer's International Symposium on Computational Heat Transfer to be held in Australia next year. A journal paper is under review by my NASA's contact. A conference paper for the ASME National Heat Transfer conference is under preparation. In summary, a total of four (4) papers (two journal and two conference) have been published, accepted or are under preparation. There are two (2) to three (3) more papers to be written for the project. In addition to the above publications, one book chapter, one journal paper and six conference papers have been published as a result of this project. Over the last contract year, the research project resulted in one Ph.D. thesis and partially supported another Ph.D. student. My NASA contact and myself have formulated radiation heat transfer procedures for materials with different indices of refraction and for combined conduction-radiation heat transfer. We are trying to find other applications for the procedures developed under this grant.

Chai, John C.↗

Forest Potential Productivity Mapping by Linking Remote-Sensing-Derived Metrics to Site Variables

A fine-resolution region-wide map of forest site productivity is an essential need for effective large-scale forestry planning and management. In this study, we incorporated Sentinel-2 satellite data into an increment-based measure of forest productivity (biomass growth index (BGI)) derived from climate, lithology, soils, and topographic metrics to map improved BGI (iBGI) in parts of North American Acadian regions. Initially, several Sentinel-2 variables including nine single spectral bands and 12 spectral vegetation indices (SVIs) were used in combination with forest management variables to predict tree volume/ha and height using Random Forest. The results showed a 10–12 % increase in out of bag (OOB) r2when Sentinel-2 variables were included in the prediction of both volume and height together with BGI. Later, selected Sentinel-2 variables were used for biomass growth prediction in Maine, USA and New Brunswick, Canada using data from 7738 provincial permanent sample plots. The Sentinel-2 red-edge position (S2REP) index was identified as the most important variable over others to have known influence on site productivity. While a slight improvement in the iBGI accuracy occurred compared to the base BGI model (~2%), substantial changes to coefficients of other variables were evident and some site variables became less important when S2REP was included.

site productivity↗

Sea-Level Static Testing of the Penn State Two-Dimensional Rocket-Based Combined Cycle (RBCC) Testbed

The present results indicated that: 1.Significant RBCC ejector mode database has been generated for single and twin thruster configuration and for global and local measurements. 2. Ongoing analysis and correlation effort for MSFC CFD modeling and turbulent shear layer analysis was completed. 3. The potential follow-on activities are: detailed measurements of air flow static pressure and velocity profiles; investigation other thruster spacing configurations; performing fundamental shear layer mixing study; and demonstrating single-shot Raman measurements.

Cramer, J. M.↗

Demonstration of Oxygen and Carbon Monoxide Propellants for Mars Missions

Currently, proposed planetary exploration missions must be small, with low costs and a short development time. Relatively high-risk technologies are being accepted for such missions if they meet these guidelines. For a Mars sample-return mission, one of the higher risk technologies is the use of return propellants produced from indigenous materials such as the Martian atmosphere. This consists of 96 percent carbon dioxide, which can be processed into oxygen and carbon monoxide. This year, the NASA Lewis Research Center completed the experimental evaluation and subscale technology development of an oxygen/carbon monoxide propellant combination. Previous research included ignition characterization, combustion performance, and heat transfer characterization with gaseous propellants at room temperature. In this year s tests, we studied the ignition characteristics and combustion of oxygen and carbon monoxide at near liquid temperatures. The mixture ratio boundaries for oxygen and carbon monoxide were determined as a function of propellant temperature in a spark torch igniter. With both propellants at room temperature, the ignition range was between 0.50 and 1.44; and with both propellants chilled to near-liquid temperatures, it was between 2.4 and 3.1. Statistical analysis of the mean value of the ignition boundaries provided models that describe the combination of oxygen temperature, carbon monoxide temperature, and mixture ratio that resulted in ignition. This range is the larger boxed area shown in the figure. The smaller boxed area indicates the range at which there is a 90-percent confidence that ignition will occur. The relatively small range at only 90-percent confidence indicates that using the oxygen/carbon monoxide combination as its own ignition source may not be the best design for a remote engine operating on Mars. Tests also were performed in a simulated small rocket engine that used oxygen/hydrogen combustion gases as the ignition source for oxygen/carbon monoxide. In these experiments, the oxygen/carbon monoxide was successfully ignited in eight of eight tests at a mixture ratio of 0.52. In addition, the oxygen/carbon monoxide maintained steady combustion after the oxygen/hydrogen ignition source was removed, verifying that the oxygen/carbon monoxide rocket engine should continue to be included in mission plans as return propulsion from Mars.

Linne, Diane L.↗

Silicon-Lattice-Matched Boron-Doped Gallium Phosphide: A Scalable Acousto-Optic Platform

The compact size, scalability, and strongly confined fields in integrated photonic devices enable new functionalities in photonic networking and information processing, both classical and quantum. Gallium phosphide (GaP) is a promising material for active integrated photonics due to its high refractive index, wide bandgap, strong nonlinear properties, and large acousto-optic figure of merit. Here, this study demonstrates that silicon-lattice-matched boron-doped GaP (BGaP), grown at the 12-inch wafer scale, provides similar functionalities as GaP. BGaP optical resonators exhibit intrinsic quality factors exceeding 25,000 and 200,000 at visible and telecom wavelengths, respectively. It further demonstrates the electromechanical generation of low-loss acoustic waves and an integrated acousto-optic (AO) modulator. High-resolution spatial and compositional mapping, combined with ab initio calculations, indicate two candidates for the excess optical loss in the visible band: the silicon-GaP interface and boron dimers. These results demonstrate the promise of the BGaP material platform for the development of scalable AO technologies at telecom and provide potential pathways toward higher performance at shorter wavelengths.

36 MATERIALS SCIENCE↗

Systems Analysis and Optimization of Circular PET Packaging Supply Chains in the United States: Environmental and Socioeconomic Impacts

Many actions are underway at global, national, and local levels to address the plastic waste problem and transition toward a circular economy of plastics. Studies evaluating environmental and socioeconomic impacts of such a transition are lacking. Here, the purpose of this study is to conduct a national systems analysis of polyethylene terephthalate (PET) packaging supply chains in the United States. Material flow data was combined with environmental and socioeconomic indicators to evaluate and compare the sustainability of the linear PET packaging supply chain, current (2019) supply chain, and possible future circular supply chain options in the United States. Environmentally optimal circular US PET packaging material flows showed 31% and 38% savings of GHG emissions and energy demand, respectively, with a circularity of 77% when compared with a linear supply chain. Additionally, the environmentally optimal system showed higher employment (29%) and wages (31%) than a linear system, but with a 5% decrease in revenue generation. A socioeconomically optimal circular PET supply chain showed increased employment (by 52%), wages (by 67%), and revenues (by 1%), with a circularity of 59% when compared with the linear system. However, it showed 14% higher GHG emissions than a linear system, indicating a trade-off between environmentally and socioeconomically optimal circular PET packaging systems. Overall, linear-to-circular material flow transition may not necessarily lead to increased revenues and decreased environmental impacts of the entire system, but it does benefit society due to increased employment and wages. Future systems analysis work should focus on improving data quality for environmental and socioeconomic dimensions.

PET↗

A Rapid Population Assessment Method for Wild Pigs Using Baited Cameras at 3 Study Sites

ABSTRACT Reliable and efficient population estimates are a critical need for effective management of invasive wild pigs ( Sus scrofa ). We evaluated the use of 10‐day camera grids for rapid population assessment (RPA) of wild pigs at 3 study sites that varied in vegetation communities and wild pig densities. Study areas included Buck Island Ranch, Florida; Tejon Ranch, California; and the Savannah River Site, South Carolina, USA, during 2016–2018. Rapid population assessments grids were composed of baited camera traps spaced approximately 500 or 750 m apart. Two RPA grids were deployed per study site and each grid was deployed twice (4–6 months apart) to assess changes in response to season or population control efforts. We assessed the ability of RPA grids to track population trends, how camera number influenced estimate precision, and how relative abundance indices related to density estimates. We detected changes in occupancy probability, detection probability, and N‐mixture estimates following removal operations and between seasons, but the ability of RPA grids to track population trends was dependent on the statistical method used and number of cameras traps. Increasing the number of cameras traps used in RPA grids increased precision, and these results can be used in determining survey design and estimate choice. We found that estimates of occupancy probability, detection probability, and N‐mixture estimates were positively correlated with spatially explicit capture–recapture density estimates. Thus, these less labor‐intensive estimates from RPA grids showed potential to index the relative abundance of wild pigs in some systems. Our evaluation of RPAs indicates that using study‐specific combinations of statistical method and number of cameras can provide a useful tool for monitoring wild pig presence, tracking population trends, and evaluating the effectiveness of management actions. © 2020 The Wildlife Society.

Schlichting, Peter E.↗

Coupled solute effects enable anomalous high-temperature strength and stability in nanotwinned Al alloys

Nanoprecipitates or grain refinement can effectively enhance the mechanical strength of Al alloys, but the room-temperature strengths of precipitation hardened and nanocrystalline Al alloys often fall below 1 GPa. Furthermore, they are largely plagued by precipitous mechanical softening at elevated temperature below 300°C, mostly due to degraded microstructural stability. Here, we report a mechanism of coupled solute effect in nanotwinned Al-Fe-Ti alloys that enables stability of nanograins up to 400°C and an unprecedented high-temperature flow stress of ~ 1.7 GPa at 300°C. The supersaturated Fe solutes in Al act as effective grain refiner, forming superstrong solid solution alloys. More importantly, empirical evidence combined with first principle calculations indicate that the Ti solutes delay the agglomeration of Fe solutes, thereby remarkably extending the temperature window for the stability of nanograins in nanotwinned Al alloys. Finally, this study highlights the opportunity to design ultrastrong and stable nanostructured alloys for potential high temperature applications via a coupled solute effect.

36 MATERIALS SCIENCE↗

Improving intrinsic oxygen reduction activity and stability: Atomic layer deposition preparation of platinum-titanium alloy catalysts

Improved activity and stability Pt-based catalysts for the oxygen reduction reaction (ORR) are needed to perpetuate the deployment of polymer electrolyte fuel cells (PEFCs) in the transportation sector. In this study, we use atomic layer deposition of TiO 2 and Pt coupled with thermal reductive annealing to prepare Pt 3 Ti electrocatalysts. The atomic level synthetic control resulted in Pt 3 Ti nanoparticles with high ORR performance, including a mass activity of 1.84 A/mg Pt and excellent electrochemical stability. The Pt 3 Ti nanoparticles show excellent specific activity — 5.3-fold higher than commercial Pt/C and 3-fold higher than polycrystalline Pt, exceeding the performance of any PtTi catalysts reported to date. Combined experimental and computational efforts indicate that Pt enrichment on the Pt 3 Ti enhances the activity, and the intrinsic stability of the Pt 3 Ti phase provides durability. Subsequently, this knowledge, along with the facile fabrication of alloys by atomic layer deposition, can be leveraged to designed improved performance catalysts.

25 ENERGY STORAGE↗

Mechanical separations of corn stover anatomical fractions in an integrated feedstock preprocessing system: An experimental and data-driven modeling study

High variabilities of material attributes in lignocellulosic biomass present risks for biofuel and biochemical productions and must be mitigated via preprocessing. Since almost no mechanical device is originally designed for processing biomass, how to operate existing apparatuses with efficient performance has not been investigated extensively. This work presents a study on an integrated screening and air classification to separate cobs and stalks from husks and leaves in corn stover. Prototype machine learning models were developed to assess the feasibility of predicting the process outcome based on the measurable parameters. The models trained upon limited experimental data rendered decent predictive accuracy of yield and purity. The experimental data and modeling results collectively suggest decreasing throughput leads to a higher purity. To the contrary, if throughput increases, a lower purity is likely. A possible trade-off between yield and purity of the separated streams indicates the need for optimal combinations of feedstock size, moisture, and throughput to achieve optimized separations. The results of this study also suggest the need to further improve model predictability by developing more accurate formulations for physics governing the integrated unit operations. To accomplish this, additional experimental data needs to be generated for model training.

09 - BIOMASS FUELS↗