Search NASA⌕ Search

SEARCH · Search NASA

Results for “population balances”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Predicting biomass comminution: Physical experiment, population balance model, and deep learning

An extended population balance model (PBM) and a deep learning-based enhanced deep neural operator (DNO+) model are introduced for predicting particle size distribution (PSD) of comminuted biomass through a large knife mill. Experimental tests using corn stalks with varied moisture contents, mill blade speeds, and discharge screen sizes are conducted to support model development. A novel mechanism in the extended PBM allows for including additional input parameters such as moisture content, which is not possible in the original PBM. The DNO+ model can include influencing factors of different data types such as moisture content and discharge screen size, which significantly extends the engineering applicability of the standard DNO model that only admits feed PSD and outcome PSD. Test results show that both models are remarkably accurate in the calibration or training parameter space and can be used as surrogate models to provide effective guidance for biomass preprocessing design.

09 BIOMASS FUELS↗
DOE OSTI2022

Population balance models for polymer upcycling: signatures of the mechanism in the molecular weight evolution

Chemical and catalytic upcycling processes could help realize a circular plastics economy, but current models for testing mechanistic hypotheses and designing catalysts remain primitive. This work shows how proposed catalytic mechanisms can be incorporated into population balance models to predict the time evolution of molecular weight distributions. We develop models for homogeneous and heterogeneous catalysts, including catalysts that cut at chain ends and catalysts that cut at random locations. For heterogeneous catalysts, we illustrate the effect of adsorption constants that depend on polymer chain length. Here, we discuss ongoing efforts and challenges in measuring and modeling the time evolving molecular weight distributions in polymer upcycling processes.

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

Population Balance Models for Catalytic Depolymerization: From Elementary Steps to Multiphase Reactors

Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.

