Search NASA⌕ Search

SEARCH · Search NASA

Results for “grouping optimization”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Influence of Demographic Variables on the Pooled Rideshare Acceptance Model Multigroup Analyses (PRAMMA)

Building on our prior research with a national survey sample of 5385 US participants, the Pooled Rideshare Acceptance Model (PRAM) was built upon two factor analyses. This exploratory study extends the PRAM framework using the Pooled Rideshare Acceptance Model Multigroup Analyses (PRAMMA) to examine how 16 demographic variables influence and interact with the acceptance of Pooled Rideshare (PR), filling a gap in understanding user segmentation and personalization. Using a national sample of 5385 US participants, this methodological approach allowed for the evaluation of how PRAM variables such as safety, privacy, service experience, and environmental impact vary across diverse groups, including gender, generation, driver’s license, rideshare experience, education level, employment status, household size, number of children, income, vehicle ownership, and typical commuting practices. Factors such as convenience, comfort, and passenger safety did not show significant differences across the moderators, suggesting their universal importance across all demographics. Furthermore, geographical differences did not significantly impact the relationships within the model, suggesting consistent relationships across different regions. The findings highlight the need to move beyond a “one size fits all” approach, demonstrating that tailored strategies may be crucial for enhancing the adoption and satisfaction of PR services among various demographic groups. The analyses provide valuable insight for policymakers and rideshare companies looking to optimize their services and increase user engagement in PR.

moderator↗
DOE OSTI2025

Performance evaluation of the SPES FEBIAD ion source

SPES (Selective Production of Exotic Species) is the Radioactive Ion Beam (RIB) facility that is currently in the installation phase at the Istituto Nazionale di Fisica Nucleare in the Laboratori Nazionali di Legnaro (LNL-INFN). The main purpose of the facility is to deliver neutron-rich radioactive beams with a range of mass between 60 and 160 amu produced as fission products originated within a Uranium Carbide (UC x ) target impinged by high energy protons (40–70 MeV). In the last few years, the SPES group and their collaborators have focused on multiple aspects: development of various ISOL targets, ion sources, a high-power beam dump, versatile Front-End designs, primary proton beam and RIB diagnostics, design of beam optic components, and remote target handling. The development and optimization of ionization sources represents one of the most important R&D activities. The collaboration between SPES and ISOLDE (CERN) groups has enabled the evaluation of the performances of a novel FEBIAD ion source taking advantage of the ISOLDE off-line Front-End. In this paper, the results of efficiency and transversal emittance measurements of the new SPES FEBIAD ion source are presented and discussed. Particular attention is dedicated to the description of the experimental setup and to the comparison of experimental data with simulation to promote the reproducibility of the experimental tests with other types of sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Structure–Activity Relationships of Inactive‐Conformation Binding EGFR Inhibitors: Linking the ATP and Allosteric Pockets

ABSTRACT The epidermal growth factor receptor (EGFR) tyrosine kinase is an important therapeutic target in non‐small cell lung cancer (NSCLC). However, the continual emergence of resistance mutations in the treatment of EGFR mutation‐positive NSCLC with currently approved tyrosine kinase inhibitors warrants the development of next‐generation inhibitors. Since research for ATP‐competitive EGFR tyrosine kinase inhibitors (TKIs) that extend into the back pocket has been neglected in the recent past, we survey the extent to which such binding functional groups can be incorporated into an ATP‐site imidazole scaffold. We find thatmeta‐substituted amide linkers derivatized with fluorine in 2,6‐positions and/or a hydroxy group in 3‐position of the back pocket phenyl exhibit the highest potency. Structural insights into how the back pocket groups are bound through points of connection provide new directions for the discovery and optimization of inactive conformation targeting agents in EGFR and other kinases.

Chemistry↗
DOE OSTI2024

Crystal structure of perfluorononanoic acid, C 9 HF 17 O 2

