Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Safety parameters”

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 145 records · Page 8

STAT7 v1.2 User Guide: The STAT7 Code for Statistical Propagation of Uncertainties in Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for thermal hydraulics analyses of research and test reactors. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to the fit functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one latera stripe of a fuel element. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads

Abstract. Floating wind turbines must withstand a unique and challenging set of loads from the wind and ocean environment. To de-risk development, accurate predictions of these loads are necessary. Uncertainty in modeling predictions leads to larger required safety factors, increasing project costs and the levelized cost of energy. Complex aero-hydro-elastic modeling tools use many input parameters to represent the wind, waves, current, aerodynamic loads, hydrodynamic loads, and structural properties. It is helpful to understand which of these parameters ultimately drives a design. In this work, an ultimate and fatigue-proxy load sensitivity analysis was performed with 35 different input parameters, using an elementary effects approach to identify the most influential parameters for a case study involving the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine atop the OC4-DeepCwind semisubmersible during normal operation. The importance of each parameter was evaluated using 14 response quantities of interest across three operational wind speed conditions. The study concludes that turbulent wind velocity standard deviation is the parameter with the strongest sensitivity; this value is important not just for turbine loads, but also for the global system response. The system center of mass in the wind direction is found to have the highest impact on the system rotation and tower loads. The current velocity is found to be the most dominating parameter for the system global motion and consequently the mooring loads. All tested wind turbulence parameters in addition to the standard deviation are also found to be influential. Wave characteristics are influential for some fatigue-proxy loading but do not significantly impact the extreme ultimate loads in these operational load cases. The required number of random seeds for stochastic environmental conditions is considered to ensure that the sensitivities are due to the input parameters and not due to the seed. The required number of analysis points in the parameter space is identified so that the conclusions represent a global sensitivity. The results are specific to the platform, turbine, and choice of parameter ranges, but the demonstrated approach can be applied widely to guide focus in parameter uncertainty.

17 WIND ENERGY↗

A review of technical and regulatory limits for hydrogen blending in natural gas pipelines

We report there is rising interest globally in the use of hydrogen for the provision of electricity or heat to industry, transport, and other applications in low-carbon energy systems. While there is attention to build out dedicated hydrogen infrastructure in the long-term, blending hydrogen into the existing natural gas pipeline network is also thought to be a promising strategy for incorporating hydrogen in the near-term. However, hydrogen injection into the existing gas grid poses additional challenges and considerations related to the ability of current gas infrastructure to operate with blended hydrogen levels. This review paper focuses on analyzing the current understanding of how much hydrogen can be integrated into the gas grid from an operational perspective and identifies areas where more research is needed. The review discusses the technical limits in hydrogen blending for both transmission and distribution networks; facilities in both systems are analyzed with respect to critical operational parameters, such as decrease in energy density, increased flow speed and pressure losses. Safety related challenges such as, embrittlement, leakage and combustion are also discussed. The review also summarizes current regulatory limits to hydrogen blending in different countries, including ongoing or proposed pilot hydrogen blending projects.

08 HYDROGEN↗

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters

Accurately predicting the criticality of an experiment before interacting with the experimental components is very important for criticality safety. Radiation transport software can be utilized to calculate the effective neutron multiplication factor of a nuclear system. Because of the integral nature of the effective neutron multiplication factor, the value calculated contains various sources of nuclear-data induced uncertainty. The sensitivity analysis and data assimilation technique presented in this paper exhibit one possible method of identifying and reducing the effective neutron multiplication factor nuclear-data induced uncertainty. The results presented in this work show that it is possible to use relative sensitivity coefficients of the prompt neutron decay constant and the effective delayed neutron fraction to 239 Pu nuclear data to reduce nuclear-data induced uncertainties in the effective neutron multiplication factor. This work has been utilized by members of the Los Alamos National Laboratory project EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) for optimally designing a new experiment, which will be used to reduce compensating errors in 239 Pu nuclear data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High Flux Isotope Reactor Low Enriched Uranium U-10Mo Fuel Design Parameters

