Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Descriptor: High Temporal Resolution Meteorological Data at Oak Ridge Reservation (ORR-HiResMet)

Access to continuous, quality assessed meteorological data is critical for understanding the climatology and atmospheric dynamics of a region. Research facilities like Oak Ridge National Laboratory (ORNL) rely on such data to assess site-specific climatology, model potential emissions, establish safety baselines, and prepare for emergency scenarios. To meet these needs, on-site towers at ORNL collect meteorological data at 15-minute and hourly intervals. However, data measurements from meteorological towers are affected by sensor sensitivity, degradation, lightning strikes, power fluctuations, glitching, and sensor failures, all of which can affect data quality. To address these challenges, we conducted a comprehensive quality assessment and processing of five years of meteorological data collected from ORNL at 15-minute intervals, including measurements of temperature, pressure, humidity, wind, and solar radiation. The time series of each variable was pre-processed and gap-filled using established meteorological data collection and cleaning techniques, i.e., the time series were subjected to structural standardization, data integrity testing, automated and manual outlier detection, and gap-filling. The data product and highly generalizable processing workflow developed in Python Jupyter notebooks are publicly accessible online. As a key contribution of this study, the evaluated 5-year data will be used to train atmospheric dispersion models that simulate dispersion dynamics across the complex ridge-and-valley topography of the Oak Ridge Reservation in East Tennessee.

