Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compartmental Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Staged progression epidemic models for the transmission of invasive nontyphoidal Salmonella (iNTS) with treatment

We develop and analyze a stage-progression compartmental model to study the emerging invasive nontyphoidal Salmonella (iNTS) epidemic in sub-Saharan Africa. iNTS bloodstream infections are often fatal, and the diverse and non-specific clinical features of iNTS make it difficult to diagnose. We focus our study on identifying approaches that can reduce the incidence of new infections. In sub-Saharan Africa, transmission and mortality are correlated with the ongoing HIV epidemic and severe malnutrition. We use our model to quantify the impact that increasing antiretroviral therapy (ART) for HIV infected adults and reducing malnutrition in children would have on mortality from iNTS in the population. We consider immunocompromised subpopulations in the region with major risk factors for mortality, such as malaria and malnutrition among children and HIV infection and ART coverage in both children and adults. We parameterize the progression rates between infection stages using the branching probabilities and estimated time spent at each stage. We interpret the basic reproduction number R 0 as the total contribution from an infinite infection loop produced by the asymptomatic carriers in the infection chain. The results indicate that the asymptomatic HIV+ adults without ART serve as the driving force of infection for the iNTS epidemic. We conclude that the worst disease outcome is among the pediatric population, which has the highest infection rates and death counts. Our sensitivity analysis indicates that the most effective strategies to reduce iNTS mortality in the studied population are to improve the ART coverage among high-risk HIV+ adults and reduce malnutrition among children.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2014

A Population Pharmacokinetic Model for Disposition in Plasma, Saliva and Urine of Scopolamine after Intranasal Administration to Healthy Human Subjects

An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials with an Investigative New Drug (IND) protocol. The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trials with INSCOP. Methods: Twelve healthy human subjects were administered three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min and 24 h after dosing and scopolamine concentrations were measured by using a validated LC‐MS‐MS assay. Pharmacokinetic Compartmental models, using actual dosing and sampling times, were built using Phoenix (version 1.2). Model selection was based on the likelihood ratio test on the difference of criteria (‐2LL) and comparison of the quality of fit plots. Results: The best structural model for INSCOP (minimal ‐2LL= 502.8) was established. It consisted of one compartment each for plasma, saliva and urine, respectively, which were connected with linear transport processes except the nonlinear PK process from plasma to saliva compartment. The best‐fit estimates of PK parameters from individual PK compartmental analysis and Population PK model analysis were shown in Tables 1 and 2, respectively. Conclusion: A population PK model that could predict population and individual PK of scopolamine in plasma, saliva and urine after dosing was developed and validated. Incorporating a non‐linear transfer from plasma to saliva compartments resulted in a significantly improved model fitting. The model could be used to predict scopolamine plasma concentrations from salivary and urinary drug levels, allowing non‐invasive therapeutic monitoring of scopolamine in space and other remote environments.

Wu, L.↗
DOE OSTI2023

Adaptive changes in sexual behavior in the high–risk population in response to human monkeypox transmission in Canada can help control the outbreak: Insights from a two–group, two–route epidemic model

Monkeypox, a zoonotic disease, is emerging as a potential sexually transmitted infection/disease, with underlying transmission mechanisms still unclear. We devised a risk-structured, compartmental model, incorporating sexual behavior dynamics. We compared different strategies targeting the high-risk population: a scenario of control policies geared toward the use of condoms and/or sexual abstinence (robust control strategy) with risk compensation behavior change, and a scenario of control strategies with behavior change in response to the doubling rate (adaptive control strategy). Monkeypox's basic reproduction number is 1.464, 0.0066, and 1.461 in the high-risk, low-risk, and total populations, respectively, with the high-risk group being the major driver of monkeypox spread. Policies imposing condom use or sexual abstinence need to achieve a 35% minimum compliance rate to stop further transmission, while a combination of both can curb the spread with 10% compliance to abstinence and 25% to condom use. With risk compensation, the only option is to impose sexual abstinence by at least 35%. Adaptive control is more effective than robust control where the daily sexual contact number is reduced proportionally and remains constant thereafter, shortening the time to epidemic peak, lowering its size, facilitating disease attenuation, and playing a key role in controlling the current outbreak.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗
DOE OSTI2024