The crystal structure of perfluorononanoic acid (PFNA) was solved via parallel tempering using synchrotron powder diffraction data obtained from the Brockhouse X-ray Diffraction and Scattering (BXDS) Wiggler Lower Energy (WLE) beamline at the Canadian Light Source. PFNA crystallizes in monoclinic space group P2 1 /c (#14) with lattice parameters a = 26.172(1) Å, b = 5.6345(2) Å, c = 10.9501(4) Å, and β = 98.752(2)°. The crystal structure is composed of dimers, with pairs of PFNA molecules connected by hydrogen bonds via the carboxylic acid functional groups. The Rietveld-refined structure was compared to a density functional theory-optimized structure, and the root-mean-square Cartesian difference was larger than normally observed for correct powder structures. The powder data likely exhibited evidence of disorder which was not successfully modeled.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Strain Effects in SrHfO 3 Films Grown by Hybrid Molecular Beam Epitaxy

Perovskite oxide heterostructures host a large number of interesting phenomena such as ferroelectricity, which are often driven by octahedral distortions in the crystal that may induce polarization. SrHfO 3 (SHO) is a perovskite oxide with a pseudocubic lattice parameter of 4.08 Å that previous density functional theory (DFT) calculations suggest can be stabilized in a ferroelectric P4mm phase when stabilized with sufficient compressive strain. Additionally, it is insulating and possesses a large band gap and a high dielectric constant, making it an ideal candidate for oxide electronic devices. Here, to test the viability of epitaxial strain as a driver of ferroic phase transitions, SHO films were grown by hybrid molecular beam epitaxy (hMBE) with a tetrakis(ethylmethylamino)hafnium(IV) source on GdScO 3 and TbScO 3 substrates. Strained SHO phases were characterized using X-ray diffraction, X-ray absorption spectroscopy, and scanning transmission electron microscopy to determine the space group of the strained films, with the results compared to those of DFT-optimized models of phase stability versus strain. Contrary to past reports, we find that compressively strained SrHfO 3 undergoes octahedral tilt distortions without associated ferroelectric polarization and most likely takes on the I4/mcm phase with the a 0 a 0 c – tilt pattern.

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

Synergistic Solvent-Surface Interactions Enable Alkyne Semihydrogenation at Palladium

Enabling higher yield and better selectivity for fine-chemical synthesis through heterogeneous catalysis is intricately linked to the interplay of active sites, reaction conditions, and mass transfer influence provided by the catalyst. Alkyne semihydrogenation is ubiquitous in the production of bulk chemicals in the pharmaceutical, polymer, or fine-chemical industries, but product selectivity remains a major challenge. Here, in this study, we demonstrate that the design of catalysts encompassing nickel (Ni) foams as contiguous monolith supports, decorated with ultralow loading of Pd/PdO x nanoparticles on a carbonized polydopamine interface and tuned with a thin layer of Al 2 O 3 , in conjunction with an optimized reaction environment leads to highly selective alkyne semihydrogenation. The reactions demonstrate good functional group tolerance and applicability to flow reactor systems. Combined computational and experimental studies are presented to describe the synergistic effect between the solvent-surface interaction and the degree of Pd surface reduction that are necessary to promote this selectivity. The system highlights the opportunity for catalyst-solvent codesign as a benign alternative to more complex reactants featuring extrinsic poisons or less-favored dopants.

atomic layer deposition↗
DOE OSTI2024

Faster Tensor Network Decoding for Topological Quantum Codes

We present a fast and Bayes-optimal-approximating tensor network decoder for planar quantum LDPC codes based on the tensor renormalization group algorithm, originally proposed by Levin, and Nave. By precomputing the renormalization group flow for the null syndrome, we need only recompute tensor contractions in the causal cone of the measured syndrome at the time of decoding. This allows us to achieve an overall runtime complexity of ($pnχ^6$) where p is the depolarizing noise rate, and χ is the cutoff value used to control singular value decomposition approximations used in the algorithm. We apply our decoder to the surface code in the code capacity noise model and compare its performance to the original matrix product state (MPS) tensor network decoder introduced by Bravyi, Suchara, and Vargo. The MPS decoder has a p-independent runtime complexity of $\mathcal{O}(nχ^3)$ resulting in significantly slower decoding times compared to our algorithm in the low-p regime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Relationship Between Filler Type, Thermomechanical Properties, and Aging of RTV Silicone Foams

Room-temperature vulcanizing (RTV) silicone foams are used in many industrial applications that require the material to perform over long time periods. However, mechanical properties tend to deteriorate when these foams age under a compressive load. The chemical aging is attributed to the presence of unreacted functional groups of the prepolymers, residues from acid, and catalytically active tin (II) species. Here, an optimized thermal treatment of an RTV foam that achieves completion of curing reactions and deactivation of reactive species is proposed. Foams that were thermally aged for three months under compressive load showed no signs of compression set, indicative of the effectiveness of the implemented post-curing approach. In addition, the effects of fillers (diatomaceous earth, fumed silica, and carbon nanofibers) on thermomechanical properties were investigated. Tensile strength, tear strength, and thermal conductivity increased when these fillers were added to the unfilled RTV formulation, with carbon nanofibers (CNFs) being the most effective filler. Rheological studies of RTV formulations indicated that 2.5 wt.% of CNFs is the upper limit that can be added to the RTV formulation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Rank-Limiting Strategies for Optimizing Tensor-Train Finite-Difference Time-Domain Simulations

We introduce rank-limiting strategies to optimize tensor-train decompositions for three-dimensional finite-difference time-domain simulations using the relationship between the tensors and their specific dimensionality. These include the use of hard caps on the inner ranks of the tensor train decomposition and the use of a group rounding algorithm taking into account all field components simultaneously. Here, several numerical examples are considered to verify the efficacy of the proposed optimization strategies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structure-guided utilization of lignocellulose for catalysis, energy, and biomaterials

As a complex composite of cellulose, hemicellulose, and lignin, plant lignocellulose has long served as a major resource for biomass conversion, materials engineering, and bio-based product development. High-resolution structural insights enabled by solid-state nuclear magnetic resonance (ssNMR) now allow the mapping of polymer interfaces, identification of functional group accessibility, and tracking of molecular organization during processing, all of which are critical factors for optimizing catalytic strategies. These insights could drive transformative progress in lignocellulose-based applications, including selective depolymerization, improved pretreatment design, and efficient upcycling of lignin into resins, plastics, and biomedical materials. In industry-relevant contexts, such as biofuel generation and renewable material manufacturing, understanding the hydration dynamics, cross-linking patterns, and structural heterogeneity is also essential. The ability to visualize these features in native biomass presents a unique opportunity to develop new strategies for sustainability and performance. As the structural toolbox continues to expand, it is becoming a central enabler for innovations in renewable energy, green chemistry, and advanced bioproducts.

bioproduct↗
DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Sampling Size Optimization for Bioburden Density Estimation in Planetary Protection

Planetary protection (PP) is a discipline that focuses on minimizing the biological contamination of spacecraft to ensure compliance with international policy. Precise estimation of bioburden - the total number of microbes in or on spacecraft hardware – and the bioburden density are of utmost importance for PP. Such estimation is the way concordance with requirements is demonstrated, and it is critical for quantifying the potential risk of inadvertently contaminating other planetary bodies. Although a suite of molecular techniques have been used to thoroughly characterize and profile the microbiome of various cleanroom environments and spacecraft, the gold standard remains the physical enumeration of microbes via culturing of samples directly taken from spacecraft and associated surfaces. However, due to technical, budgetary, and programmatic constraints, only a manageable portion (around 10%) of the entire spacecraft surface is directly sampled with cotton swabs or wipes. To generate the bioburden current best estimate (CBE) for components not directly verifiable, the accepted approach is to apply a NASA-defined bioburden estimate based on the components’ manufacturing or assembly environment. This approach utilizes a prespecified bioburden density estimation that applies a maximum value across the total surface area of the specified component. For hardware components that underwent similar assembly processes, an implied bioburden is adopted for all components, based on a direct verification of a representative component within the same lot. Once all components have a CBE, the bioburden estimates are generated. In previous publication [ 1], we have shown that statistical risks quantifying the accuracy of the estimates for sampled, prespecified, and implied components can be derived and ranked. For mean squared error (MSE) function, the risks are available analytically and hence a cost function can be obtained to optimize the risks with respect to the sampling area and sampling cost. Since the sampling area and sampling cost are two complimentary variables, their sum will have a well-defined minimum. This paper presents the multivariate optimization of the integrated risk of an empirical Bayes estimator to determine the optimal sampling schedule for a given number of components. It is assumed that given a number of components, N, the bioburden density for each component can either be sampled, implied, or prespecified. The multivariate optimization searches through different options to sample, imply or prespecify the bioburden density for a component, and account for the component’s surface area and cost of sampling. The idea of the optimization is based on the observation that the statistical risk of using an estimator is a monotonically decreasing function of the sampled area. The larger the sampled area, the lower the risk of using the estimator as the estimator becomes more and more accurate as the sampling area increases. On the other hand, the cost of sampling is monotonically increasing as the sampled surface grows. This makes the risk and total cost of sampling complimentary variables which can be counterbalanced to achieve an optimal overall value with respect to the sampled surface. In this paper, the integrated risk has been used to quantify the accuracy of the estimator. This risk has been selected because it depends on neither the true value of the parameter nor on the collected data. The cost of each sample was also available to obtain the total cost of sampling of N components. The paper will present the results based on computer-simulated data as well as the data collected during the InSight mission. The computer-simulated data have N components with randomly generated total areas and each component assigned to one of the three categories according to the method of estimating of bioburden density: sampled, implied, or prespecified. The cost of sampling is also available. The cost of sampling is estimated based on a cost model provided by the planetary protection group at JPL. For this paper, the overall cost was assumed to be a linear function of exposure. The optimization process finds the allocation of the components to the three categories that minimizes the tradeoff between integrated risk and total cost. For the InSight data, a set of components is selected representing all three categories, and optimization is performed to determine if the performed allocation was optimal or if a better allocation could have been obtained. To the best of our knowledge, this work is the first attempt not only perform an accurate estimation of bioburden density but also do it in an optimal way.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Insulating moiré homobilayers lack a threefold symmetric second-harmonic generation

Atoms within moiré bilayers relax in plane to minimize elastic energy; such relaxation brings their space group symmetries down to P1. Here, the ab initio second harmonic generation (SHG) of twisted and atomistically optimized hBN bilayers was determined at four twist angles (θ = 38.21°, 60.00°, 73.17°, and 98.21°) and for three displacements τ measured away from the ground state AA' configuration. All moiré bilayers have a P1 space symmetry after structural optimization. This situation is quite different to monolayers with hexagonal lattices, which retain a threefold symmetry. We point out that the actual symmetries of the SHG reported for hBN bilayers on two experimental works do not coincide with the sixfold symmetric theoretical profiles they provide [either sin 2 ⁡(3⁢Φ) or cos 2 ⁡(3⁢Φ)], and show that the intrinsic low structural symmetry of (atomically optimized) hBN bilayer moirés can in fact be read out from experimental SHG intensity profiles—which are tunable by θ and by the frequency of light ω: The SHG is most definitely not sixfold symmetric because moirés do not retain a threefold symmetry. Furthermore, an extrinsic twofold symmetry of the SHG emission is realized by tilting the pump by an angle α away from the 2D material's normal, regardless of θ and ω. Furthermore, the design of in-plane and ultrathin sources of SHG with low symmetry could be useful for the eventual creation of entanglement sources from 2D materials.

2-dimensional systems↗
DOE OSTI2024

Improving I/O-aware Workflow Scheduling via Data Flow Characterization and trade-off Analysis

The scientific computing paradigm has transitioned from compute-intensive to I/O-intensive and memory-intensive in the past decade, especially when data-driven science has become common practice. Numerous empirical I/O-aware scheduling optimizations have been developed by incorporating I/O capacity and bandwidth as constraints into scheduling. Unfortunately, there is a lack of data flow (I/O) characterization tool and an understanding of trade-offs between concurrency, locality, and I/O bandwidth. To bridge the gap, this work 1) presents a set of descriptors to characterize, organize, and visualize I/O profiles, including flow size, I/O bandwidth, and operation count, which group data flows by I/O types, tasks, and files; 2) proposes an I/O Roofline model-based trade-off analysis to find the optimal trade-off between flow operational intensity, concurrency, and flow performance. The I/O descriptors generate useful insights into complicated I/O behaviors, suggesting distinct concurrency, storage, and scheduling to be used by types, tasks, and files. The proposed trade-off analysis guides scheduling decisions that generate resource assignment with the best flow parallelism. We evaluate our I/O-aware scheduling methodology on a highly I/O-intensive workflow–1000 Genomes. The experimental results demonstrate speedups of up to 2.4× compared to the state-of-the- art methods.

