Search NASA⌕ Search

SEARCH · Search NASA

Results for “mathematical model”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Leveraging public AI tools to explore systems biology resources in mathematical modeling

Predictive mathematical modeling is an essential part of systems biology and is interconnected with information management. Systems biology information is often stored in specialized formats to facilitate data storage and analysis. These formats are not designed for easy human readability and thus require specialized software to visualize and interpret results. Therefore, comprehending modeling and underlying networks and pathways is contingent on mastering systems biology tools, which is particularly challenging for users with no or little background in data science or system biology. To address this challenge, we investigated the usage of public Artificial Intelligence (AI) tools in exploring systems biology resources in mathematical modeling. We tested public AI’s understanding of mathematics in models, related systems biology data, and the complexity of model structures. Our approach can enhance the accessibility of systems biology for non-system biologists and help them understand systems biology without a deep learning curve.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Bioelectrochemical reactor to manage anthropogenic sulfate pollution for freshwater ecosystems: Mathematical modeling and experimental validation

Anthropogenic sulfate loading into otherwise low-sulfate freshwater systems can cause significant ecological consequences as a biogeochemical stressor. To address this challenge, in situ bioremediation technologies have been developed to leverage naturally occurring microorganisms that transform sulfate into sulfide rather than implementing resource-intensive physio-chemical processes. However, bioremediation technologies often require the supply of electron donors to facilitate biological sulfate reduction. Bioelectrochemical systems (BES) can be an alternative approach for supplying molecular hydrogen as an electron donor for sulfate-reducing bacteria through water electrolysis. Although the fundamental mechanisms behind BESs have been studied, limited research has evaluated the design and operational parameters of treatment systems when developing BESs on a scale relevant to environmental systems. Further, this study aimed to develop an application-based mathematical model to evaluate the performance of BESs across a range of reactor configurations and operational modes. The model was based on sulfate transformation by hydrogenotrophic sulfate-reducing bacteria coupled with the recovery of solid iron sulfide species formed by the oxidative dissolution of dissolved ferrous iron from a stainless steel anode. Sulfate removal closely corresponded to the rate of electrolytic hydrogen production and hydraulic residence time but was less sensitive to specific microbial rate constants. The mathematical model results were compared to experimental data from a pilot-scale BES tested with nonacidic mine drainage as a case study. The close agreement between the mathematical model and the pilot-scale BES experiment highlights the efficacy of using a mathematical model as a tool to develop a conceptual design of a scaled-up treatment system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanistic within-host mathematical model of inhalational anthrax

We present a mathematical model of the dynamics of Bacillus anthracis bacteria within the lymph nodes and blood of a host, following inhalation of an initial dose of spores. We also incorporate the dynamics of protective antigen, which is the binding component of the anthrax toxin produced by the bacteria. The model offers a mechanistic description of the early infection dynamics of inhalational anthrax, while its stochastic nature allows us to study the probabilities of different outcomes (for example, how likely it is that the infection will be cleared for a given inhaled dose of spores) in order to explain dose-response data for inhalational anthrax. The model is calibrated via a Bayesian approach, using in vivo data from New Zealand white rabbit and guinea pig infection studies, enabling within-host parameters to be estimated. We also leverage incubation-period data from the Sverdlovsk 1979 anthrax outbreak to show that the model can accurately describe human time-to-symptoms data under reasonable parameter regimes. Finally, we derive a simple approximate formula for the probability of symptom onset before time t, assuming that the number of inhaled spores has a Poisson distribution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Mathematical Modeling of the Potential and Time Dependence of Ir Dissolution from Hydrous Ir Oxide Oxygen Evolution Catalysts

One of the main degradation mechanisms of hydrous iridium oxide acidic oxygen evolution reaction (OER) catalysts is dissolution and loss into the acidic membrane. While degradation models have been proposed, there is a gap in understanding the potential and time dependence of the iridium dissolution reaction and its mechanistic underpinnings. In this work, Ir dissolution rates measured as a function of time and potential via time-resolved inductively-coupled plasma mass spectrometry (ICP-MS) in aqueous acidic electrolyte are used to establish a mathematical model for Ir dissolution. The mathematical model is generated using proposed formation and dissolution reactions for Ir species. Through comparison with the ICP-MS data and existing information on the potential-dependent Ir phase, we find that the potential and time-dependence of dissolution can be modeled as dissolution of an oxide phase, here represented as IrO 2 , with potential dependent kinetics and formation of a passivating species, a process with a rate-limiting step that is not potential dependent. This understanding of the potential dependence of dissolution and passivation kinetics using aqueous electrolyte half-cell measurements can be used to predict the degradation of Ir oxide in operating energy conversion devices relying on the OER, such as proton-exchange membrane water electrolyzers.

