Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Correlation coefficient”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Estimation of micro-Hall-Petch coefficients for prismatic slip system in Mg-4Al as a function of grain boundary parameters

Grain size strengthening, referred to as the Hall-Petch effect, is a common strategy to improve the yield strength of magnesium (Mg) alloys. Several experimental studies have reported that the Hall-Petch slope strongly depends on the texture of the alloy. This effect arises from altering grain boundaries (GBs) resistance to different slip systems to transfer across adjacent grains. The grain boundary barrier strength of certain grain boundaries to basal slip, referred to as basal micro-Hall-Petch, was investigated in the previous work. Here in this study, the micro-Hall-Petch coefficient values for the prismatic slip $(k^{prismatic}_{μ})$ in Mg-4Al and their correlation with the grain boundary parameters were investigated. An experimental method was developed to initiate the prismatic slip band at low-stress levels. High-resolution electron backscatter diffraction (HR-EBSD) was used to measure the residual stress tensor, from which the resolved shear stress ahead of blocked prismatic slip bands was computed for seven different grain boundaries. $(k^{prismatic}_{μ})$ values for each individual GB were calculated by coupling the stress profile information with a continuum dislocation pile-up model. The $(k^{prismatic}_{μ})$ values vary from 0.138 MPa.m 1/2 to 0.685 MPa. m 1/2 which are almost three times larger than the calculated values for the basal micro-Hall-Petch. The $(k^{prismatic}_{μ})$ values were correlated with the GB parameters, and a functional relationship depending on the two most effective angles, the angle between the traces of the slip planes on the GB plane ($\phi$) and the angle between incoming and outgoing slip directions (κ), was proposed to estimate the Hall-Petch barrier for prismatic slip system. The work provides coefficients that can be supplied as input to crystal plasticity models to couple the effect of texture and grain size effectively.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ab Initio Insights to Metal Passivation: Diffusion and Defect Formation in Amorphous Zirconia and Alumina

Metal passivation refers to the formation of protective oxide films on metals, which shield them from further corrosion and oxidation, playing a crucial role in maintaining their stability. The mesoscopic Point Defect Model has successfully predicted passivity as a steady state process where oxide growth from oxygen vacancies at the metal/film interface competes with oxide dissolution at the film/environment interface. In this work, informed by the Point Defect Model parameters, we use first-principles calculations to calculate defect formation and atomic diffusion in amorphous materials and correlate these descriptors with the behavior and growth of the oxide film. Focusing on amorphous zirconia and alumina, we demonstrate that defect formation energies exhibit significant variability in amorphous systems. In alumina, vacancies dominate, with cation and anion vacancies occurring at comparable concentrations. Diffusion calculations for stoichiometric amorphous alumina and zirconia, as well as oxygen-deficient zirconia, reveal faster diffusion in the oxygen-deficient case, highlighting the impact of defects on transport. Comparison of calculated self-diffusion coefficients for the dominant defect species with experimentally measured oxide thicknesses shows a clear correlation, suggesting that first-principles-derived diffusivity information can serve as a key descriptor for surface passivation film growth.

