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 55 records · Page 3

Cellular and Network Mechanisms Underlying Information Processing in a Simple Sensory System

Realistic, biophysically-based compartmental models were constructed of several primary sensory interneurons in the cricket cercal sensory system. A dynamic atlas of the afferent input to these cells was used to set spatio-temporal parameters for the simulated stimulus-dependent synaptic inputs. We examined the roles of dendritic morphology, passive membrane properties, and active conductances on the frequency tuning of the neurons. The sensitivity of narrow-band low pass interneurons could be explained entirely by the electronic structure of the dendritic arbors and the dynamic sensitivity of the SIZ. The dynamic characteristics of interneurons with higher frequency sensitivity required models with voltage-dependent dendritic conductances.

Jacobs, Gwen↗

Calcium kinetics with microgram stable isotope doses and saliva sampling

Studies of calcium kinetics require administration of tracer doses of calcium and subsequent repeated sampling of biological fluids. This study was designed to develop techniques that would allow estimation of calcium kinetics by using small (micrograms) doses of isotopes instead of the more common large (mg) doses to minimize tracer perturbation of the system and reduce cost, and to explore the use of saliva sampling as an alternative to blood sampling. Subjects received an oral dose (133 micrograms) of 43Ca and an i.v. dose (7.7 micrograms) of 46Ca. Isotopic enrichment in blood, urine, saliva and feces was well above thermal ionization mass spectrometry measurement precision up to 170 h after dosing. Fractional calcium absorptions determined from isotopic ratios in blood, urine and saliva were similar. Compartmental modeling revealed that kinetic parameters determined from serum or saliva data were similar, decreasing the necessity for blood samples. It is concluded from these results that calcium kinetics can be assessed with micrograms doses of stable isotopes, thereby reducing tracer costs and with saliva samples, thereby reducing the amount of blood needed.

Clinical Trial↗

Evidence of a sensory processing unit in the mammalian macula

We cut serial sections through the medial part of the rat vestibular macula for transmission electron microscopic (TEM) examination, computer-assisted 3-D reconstruction, and compartmental modeling. The ultrastructural research showed that many primary vestibular neurons have an unmyelinated segment, often branched, that extends between the heminode (putative site of the spike initiation zone) and the expanded terminal(s) (calyx, calyces). These segments, termed the neuron branches, and the calyces frequently have spine-like processes of various dimensions with bouton endings that morphologically are afferent, efferent, or reciprocal to other macular neural elements. The major questions posed by this study were whether small details of morphology, such as the size and location of neuronal processes or synapses, could influence the output of a vestibular afferent, and whether a knowledge of morphological details could guide the selection of values for simulation parameters. The conclusions from our simulations are (1) values of 5.0 k omega cm2 for membrane resistivity and 1.0 nS for synaptic conductance yield simulations that best match published physiological results; (2) process morphology has little effect on orthodromic spread of depolarization from the head (bouton) to the spike initiation zone (SIZ); (3) process morphology has no effect on antidromic spread of depolarization to the process head; (4) synapses do not sum linearly; (5) synapses are electrically close to the SIZ; and (6) all whole-cell simulations should be run with an active SIZ.

NASA Center ARC↗

Physics-Informed Machine Learning for Epidemiological Models

One challenge of using compartmental SEIR models for public health planning is the difficulty in manually tuning parameters to capture behavior reflected in the real-world data. This team conducted initial, exploratory analysis of a novel technique to use physics-informed machine learning tools to rapidly develop data-driven models for physical systems. This machine learning approach may be used to perform data assimilation of compartment models which account for unknown interactions between geospatial domains (i.e. diffusion processes coupling across neighborhoods/counties/states/etc.). Results presented here are early, proof-of-concept ideas that demonstrate initial success in using a physically informed neural network (PINN) model to assimilate data in a compartmental epidemiology model. The results demonstrate initial success and warrant further research and development.

60 APPLIED LIFE SCIENCES↗

