Search NASA⌕ Search

SEARCH · Search NASA

Results for “Population modeling”

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 199 records · Page 11

Nested active learning for efficient model contextualization and parameterization: pathway to generating simulated populations using multi-scale computational models

There is increasing interest in the use of mechanism-based multi-scale computational models (such as agent-based models (ABMs)) to generate simulated clinical populations in order to discover and evaluate potential diagnostic and therapeutic modalities. The description of the environment in which a biomedical simulation operates (model context) and parameterization of internal model rules (model content) requires the optimization of a large number of free parameters. In this work, we utilize a nested active learning (AL) workflow to efficiently parameterize and contextualize an ABM of systemic inflammation used to examine sepsis. Contextual parameter space was examined using four parameters external to the model’s rule set. The model’s internal parameterization, which represents gene expression and associated cellular behaviors, was explored through the augmentation or inhibition of signaling pathways for 12 signaling mediators associated with inflammation and wound healing. We have implemented a nested AL approach in which the clinically relevant (CR) model environment space for a given internal model parameterization is mapped using a small Artificial Neural Network (ANN). The outer AL level workflow is a larger ANN that uses AL to efficiently regress the volume and centroid location of the CR space given by a single internal parameterization. We have reduced the number of simulations required to efficiently map the CR parameter space of this model by approximately 99%. In addition, we have shown that more complex models with a larger number of variables may expect further improvements in efficiency.

97 MATHEMATICS AND COMPUTING↗

Utility of Physiologically Based Pharmacokinetic Modeling to Investigate the Impact of Physiological Changes of Pregnancy and Cancer on Oncology Drug Pharmacokinetics

The treatment of cancer during pregnancy remains challenging with knowledge gaps in drug dosage, safety, and efficacy due to the under-representation of this population in clinical trials. Our aim was to investigate physiological changes reported in both pregnancy and cancer populations into a PBPK modeling framework that allows for a more accurate estimation of PK changes in pregnant patients with cancer. Methods: Paclitaxel and docetaxel were selected to validate a population model using clinical data from pregnant patients with cancer. The validated population model was subsequently used to predict the PK of acalabrutinib in pregnant patients with cancer. Results: The Simcyp pregnancy population model reasonably predicted the PK of docetaxel in pregnant patients with cancer, while a modified model that included a 2.5-fold increase in CYP2C8 abundance, consistent with the increased expression during pregnancy, was needed to reasonably predict the PK of paclitaxel in pregnant patients with cancer. Changes in protein binding levels of patients with cancer had a minimal impact on the predicted clearance of paclitaxel and docetaxel. PBPK modeling predicted approximately 60% lower AUC and Cmax for acalabrutinib in pregnant versus non-pregnant patients with cancer. Conclusions: Our results suggest that PBPK modeling is a promising approach to investigate the effects of pregnancy and cancer on the PK of oncology drugs and potentially inform dosing for pregnant patients with cancer. Further evaluation and refinement of the population model are needed for pregnant patients with cancer with additional compounds and clinical PK data.

60 APPLIED LIFE SCIENCES↗

Numerical investigation of nonequilibrium electron effects on the collisional ionization rate in the collisional-radiative model

In this study, the interplay of kinetic electron physics and atomic processes in ultrashort laser-plasma interactions provides a comprehensive understanding of the impact of the electron energy distribution on plasma properties. Notably, nonequilibrium electrons play a vital role in collisional ionization, influencing ionization degrees and spectra. This paper introduces a computational model that integrates the physics of kinetic electrons and atomic processes, utilizing a Boltzmann equation for nonequilibrium electrons and a collisional-radiative model for atomic state populations. The model is used to investigate the influence of nonequilibrium electrons on collisional ionization rates and its effect on the population distribution, as observed in a widely known experiment. The study reveals a significant nonequilibrium electron presence during XFEL-matter interactions, profoundly affecting collisional ionization rates in the gas plasma, thereby necessitating careful consideration of the Collisional-Radiative model applied to such systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Distinct critical behaviors from the same state in quantum spin and population dynamics perspectives

There is a deep connection between the ground states of transverse-field spin systems and the late-time distributions of evolving viral populations—within simple models, both are obtained from the principal eigenvector of the same matrix. However, that vector is the wave-function amplitude in the quantum spin model, whereas it is the probability itself in the population model. We show that this seemingly minor difference has significant consequences: Phase transitions that are discontinuous in the spin system become continuous when viewed through the population perspective, and transitions that are continuous become governed by new critical exponents. We introduce a more general class of models that encompasses both cases and that can be solved exactly in a mean-field limit. Numerical results are also presented for a number of one-dimensional chains with power-law interactions. We see that well-worn spin models of quantum statistical mechanics can contain unexpected new physics and insights when treated as population-dynamical models and beyond, motivating further studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗

Population Viability Analysis (PVA) as a Platform for Predicting Outcomes of Management Options for the Florida Scrub-Jay in Brevard County

The Florida scrub-jay (FSJ) is a species in decline because of extinction debt caused by habitat fragmentation and degradation. Understanding and managing the species and its habitats is challenging due to the complex interactions among the social system, a dynamic mosaic of scrub habitat, and active management. To provide insights into the possible fates of the FSJ populations of mainland, cape, and island sections of Brevard County, we modeled the population dynamics using the Vortex population viability analysis (PVA) software. Vortex is an individual-based simulation that allowed us to include such factors as demographic rates dependent on habitat state, impact of helpers on breeding success, and helper to breeder transition probabilities responding to availability of nearby vacant optimal habitat and vacancies due to the death of breeders. Detailed modeling of the FSJ was possible only because a lot of data about the species and its habitats have been gathered over decades of intensive research. We followed a phased approach to constructing population models that incorporated the best available science and data to address a variety of conservation actions. The first phase focused on constructing a model that incorporates sociobiology and source-sink habitat dynamics. This model allowed us to address some of the most important questions about population size and habitat quality. Once the model framework was in place, we considered the real landscapes and actual local populations, rather than just generic representations of typical FSJ dynamics. After examining the viability of the metapopulations under current conditions, we explored the likely consequences of various management actions that might slow or reverse population declines.

Population Viability Analysis↗

Comparing Two Statistical Models for Low Boom Dose-response Relationships with Correlated Responses

This study compares two statistical modeling approaches to construct summary dose-response curves for low boom community noise surveys. NASA field survey data were used that consist of multiple responses from a survey participant. These data require an approach that accounts for the correlation among repeated annoyance observations from the same participant. A multilevel model accounts for the correlation by allowing estimated parameters to vary with each survey participant. On the other hand, a population average model utilizes generalized estimating equations and accounts for the correlation via a user-specified within-subject correlation structure. A visual comparison of the dose-response curves for these two methods reveals similar results. When comparing specific points along the summary curves, the multilevel model yields more precise confidence bounds than the population average model. The similarity between the summary curves derived from each model lends validity to both approaches for approximating a population representative summary curve.

X-59↗

Spectroscopy of candidate young globular clusters in NGC 1275

We present spectra of the brightest member of the population of compact blue objects discovered in the peculiar galaxy NGC 1275 by Holtzman et al. (1992) using Hubble Space Telescope images. These spectra show strong Balmer absorption lines like those observed in A-type stars, as expected if the object is a young globular cluster. The age estimated from the strength of the Balmer lines is about 0.5 Gyr, although ages ranging from 0.1 Gyr to 0.9 Gyr cannot be confidently excluded given current models of stellar populations. If these estimated ages are adopted for the young cluster population of NGC 1275 as a whole, the fading predicted by stellar populations models gives a luminosity function which is consistent with that of the Galactic globular cluster system convolved with the observational selection function for the NGC 1275 system. We also use the equivalent widths of the Mg b and Fe 5270 features to constrain the metallicity of the young cluster. Combining these absorption-line widths with the age estimates from the Balmer lines and stellar population models, we find a metallicity of roughly solar, based on the Mg b index, and somewhat higher for the Fe 5270 index. The radial velocity of the absorption lines of the cluster spectrum is offset from the emission lines of the galaxy spectrum at the same position by - 130 km/s, providing further evidence for the identification of the object as a global cluster and opening up the future possibility of studying the kinematics of young cluster systems. The discovery of objects with the characteristics of young globular clusters in NGC 1275, which shows evidence of a recent interaction or merger, supports the hypothesis that galaxy interactions and mergers are favorable sites for the formation of globular clusters.

Zepf, Stephen E.↗

Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 10$^{17}$ eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5 × 10$^{18}$ eV (the so-called “ankle” feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources.To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).

79 ASTRONOMY AND ASTROPHYSICS↗

Improved biomass feedstock materials handling and feeding engineering data sets, design methods, and modeling/simulation tools

Forest Concepts led a project to develop a set of tools which enable feedstock handing equipment designers to better understand and model the flowability of bulk particulate biomass materials. The work products included improved flowability mathematical models, data sets used to populate the models, and two new laboratory devices – a biomass-scale true cubical triaxial tester and a biomass-scale gas pycnometer. This work was funded in part by the US Department of Energy under contact DE-EE0008254.

09 BIOMASS FUELS↗

A strong blend in the morning: studying the circumgalactic medium before cosmic noon with strong, blended Ly α forest systems