Farnell, Mackinzie S. [University of California, B↗
DOE OSTI2020

Water Energy Resource Innovation on the Cavitation Characteristics

The main purpose of this study is to numerically correlate the amount of generated vapor over a hydrofoil to the lift and drag coefficients acting on it. As such, cavitation characteristics were investigated of a hydrofoil in the cavitating, sub-cavitating, and non-cavitating flows for different angles of attacks (AoA) with the high upstream flow velocity. The hydrofoil was tested in a square water tunnel with water entering the tunnel at various velocities for each AoA ranges from 9.1 m/s to 12.2 m/s. It was found that lift and drag forces acting on the hydrofoil follow the trend of the experimental data quite closely. While the cavitation can be identified by a unique number (averaged vapor volume fraction), the work done created an inverse correlation between this number and the cavitation number at the same angle of attack. The lift force declines with the increase in the vapor content on the hydrofoil surface, meanwhile the drag force peaks at certain vapor volume fraction, and then, a huge reduction occurs with the considerable decrease in the corresponding cavitation number. A fourth-order correlation generated between the lift to drag (L/D) and the cavitation number (σ). It was found the lift-to-drag ratio decreases by the formation of the cavitation over the hydrofoil, thus causing a drop in the efficiency of the turbomachines.

42 ENGINEERING↗
DOE OSTI2023

Film condensation with high heat fluxes and scaled experiments using pure steam for reactor containment cooling

Condensation tests were performed using a newly developed test facility for scaling the passive containment cooling system (PCCS) to a small modular reactor (SMR). The PCCS of the SMR plays a pivotal role in ensuring greater safety, reliability, and compactness than what is afforded by traditional reactors. Therefore, a well-designed PCCS is essential to SMRs. However, previous studies and test data were unsuitable for scaling, due to high variation in the test geometry and operating conditions. This study intends to close this research gap by using a novel designed scaled test facility consisting of vertical condensing test sections featuring 1-, 2-, and 4-inch-diameter condensing tubes with annular water cooling, and by applying superheated and saturated steam with different steam mass flow ranges of 5–25 g/s. Further, the primary test data, including axial temperatures, mass flow rates, and pressures, were used in conjunction with a standard data reduction method to estimate critical parameters such as heat fluxes, heat transfer coefficients, and condensation rates. These scaled test data would support improving empirical correlations and validating condensation models to identify scaling distortion for SMR PCCSs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Predicting roughness effects in additively manufactured coolant channels with helical enhancements

Additive manufacturing (AM) is a promising technique for fabrication of complex geometries such as those expected to be utilized in the blanket, first wall, and divertor. In the case of cooling, metallic AM may be exploited to embed geometric enhancements (ribs, rifling, etc.) to improve cooling performance. However, due to the roughness of these unfinished internal AM surfaces, prediction of thermal hydraulic performance in such channels is difficult. In this work, we consider a methodology for predicting pressure drop and heat transfer in AM channels containing helical enhancements (e.g. helical ribs, twisted tapes) that allows the incorporation of roughness data through conventional pipe flow correlations. This methodology is tested using experimental friction factor and heat transfer coefficient data from high-pressure helium coolant flow measurements in AM stainless steel tubes fabricated by laser powder bed fusion. Both a featureless AM tube and one containing helical ribs were considered alongside a conventionally manufactured smooth tube. The AM surface roughness is obtained by profilometry and used to predict an equivalent sand-grain roughness, with this equivalent roughness confirmed through AM featureless tube measurements. Under the proposed methodology, this roughness information is incorporated into predictions of friction factor and Nusselt number for the rifled tube. Furthermore, these predictions agree well with experimental data across a large range of Reynolds numbers, encouraging the use of this methodology for thermal hydraulic analysis of similar systems and design of future coolant channel geometries.

Additive manufacturing↗
DOE OSTI2024

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗
DOE OSTI2022

Characterizing how size distribution and concentration affect echogenicity of ultrasound contrast agents

Here, we investigated the effects of UCA gas bubble size distribution and concentration on the generated ultrasound echogenicity signal. Gas bubble size characterization using Coulter Counter and cryogenic-SEM revealed the hollow structure and rare presence of microbubbles >10 µm in a commercial UCA product, Lumason™. Volume-weighed size and concentration were observed to be more sensitive to changes in UCA bubble stability than number-weighted size and concentration. Size distribution measurements showed that the force (e.g., shaking/agitation energy) used to redisperse the sample did not affect the size distribution, concentration, or echogenicity of the UCA sample. The ultrasound backscattering coefficient (BSC) of size fractionated and serial diluted microbubbles showed that the echogenicity signal correlates most with UCA bubble concentration, especially volume-weighted concentration. Findings from this study may be used to support demonstrating the equivalence of a generic UCA product to the reference listed drug.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Measurement of Convective Heat Transfer Coefficients With Supercritical CO 2 Using the Wilson-Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in a smooth tube and compares the results with published correlations for validation. The paper also describes the Heat Exchange and Experimental Testing (HEET) rig recently designed and built at the U.S. Department of Energy’s (DoE’s) National Energy Technology Laboratory (NETL) in Morgantown, WV. Here, the Wilson-plot technique used for measuring the heat transfer coefficients is described along with the data reduction process. The Wilson-plot technique was chosen as the basis for the design of NETL’s HEET rig. Advantages of the Wilson-plot technique include the (1) ability to measure high convective heat transfer coefficients accurately, (2) ability to measure average heat transfer coefficient for complicated heat exchange geometries like those produced using additive manufacturing, (3) ability to measure heat transfer coefficients on both sides of a heat exchanger independently, and (4) simplicity of experimental setup. Capabilities of the HEET rig include pressure to 24 MPa (3500 psig), temperature to 538 °C (1000 °F), mass flow rate to 1.5 kg/s (3 lb/s), and Re to 500,000. The rig is designed to operate with pure CO 2 or a mixture of CO 2 and up to 10% N2 by volume to study the impact of gas mixtures typical of direct-fired sCO 2 power cycles on the convective heat transfer and pressure drop. Preliminary tests in the HEET rig were performed with smooth stainless-steel tube and pure CO 2 , and the results were compared with published correlations for Nusselt number (Nu) and friction factor. Over a Reynolds number (Re) range from 58,000 to 228,000, measured Nu was compared to predictions using the Dittus and Boelter equation (Kreith and Bohn, 1993, “Principles of Heat Transfer, West Publishing Company”) within 5% and measured friction factors were compared to predictions using the McAdams correlation (“McAdams, 1954, “Heat Transmission,” 3rd ed., McGraw Hill, New York)” for smooth tube to be within 5%.

42 ENGINEERING↗
DOE OSTI2024

Dihadron correlations in small- x DIS at NLO: transverse momentum dependent fragmentation

We compute the inclusive dihadron cross-section in Deep Inelastic Scattering at next-to-leading order (NLO) and small x in the Color Glass Condensate. We focus on the kinematic limit where the hadrons are produced at forward rapidities (in the direction of the virtual photon) and back-to-back in the transverse plane. Our calculation demonstrates that the coefficient of the Sudakov double logarithm for this process is –$\frac{αs}{2π}$ [C F + $\frac{Nc}{2}$] instead of –$\frac{αsNc}{4π}$ when back-to-back jets are measured in the final state. To preserve the universality of the Sudakov soft factor associated with the Weizsäcker-Williams transverse momentum dependent (TMD) gluon distribution, we promote the collinear fragmentation functions into TMD fragmentation functions. We then perform the resummation of the Sudakov logarithms through Collins-Soper-Sterman evolution of the TMD fragmentation functions and the Weizsäcker-Williams TMD gluon distribution. Finally, analytic expressions are obtained for the NLO coefficient functions in the $\overline{MS}$-scheme. These results pave the way towards numerically calculating dihadron correlations at small x at the future Electron-Ion Collider with full NLO accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Molecular origin of negative lithium transference in electrolytes with star-shaped multivalent anions

Large-scale molecular dynamics simulations illustrate that highly correlated cation–anion motion leads to negative t 0+ on the order of 1 in lithium electrolytes with star-shaped multivalent anions. Large multivalent anions have gained increasing attention for their potential to improve lithium transference in electrolytes. We employ large-scale molecular dynamics simulations based on the Onsager transport framework to investigate ion transport in a lithium electrolyte with star-shaped multivalent anions. The simulations show that t 0+, the cation transference number with respect to solvent velocity, is negative over a wide range of concentration. This is consistent with experimental data reported previously. The simulation-based Onsager transport coefficients reveal that the magnitudes of the cation–cation, anion–anion, and cation–anion correlations are comparable, a signature of highly correlated motion in the electrolyte. Examination of the cation solvation environment indicates the presence of strong cation–anion association across the entire concentration range, which leads to negative t 0+ on the order of −1. Both simulation and experiment also show that the maximum value of t 0+ reaches 0 when the cation concentration is c + = 0.4 M. This is the concentration at which the anions begin to spatially overlap, and lithium ions serve as dynamic linkers to balance cation–cation and cation–anion correlations. Our results provide molecular-level insights into the origin of transference in multivalent electrolytes.

Fang, Chao↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A physical basis for cosmological correlators from cuts

Significant progress has been made in our understanding of the analytic structure of FRW wavefunction coefficients, facilitated by the development of efficient algorithms to derive the differential equations they satisfy. Moreover, recent findings indicate that the twisted cohomology of the associated hyperplane arrangement defining FRW integrals overestimates the number of integrals required to define differential equations for the wave-function coefficient. We demonstrate that the associated dual cohomology is automatically organized in a way that is ideal for understanding and exploiting the cut/residue structure of FRW integrals. Utilizing this understanding, we develop a systematic approach to organize compatible sequential residues, which dictates the physical subspace of FRW integrals for any n -site, ℓ-loop graph. In particular, the physical subspace of tree-level FRW wavefunction coefficients is populated by differential forms associated to cuts/residues that factorize the integrand of the wavefunction coefficient into only flat space amplitudes. After demonstrating the validity of our construction using intersection theory, we develop simple graphical rules for cut tubings that enumerate the space of physical cuts and, consequently, differential forms without any calculation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Two-particle azimuthal correlations in γp interactions using pPb collisions at s NN = 8.16 TeV

The first measurements of the Fourier coefficients ( V n Δ ) of the azimuthal distributions of charged hadrons emitted from photon-proton (γp) interactions are presented. The data are extracted from 68.8 nb −1 of ultra-peripheral proton-lead (pPb) collisions at s NN = 8.16 TeV using the CMS detector. The high energy lead ions produce a flux of photons that can interact with the oncoming proton. This γp system provides a set of unique initial conditions with multiplicity lower than in photon-lead collisions but comparable to recent electron-positron and electron-proton data. The V n Δ coefficients are presented in ranges of event multiplicity and transverse momentum ( p T ) and are compared to corresponding hadronic minimum bias pPb results. For a given multiplicity range, the mean p T of charged particles is smaller in γp than in pPb collisions. For both the γp and pPb samples, V 1 Δ is negative, V 2 Δ is positive, and V 3 Δ consistent with 0. For each multiplicity and p T range, V 2 Δ is larger for γp events. The γp data are consistent with model predictions that have no collective effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The Cosmological Bootstrap: Spinning correlators from symmetries and factorization

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model

The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.

DIII-D↗
DOE OSTI2022

Predictive models enhance feedstock quality of corn stover via air classification

Feedstock heterogeneity is a fundamental obstacle to cost-competitive biobased products. Agricultural products like corn stover have anatomical components that vary in their chemical composition, mechanical properties, structure, and response to chemical and biological treatments. A technique that can enrich streams in select anatomical fractions would allow a tailored deconstruction approach to increase overall process efficiency. Air classification can be leveraged for such refining, however, fundamental characterization and understanding of the particle properties that underly the physics of air classification are only modestly documented. Here, we determine fundamental particle properties including mass-to-area ratio, drag coefficient, and partition velocity that describe how anatomical tissues of corn stover behave during air classification. In this work, mass-to-area ratios of anatomical tissues vary by nearly two orders of magnitude from 2.3 mg/mm 2 for cob to 0.04 mg/mm 2 for leaf. Drag coefficients of longer, fibrous materials (i.e., rind, husk, and sheath) are shown to correlate with particle area (p-value < 0.001) whereas granular tissues (i.e., cob, pith, and leaf) correlate better with mass-to-area ratio (p-values < 0.001). When compared to experimental observations, a simulated two-stage air classification and size reduction scenario predicts the overall partitioning of anatomical tissues within 15% for pith, husk, rind, and cob tissues. The model predicts an air-classified fraction preferentially enriched in cob (purity = 20%), rind (purity = 74%), and pith (purity = 4.5%) with a mass yield of 47%. Empirical relations for these properties can be used to predict the partitioning of corn stover during air classification based on anatomical type and size.

09 BIOMASS FUELS↗