Steckler, Morgan R. [Oak Ridge National Laboratory↗
DOE OSTI2022

Artificial neural network approach for multiphase segmentation of battery electrode nano-CT images

The segmentation of tomographic images of the battery electrode is a crucial processing step, which will have an additional impact on the results of material characterization and electrochemical simulation. However, manually labeling X-ray CT images (XCT) is time-consuming, and these XCT images are generally difficult to segment with histographical methods. We propose a deep learning approach with an asymmetrical depth encode-decoder convolutional neural network (CNN) for real-world battery material datasets. This network achieves high accuracy while requiring small amounts of labeled data and predicts a volume of billions voxel within few minutes. While applying supervised machine learning for segmenting real-world data, the ground truth is often absent. The results of segmentation are usually qualitatively justified by visual judgement. We try to unravel this fuzzy definition of segmentation quality by identifying the uncertainty due to the human bias diluted in the training data. Further CNN trainings using synthetic data show quantitative impact of such uncertainty on the determination of material’s properties. Nano-XCT datasets of various battery materials have been successfully segmented by training this neural network from scratch. We will also show that applying the transfer learning, which consists of reusing a well-trained network, can improve the accuracy of a similar dataset.

25 ENERGY STORAGE↗
DOE OSTI2021

Reactive boride infusion stabilizes Ni-rich cathodes for lithium-ion batteries

Engineered polycrystalline electrodes are critical to the cycling stability and safety of lithium-ion batteries, yet it is challenging to construct high-quality coatings at both the primary- and secondary-particle levels. Here, we present a room-temperature synthesis route to achieve full surface coverage of secondary particles and facile infusion into grain boundaries, thus offering a complete “coating-plus-infusion” strategy. Cobalt boride metallic glass is successfully applied to Ni-rich layered cathode LiNi 0.8 Co 0.1 Mn 0.1 O 2 . Here, it dramatically improves the rate capability and cycling stability, including under high-discharge-rate and elevated-temperature conditions and in pouch full cells. The superior performance originates from simultaneous suppression of microstructural degradation of intergranular cracking and side reactions with electrolyte. Atomistic simulations identified the critical role of strong selective interfacial bonding, which offers not only a large chemical driving force to ensure uniform reactive wetting and facile infusion but also lowered the surface/interface oxygen activity, contributing to the exceptional mechanical and electrochemical stabilities of the infused electrode.

25 ENERGY STORAGE↗
DOE OSTI2025

Establishing model credibility for process-microstructure-property relationships in additive manufacturing using exascale computing

Additive Manufacturing (AM) of alloys holds significant promise as a disruptive technology in various industries, yet its adoption is often hindered by challenges in achieving consistent part quality. These issues are primarily due to the complex process-microstructure-property (PSP) relationships inherent to AM. Computational models can greatly aid in understanding these relationships, but their widespread impact and adoption has been limited by a lack of validated, open-source, and computationally efficient PSP modeling frameworks and hardware limitations. Here, this study leverages the ExaAM software suite and data from the AMBench-2018 series of laser powder bed fusion (LPBF) benchmark experiments to perform a comprehensive model assessment, including verification, validation, sensitivity analysis, and uncertainty quantification. The RADICAL-EnTK workflow manager was used to perform an ensemble of heat transport, solidification, and mechanical response simulations on the exascale computer Frontier, considering uncertainties in critical model inputs such as laser spot size and nucleation parameters, and consisting of 125 explicit grain structure simulations and 7875 crystal plasticity simulations. For a selected location within the Inconel 625 AMBench-2018 test artifact, sensitivity analysis and uncertainty quantification were performed using the predicted distributions of grain structure and mechanical properties. Qualitative agreement was found between the predicted grain size and texture and the observed AMBench-2018 microstructure, the mean predicted yield stress was within 5% of the experimental measurement mean, and the mean predicted engineering stress at 5% strain was within 10% of the experimental measurement mean. The insights gained from development and validation of the ExaAM PSP modeling framework will help guide future directions for enhancing the credibility and reliability of PSP models in AM, thereby accelerating the adoption of AM technologies in various industries.

Additive manufacturing↗
DOE OSTI2020

Hot-spot mix in large-scale HDC implosions at NIF

Mix of high-Z material from the capsule into the fuel can severely degrade the performance of inertial fusion implosions. On the Hybrid B campaign, testing the largest high-density-carbon capsules yet fielded at the National Ignition Facility, several shots show signatures of high levels of hot-spot mix. We attribute a ~40% yield degradation on these shots to the hot-spot mix, comparable to the level of degradation from large $P_2$ asymmetries observed on some shots. Additionally, a range of instability growth factors and diamond crystallinity were tested and they do not determine the level of mix for these implosions, which is instead set by the capsule quality.

Fusion experiments↗
DOE OSTI2024

Modeling ablator defects as a source of mix in high-performance implosions at the National Ignition Facility

Recent indirect drive inertial confinement fusion implosions on the National Ignition Facility (NIF) [Spaeth et al., Fusion Sci. Technol. 69, 25 (2016)] have crossed the threshold of ignition. However, performance has been variable due to several factors. One of the leading sources of variability is the quality of the high-density carbon (HDC) shells used as ablators in these experiments. In particular, these shells can have a number of defects that have been found to correlate with the appearance of ablator mix into the hot spot and a degradation in nuclear yield. These defects include pits on the ablator surface, voids in the ablator bulk, high-Z debris from the Hohlraum wall that adheres to the capsule surface, and finally the inherent granular micro-structure of the crystalline HDC itself. This paper summarizes high-resolution modeling of each of these mix sources in two recent high-performance NIF implosion experiments. The simulated impact from a range of individual capsule defects is found to be broadly consistent with the trends seen in experiment, lending credence to the modeling results and the details of the mixing process that they reveal. Interestingly, modeling of the micro-structure inherent to HDC shows that this perturbation source results in considerable mixing of the deuterium–tritium fuel with ablator material during the implosion. The reduction in fuel compression from this mix results in an approximately factor of two reduction in neutron yield in current implosions and emphasizes the importance of mitigating this significant performance degradation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Efficient Optimization of Plasma Radiation Detector Configurations using Imperfect Inference Models

The configurations of instruments fielded on an experiment affect the amount of information captured and the quality of subsequent inference. Here, we investigate the problem of optimizing plasma x-ray radiation detectors in a magneto-inertial fusion experiment at Sandia National Laboratories. It is impossible to directly measure properties such as the temperature of the thermonuclear fusion plasma produced in these experiments because of the extreme environment and destructive nature of the experiment. Among other diagnostics, several detectors are placed with significant standoff from the fusion target to capture the x-rays emitted by the fusion plasma, which can be used to infer some of its properties. To optimize the configuration of these detectors, a high-fidelity model (HFM) is used for simulating outputs and a low-fidelity model (LFM) is used for inference. We develop methods based on A- and L-optimality criteria that are efficient to compute while explicitly accounting for the discrepancy between the HFM and the LFM. The method allows us to find detector configurations that perform similarly to or better than the configuration obtained using an existing sampling-based optimization method while decreasing computational time by a factor of 50. Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work.

Bayesian optimization↗
DOE OSTI2024

Understanding the Cascade: Removing GCM Biases Improves Dynamically Downscaled Climate Projections

Polarization surrounding bias correction (BC) in creating climate projections arises from its lack of physicality. Here, we perform and analyze 18 dynamical downscaling simulations (with and without BC) to better understand the physical impacts of BC, applied before downscaling, on regional climate output across the western United States. Without BC, downscaled precipitation is systematically and unrealistically wet biased compared to a hierarchy of observationally based datasets over the 1980–2014 period due to cascading mean–state Global Climate Model (GCM) biases: (a) overly strong lower–tropospheric lapse rates (5 K/km), (b) overly cold (2 K) tropospheric temperatures, and (c) anomalous mid–tropospheric cyclonic vorticity advection. With BC, downscaled precipitation (snow) biases are virtually eliminated (halved). Identified GCM biases are common to the broader Coupled Model Intercomparison Project ensemble. Physical effects of BC on the quality of the regionalized projections, pending an evaluation of BC's distortion of the downscaled climate response, may motivate its broader application by dynamical downscalers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation of Plot File Testing in the DYNA3D/ParaDyn Software Quality Assurance Suite

Automated testing of DYNA3D/ParaDyn plot files was added to the DYNA3D/ParaDyn software quality assurance (SQA) test suite. The new capability extracts select data from the plot files generated during each verification run and compares it to the same baseline answers used to verify the problem. Deviations between baseline answers and plot file values are reported in the same manner as solution discrepancies, and differences in precision levels between the baseline answers and plot file results are accounted for. The new testing leverages the existing SQA test suite framework and test problems and the Python Mili reader and minimally increases the overall run time (< 5%) of the SQA test suite. This new capability provides incremental end-toend testing of the most common DYNA3D/ParaDyn simulation workflows.

42 ENGINEERING↗
DOE OSTI2025

Observed declines in leaf nitrogen explained by photosynthetic acclimation to CO 2

Widespread evidence of decreasing leaf nutrients has raised concerns about ecosystem productivity under global change. Interpreting trends in leaf nutrients has important implications for the fate of ecosystem services, particularly the role of forests in mitigating climate change and sustaining quality food sources. Here, we challenge the common interpretation that decreasing leaf nitrogen concentration (LNC) is evidence of increasing nutrient limitations on ecosystem primary productivity. Instead, we show that declines in LNC (4% decrease per 50 ppm CO 2 increase), observed across 409 European forest plots over 22 y, can be explained by reduced photosynthetic nitrogen demand. This regional trend is consistent with leaf acclimation to increasing atmospheric CO 2 according to optimality theory. This finding suggests that enhanced photosynthetic nitrogen use efficiency due to CO 2 fertilization may lead to less nitrogen uptake and/or reallocation of nitrogen for plant growth and other functions. Our results have large implications for understanding and simulating interactions between ecosystem nitrogen and carbon cycles and suggest nitrogen requirements for terrestrial carbon uptake under elevated CO 2 may be lower than previously thought.

CO2 fertilization↗
DOE OSTI2024

Development of APCVD BSG and POCl 3 Codiffusion Process for Double-Side TOPCon Solar Cell Precursor Fabrication

This paper presents a commercially viable process for fabricating a high-quality double-side tunnel oxide passivating contact (DS-TOPCon) cell precursor using APCVD-deposited boron silicate glass and ex-situ POCl 3 diffusion in a single high-temperature step, eliminating the need for additional masking and diffusion processes. A two-tier temperature profile was developed, involving a pre-annealing at above 900°C in nitrogen (N2) ambient followed by POCl 3 diffusion at 840°C. We investigated the effect of varying pre-annealing temperatures, ranging from 875°C to 950°C, on the passivation quality and metal-Si contact properties of both n-TOPCon and p-TOPCon layers. The resultant DS-TOPCon cell precursor after silicon nitride (SiNX) passivation exhibited an excellent iV OC of close to 730 mV. In addition, a rapid asymmetric poly-Si thinning technique, developed in this work, enabled adjustment of the front n + poly-Si thickness while maintaining the rear p + poly-Si thickness. Two types of DS-TOPCon cell architectures can be fabricated: i) full-area thin (≈40nm) n-TOPCon layer on the front and ii) selective-area thick (≈200nm) n-TOPCon fingers underneath the metal grid. Device simulations suggest that full-area DS-TOPCon cell with 40 nm n + poly-Si and selective-area DS-TOPCon cell with 200 nm n + poly fingers on the front, fabricated from our current DS-TOPCon cell precursor, can achieve cell efficiencies of 22.1 and 23.5%, respectively. Detailed power loss analysis and device simulation reveal that further improvements in material and device parameters have potential to push the cell efficiencies of DS-TOPCon cell structure beyond 25%, making it a promising alternative to fabricate a high-efficiency next-generation solar cells at low cost.