Adaptive Recovery Model: Designing Systems for Testing Tracing and Vaccination to Support COVID-19 Recovery Planning.

This report documents a new approach to designing disease control policies that allocate scarce testing, contact tracing, and vaccination resources to better control community transmission of COVID19 or similar diseases. The Adaptive Recovery Model (ARM) combines a deterministic compartmental disease model with a stochastic network disease propagation model to enable us to simulate COVID-19 community spread through the lens of two complementary modeling motifs. ARM contact networks are derived from cell-phone location data that have been anonymized and interpreted as individual arrivals to specic public locations. Modeling disease spread over these networks allows us to identify locations within communities conducive to rapid disease spread. ARM applies this model- and data-derived abstractions of community transmission to evaluate the effectiveness of disease control measures including targeted social distancing, contact tracing, testing and vaccination. The architecture of ARM provides a unique capacity to help decision makers understand how best to deploy scarce testing, tracing and vaccination resources to minimize disease-spread potential in a community. This document details the novel mathematical formulations underlying ARM, presents a dynamical stability analysis of the deterministic model components, a sensitivity analysis of control parameters and network structure, and summarizes a process for deriving contact networks from cell-phone location data. An example use case steps through applying ARM to evaluate three targeted social distancing policies using Bernalillo County, New Mexico as an exemplar test locale. This step-by-step analysis demonstrates how ARM can be used to measure the relative performance of competing public health policies. Initial scenario tests of ARM shows that ARMs design focus on resource utilization rather than simple incidence prediction can provide decision makers with additional quantitative guidance for managing ongoing public health emergencies and planning future responses.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Modeling efficient and equitable distribution of COVID-19 vaccines

Producing and distributing COVID-19 vaccine during the pandemic is a major logistical challenge requiring careful planning and efficient execution. This report presents information on logistical, policy and technical issues relevant to rapidly fielding a COVID-19 vaccination program. For this study we (a) conducted literature review and subject matter expert elicitation to understand current vaccine manufacturing and distribution capabilities and vaccine allocation strategies, (b) designed a baseline vaccine distribution strategy and modeling strategy to provide insight into the potential for targeted distribution of limited initial vaccine supplies, and (c) developed parametric interfaces to enable vaccine distribution scenarios to be analyzed in depth with Sandias Adaptive Recovery Model that will allow us evaluate the additional sub- populations and alternative distribution scenarios from a public health benefit and associated economic disruption Principal issues, challenges, and complexities that complicate COVID-19 vaccine delivery identified in our literature and subject matter expert investigation include these items: The United States has not mounted an urgent nationwide vaccination campaign in recent history. The existing global manufacturing and distribution infrastructure are not able to produce enough vaccine for the population immediately. Vaccines, once available will be scarce resources. Prioritization for vaccine allocation will be built on existing distribution networks. Vaccine distribution may not have a universal impact on disease transmission and morbidity because of scarcity, priority population demographics, and underlying disease transmission rates. Considerations for designing a vaccine distribution strategy are discussed. A baseline distribution strategy is designed and tested using the Adaptive Recovery Model, which couples a deterministic compartmental epidemiological model and a stochastic network model. We show the impact of this vaccine distribution strategy on hospitalizations, mortality, and contact tracing requirements. This model can be used to quantitatively evaluate alternative distribution scenarios, guiding policy decisions as vaccine candidates are narrowed down.

59 BASIC BIOLOGICAL SCIENCES↗

Community Modeling and Flux Balance Analysis of Synthetic Bacterial Community in MOPS Medium