Activities to convert the HFIR from HEU to LEU are ongoing as part of the US Department of Energy (DOE) National Nuclear Security Administration (NNSA) nuclear nonproliferation mission. Design activities to study the conversion of HFIR from HEU to LEU fuel explored different fuel design features and shapes with a uranium-molybdenum (U-10Mo) monolithic alloy fuel. This high-density alloy contains 90 wt % uranium and 10 wt % molybdenum and has a uranium density of 15.318gU/cm 3 . The goal of these studies is to generate several candidate HFIR LEU fuel designs of varying fuel fabrication complexity that meet the current HEU performance metrics and safety requirements. Recent advancements in modeling and simulation tools and design methods enabled a thorough analysis of the available design space with U-10Mo fuel. A surrogate model used this analysis as training data to quickly determine the performance of a design given specific design parameters. An optimization module used this surrogate model to quickly search this multidimensional search space given specific desired performance characteristics. This approach was made possible by the large available design space with U-10Mo fuel. Shift, a Monte Carlo tool optimized for high-performance computing (HPC) architectures, was used for faster calculation and better data management for reactor physics simulations. Once most of these design studies were complete, a new suite called the Python HFIR Analysis and Measurement Engine (PHAME) was developed to connect all fuel design analysis steps, making design studies more efficient and reproducible. The post-processing capabilities of these new tools are leveraged for the information provided herein. Leveraging these tools, several candidate fuel designs were selected with varying levels of feature complexity and reactor performance. This report provides design feature details for four selected HFIR LEU U-10Mo fuel designs and their corresponding performance and safety metrics. Nominal best-estimate design parameters and irradiation conditions, including fission rate densities, power densities, heat fluxes, and cumulative fission densities, are provided. Simulations show that the high uranium density of U-10Mo fuel provides a large potential design space that enables various LEU designs to meet HEU core performance metrics and safety requirements with a power increase from 85 MW (HEU) to 95 MW or 100 MW (LEU).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

CORE DESIGN AND NEUTRONIC ANALYSIS OF THE EUROPEAN SODIUM FAST REACTOR WITH METALLIC FUEL

The current ESFR (European Sodium Fast Reactor) design was proposed and in-depth evaluated in the frame of the past ESFR-SMART project. As a follow-up project, the ESFR-SIMPLE has been launched with the aim of challenging the current commercial-size ESFR design in terms of safety features and economic performance. Among the new safety measures to be developed and assessed in ESFR-SIMPLE, the current oxide fuel ESFR design will be challenged by a modified version of the core with metallic fuel. This intends to conclude on what types of benefits can be obtained with high-density fuel, under similar safety and design constraints. In this paper, the designing approach for enabling the use of metallic fuel in the current ESFR core is described and a preliminary neutronic evaluation is carried out. The optimal configuration is established through the optimization of key neutronic parameters aiming at the potential reduction of the plutonium inventory. The resulting core configuration serves as a basis for further safety assessment analyses, which will provide insight into the advantages and drawbacks of the two types of fuels.

Jiménez-Carrascosa, Antonio↗

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗

Hand Calculation Methods for Nuclear Criticality Safety

This primer provides an overview of the most common hand calculation methods used for criticality safety calculations. The most widely used tools available to a nuclear criticality safety (NCS) practitioner are probably the common Monte Carlo or deterministic criticality safety codes, which can be used to model very complex systems. However, use of these codes can obscure the parameters to which a particular fissile system may be sensitive, whereas the hand calculation methods can be used to delve into the ways each parameter may affect the reactivity of a fissile material system. Furthermore, practitioners must avoid using computer codes as devices that take inputs and simply provide outputs (i.e., a “black box”). Many years ago, pioneers such as Joe Thomas, David Smith, and Hugh Paxton, among others in the field of nuclear criticality safety, took the time before the advent of high-speed desktop computers to create simple hand methods for criticality safety analyses. Some of the methods can be used for single fissile units; others are applicable to fissile units arranged into simple array configurations. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the various results. The NCS practitioner will need to spend time to master the methods that could be most useful; however, they can provide the practitioner with fast and accurate answers to criticality safety problems if they are used correctly and if critical data exist for the problem at hand. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information more quickly than using a criticality code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sensitivity Analyses of Natural Convection in the AP1000 Passive Containment Cooling System Following LBLOCA Using CFD

The AP1000 power plant has implemented multiple layers of passive safety systems to enhance reactor safety and ensure system integrity during postulated scenarios, such as loss of coolant accident (LOCA) and main steam line break (MSLB). Among these, the passive containment cooling system (PCCS) is a safety-related system designed to prevent the containment from exceeding the design limits of both pressure and temperature following a postulated design basis accident, by transferring heat from the steel containment vessel to the atmosphere. During LOCA, natural air convection plays a vital role in removing heat from the steel containment vessel to the atmosphere. This is integrated with the condensate film from the Passive Containment Cooling Water Storage Tank (PCCWST). Air natural convection is also the only heat removal mechanism after the dry-out of the PCCWST and during the loss of shutdown decay heat removal (LOSDHR) event. Therefore, it is essential to investigate the thermal hydraulics behavior of the containment under transient conditions to ensure its safety and integrity. This paper presents a simplified CFD/ANSYS FLUENT model developed to analyze the impact of different parameters, such as air relative humidity, air temperature, and steel containment temperature, on the natural convection capability to remove the decay heat following LOCA. The model aims to examine the thermal behavior of the containment and provides a comprehensive understanding of the system's natural convection capability during postulated accidents. The results obtained from the model can be used to improve the safety and integrity of the containment system and provide valuable insights for future research in the field.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

