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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Microgravity↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Crystal Growth↗

Use of MODIS Terra Imagery to Estimate Surface Water Quality Standards, Using Lake Thonotosassa, Florida, as a Case Study

Lake Thonotosassa is a highly eutrophied lake located in an area with rapidly growing population in the Tampa Bay watershed, Florida. The Florida Administrative Code has designated its use for "recreation, propagation and maintenance of a healthy, well-balanced population of fish and wildlife." Although this lake has been the subject of efforts to improve water quality since 1970, overall water quality has remained below the acceptable state standards, and has a high concentration of nutrients. This condition is of great concern to public health since it has favored episodic blooms of Cyanobacteria. Some Cyanobacterial species release toxins that can reach humans through drinking water, fish consumption, and direct contact with contaminated water. The lake has been historically popular for fishing and water sports, and its overflow water drains into the Hillsborough River, the main supply of municipal water for the City of Tampa, this explains why it has being constantly monitored in situ for water quality by the Environmental Protection Commission of Hillsborough County (EPC). Advances in remote sensing technology, however, open the possibility of facilitating similar types of monitoring in this and similar lakes, further contributing to the implementation of surveillance systems that would benefit not just public health, but also tourism and ecosystems. Although traditional application of this technology to water quality has been focused on much larger coastal water bodies like bays and estuaries, this study evaluates the feasibility of its application on a 46.6 km2 freshwater lake. Using surface reflectance products from Moderate-Resolution Imaging Spectroradiometer (MODIS) Terra, this study evaluates associations between remotely sensed data and in situ data from the EPC. The parameters analyzed are the surface water quality standards used by the State of Florida and general indicators of trophic status.

Moreno, Max J.↗

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↗

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↗

Estimating the Rate of Occurrence of Renal Stones in Astronauts

Changes in urine chemistry, during and post flight, potentially increases the risk of renal stones in astronauts. Although much is known about the effects of space flight on urine chemistry, no inflight incidence of renal stones in US astronauts exists and the question "How much does this risk change with space flight?" remains difficult to accurately quantify. In this discussion, we tackle this question utilizing a combination of deterministic and probabilistic modeling that implements the physics behind free stone growth and agglomeration, speciation of urine chemistry and published observations of population renal stone incidences to estimate changes in the rate of renal stone presentation. The modeling process utilizes a Population Balance Equation based model developed in the companion IWS abstract by Kassemi et al. (2016) to evaluate the maximum growth and agglomeration potential from a specified set of urine chemistry values. Changes in renal stone occurrence rates are obtained from this model in a probabilistic simulation that interrogates the range of possible urine chemistries using Monte Carlo techniques. Subsequently, each randomly sampled urine chemistry undergoes speciation analysis using the well-established Joint Expert Speciation System (JESS) code to calculate critical values, such as ionic strength and relative supersaturation. The Kassemi model utilizes this information to predict the mean and maximum stone size. We close the assessment loop by using a transfer function that estimates the rate of stone formation from combining the relative supersaturation and both the mean and maximum free stone growth sizes. The transfer function is established by a simulation analysis which combines population stone formation rates and Poisson regression. Training this transfer function requires using the output of the aforementioned assessment steps with inputs from known non-stone-former and known stone-former urine chemistries. Established in a Monte Carlo system, the entire renal stone analysis model produces a probability distribution of the stone formation rate and an expected uncertainty in the estimate. The utility of this analysis will be demonstrated by showing the change in renal stone occurrence predicted by this method using urine chemistry distributions published in Whitson et al. 2009. A comparison to the model predictions to previous assessments of renal stone risk will be used to illustrate initial validation of the model.

gravitational physiology↗

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↗