ABSTRACT We study of the properties of a new class of circumgalactic medium absorbers identified in the Ly α forest: ‘Strong, Blended Lyman-α’ (or SBLA) absorption systems. We study SBLAs at 2.4 < z < 3.1 in SDSS-IV/eBOSS spectra by their strong extended Ly α absorption complexes covering 138 $\, \, {\rm km}\, {\rm s}^{-1}$ with an integrated $\log (N_{\rm H\, {\small I}}/\mathrm{cm}^{-2}) =16.04$$\substack{+0.05 \\ -0.06}$ and Doppler parameter b = 18.1$\substack{+0.7 \\ -0.4}$$\, \, {\rm km}\, {\rm s}^{-1}$. Clustering with the Ly α forest provides a large-scale structure bias of b = 2.34 ± 0.06 and halo mass estimate of $M_h \approx 10^{12}\, h^{-1}\, {\rm M_{\odot }}$ for our SBLA sample. We measure the ensemble mean column densities of 22 metal features in the SBLA composite spectrum and find that no single-population multiphase model for them is viable. We therefore explore the underlying SBLA population by forward modelling the SBLA absorption distribution. Based on covariance measurements and favoured populations we find that ≈25 per cent of our SBLAs have stronger metals. Using silicon only we find that our strong metal SBLAs trace gas with a log (nH/cm−3) > −2.40 for T = 103.5 K and show gas clumping on <210 parsec scales. We fit multiphase models to this strong subpopulation and find a low ionization phase with nH = 1 cm−3, T = 103.5 K, and [X/H] = 0.8, an intermediate ionization phase with log (nH/cm−3) = −3.05, T = 103.5 K and [X/H] = −0.8, and a poorly constrained higher ionization phase. We find that the low ionization phase favours cold, dense super-solar metallicity gas with a clumping scale of just 0.009 parsecs.

Morrison, Sean (ORCID:0000000267702627)↗

LandScan Global 2023: Silver Edition

For a quarter of a century, the LandScan Global (LSG) project has annually released a global, high-resolution gridded population dataset representing the ambient or unwarned population at a 30 arcsecond resolution. LSG supports a range of applications such as emergency management, disaster response, and human health and security for understanding populations at risk. The 2023 release of LSG, the LandScan Silver Edition, represents a major methodological leap forward while also leveraging previous knowledge—the previous year was the baseline for the current annual update carrying forward valuable knowledge of the built environment for the past quarter century—to train the machine learning models. Compared with annual releases over the past 24years, multiple advancements were made to different aspects of the methodology to achieve reproducibility, transparency, and consistent global propagation of solutions to modeling or population distribution issues identified during the review process. These novel changes include incorporation of the latest available geospatial inputs across the globe, machine learning models instead of manual modifications, population feature importance analysis, open-source solutions vs. proprietary software, generation of multiple global versions, analytic validations, and human-in-the-loop revisions to produce the final version. Additionally, algorithms—such as anomaly detection—were introduced to quickly identify areas of focus to develop a new and robust systematic review. Significant changes in modeled population distributions were observed between the 2022 and 2023 releases, largely attributable to improvements in data and methods and discussed thoroughly within this report. In summation, the LandScan Silver Edition leverages the best of the past quarter century of LSG legacy knowledge and continues a tradition of applying cutting-edge enhancements to serve as a new benchmark for accurate, actionable gridded population data

Lebakula, Viswadeep↗

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗

Modeling of the Orbital Debris Population of RORSAT Sodium-Potassium Droplets

A large population resident in the orbital debris environment is composed of eutectic sodium-potassium (NaK) droplets, released during the reactor core ejection of 16 nuclear-powered Radar Ocean Reconnaissance Satellites (RORSATs) launched in the 1980s by the former Soviet Union. These electrically conducting RORSAT debris objects are spherical in shape, generating highly polarized radar returns. Their diameters are mostly in the centimeter and millimeter size regimes. Since the Space Surveillance Network catalog is limited to objects greater than 5 cm in low Earth orbit, our current knowledge about this special class of orbital debris relies largely on the analysis of Haystack radar data. This paper elaborates the simulation of the RORSAT debris populations in the new NASA Orbital Debris Engineering Model ORDEM2010, which replaces ORDEM2000. The estimation of the NaK populations uses the NASA NaK-module as a benchmark. It follows the general statistical approach to developing all other ORDEM2010-required LEO populations (for various types of debris and across a wide range of object sizes). This paper describes, in detail, each major step in the NaK-population derivation, including a specific discussion on the conversion between Haystack-measured radar-cross-sections and object-size distribution for the NaK droplets. Modeling results show that the RORSAT debris population is stable for the time period under study and that Haystack data sets are fairly consistent over the observations of multiple years.

Xu, Y.-L.↗