Manis, Lela K. [University of Illinois at Urbana-C↗
DOE OSTI2023

Nanoparticle Catalyst Formation and Subsequent Sintering: Size-Distribution Fundamental Understanding and Catalyst Control via Mechanism-Enabled Population Balance Modeling (Final Report)

The challenge of elucidating the chemical mechanisms by which nano- and other particles form, and then using those mechanisms to understand, analyze, and ideally predict particle-size and particle-size distributions, was undertaken. Thanks to multi-year, longer-term DOE support each of these challenges was successfully met. The community is now able to both predict, and control in cases where the underlying mechanism of particle formation has been determined, the key PS and PSD that are, in turn, critical properties that strongly influence the catalytic and other physical properties of nano- and larger particles across nature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Modeling simultaneous particle shrinkage, dissolution and breakage using the modified moving grid technique

Simultaneous shrinkage, dissolution and breakage are important particle size reduction phenomena that characterize processes like the reactive degradation of solid chemicals. The dynamics of the particle size distribution (PSD) for such processes are non-trivial to model due to the number expending processes brought about by the eventual dissolution of particles. To this end, Population Balance Model (PBM) resolved through the sectional techniques is the natural approach. Here, we introduce a modified Moving Grid technique (m-MGT) to accurately resolve the particle size reduction phenomena. Our technique mimics the perpetual particle shrinkage through a continuously left-moving size grid and incorporates a strategic grid removal routine to capture the disappearance of particles. Coupled with the Fixed Pivot (FP) discretization for breakage, our m-MGT not only preserves the moment-related properties, but also benchmarked very well against the analytical number densities and exhibited a minimum of first-order convergence in all assessed case studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mathematical Modeling of Multi-Phenomena Anisotropic Systems: Ejection of Primary Aerosols during the Fast Pyrolysis of Biomass

This study introduces a novel particle model for biomass fast pyrolysis, incorporating an anisotropic cylindrical particle to address mass and energy transport coupled with aerosol ejection, which previous models have overlooked. The main contribution lies in developing a model that considers aerosol generation in anisotropic cylindrical particles for the first time, addressing bubbling dynamics and bursting within the liquid phase. The population balance equation describes bubble dynamics and aerosol formation, capturing phenomena like nucleation, growth, coalescence, and bursting. The model employs the method of moments with bubble volume as an internal variable, substantially reducing computational costs by eliminating dependence on this variable. Results highlight the significant impact of anisotropy and particle size on aerosol ejection: smaller, less elongated particles experience faster heating, quicker conversion, and the increased accumulation of the liquid intermediate phase. Specifically, 1 mm diameter particles yield higher concentrations of metaplast and bio-oil aerosols, exceeding 15%, compared to concentrations below 11% for 3 mm particles. This model provides insights into aerosol structure (volume, surface area), aiding in understanding aerosol reactivity at the reactor scale.

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

Modelling the ejection of primary aerosols during the fast pyrolysis of biomass anisotropic particles

A model for the fast pyrolysis of anisotropic biomass particles is presented which considers bubbling dynamics within the liquid intermediate phase (metaplast) and aerosol ejection from this phase. The model employs the population balance equation and the method of moments to estimate the production rate and resultant size distribution of aerosol ejections, incorporating a detailed CRECK reaction mechanism, and considers the effect of anisotropic biomass microstructure on the intraparticle transport of mass and energy. Here, this study investigates the impact of particle size, heating rate (heat transfer coefficient), and lignocellulosic composition on aerosol ejection. The model predicts that, at high heating rates (convective heat transfer coefficient of 359 W/m 2 .K), aerosols can contribute over 20% to the heavy fraction yield in bio-oil for small particles (1 mm diameter, 4 mm length). The model can predict aerosol size distribution and surface area, indicating an average size of 20 μm for bubbles and 5 μm for aerosols during increased bubble production and aerosol ejection rates. These findings are consistent with prior experimental results and provide essential information for future modeling of extra-particle reactions of the aerosols as they progress through the reactor.

09 BIOMASS FUELS↗
DOE OSTI2023

Smart Preprocessing & Robust Integration Emulator

To achieve the desired particle size of biomass feedstocks during preprocessing for trouble-free handling and conversion to produce biofuels and bioproducts, the raw materials must undergo a crucial milling process. The particle size of biomass plays a critical role in subsequent biofuel manufacturing, where a larger area-to-volume ratio facilitates efficient synthesis while balancing the impact of moisture on biomass storage. To optimize biofuel production efficiency and overcome these challenges, it is imperative to accurately predict the particle size distribution (PSD) of the biomass in the design of efficient preprocessing systems. The population balance model (PBM), upon empirical calibration and validation, can provide rapid prediction of post-milling PSD of granular biomass. However, PSD has limitations related to mass conservation and the absence of moisture considerations. To overcome these drawbacks, a deep learning model called the enhanced deep neural operator (DNO+) is implemented in the code. This model not only retains the capabilities of the PBM in handling complex mapping functions but also incorporates additional factors influencing the system. By considering various experimental conditions such as sieve size and moisture content, the trained DNO+ model can effectively predict the PSD after milling for any given feed PSD. To further reduce the reliance on experimental data, the PBM is integrated into the DNO+ model, resulting in a physics-informed DNO+ (PIDNO+). The PIDNO+ model addresses the non-conservation of quality exhibited by the PBM while inheriting the advantages of the DNO+ model in considering multiple influencing factors. Moreover, the PIDNO+ model significantly reduces the amount of data required for model training. Both deep learning models, i.e., DNO+ and PIDNO+, are excellent in predictive performance, offering swift and accurate machine learning-based predictions. The use of this code that contains these models will assist in guiding the proper milling equipment selection and operational conditions to achieve the desired biomass particle sizes, ensuring the efficiency of subsequent biofuel and bioproduct production processes.

Xia, Yidong [Idaho National Laboratory (INL), Idah↗
DOE OSTI2024

Reducing Sample Size While Improving Equity in Vaccine Clinical Trials: A Machine Learning-Based Recruitment Methodology with Application to Improving Trials of Hepatitis C Virus Vaccines in People Who Inject Drugs

Despite the availability of direct-acting antivirals that cure individuals infected with the hepatitis C virus (HCV), developing a vaccine is critically needed in achieving HCV elimination. HCV vaccine trials have been performed in populations with high incidence of new HCV infection such as people who inject drugs (PWID). Developing strategies of optimal recruitment of PWID for HCV vaccine trials could reduce sample size, follow-up costs and disparities in enrollment. We investigate trial recruitment informed by machine learning and evaluate a strategy for HCV vaccine trials termed PREDICTEE—Predictive Recruitment and Enrichment method balancing Demographics and Incidence for Clinical Trial Equity and Efficiency. PREDICTEE utilizes a survival analysis model applied to trial candidates, considering their demographic and injection characteristics to predict the candidate’s probability of HCV infection during the trial. The decision to recruit considers both the candidate’s predicted incidence and demographic characteristics such as age, sex, and race. We evaluated PREDICTEE using in silico methods, in which we first generated a synthetic candidate pool and their respective HCV infection events using HepCEP, a validated agent-based simulation model of HCV transmission among PWID in metropolitan Chicago. We then compared PREDICTEE to conventional recruitment of high-risk PWID who share drugs or injection equipment in terms of sample size and recruitment equity, with the latter measured by participation-to-prevalence ratio (PPR) across age, sex, and race. Comparing conventional recruitment to PREDICTEE found a reduction in sample size from 802 (95%: 642–1010) to 278 (95%: 264–294) with PREDICTEE, while also reducing screening requirements by 30%. Simultaneously, PPR increased from 0.475 (95%: 0.356–0.568) to 0.754 (95%: 0.685–0.834). Even when targeting a dissimilar maximally balanced population in which achieving recruitment equity would be more difficult, PREDICTEE is able to reduce sample size from 802 (95%: 642–1010) to 304 (95%: 288–322) while improving PPR to 0.807 (95%: 0.792–0.821). PREDICTEE presents a promising strategy for HCV clinical trial recruitment, achieving sample size reduction while improving recruitment equity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗
DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2022

Particle scale impacts on deconstruction energy of pine residues

The goal of this Case Study was to quantify the impacts of variable moisture and ash on hammer mill throughput and energy consumption and on generation of fines that are not able to be fed to conversion, as compared to a status quo Base Case system. Also considered was convertible carbon content (minimum carbon specification) and maximum ash content and the delivered feedstock cost impacts of not being able to feed residue not meeting both specifications to the conversion reactor. Laboratory data on the impacts of input particle size and moisture content on the exit particle size were received from FCIC Subtask 5.2 from their single particle impact population balance modeling study. Additional throughput and energy consumption data were obtained from FCIC Subtask 5.2 for the same grinder with a 6 mm screen in place. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the ash contents in the separated fines had not been analyzed in the laboratory at the time of the model runs, we chose to assume that the ash distributed proportionally with total mass into the overs and unders in the disk screen following grinding. Key takeaways from this Case Study are that it is significantly more cost effective to hammer mill the residue prior to drying, even though the grinder throughput is lower and energy consumption is higher versus drying first before grinding. An effect of dry grinding versus high moisture grinding is the production of higher amounts of fines during dry grinding, leading to significantly more of the ground feedstock being rejected by conversion for being below a minimum particle size. With wet grinding the system is still able to produce more preprocessed feedstock meeting the minimum particle size specification even though the instantaneous throughput is lower than for the case of grinding dry feedstock. Additionally, even without the higher fines production from dry grinding, the status quo would still be more costly than wet grinding because the material is rejected after the drying energy has already been input for the dry grinding case. Finally, significant reductions in drying energy are obtained by drying after grinding, and those reductions are of far greater magnitude than the grinding energy increase.

09 BIOMASS FUELS↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Charging, aggregation, and electrostatic dispersion of radioactive and nonradioactive particles in the atmosphere

Electrostatic dispersion can significantly impact the microphysical behavior of charged particles and ions until reaching zero space charge. However, although radioactive particles can be strongly charged in air, the influence of electrostatic dispersion has been neglected in understanding their behavior. This study is aimed at investigating the time evolution of the charge and size distributions of radioactive and nonradioactive particles in air and developing simple approaches for applications. With processes involving charging, aggregation, and electrostatic dispersion, a comprehensive population balance model (PBM) has been developed to examine particle charge/size distribution dynamics. It is shown that compared to nonradioactive particles, the charge and size distributions of radioactive particles may evolve differently with time because radioactivity and electrostatic dispersion can significantly affect the charging and aggregation kinetics of the particles. It is found that, after the Fukushima accident, background aerosols in the pathway of radioactive plumes might be highly charged due to ionizing radiation, suggesting that radiation fields may strongly influence in situ measurements of charged atmospheric particles. The comprehensive PBM is simplified, and then the verification and application of the simplified PBMs are discussed. This study provides useful insight into how radioactivity can affect the dynamic behavior of particles in atmospheric systems including radiation sources.

54 ENVIRONMENTAL SCIENCES↗