GP-BayesOpInf

SAND2025-01851O GP-BayesOpInf is a software tool that uses algorithms to combine Gaussian process regression, principal component analysis, and linear Bayesian inference to produce a probabilistic reduced-order model for time-dependent systems. Numerical examples include the compressible Euler equations for an ideal gas, a heat diffusion process with a nonlinear reaction term, and a set of ordinary differential equations describing a compartmental model in epidemiology. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2025

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗
DOE OSTI2020

Radon kinetics in a natural indoor radon chamber

An unusual 180 m 3 storage room in the basement of a two-story laboratory building is unventilated, and separated from occupiable rooms by double steel doors. The space completely borders on soil through the concrete floor and two of its concrete walls. The room also contains a separate inner chamber with 1 m thick concrete walls designed to damp vibrations in the room above it. The space boasts a relatively high radon level, 1083 Bq m –3 , which varies with local outdoor environmental conditions. Measurements were made of radon concentrations at various locations and heights within the facility. More than a year of continuous radon concentration data corresponding to a single location are also available, along with measurements of indoor and outdoor pressure, temperature, and humidity. Data were also collected with as many as five fans placed in different locations and cycled on for variable time periods. First order linear kinetic models were created to explain the observed approaches to steady state due to changing conditions and wash-out resulting from intentional ventilation. Results demonstrate a good fit between changes in the radon concentration level and the developed compartmental models. Furthermore, no significant differences were observed between radon concentration at different locations or heights in the chamber.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Compartmentalized Storage Reservoir Model and Simulation Results

A hypothetical closed (or compartmentalized) reservoir was created based on the numerical mesh used for Wainwright et al. (2013). The reservoir extent is smaller than the original model and limited between the two faults. The caprock properties are assumed in the rest of the storage formation.

Carbon Storage↗
DOE OSTI2024

Modeling with uncertainty quantification reveals the essentials of a non-canonical algal carbon-concentrating mechanism

The thermoacidophilic red alga Cyanidioschyzon merolae survives its challenging environment likely in part by operating a carbon-concentrating mechanism (CCM). Here, we demonstrated that C. merolae 's cellular affinity for CO 2 is stronger than the affinity of its rubisco for CO 2 . This finding provided additional evidence that C. merolae operates a CCM while lacking the structures and functions characteristic of CCMs in other organisms. To test how such a CCM could function, we created a mathematical compartmental model of a simple CCM, distinct from those we have seen previously described in detail. The results of our modeling supported the feasibility of this proposed minimal and non-canonical CCM in C. merolae . To facilitate the robust modeling of this process, we measured and incorporated physiological and enzymatic parameters into the model. Additionally, we trained a surrogate machine-learning model to emulate the mechanistic model and characterized the effects of model parameters on key outputs. This parameter exploration enabled us to identify model features that influenced whether the model met the experimentally derived criteria for functional carbon concentration and efficient energy usage. Such parameters included cytosolic pH, bicarbonate pumping cost and kinetics, cell radius, carboxylation velocity, number of thylakoid membranes, and CO 2 membrane permeability. Our exploration thus suggested that a non-canonical CCM could exist in C. merolae and illuminated the essential features generally necessary for CCMs to function.

