Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete element method (DEM)”

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

Modeling powder spreadability in powder-based processes using the discrete element method

Powder-bed fusion (PBF) processes refer to a subset of Additive Manufacturing (AM) techniques where powder is spread on the build-plate before melting (by a laser or electron beam). While PBF processes are attractive due to their ability for realizing complex structures that are either difficult or impossible to create through conventional means, the parts fabricated with these techniques can exhibit defects such as pores, inclusions, and excessive surface roughness. To minimize these defects, much research has been dedicated towards process maturation by optimizing laser or electron beam parameters. However, these developmental efforts typically do not address the recoating process where achieving dense and uniform layers of powder is a necessity for ensuring process repeatability and part quality. While the recoating process can be studied through experimentation, the dynamics of particle movement are difficult to analyze experimentally. Therefore, here, in this study, powder spreading in PBF was simulated through the Discrete Element Method (DEM) to elucidate the mechanisms that control powder-bed quality. Utilizing the Buckingham Pi theorem, a dimensionless metric referred to as the spreading index is developed that combines powder-bed density, roughness, and particle size to assess the quality of powder layers. The formulated spreading index is then related to several dimensionless quantities that provide insight into the mechanisms dominating powder spreading in PBF. The DEM simulations conducted in this work focused on the scenario where powder is spread onto an existing powder bed and revealed that a reduction in the recoating velocity causes an increase in the spreading index while little to no impact on the spreading index was observed when varying layer thickness from 30 μm to 75 μm.Particle size effects on the powder-bed quality were also investigated.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
NASA NTRSDate not supplied

Monodisperse Single-Material Granular Tribocharging Modeling and Experimental Validation

Charge transfer between insulating grains has been a topic of interest for many years as this phenomenon is extremely important to many areas in industry. Of particular interest to NASA is the behavior of electrically insulative dust grains such as those found on the lunar surface. Whether poured from a scoop during sample collection activities, agitated inside a drum via mining robots, or fluidized by a gas plume, particle-particle interactions between similar granular materials will be widespread on the Moon. Designers hope to have a predictive model for how dust grains are charged, transported, and deposited on flight systems to better understand issues electrostatics may cause in future missions. While difference in work function is often seen as a driver for charge transfer between two materials, this is an incomplete representation for insulative particles. Other asymmetries in the system can also lead to charge being transferred: one such being the difference between static vs dynamic particles. Experiments performed under vacuum in the Electrostatics and Surface Physics Laboratory (ESPL) clearly show that particles having more contacts with other particles tend to charge positively and particles having fewer contacts tend to charge negatively, despite being made from the same material. During these experiments, monodisperse (population of single diameter) spherical particles were held in a reservoir which was then elevated at one end to slope into a Faraday cup used to measure bulk charge of the grains. Care was taken to limit the number of spheres bouncing out of the cup, to isolate the beads from the container itself so that only grain-to-grain interactions occurred, and to assure neutrality of the beads before the experiment began. A discrete element method (DEM) modeling package already incorporating many granular mechanics interactions was augmented to include the electrostatic behaviors seen in experiments. The model can reproduce the experimental results above by tracking the number of contacts between pairs of individual grains and applying a charge transfer condition related to this contact parameter. The experimental results as well as their incorporation into the improved modeling suite will be discussed. Future work includes adding additional asymmetries into the experiments/model such as polydisperse populations of spherical particles and irregularly shaped grains.

granular↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗
DOE OSTI2022

Very small-scale, segregating-fluidized-bed experiments: A dataset for CFD-DEM validation and uncertainty quantification

IWe report fluidization experiments were conducted on a small scale and with a rapid response (short duration) to enable corresponding simulations at low-computational cost. Rise times are reported for four or fewer polyethylene particles (intruders) in an air-fluidized bed of ~5000 group D glass beads. Experimental inputs were completely characterized—particle properties, system dimensions and operating conditions—which is necessary for validating computational fluid mechanics (CFD)-discrete element method (DEM) including a comprehensive uncertainty quantification (UQ) analysis. Input uncertainties are reported as bounds or cumulative distribution functions of measured values. The staggering number of simulations required to complete a UQ analysis (~O[10 4 ] simulations corresponding to ~5 uncertain inputs) motivates this study. These segregating-bed experiments are designed to permit analogous CFD-DEM simulations to complete in less than a day on a single (~2.5 GHz) computational processor unit (CPU). Segregation times are reported for several operating conditions, intruder sizes, and initial configurations, providing a rich dataset for numerical model testing, validation and UQ.

42 ENGINEERING↗
NASA NTRSDate not supplied

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗
DOE OSTI2025

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗
NASA NTRS2022

Investigating Particle-Particle Electrostatic Effects on Charged Lunar Dust Transport via Discrete Element Modeling

NASA surface exploration missions have always seen negative effects of dust including the Apollo missions. The astronaut-witnessed unusual behavior of the dust particles that surround the vehicle after engine cutoff has the potential to have more of an influence on surface systems dust loading than the high velocity lunar rocket plume ejecta in the landing process. The levitation and transport of the fine components of regolith on lunar surface has been linked to electrostatic effects and electric field, but so far there is no accurate model considering the inter-particle electrostatic interactions, especially when the particles are charged by rocket plume or other mechanical interactions due to exploration activities. This study is proposed to investigate the dynamics of charged lunar regolith with a discrete element method (DEM) approach focusing on the inter-particle interactions and contact charge transfer. The grain dynamics is coupled with mechanical and electrical particle interactions, and both short- and long-range interactions between spherical particles are incorporated. A tribo-charging model based on instantaneous collisions between particles is adopted and validated by comparing the simulation results to existing experimental data. Sensitivity analysis is conducted to quantify the effects of initial charge, tribo-charging, and E-field on transport of lunar dust based on JSC-1 simulants with a radius of 50 lm. DEM simulations are also conducted in a near realistic lunar environment with the estimations of initial conditions that shows the difference in position and velocity distributions between charged particles and uncharged particles. The results indicate that the charged dust particles have higher dispersion of position and velocity by several orders of magnitude due to electrostatic effects. This provides a potential explanation for the phenomena of the approximately 30 s dust lofting following Apollo Lunar Module landing.