14 SOLAR ENERGY↗
DOE OSTI2025

Covariance-Free Bifidelity Control Variates Importance Sampling for Rare Event Reliability Analysis

Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.

Markov chain Monte Carlo↗
DOE OSTI2020

Theoretical formulation of phase space microbunching instability in the presence of intrabeam scattering for single-pass or recirculation accelerators

Microbunching instability (MBI) has been one of the most challenging issues in designs of high-brightness beam transport lines for single-pass or recirculating accelerators. Although the intrabeam scattering (IBS) has long been studied in lepton or hadron storage rings as a slow diffusion process or in high-intensity proton linear accelerators as one mechanism for the beam halo, the effects of IBS on single-pass or recirculating electron accelerators have drawn attention only in the recent two decades due to emergence of linac-based or energy-recovery-linac-based fourth-generation light sources, which require high-quality electron beams during the beam transport. In this paper we develop a theoretical formulation of microbunching instability in the presence of IBS for single-pass or recirculation accelerators. To quantify MBI with inclusion of IBS, we start from the Vlasov-Fokker-Planck (VFP) equation, combining both collective interactions and incoherent IBS effects. The linearized VFP equation and the corresponding friction and diffusion coefficients are derived. The evolutions of the resultant density and energy modulations are formulated as a set of coupled integral equations. The theoretical formulation is then applied to a recirculating beamline design. The results from the semianalytical calculation are compared and show good agreement with massive particle tracking simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