Steensma, Anne K. [Michigan State Univ., East Lans↗
NASA NTRS2014

Pharmacokinetic Modeling of Intranasal Scopolamine in Plasma Saliva and Urine

An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials for an Investigative New Drug (IND). The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial with INSCOP. METHODS: Twelve healthy human subjects were administered three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min to 24 h after dosing and scopolamine concentrations measured by using a validated LC-MS-MS assay. Pharmacokinetic Compartmental models, using actual dosing and sampling times, were built using Phoenix (version 1.2). Model discrimination was performed, by minimizing the Akaike Information Criteria (AIC), maximizing the coefficient of determination (r²) and by comparison of the quality of fit plots. RESULTS: The best structural model to describe scopolamine disposition after INSCOP administration (minimal AIC =907.2) consisted of one compartment for plasma, saliva and urine respectively that were inter-connected with different rate constants. The estimated values of PK parameters were compiled in Table 1. The model fitting exercises revealed a nonlinear PK for scopolamine between plasma and saliva compartments for K21, Vmax and Km. CONCLUSION: PK model for INSCOP was developed and for the first time it satisfactorily predicted the PK of scopolamine in plasma, saliva and urine after INSCOP administration. Using non-linear PK yielded the best structural model to describe scopolamine disposition between plasma and saliva compartments, and inclusion of non-linear PK resulted in a significant improved model fitting. The model can be utilized to predict scopolamine plasma concentration using saliva and/or urine data that allows non-invasive assessment of pharmacotherapeutics of scopolamine in space and other remote environments without requiring blood sampling.

Wu, L.↗
NASA NTRS2015

Pharmacokinetic Modeling of Intranasal Scopolamine in Plasma Saliva and Urine

An intranasal gel dosage formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness (SMS). The bioavailability and pharmacokinetics (PK) were evaluated under IND (Investigational New Drug) guidelines. The aim of the project was to develop a PK model that can predict the relationships among plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial protocol with INSCOP. Twelve healthy human subjects were administered at three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min to 24 h after dosing and scopolamine concentrations were measured by using a validated LC‐MS‐MS assay. PK compartmental models, using actual dosing and sampling time, were established using Phoenix (version 1.2). Model selection was based on a likelihood ratio test on the difference of criteria (‐2LL (i.e. log-likelihood ratio test)) and comparison of the quality of fit plots. The results: Predictable correlations among scopolamine concentrations in compartments of plasma, saliva and urine were established, and for the first time the model satisfactorily predicted the population and individual PK of INSCOP in plasma, saliva and urine. The model can be utilized to predict the INSCOP plasma concentration by saliva and urine data, and it will be useful for monitoring the PK of scopolamine in space and other remote environments using non‐invasive sampling of saliva and/or urine.

Wu, L.↗
NASA NTRS1999

Modeling Calcium Loss from Bones During Space Flight

Calcium loss from bones during space flight creates a risk for astronauts who travel into space, and may prohibit space flights to other planets. The problem of calcium loss during space flight has been studied using animal models, bed rest (as a ground-based model), and humans in-flight. In-flight studies have typically documented bone loss by comparing bone mass before and after flight. To identify changes in metabolism leading to bone loss, we have performed kinetic studies using stable isotopes of calcium. Oral (Ca-43) and intravenous (Ca-46) tracers were administered to subjects (n=3), three-times before flight, once in-flight (after 110 days), and three times post-flight (on landing day, and 9 days and 3 months after flight). Samples of blood, saliva, urine, and feces were collected for up to 5 days after isotope administration, and were analyzed for tracer enrichment. Tracer data in tissues were analyzed using a compartmental model for calcium metabolism and the WinSAAM software. The model was used to: account for carryover of tracer between studies, fit data for all studies using the minimal number of changes between studies, and calculate calcium absorption, excretion, bone calcium deposition and bone calcium resorption. Results showed that fractional absorption decreased by 50% during flight and that bone resorption and urinary excretion increased by 50%. Results were supported by changes in biochemical markers of bone metabolism. Inflight bone loss of approximately 250 mg Ca/d resulted from decreased calcium absorption combined with increased bone resorption and excretion. Further studies will assess the time course of these changes during flight, and the effectiveness of countermeasures to mitigate flight-induced bone loss. The overall goal is to enable human travel beyond low-Earth orbit, and to allow for better understanding and treatment of bone diseases on Earth.

Wastney, Meryl E.↗
DOE OSTI2024

Impacts of Vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 Variants Alpha and Delta on Coronavirus Disease 2019 Transmission Dynamics in Four Metropolitan Areas of the United States

To characterize Coronavirus Disease 2019 (COVID-19) transmission dynamics in each of the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix in 2020 and 2021, we extended a previously reported compartmental model accounting for effects of multiple distinct periods of non-pharmaceutical interventions by adding consideration of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2). For each MSA, we found region-specific parameterizations of the model using daily reports of new COVID-19 cases available from January 21, 2020 to October 31, 2021. In the process, we obtained estimates of the relative infectiousness of Alpha and Delta as well as their takeoff times in each MSA (the times at which sustained transmission began). The estimated infectiousness of Alpha ranged from 1.1x to 1.4x that of viral strains circulating in 2020 and early 2021. The estimated relative infectiousness of Delta was higher in all cases, ranging from 1.6x to 2.1x. The estimated Alpha takeoff times ranged from February 1 to February 28, 2021. The estimated Delta takeoff times ranged from June 2 to June 26, 2021. In conclusion, estimated takeoff times are consistent with genomic surveillance data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Differential contagiousness of respiratory disease across the United States

