Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter compression”

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

Diagnostics of Magnetohydrodynamic Modes in the Interstellar Medium through Synchrotron Polarization Statistics

One of the biggest challenges in understanding magnetohydrodynamic (MHD) turbulence is identifying the plasma mode components from observational data. Previous studies on synchrotron polarization from the interstellar medium (ISM) suggest that the dominant MHD modes can be identified via statistics of Stokes parameters, which would be crucial for studying various ISM processes such as the scattering and acceleration of cosmic rays, star formation, and dynamo. In this paper, we present a numerical study of the synchrotron polarization analysis (SPA) method through systematic investigation of the statistical properties of the Stokes parameters. We derive the theoretical basis for our method from the fundamental statistics of MHD turbulence, recognizing that the projection of the MHD modes allows us to identify the modes dominating the energy fraction from synchrotron observations. Based on the discovery, we revise the SPA method using synthetic synchrotron polarization observations obtained from 3D ideal MHD simulations with a wide range of plasma parameters and driving mechanisms, and present a modified recipe for mode identification. We propose a classification criterion based on a new SPA+ fitting procedure, which allows us to distinguish between Alfvén mode and compressible/slow mode dominated turbulence. We further propose a new method to identify fast modes by analyzing the asymmetry of the SPA+ signature and establish a new asymmetry parameter to detect the presence of fast mode turbulence. Additionally, we confirm through numerical tests that the identification of the compressible and fast modes is not affected by Faraday rotation in both the emitting plasma and the foreground.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Anelasticity and Phase Transition During Ramp-Release in Tin

This study examines the qualitative features of an anelasticity model associated with the bowing of dislocations in the presence of phase transition. A simple physically plausible mechanism is introduced to describe the interaction of anelasticity and the transformation. Varying the anelastic parameters results in strong differences in the deviatoric stress response. The model is applied to study the behavior of tin (Sn) and compared to data from ramp driven compression-release experiments. Tin exhibits a complex phase diagram within a relatively accessible range of temperature and pressures and the characterization of its phases is considered an open problem with significant scientific merit. The coupling between anelasticity, plasticity, and phase transformation contributes to release wave features traditionally associated with the phase transition effect alone suggesting the importance of accounting for the effects jointly. Posterior distributions of the plastic and anelastic parameters are computed using Bayesian-inference-based methods, further highlighting the importance of anelasticity in this regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-stage heat release of multi-component fuels: Insights and implications for advanced engine operation

Multi-stage heat release (MSHR) is a unique phenomenon typically seen in lean/diluted and low- to intermediate-temperature combustion. Despite its relevance to advanced engine operation, the MSHR of multi-component fuels has barely been quantified. This study aims to characterize the MSHR of multi-component fuels in a rapid compression machine (RCM) at conditions representative of advanced combustion engines. New experimental data are first reported in the RCM at an equivalence ratio of 0.4, pressure of 60 bar and temperatures from 702 to 795 K for a research grade, multi-component gasoline surrogate, termed PACE-20. Here, experiments confirm the existence of MSHR for PACE-20 at all temperatures, and reveal the strong inhibiting effect of temperature on MSHR, where increasing temperature inhibits MSHR by suppressing first stage heat release and promoting second and third stage heat release. A detailed chemical kinetic model is also adopted to model the experiments, with good agreement observed. The response of MSHR characteristics to changes in different engine operating parameters (i.e., temperature, pressure, equivalence ratio and CO 2 dilution level) and fuel compositions (i.e., mole fraction of n-pentane, n-heptane, isooctane, cyclopentane, 1-hexene, ethanol and aromatics in PACE-20) is further evaluated via statistical analysis coupling quasi-random sampling, extensive computer experiments and high-dimensional model representation. The change in MSHR is characterized through 8 quantities of interest (QOIs), i.e., duration and extent of each heat release stage, and their standard deviations. The analysis highlights the dependence of the QOIs on the individual parameters as well as their interplays, with temperature exhibiting the strongest impact among all parameters. Furthermore, it is demonstrated that if MSHR is to be facilitated and combustion phasing is to be extended to enable engine operation at higher compression ratios, it is recommended to operate the engine at low intake temperatures, boosted intake pressures and high CO 2 dilutions (i.e., high EGR levels) with gasoline fuels containing more n-pentane and less aromatics.

42 ENGINEERING↗
DOE OSTI2024

Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales

The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.

58 GEOSCIENCES↗
DOE OSTI2020