Guo, Luanzheng [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2025

Choline Chloride-Based Water-in-Salt Electrolyte for Efficient Iron Electrodeposition

Electrochemical production of iron is a promising low-cost and modular approach to replace the traditional blast furnace. Aqueous electrolytes for iron electrolysis are advantageous as they can be operated at near-ambient temperatures, but they suffer from inefficiencies due to the parasitic hydrogen evolution reaction. In this work, we identify a new water-in-salt electrolyte (WiSE) based on choline chloride (ChCl) for high coulombic efficiency (>85%) iron deposition. Electrochemical analysis of the partial current densities of iron plating and hydrogen co-evolution revealed that, at optimal WiSE compositions, water reduction is kinetically suppressed resulting in an increase in the Fe plating efficiency. Decreased coordination of water and increased coordination of choline’s alcohol group with the Fe 2+ ion were observed through 1 H NMR providing evidence that water reduction is kinetically suppressed in WiSE. Additionally, Raman spectroscopy revealed that complexation effects (with Cl – ) reduce Fe 2+ diffusion coefficients and corresponding limiting currents as the ChCl concentration is increased. This results in an optimal WiSE composition (4 M ChCl + 1 M FeCl 2 ) that provides kinetic suppression of HER but also low transport resistance to Fe plating yielding 85% coulombic efficiency at high current densities.

Sinclair, Nicholas Scott [Case Western Reserve Uni↗
DOE OSTI2024

Carbon‐13 Hyperpolarization of α‐Ketocarboxylates with Parahydrogen in Reversible Exchange

Abstract Signal Amplification by Reversible Exchange (SABRE) is a relatively simple and fast hyperpolarization technique that has been used to hyperpolarize the α‐ketocarboxylate pyruvate, a central metabolite and the leading hyperpolarized MRI contrast agent. In this work, we show that SABRE can readily be extended to hyperpolarize 13 C nuclei at natural abundance on many other α‐ketocarboxylates. Hyperpolarization is observed and optimized on pyruvate (P 13C =17 %) and 2‐oxobutyrate (P 13C =25 %) with alkyl chains in the R‐group, oxaloacetate (P 13C =11 %) and alpha‐ketoglutarate (P 13C =13 %) with carboxylate moieties in the R group, and phenylpyruvate (P 13C =2 %) and phenylglyoxylate (P 13C =2 %) with phenyl rings in the R‐group. New catalytically active SABRE binding motifs of the substrates to the hyperpolarization transfer catalyst – particularly for oxaloacetate – are observed. We experimentally explore the connection between temperature and exchange rates for all of these SABRE systems and develop a theoretical kinetic model, which is used to fit the hyperpolarization build‐up and decay during SABRE activity.

McBride, Stephen J. [Department of Chemistry North↗
DOE OSTI2025

Optimized NiTe 2 /Ni back contacts via chemical deposition for high performance CdTe solar cells

The development of efficient and stable back contacts remains a major challenge in achieving high performance and long-term stability of CdTe thin-film solar cells. Here, this work revisits the formation of NiTe 2 by chemical bath deposition (CBD) as a back contact for CdTe devices. An optimized CBD recipe, based on high-purity precursors and the addition of copper chloride directly into the bath, was developed and applied to fabricate Cu-doped CdTe solar cells. A modified Cu-free methodology was also applied to Group V doped absorbers. The process included pinhole filling, ion milling, CBD, annealing, and sputtering to form a low-barrier back contact. Devices fabricated using this method achieved consistent open-circuit voltages (V oc ) above 800 mV and fill factors (FF) exceeding 70%. The best Cu-doped devices reached power conversion efficiencies (PCE) above 18 %, and preliminary results with Group V-doped material demonstrated compatibility of the method with high-efficiency, state-of-the-art CdTe devices. This study shows that NiTe 2 /Ni back contacts, formed via an optimized chemical process followed by sputtering of Ni, represent a promising pathway for achieving low-barrier and potentially stable back contacts in modern CdTe photovoltaics.

14 SOLAR ENERGY↗
DOE OSTI2025

Characterization of Kariya ( Hildegardia barteri (Mast.) Kosterm) Seed Oil Fatty Acid Methyl Ester Prepared from Basic Catalytic Transesterification

The rising global energy demand, alongside concerns regarding environmental deterioration due to the use of fossil fuels, has spurred extensive investigation into renewable energy alternatives. Biomass-derived biodiesel, especially from lesser-known oil sources, emerges as a promising option. This research focuses on analyzing the fatty acid methyl esters (FAMEs) derived from Kariya (Hildegardia barteri (Mast.) Kosterm) seed oil through basic catalytic transesterification using gas chromatography–flame ionization detector (GC–FID) analysis, assessing its potential as a biodiesel feedstock. Oil extraction from Kariya seeds was carried out using three solvents (n-hexane, ethanol, and a 1:1 blend of hexane and ethanol), followed by transesterification with methanol. Gas chromatography–mass spectrometry (GC–MS) and GC–FID analyses were utilized to identify and quantify FAMEs in the resulting biodiesel. The results revealed various FAMEs, including methyl myristate, methyl palmitate, methyl stearate, methyl oleate, methyl linoleate, and methyl linolenate. Significant differences in FAME composition were observed among the samples, with hexane–ethanol Kariya oil biodiesel (HE-KOB) showing the highest FAME content (76.1%). This combination of solvents exhibited synergistic effects on the composition of HE-KOB, suggesting potential optimization strategies for biodiesel production. Fourier transform infrared spectroscopy (FTIR) provided additional insights into the molecular composition of the biodiesel samples, confirming their biodiesel nature through the identified functional groups such as methyl, methylene, hydrocarbon, ester, aldehyde, and alkene. Thermogravimetric analysis (TGA) for thermal decomposition also gave an insight into FAME composition and its contribution to the degree of conversion of biodiesel to energy. These findings highlight the feasibility of utilizing Kariya seed oil as a biodiesel feedstock, emphasizing the importance of solvent selection and transesterification conditions in optimizing FAME yield and composition. This research contributes to the exploration of underutilized oil sources for sustainable biodiesel production, aligning with the global shift towards cleaner and renewable energy sources.

09 BIOMASS FUELS↗