$\overline{TKE}$ Parameterization and $\bar{v}$ Uncertainty Analysis for CGMF

Previous work was performed on tuning CGMF parameters for 235 U, 238 U, and Plutonium isotopes. Now work is being done to tune minor uranium isotopes. However, uranium isotopes like 232 U and 236 U have almost no experimental data. We are applying cross-isotope models to extrapolate and tune CGMF on isotopes that lack experimental data. There exist several internal CGMF physics quantities that affect the output of CGMF—multi-chance fission probability, excitation energy sharing, spin-cutoff factor, spin scaling, and fragment total kinetic energy to name a few. The mean fragment total kinetic energy, $\overline{TKE}$, is particularly interesting because of its strong anti-correlation with $\bar{v}$. We are most interested in the mean fragment total kinetic energy before neutron emissions. $\overline{TKE}$ is assumed to be pre-neutron emission unless otherwise stated. Currently in CGMF, the $\overline{TKE}$ model for 233,234,235,238 U are tuned independently to reproduce ν for the associated isotopes. In this report, we will tune a cross-isotope $\overline{TKE}$ model to experimental $\overline{TKE}$ data for 232,233,234,235,236,238 U. Because of the unreliable and sparse nature of $\overline{TKE}$ experimental data, future work will use more reliable experimental $\bar{v}$ data to infer the $\overline{TKE}$ model (and likely other internal CGMF parameters) for uranium isotopes. Such work has been performed previously using a sensitivity analysis and Kalman filter methods.

07 ISOTOPE AND RADIATION SOURCES↗

Restructuring and Optimizing Reactor Building Lifting Processes & Designing a Solution for an Overhead Door Handling Forklift Attachment

This poster presents two innovative projects aimed at enhancing operational efficiency and safety at the Advanced Test Reactor (ATR) Complex. The first project focuses on restructuring the existing lift book used for hoisting operations within the ATR Reactor Building. The current lift book, a cumbersome 150+ page document, is being transformed into a more efficient format using charts that allow for quick identification of maximum lift heights based on object footprint and weight. This new method also considers various parameters such as floor integrity and reactor status, dividing entries into four distinct charts to improve usability and safety during lifts. The second project involves designing a specialized forklift attachment for handling an overhead door at the ATR Reactor Building. The attachment is engineered to securely lift and lower a 1200lb, 14ft long door, ensuring safe replacement operations. The design process included research on forklift attachment codes, modeling in Autodesk Inventor, and testing through Finite Element Analysis (FEA) and hand calculations to confirm the attachment's structural integrity. Future work includes completing detailed drawings and an Engineering Calculations Analysis and Review (ECAR) document to proceed with manufacturing and assembly. Together, these projects demonstrate a commitment to optimizing reactor building operations through innovative engineering solutions and safety considerations.

42 - ENGINEERING↗

Dose Consequence and Probabilistic Risk Assessment Integration into Digital Documented Safety Analysis

This is an intern poster presentation. The current reactor authorization process is complex and error prone. The development of a digital Documented Safety Analysis has been proposed to provide an automated and integrated solution to enhance the design and authorization process. The digital DSA will consist of interlinked models, analyses, and reports, all of which will be updated using automated workflows when design changes are made. This poster examines the integration of the probabilistic risk assessment (PRA) with the transient and dose consequence analyses. A PRA for a generic high temperature gas-cooled reactor (HTGR) has been constructed which will drive the input parameters for a transient analysis model currently being constructed. Dose consequence will be determined using the results of the transient analysis and the Radiation Safety Analysis Computer (RSAC) code. Dose consequence data will then be input back into the PRA to drive design parameters.

42 ENGINEERING↗

Pilot-Scale Validation of Distributed Optical Fiber Sensors for Underground Pipeline Monitoring