Lunar dust↗
DOE OSTI2020

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

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

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CFD-DEM and PR-DNS studies of low-temperature densely packed beds

Over the past few decades, granular media is gaining attention as a viable option for heat transfer fluids (HTFs). Several research efforts are studying the use of particle-based heat transfer fluids in a wide variety of applications. With this motivation, the current work focusses on analyzing the different heat transfer mechanisms in low-temperature mono-sized densely packed granular media. To study the heat transfer behavior of granular media at different scales, the current work employs a two-way coupled computational strategy. The motion of particles is solved using the Discrete Element Method (DEM) and the interstitial air is solved using a Finite-Volume (CFD) approach. The Open-Source library CFDEM Coupling® is used in the current study to join the Finite Volume PISO solver of OpenFOAM® and the DEM solver of LIGGGHTS®. Typically, particle-particle contact conduction and particle-air convection are the most popular closure models. But recent research identified a different heat transfer phenomenon in packed beds that cannot be identified by conduction or convection models. While closure models were developed to implement this on a CFD-DEM framework, they did not capture the effect of intra-particulate thermal gradients on this phenomenon. Hence the current work also employs Particle-Resolved Direct Numerical Simulations (PR-DNS) to gain valuable insights allowing for the modification of existing models. A new closure model is then proposed here and is implemented in the CFD-DEM framework. This model provides key insights into the different heat transfer mechanism of packed beds.

42 ENGINEERING↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
NASA NTRS2010

Charge-Spot Model for Electrostatic Forces in Simulation of Fine Particulates

The charge-spot technique for modeling the static electric forces acting between charged fine particles entails treating electric charges on individual particles as small sets of discrete point charges, located near their surfaces. This is in contrast to existing models, which assume a single charge per particle. The charge-spot technique more accurately describes the forces, torques, and moments that act on triboelectrically charged particles, especially image-charge forces acting near conducting surfaces. The discrete element method (DEM) simulation uses a truncation range to limit the number of near-neighbor charge spots via a shifted and truncated potential Coulomb interaction. The model can be readily adapted to account for induced dipoles in uncharged particles (and thus dielectrophoretic forces) by allowing two charge spots of opposite signs to be created in response to an external electric field. To account for virtual overlap during contacts, the model can be set to automatically scale down the effective charge in proportion to the amount of virtual overlap of the charge spots. This can be accomplished by mimicking the behavior of two real overlapping spherical charge clouds, or with other approximate forms. The charge-spot method much more closely resembles real non-uniform surface charge distributions that result from tribocharging than simpler approaches, which just assign a single total charge to a particle. With the charge-spot model, a single particle may have a zero net charge, but still have both positive and negative charge spots, which could produce substantial forces on the particle when it is close to other charges, when it is in an external electric field, or when near a conducting surface. Since the charge-spot model can contain any number of charges per particle, can be used with only one or two charge spots per particle for simulating charging from solar wind bombardment, or with several charge spots for simulating triboelectric charging. Adhesive image-charge forces acting on charged particles touching conducting surfaces can be up to 50 times stronger if the charge is located in discrete spots on the particle surface instead of being distributed uniformly over the surface of the particle, as is assumed by most other models. Besides being useful in modeling particulates in space and distant objects, this modeling technique is useful for electrophotography (used in copiers) and in simulating the effects of static charge in the pulmonary delivery of fine dry powders.

Walton, Otis R.↗
DOE OSTI2021

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗
DOE OSTI2024

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. In conclusion, we present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

cohesive powders↗
DOE OSTI2023

An experiment-informed discrete element modelling study of knife milling for flexural biomass feedstocks

A discrete element method (DEM) based approach is used to study the relationships between material attributes (MAs), processing parameters (PPs), and quality attributes (QAs) for the knife milling of maize stalks. An approximate DEM shape model was conceptualized based on real maize stalks and calibrated based on experimental bending test data for flexural properties (elastic bending stiffness, elastic bending angle limit, elastoplastic ratio, etc.). DEM simulations of maize stalk comminution in a Jordan Reduction Solutions (“JRS”) knife mill were performed to investigate the relationships between the MAs (maize stalk size and breakage stress limit), PPs (impeller rotational speed), and QAs (mass throughput and output particle size distribution (PSD)). The DEM results suggest that stalk length has little influence on mass throughput and PSD, whilst stalks with larger cross sections tend to generate larger sizes of milled particles given the same breakage stress limit. Both the DEM and experimental results show that faster impeller rotation (or higher power) does not necessarily generate higher throughput or smaller output PSD, especially for maize stalks of higher breakage stress limit. The correlations between these MAs, PPs and QAs are found highly stochastic, though breakage stress limit dictates mass throughput, regardless of stalk size. The DEM-predicted output particle size tended to match the experimental data with coarse PSDs based on sieve size but showed weakened fidelity with finer material, indicating the potential for further model improvement.

09 BIOMASS FUELS↗