This narrative was used for the 4-member community modeling using MOPS medium for the bacteria isolated from Populus deltoides (PD). 10 bacterial strains isolated from PD were co-cultured as a synthetic microbial community in MOPS minimal medium. After 15 passages there were 4 members survived in MOPS medium. The genome of each community member was sequenced and uploaded into KBase for genome annotation and metabolic modeling. To simulate the metabolic interactions among the finally surviving community members, the individual metabolic models of each strain were merged into a compartmentalized community model, and then the community model was gapfilled using MOPS medium. The FBA modeling was performed and the metabolites exchange among community members was proposed. The publication by Jia Wang, Dana L. Carper, Leah H. Burdick, Him Shrestha, Manasa Appidi, Paul E. Abraham, Collin M. Timm, Robert L. Hettich, Dale A. Pelletier, Mitchel J. Doktycz can be found here: https://doi.org/10.​1016/​j.​csbj.​2021.​03.​034

59 BASIC BIOLOGICAL SCIENCES↗

Community Modeling and Flux Balance Analysis of Synthetic Bacterial Community in R2A Medium

This narrative was used for the 3-member community modeling using R2A medium for the bacteria isolated from Populus deltoides (PD). 10 bacteiral strains isolated from PD were co-cultured as a synthetic microbial community in R2A complex medium. After 15 passages, there were 3 members survived in R2A medium. The genome of each community member was sequenced and uploaded into KBase for genome annotation and metabolic modeling. To simulate the metabolic interactions among the finally surviving community members, the individual metabolic models of each strain were merged into a compartmentalized community model, and then the community model was gapfilled using R2A medium. The FBA modeling was performed and the metabolites exchange among community members was proposed. The publication by Jia Wang, Dana L. Carper, Leah H. Burdick, Him Shrestha, Manasa Appidi, Paul E. Abraham, Collin M. Timm, Robert L. Hettich, Dale A. Pelletier, Mitchel J. Doktycz can be found here: https://doi.org/10.​1016/​j.​csbj.​2021.​03.​034

59 BASIC BIOLOGICAL SCIENCES↗

Distributed Delay Differential Equation Representations of Cyclic Differential Equations

Compartmental ordinary differential equation (ODE) models are used extensively in mathematical biology. When transit between compartments occurs at a constant rate, the well-known linear chain trick can be used to show that the ODE model is equivalent to an Erlang distributed delay differential equation (DDE). Here, we demonstrate that compartmental models with nonlinear transit rates and possibly delayed arguments are also equivalent to a scalar distributed DDE. To illustrate the utility of these equivalences, we calculate the equilibria of the scalar DDE, and compute the characteristic function---without calculating a determinant. Further, we derive the equivalent scalar DDE for two examples of models in mathematical biology and use the DDE formulation to identify physiological processes that were otherwise hidden by the compartmental structure of the ODE model.

97 MATHEMATICS AND COMPUTING↗

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Design, Synthesis, and Validation: Genome Scale Optimization of Energy Flux through Compartmentalized Metabolic Networks in a Model Photosynthetic Eukaryotic Microbe (Final Report)

Photosynthetic organisms have recently gained considerable attention for a role in development of renewable energy sources. Genome-enabled systems biology methodology and modeling, coupled with high throughput genome engineering strategies, present opportunities to develop sustainable and economical applications such as fuel production within the next 10 to 15 years. However, optimization of light-driven metabolism for biomass or biofuel production will require significant advances in methodological throughput as well as improvements in detailed systems biology understanding of photosynthetic processes and cellular metabolism. The ability of diatoms to thrive in upwelling-induced, periodically nutrient-rich conditions makes them the base for the world’s shortest and most energy-efficient food webs. Diatom photosynthesis is estimated to account for between 25% and 40% of the 45-50 billion tons of organic carbon fixed annually in the sea.

60 APPLIED LIFE SCIENCES↗

Comparing modelling approaches for a generic nuclear waste repository in salt

This paper contains a comparison of five modelling approaches for a simplified nuclear waste repository in a domal salt formation. It is the result of a four-year collaboration between five international teams on Task F of the DECOVALEX-2023 project on performance assessment modelling. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. This work demonstrates how these objectives are accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Participating teams made a wide range of model assumptions, ranging from compartmentalized networks to full 3D models of the salt formation and repository. Despite differences in the modelling strategies, all models indicate that salt compaction and diffusion of radionuclides in brine are key processes in the repository. For the isothermal spent nuclear fuel and vitrified waste scenario with multiple early failures considered, all models indicate little of the disposed radionuclides will migrate beyond the repository seal over the 100,000-year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. Disparities between the models are believed to be due to differing simplifications from the conceptual model.