Distributed fiber optic sensing is a cutting-edge technology that has found extensive applications in the monitoring of Ensuring the safety, integrity, and operational efficiency of underground product pipelines is vital for maintaining the nation’s critical infrastructure. Monitoring parameters such as hoop strain, pressure, and acoustic vibrations is key to detecting potential leaks, intrusions, or structural issues. Distributed optical fiber sensor (DOFS) systems provide a compelling solution for continuous, real-time monitoring over long distances. This paper details the development and pilot-scale implementation of DOFS systems for underground pipeline monitoring, evolving from a proof-of-concept stage. Multiple custom-designed DOFS interrogator units—such as optical frequency-domain reflectometry (OFDR), Brillouin optical time-domain analysis (BOTDA), and multimodal interferometer-based fiber acoustic sensors—were employed to measure key parameters like hoop strain, pressure, and acoustic vibrations. The underground product pipeline's outer diameter is 30 inches, the wall thickness is 1.28 inches, and the 3-foot depth. The fiber deployment strategies, and sensing data acquisition methods for these systems are discussed. The results demonstrate the effectiveness of DOFS in detecting hoop strain, temperature changes, and acoustic vibrations, showcasing their potential for real-time monitoring and enhancing pipeline safety.

distributed fiber sensing↗

A review on the sizing and selection of control valves for thermal hydraulics for reactor system applications

This study presents an overview regarding how to model/size, choose, and maintain control valves (CVs) for a reactor system thermal hydraulics experimentation—specifically simulating pipe break(s) and a loss of coolant accident (LOCA) analysis. In a thermal hydraulics test loop, CVs modulate fluid flow according to the control signal generated by the controller to control the process parameter. The selection of a CV in industrial applications involves many factors, including process diversity, safety, and reliability. Further, the key issues and challenges for performing proper mathematical modeling, installation, servicing, and calibration are also discussed for the benefit of operators and end-users. Additionally, this research includes case studies illustrating major plant-level accidents associated with the malfunctioning of valves. Therefore, the study will serve as a knowledge base or reference guide for young professionals and operators who wish to better understand process control applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

HFIR LEU High Density Silicide Dispersion Optimized Design Steady-State Heat Transfer Analyses

Steady-state heat transfer simulations of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design were performed to support comprehensive performance and safety metric studies concerning this design. The LEU Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium (HEU) core performance level at 85 MW. Full cycle Mode 1 full flow Case 1 (inlet temperature), Case 2 (flux-to-flow), and Case 3 (inlet pressure) safety limit analyses were performed to assess the margins to critical heat flux. Under the prescribed conditions, this LEU design meets the safety limit and limiting control setting requirements outlined in HFIR’s documented safety analysis; however, the safety margins are less than those for the 85 MW HEU core, and several assumptions were made where fuel fabrication and qualification data are currently lacking for the silicide fuel design. Effects of changes to pertinent fuel fabrication assumptions and uncertainty factors on thermal safety margins were also evaluated, showing that the margins are sensitive to many of these parameters. Power and pressure perturbations were also performed, indicating that significant steady-state thermal margins could be gained by increasing the coolant inlet pressure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

Empirical Bounding Analysis and User Recommendations for a Neutron Multiplicity Detector

Neutron multiplicity detectors are useful for a variety of applications including nuclear emergency response, nuclear nonproliferation, safeguards, and criticality safety. When measuring black-box problems (i.e., when the system being measured is completely unknown) with systems that have relatively low detection efficiency, expert analysis is frequently required to determine appropriate bounds for system parameters such as neutron multiplication and neutron leakage. This is because the detection efficiency can vary wildly, which also means that the resulting system parameters can have large variations. This work applies an empirical approach to provide bounds on system parameters (such as system multiplication, neutron leakage, and detector efficiency). The results of this bounding analysis are then used as inputs to provide recommendations to users regarding criticality safety. In addition to describing the method used, this work provides sample results from a measurement campaign performed at the National Criticality Experiments Research Center (NCERC).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analyzing Potential Failures and Effects in a Pilot-Scale Biomass Preprocessing Facility for Improved Reliability

This study demonstrates a failure identification methodology applied to a preprocessing facility generating conversion-ready feedstocks from biomass meeting conversion process critical quality attribute (CQA) specifications. Failure Modes and Effects Analysis (FMEA) was used as an industrially relevant risk analysis approach to evaluate a logging residue preprocessing system to prepare feedstock for pyrolysis conversion. Risk evaluations considered both system-level and operation unit-level assessments considering process efficiency, product quality, cost, sustainability, and safety. Key outputs included estimations of semi-quantitative risk scores for each failure, identification of the failure impacts, identification of failure causes associated with material attributes and process parameters, ranking success rates of failure detection methods, and speculation of potential mitigation strategies for decreasing failure risk scores. Results showed that deviations from moisture specifications had cascading consequences for other CQAs along with process safety implications. Failures linked to fixed carbon specifications carried the highest risk scores for product quality and process efficiency impacts. As increased throughput can be inversely related to meeting product quality specifications; achieving throughput and other material-based CQAs simultaneously will likely require system optimization or prioritization based on system economics. Ultimately, this work successfully demonstrates FMEA as a risk analysis approach for other bioenergy process systems.

09 BIOMASS FUELS↗