SAM Co-Development to Support Fluoride-salt-cooled High-temperature Reactor Design and Licensing (CRADA 2018-18177 Final Report)

This project aims to raise the SAM code's technical and commercial maturity level to enable the Kairos Power ("Kairos") to use SAM to support its KP-FHR design analysis and licensing application. Argonne National Laboratory ("Argonne") has actively developed and maintained SAM, a modern system-level analysis tool for advanced nonlight water reactors safety analysis. Kairos is currently actively pursuing development of an FHR design and associated technology, and consequently requires a pedigreed safety analysis tool. The purpose of this project is to increase the maturity of the SAM code for the modeling and simulation of the KP-FHR design, thus enabling its use in safety analyses that support licensing application. Argonne will primarily be responsible for software design, development, and testing; while Kairos will be responsible for software requirements, assessment of software capabilities, needs, gaps, and priorities, and development of proprietary models. The research and development activities for the joint SAM development project include: a series of identification and prioritization studies on design characteristics, event sequences, relevant phenomena, and software capabilities; SAM capability enhancements for specific KP-FHR systems and components; performing code verification and validations; integrating uncertainty quantification (UQ), model calibration, and sensitivity analysis (SA) techniques in safety analyses; and raising the software quality rigor level for commercial-grade applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

In-Situ Spatial Mapping of Hydrogen in Yttrium Hydrides at LANSCE (FY23 Version, Rev. 1)