The initial contagiousness of a communicable disease within a given population is quantified by the basic reproduction number, $R_0$. This number depends on both pathogen and population properties. On the basis of compartmental models that reproduce Coronavirus Disease 2019 (COVID-19) surveillance data, we used Bayesian inference and the next-generation matrix approach to estimate region-specific $R_0$ values for 280 of 384 metropolitan statistical areas (MSAs) in the United States (US), which account for 95% of the US population living in urban areas and 82% of the total population. We focused on MSA populations after finding that these populations were more uniformly impacted by COVID-19 than state populations. Our maximum a posteriori (MAP) estimates for $R_0$ range from 1.9 to 7.7 and quantify the relative susceptibilities of regional populations to spread of respiratory diseases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Radon Kinetics in a Basement Space Measured with Different Devices

Although indoor monitoring of radon and benchmarking of radon measurement devices remain important research topics, few intercomparisons of active radon measurement devices have been performed under realistic conditions, let alone dynamic ones enabling comparison of their transient behavior. Five different radon monitors were therefore placed in a poorly ventilated basement space under three different conditions: 24 h under a steady, elevated radon level, 24 h with fans turned on to produce a radon washout transient, and 9 d with fans turned off for a radon buildup transient. Resulting radon concentrations varied between ~200 and ~2,000 Bq m -3 . Accuracy of the devices were evaluated using root-mean-square error, and ventilation data were fit to first order linear compartmental models. To more accurately model behaviors such as cyclic diurnal variations, the source term corresponding to entry of radon from soil into the basement was considered to be non-constant, as it is likely to vary drastically with both the indoor-outdoor pressure differential and soil concentration variations. Here, the improved radon washout model fit very well with the measurements. Despite a wide variety in list prices, all devices performed similarly during transients and at different radon concentrations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah)

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Bayesian Inference of State-Level COVID-19 Basic Reproduction Numbers across the United States

Although many persons in the United States have acquired immunity to COVID-19, either through vaccination or infection with SARS-CoV-2, COVID-19 will pose an ongoing threat to non-immune persons so long as disease transmission continues. We can estimate when sustained disease transmission will end in a population by calculating the population-specific basic reproduction number ℛ0, the expected number of secondary cases generated by an infected person in the absence of any interventions. The value of ℛ0 relates to a herd immunity threshold (HIT), which is given by 1−1/ℛ0. When the immune fraction of a population exceeds this threshold, sustained disease transmission becomes exponentially unlikely (barring mutations allowing SARS-CoV-2 to escape immunity). Here, we report state-level ℛ0 estimates obtained using Bayesian inference. Maximum a posteriori estimates range from 7.1 for New Jersey to 2.3 for Wyoming, indicating that disease transmission varies considerably across states and that reaching herd immunity will be more difficult in some states than others. ℛ0 estimates were obtained from compartmental models via the next-generation matrix approach after each model was parameterized using regional daily confirmed case reports of COVID-19 from 21 January 2020 to 21 June 2020. Our ℛ0 estimates characterize the infectiousness of ancestral strains, but they can be used to determine HITs for a distinct, currently dominant circulating strain, such as SARS-CoV-2 variant Delta (lineage B.1.617.2), if the relative infectiousness of the strain can be ascertained. On the basis of Delta-adjusted HITs, vaccination data, and seroprevalence survey data, we found that no state had achieved herd immunity as of 20 September 2021.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1971

A program for identification of linear systems

A program has been written for the identification of parameters in certain linear systems. These systems appear in biomedical problems, particularly in compartmental models of pharmacokinetics. The method presented here assumes that some of the state variables are regularly modified by jump conditions. This simulates administration of drugs following some prescribed drug regime. Parameters are identified by a least-square fit of the linear differential system to a set of experimental observations. The method is especially suited when the interval of observation of the system is very long.

Buell, J.↗