Kariuki, Nancy N. [Argonne National Laboratory (AN↗
DOE OSTI2024

Machine learning mathematical models for incidence estimation during pandemics

Accurate estimates of the incidence of infectious diseases are key for the control of epidemics. However, healthcare systems are often unable to test the population exhaustively, especially when asymptomatic and paucisymptomatic cases are widespread; this leads to significant and systematic under-reporting of the real incidence. Here, we propose a machine learning approach to estimate the incidence of a pandemic in real-time, using reported cases and the overall test rate. In particular, we use Bayesian symbolic regression to automatically learn the closed-form mathematical models that most parsimoniously describe incidence. We develop and validate our models using COVID-19 incidence values for nine different countries, confirming their ability to accurately predict daily incidence. Remarkably, despite the differences in epidemic trajectories and dynamics across countries, we find that a single model for all countries offers a more parsimonious description and is more predictive of actual incidence compared to separate models for each country. Our results show the potential to accurately model incidence in real-time using closed-form mathematical models, providing a valuable tool for public health decision-makers.

Fajardo-Fontiveros, Oscar (ORCID:0000000207058972)↗
DOE OSTI2023

On Soft Mathematical Models of Subjective Time Acceleration with Age

It is a commonplace perception that speed of time subjectively experienced by humans significantly differs from chronological (objective) time and shows a great deal of variability. An often cited example is the phenomenon of the time acceleration with age - subjectively, the time passes faster as we get older. While the exact mechanisms behind it are not yet fully established, here we consider three 'soft' (conceptual) mathematical models that might be applicable to the speeding time phenomenon: two proportionality theories widely discussed in the past and the original model that takes into account the novelty of experience effect. Furthermore, the latter is found the most plausible, as not only it satisfactorily describes the decadal subjective time acceleration, but also offers a reasonable explanation of the human life experience accumulation with age.

acceleration↗
DOE OSTI2023

Ionization detail parameters and cluster dose: a mathematical model for selection of nanodosimetric quantities for use in treatment planning in charged particle radiotherapy

Abstract Objective . To propose a mathematical model for applying ionization detail (ID), the detailed spatial distribution of ionization along a particle track, to proton and ion beam radiotherapy treatment planning (RTP). Approach . Our model provides for selection of preferred ID parameters ( I p ) for RTP, that associate closest to biological effects. Cluster dose is proposed to bridge the large gap between nanoscopic I p and macroscopic RTP. Selection of I p is demonstrated using published cell survival measurements for protons through argon, comparing results for nineteen I p : N k , k = 2, 3, …, 10, the number of ionizations in clusters of k or more per particle, and F k , k = 1, 2, …, 10, the number of clusters of k or more per particle. We then describe application of the model to ID-based RTP and propose a path to clinical translation. Main results . The preferred I p were N 4 and F 5 for aerobic cells, N 5 and F 7 for hypoxic cells. Significant differences were found in cell survival for beams having the same LET or the preferred N k . Conversely, there was no significant difference for F 5 for aerobic cells and F 7 for hypoxic cells, regardless of ion beam atomic number or energy. Further, cells irradiated with the same cluster dose for these I p had the same cell survival. Based on these preliminary results and other compelling results in nanodosimetry, it is reasonable to assert that I p exist that are more closely associated with biological effects than current LET-based approaches and microdosimetric RBE-based models used in particle RTP. However, more biological variables such as cell line and cycle phase, as well as ion beam pulse structure and rate still need investigation. Significance . Our model provides a practical means to select preferred I p from radiobiological data, and to convert I p to the macroscopic cluster dose for particle RTP.

Engineering↗
DOE OSTI2025

A consensus mathematical model of vaccine-induced antibody dynamics for multiple vaccine platforms and pathogens

Introduction: Vaccine platforms used in successful, licensed vaccines have varied among pathogens. However, antibody level is still the main clinical correlate of protection in most approved vaccines. Decisions as to the best vaccine platform to pursue for a given pathogen may be informed through improved understanding of the process of antibody generation and its temporal dynamics, as well as the relationship between these processes and the type of vaccine. Methods: We have analyzed the dynamics of antibody generation for different vaccine platforms against diverse pathogens, and developed a consensus mathematical model that captures antibody dynamics across these diverse systems. Initially, the model was fitted to a rich dataset of antibody and immune cell concentrations in a SARS-CoV-2 vaccine experiment. We then used concepts from machine learning, such as transfer learning, to apply the same model to a variety of systems, involving different pathogens, vaccine platforms, and booster dose use/timing, fixing most parameter values relating to the dynamics of the immune system. Results: The model includes B cell proliferation and differentiation, as well as the generation of plasma cells, which secrete large amounts of antibody, and memory B cells. Overall, the model describes antibody generation in all systems tested well and shows that the main differences across platforms are related to the dynamics of antigen presentation. Discussion: This model can be used to predict antibody generation in pairs of vaccine platform/pathogen, allowing for the use of in silico results to narrow down experimental burden in vaccine development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways

Post-translational modifications (PTMs) are key regulatory mechanisms that can control protein function. Of these, phosphorylation is the most common and widely studied. Because of its importance in regulating cell signaling, precise and accurate measurements of protein phosphorylation across wide dynamic ranges are crucial to understanding how signaling pathways function. Although immunological assays are commonly used to detect phosphoproteins, their lack of sensitivity, specificity, and selectivity often make them unreliable for quantitative measurements of complex biological samples. Recent advances in Mass Spectrometry (MS)-based targeted proteomics have made it a more useful approach than immunoassays for studying the dynamics of protein phosphorylation. Selected reaction monitoring (SRM)—also known as multiple reaction monitoring (MRM)—and parallel reaction monitoring (PRM) can quantify relative and absolute abundances of protein phosphorylation in multiplexed fashions targeting specific pathways. In addition, the refinement of these tools by enrichment and fractionation strategies has improved measurement of phosphorylation of low-abundance proteins. The quantitative data generated are particularly useful for building and parameterizing mathematical models of complex phospho-signaling pathways. Potentially, these models can provide a framework for linking analytical measurements of clinical samples to better diagnosis and treatment of disease.

mathematical modeling↗
DOE OSTI2023

Mathematical modeling of plus-strand RNA virus replication to identify broad-spectrum antiviral treatment strategies

Plus-strand RNA viruses are the largest group of viruses. Many are human pathogens that inflict a socio-economic burden. Interestingly, plus-strand RNA viruses share remarkable similarities in their replication. A hallmark of plus-strand RNA viruses is the remodeling of intracellular membranes to establish replication organelles (so-called “replication factories”), which provide a protected environment for the replicase complex, consisting of the viral genome and proteins necessary for viral RNA synthesis. In the current study, we investigate pan-viral similarities and virus-specific differences in the life cycle of this highly relevant group of viruses. We first measured the kinetics of viral RNA, viral protein, and infectious virus particle production of hepatitis C virus (HCV), dengue virus (DENV), and coxsackievirus B3 (CVB3) in the immuno-compromised Huh7 cell line and thus without perturbations by an intrinsic immune response. Based on these measurements, we developed a detailed mathematical model of the replication of HCV, DENV, and CVB3 and showed that only small virus-specific changes in the model were necessary to describe the in vitro dynamics of the different viruses. Our model correctly predicted virus-specific mechanisms such as host cell translation shut off and different kinetics of replication organelles. Further, our model suggests that the ability to suppress or shut down host cell mRNA translation may be a key factor for in vitro replication efficiency, which may determine acute self-limited or chronic infection. We further analyzed potential broad-spectrum antiviral treatment options in silico and found that targeting viral RNA translation, such as polyprotein cleavage and viral RNA synthesis, may be the most promising drug targets for all plus-strand RNA viruses. Moreover, we found that targeting only the formation of replicase complexes did not stop the in vitro viral replication early in infection, while inhibiting intracellular trafficking processes may even lead to amplified viral growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Human–mouse comparison of the multistage nature of radiation carcinogenesis in a mathematical model

Abstract Mouse models are vital for assessing risk from environmental carcinogens, including ionizing radiation, yet the interspecies difference in the dose response precludes direct application of experimental evidence to humans. Herein, we take a mathematical approach to delineate the mechanism underlying the human–mouse difference in radiation‐related cancer risk. We used a multistage carcinogenesis model assuming a mutational action of radiation to analyze previous data on cancer mortality in the Japanese atomic bomb survivors and in lifespan mouse experiments. Theoretically, the model predicted that exposure will chronologically shift the age‐related increase in cancer risk forward by a period corresponding to the time in which the spontaneous mutational process generates the same mutational burden as that the exposure generates. This model appropriately fitted both human and mouse data and suggested a linear dose response for the time shift. The effect per dose decreased with increasing age at exposure similarly between humans and mice on a per‐lifespan basis (0.72‐ and 0.71‐fold, respectively, for every tenth lifetime). The time shift per dose was larger by two orders of magnitude in humans (7.8 and 0.046 years per Gy for humans and mice, respectively, when exposed at ~35% of their lifetime). The difference was mostly explained by the two orders of magnitude difference in spontaneous somatic mutation rates between the species plus the species‐independent radiation‐induced mutation rate. Thus, the findings delineate the mechanism underlying the interspecies difference in radiation‐associated cancer mortality and may lead to the use of experimental evidence for risk prediction in humans.

Oncology↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

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

Analysis of tumor-immune functional responses in a mathematical model of neoantigen cancer vaccines

Cancer neoantigen vaccines have emerged as a promising approach to stimulating the immune system to fight cancer. We propose a simple model including key elements of cancer-immune interactions and conduct a phase plane analysis to understand the immunological mechanisms of cancer neoantigen vaccines. Analytical results are obtained for two widely used functional forms that represent the killing rate of tumor cells by immune cells: the law of mass action (LMA) and the dePillis-Radunskaya Law (LPR). Using the LMA, our results reveal that a slowly growing tumor can escape the immune surveillance and that there is a unique periodic solution. The LPR offers richer dynamics, in which tumor elimination and uncontrolled tumor growth are both present. We show that tumor elimination requires sufficient number of initial activated T cells in relationship to the malignant cells, which lends support to using the neoantigen cancer vaccine as an adjuvant therapy after the primary tumor is surgically removed or treated using radiotherapy. We also derive a sufficient condition for uncontrolled tumor growth under the assumption of the LPR. Here, the juxtaposition of analyses with these two different choices for the killing rate function highlights their importance on model behavior and biological implications, by which we hope to spur further theoretical and experimental work to understand mechanisms underlying different functional forms for the killing rate.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