This report summarizes the development of neutron imaging capabilities and experimental activities performed at the Los Alamos Neutron Science Center (LANSCE) with the main goal of measuring temperature-driven hydrogen diffusion within bulk-yttrium hydride (YH x ) materials. Yttrium hydride is the leading candidate to serve as a solid neutron moderator in microreactor cores, owing to its high density of hydrogen atoms as well as its superior thermal stability compared to all other metal hydrides. The experimental results and technique developments reported herein support the U.S. Department of Energy Office of Nuclear Energy’s (DOE-NE) Microreactor Program under Technology Maturation. In particular, it addresses the critical need to experimentally validate and verify hydrogen-diffusion models of metal hydrides used in high-temperature microreactor designs by means of high-spatial-resolution neutron imaging. These capabilities were designed to apply large temperature gradients across centimeter-sized YH x pellets to simulate conditions faced in the microreactor environment. In principle, neutron imaging, combined with in-situ sample heating, enables near real-time tracking of hydrogen diffusion in YH x on the sub-millimeter scale. In this report, an overview of neutron imaging methodology and technologies are given in the context of recent spatial measures of hydrogen concentrations in similar metal hydrides. Additionally, the commissioning and operation of a custom-built compact dual-zone furnace is given along with details on three in-situ heating measurements of YH x performed over the 2020 to 2022 LANSCE operation cycles. The aims of these experiments ranged from furnace commissioning, determining sample quality, i.e., hydrogen uniformity via neutron computed tomography, and studying the effects of applied temperature-gradients on YH x pellets. Analyses and results from these neutron imaging measurements are given along with outlooks and guidelines for optimal future hydrogen diffusion measurements. Our conclusions are as follows. Image analyses indicate that centimeter-sized yttrium hydride cylindrical pellets exhibit uniform, whole-body hydrogen desorption and absorption without clear temperature dependence as reflected in the image attenuation at the opposing ends of each sample. This suggests that despite the large magnitude in temperature gradients applied by the furnace heating elements, the sample equilibrates to an unknown intermediate temperature. The origin of this result is likely the combination of short sample length (∼1cm) and use of a TZM can for containment where the latter created a thermal short across the sample. Nevertheless, the results from the most recent measurements indicate that neither significant concentration gradients of hydrogen were formed in centimeter-sized samples through the entire temperature range (25 °C to 950 °C) nor any formed due to temperature gradients on the order of 50 °C/cm up to 700 °C/cm. Furthermore, images from the FY2021 and FY2022 measurements indicate that samples of YH x , fabricated from either the direct hydride or powder metallurgy methods, are highly uniform in their hydrogen concentration to within the measurements’ spatial resolutions. The following questions arise from these latest results: 1) What is the intermediate temperature of the pellets in the TZM cans? 2) How quickly does the temperature equilibrate within the sample? and, 3) Do the observed changes in image attenuation follow known pressure-composition-temperature relations of yttrium hydride?

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A signal compensation based cascaded PI control for an industrial heat exchange system

Since the characteristics of the industrial heat exchange system (HES) are drastically affected by variations of the steam pressure, ambient temperature and water quality, it is difficult for the classical cascade PI control system to meet the desirable control objectives. Therefore, in this paper, a deterministic linear model with unmodeled dynamics is used to describe the dynamic characteristics of the HES. It should be noted that the unmodeled dynamics at the previously sampled instant can be calculated using actual measurements and that the tracking error, which includes effects of the rate of changes of unmodeled dynamics on the controlled model, can be accurately determined. An additional compensating signal is obtained with the feedforward control together with a one-step optimal control law. Consequently, the cascade PI control strategy with the signal compensation is established. Furthermore, the simulation and industrial application show that the proposed method can effectively improve the performance of the cascade PI control system, especially when the unmodeled dynamics is subjected to unknown variations.

42 ENGINEERING↗
DOE OSTI2020

Crystal Structure and Atomic Vacancy Optimized Thermoelectric Properties in Gadolinium Selenides

Thermoelectric materials enable the energy conversion of waste heat into electricity, helpful to relieve global energy crisis. Here, we report a systematic investigation on high-temperature thermoelectric gadolinium selenides, cubic Gd 3-x Se 4 (x = 0.16, 0.21 and 0.25) and orthorhombic Gd 2 Se 3-y (y = 0.02, 0.06 and 0.08). High energy synchrotron x-ray diffraction and total scattering have been used to investigate the crystallographic and local structures. Atomic-scale clusters of Gd vacancy in the cubic phase are observed by employing the reverse Monte Carlo simulation. For cubic Gd 3-x Se 4 , adjusting Gd vacancy triggers the effect of multiple conduction bands, confirmed by the increase in effective masses and theoretical calculations. A reasonable peak zT value of 0.27 is achieved at 850 K for Gd 3-x Se 4 (x = 0.16). On the other hand, tuning Se vacancy enables the optimization of electron concentration for the orthorhombic Gd 2 Se 3-y . More significantly, its low deformation potential (Ξ = 12eV) gives rise to enhanced electron mobility and higher peak zT value of 0.54 at 850 K for Gd 2 Se 3-y (y = 0.02). Intriguingly, a higher zT of 1.2 at 1200 K is reasonably predicted by quality factor analysis. Finally, this work extends the scope of high-temperature thermoelectric materials and facilitates the exploration of novel high-temperature thermoelectric materials

36 MATERIALS SCIENCE↗
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↗