Localization of deformation in fluid-saturated sandstone

Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.

42 ENGINEERING↗
DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗
DOE OSTI2024

Unraveling the structural dependency of Weyl nodes in Co 2 MnGa

Conventionally, the modulation of the intrinsic Weyl nodes in Weyl semimetals is challenging, due to topological protection. Here we report the structural dependence of the Weyl nodes in a Co2⁢MnGa Heusler thin film via a temperature-dependent tetragonal distortion. The ability to manipulate these Weyl nodes allows for the control of the intrinsic electromagnetic properties. Temperature-dependent x-ray diffraction (XRD) measurements identify a compressive tetragonal distortion with decreasing temperature from 300 to 20 K. The calculated Weyl properties can be directly compared with experimental parameters through the temperature-dependent XRD measurements which show the intrinsic correlation between Weyl properties and important magnetic parameters. Further, the microscopic momentum space properties of Weyl nodes such as the distance (d W ), solid angle (Ω W ), tilt (φ W ), and nodal point energy (E W ) directly affect the macroscopic observable properties such as exchange stiffness (A), magnetization (M), and effective anisotropy field $H^{eff}_{K}$, as shown via structure-dependent density functional theory calculations. These predictions are experimentally observed as large variations in the bulk magnetization and effective anisotropy field as a function of temperature. These results highlight a unique degree of freedom in the control of macroscopic magnetic properties via the modulation of the intrinsic properties of Weyl nodes through structural distortions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Material extrusion with integrated compression molding of NdFeB/SmFeN nylon bonded magnets using small- and large-scale pellet-based 3D-printers

High-density bonded rare-earth magnets are manufactured using pellet-fed additive manufacturing (AM)/material extrusion and an integrated additive manufacturing-compression molding (AM-CM) process. Neodymium iron boron – samarium iron nitride in polyamide 12 (NdFeB-SmFeN/PA12) of 93 % weight fraction (65 % volume fraction) are used for the study. The mechanical properties (tensile strength and modulus), magnetic properties (maximum energy density, coercivity, remanence) are reported. Manufacturing parameters such as layer height, barrel temperatures, screw speed and gantry feed rate are optimized to obtain the highest possible density of the magnets using a small-scale desktop material extrusion printer. Large scale integrated additive manufacturing-compression molding (AM-CM) is then utilized to increase the density of the magnets by reducing porosity defects common in the material extrusion process. The density of as-printed magnets was 5.2 g/cm 3 with a BH max value of 124.14 kJ/m 3 , tensile strength of 20 MPa and a modulus of 2 GPa. AM-CM increased the density of the compound by 5.5 % (5.49 g/cm 3 ). The reduction in porosity was confirmed using X-ray tomography (XCT). Improvement in mechanical strength of the material was also observed, with an increase in tensile strength of 25 % (25.09 MPa) and increase in tensile modulus of 275 % (5.49 GPa). Scanning electron microscopy showed increased particle-matrix adhesion with the integrated AM-CM process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2023

Investigation of the High-Pressure Behaviors of Amblygonite by Single-Crystal X-ray Diffraction, Raman Spectroscopy, and DFT Calculations

In the present study, we extensively explored the high-pressure behaviors and vibrational properties of amblygonite LiAlPO 4 F with elevated pressures up to 34.3 GPa based on single-crystal X-ray diffraction measurements, Raman spectroscopy, and DFT calculations. The compressibility and elastic properties of amblygonite were determined first. Specifically, the obtained isothermal bulk modulus of LiAlPO 4 F is 128(4) GPa and the triclinic phase exhibited anisotropic compression with axial compressibility β c > β a > β b with a ratio of 1.11:1.00:1.20. The Raman spectra showed no indication of phase transformation and were used to obtained mode Grüneisen parameters. The average Grüneisen parameter for PO 4 tetrahedral sites was smaller than for the LiO 4 F sites. Our results provide new insights into the phase stability and elastic properties of lithium-fluorite granites at extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Compressive properties of silicone Bouligand structures