DECOVALEX↗

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗

Chromosome compartmentalization alterations in prostate cancer cell lines model disease progression

Prostate cancer aggressiveness and metastatic potential are influenced by gene expression and genomic aberrations, features that can be influenced by the 3D structure of chromosomes inside the nucleus. Using chromosome conformation capture (Hi-C), we conducted a systematic genome architecture comparison on a cohort of cell lines that model prostate cancer progression, from normal epithelium to bone metastasis. We describe spatial compartment identity (A-open versus B-closed) changes with progression in these cell lines and their relation to gene expression changes in both cell lines and patient samples. In particular, 48 gene clusters switch from the B to the A compartment, including androgen receptor, WNT5A, and CDK14. These switches are accompanied by changes in the structure, size, and boundaries of topologically associating domains (TADs). Further, compartment changes in chromosome 21 are exacerbated with progression and may explain, in part, the genesis of the TMPRSS2-ERG translocation. These results suggest that discrete 3D genome structure changes play a deleterious role in prostate cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗

Face Recognition Oak Ridge (FaRO): A Framework for Distributed and Scalable Biometrics Applications

The facial biometrics community has seen a recent abundance of high-accuracy facial analytic models become freely available. Although these models' capabilities in facial detection, landmark detection, attribute analysis, and recognition are ever-increasing, they aren't always straightforward to deploy in a real-world environment. In reality, the use of the field's ever growing collection of models is becoming exceedingly difficult as library dependencies update and deprecate. Researchers often encounter headaches when attempting to utilize multiple models requiring different or conflicting software packages. Face Recognition Oak Ridge (FaRO) is an open-source project designed to provide a highly modular, flexible framework for unifying facial analytic models through a compartmentalized plug-and-play paradigm built on top of the gRPC (Google Remote Procedure Call) protocol. FaRO's server-client architecture and flexible portability allows easy construction of modularized and heterogeneous face analysis pipelines, distributed over many machines with differing hardware and software resources. This paper outlines FaRO's architecture and current capabilities, along with some experiments in model testing and distributed scaling through FaRO.

Bolme, David↗

Computational Modeling of Cephalad Fluid Shift for Application to Microgravity-Induced Visual Impairment

An improved understanding of spaceflight-induced ocular pathology, including the loss of visual acuity, globe flattening, optic disk edema and distension of the optic nerve and optic nerve sheath, is of keen interest to space medicine. Cephalad fluid shift causes a profoundly altered distribution of fluid within the compartments of the head and body, and may indirectly generate phenomena that are biomechanically relevant to visual function, such as choroidal engorgement, compromised drainage of blood and cerebrospinal fluid (CSF), and altered translaminar pressure gradient posterior to the eye. The experimental body of evidence with respect to the consequences of fluid shift has not yet been able to provide a definitive picture of the sequence of events. On earth, elevated intracranial pressure (ICP) is associated with idiopathic intracranial hypertension (IIH), which can produce ocular pathologies that look similar to those seen in some astronauts returning from long-duration flight. However, the clinically observable features of the Visual Impairment and Intracranial Pressure (VIIP) syndrome in space and IIH on earth are not entirely consistent. Moreover, there are at present no experimental measurements of ICP in microgravity. By its very nature, physiological measurements in spaceflight are sparse, and the space environment does not lend itself to well-controlled experiments. In the absence of such data, numerical modeling can play a role in the investigation of biomechanical causal pathways that are suspected of involvement in VIIP. In this work, we describe the conceptual framework for modeling the altered compartmental fluid distribution that represents an equilibrium fluid distribution resulting from the loss of hydrostatic pressure gradient.

Nelson, Emily S.↗