Abstract This article presents an examination of silicone, Bouligand lattices in compression. Appearing frequently in biological organisms and manufacturing design, Bouligand structures comprise layers of parallel strands or fibers oriented in a helical fashion. They can exhibit exceptional fracture resistance when composed of rigid or composite materials. The behavior of elastomeric Bouligand structures, however, is less well understood. Additively manufactured (AM) elastomeric lattices have applications in stress mitigation, medical devices, and soft robotics. This article demonstrates that Bouligand structures are a useful addition to the design space of AM elastomers. By adjusting the layer-rotation parameters, lattice stress can increase by more than 300% without altering the porosity. Additionally, we introduce path length metrics that help explain the observed relationship between layer rotation and compression response. Impact statement Additive manufacturing (AM) continues to push the boundary of manufacturable structures and enhance the ability to robustly design for specific properties and behaviors. The more we understand the design space of a novel AM microstructure, the greater its application range. In this article, we describe the mechanical behavior of helicoidal, elastomeric lattices and introduce path length metrics to help explain their stress response. We show that these structures can exhibit a large range of mechanical behaviors in compression, making them well suited for applications such as stress mitigation and impact absorption. Additionally, the path length metrics could become useful design tools and may be applicable to a larger set of cellular structures. These findings expand our ability to rapidly design materials with highly specific and customizable properties to meet the needs of modern engineering challenges. Graphical abstract

Van Meerbeek, I. M. (ORCID:000000021258132X)↗
DOE OSTI2026

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Cascade of transitions in twisted and non-twisted graphene layers within the van Hove scenario

Abstract Motivated by measurements of compressibility and STM spectra in twisted bilayer graphene, we analyze the pattern of symmetry breaking for itinerant fermions near a van Hove singularity. Making use of an approximate SU(4) symmetry of the Landau functional, we show that the structure of the spin/isospin order parameter changes with increasing filling via a cascade of transitions. We compute the feedback from different spin/isospin orders on fermions and argue that each order splits the initially 4-fold degenerate van Hove peak in a particular fashion, consistent with the STM data and compressibility measurements, providing a unified interpretation of the cascade of transitions in twisted bilayer graphene. Our results follow from a generic analysis of an SU(4)-symmetric Landau functional and are valid beyond a specific underlying fermionic model. We argue that an analogous van Hove scenario explains the cascade of phase transitions in non-twisted Bernal bilayer and rhombohedral trilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

LDA + $U_{sc}$ calculations of phase relations in FeO

Using the LDA + $U_{sc}$ method, we present calculation phase relations of iron monoxides involving five polytypes in multiple spin-state configurations. In this work, the Hubbard parameter $\textit{U}$ is determined self-consistently simultaneously with the occupation matrix and structures at arbitrary pressures. The Hubbard parameter strongly depends on pressure, structure, and spin state. Comparison with experimental structural data indicates the LDA + $U_{sc}$ can predict structure, compression curves, phase relations, and transition pressures very well for the insulating $\textit{B}$1 and $\textit{iB}$8 states. However, it requires additional calculations using the Mermin functional that includes the electronic entropic contribution to the free energy to obtain an $\textit{nB}$8 metallic state and a consistent $\textit{iB}$8 to $\textit{nB}$8 insulator to metal transition pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multigroup Thermal Radiation Transport with Tensor Trains

We investigate the application of tensor-train (TT) algorithms to multigroup thermal radiation transport (i.e., photon radiation transport). The TT framework enables simulations at discretizations that might otherwise be computationally infeasible on conventional hardware. We show that solutions to certain multigroup problems possess an intrinsic low-rank structure, which the TT representation leverages effectively. This enables us to solve problems where the discretized solution size exceeds a trillion parameters on a single node. The solver is evaluated on a range of test problems with varying levels of complexity, consistently achieving compression factors greater than 100× and speedups exceeding 2×. We also investigate alternative TT topologies by analyzing the low-rank structure of the merged spatio-spectral core to assess the potential for greater compression. This analysis suggests that compression gains could increase by factors as large as 7. Our results indicate that the low-rank structure of the merged spatio-spectral core captures the spatio-spectral complexity of the solution, largely driven by the opacity structure of the medium. Beyond identifying opportunities for improved compression, this analysis highlights the types of errors that may arise in angle-integrated quantities when exploiting this low-rank structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Mechanical Characterization of Loblolly Pine

Biorefineries are faced with various material flow challenges during feeding and handling, due to inherent variability and inhomogeneity of biomass feedstocks. This impacts the efficiency of refinery processes and unit operations. Proper characterization of bulk properties allows for optimal process design and reduction of overall costs. This work is aimed at evaluating the physical and mechanical properties that impact bulk flow of loblolly pine. Cyclic axial compression tests are performed to elucidate consolidation behavior under different stress conditions. The results are then used to develop correlations between material parameters (particle size, tissue type, moisture content) and Elastic Modulus of bulk samples.

09 BIOMASS FUELS↗