Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mathematical 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

New insights into RAS biology reinvigorate interest in mathematical modeling of RAS signaling

RAS is the most frequently mutated gene across human cancers, but developing inhibitors of mutant RAS has proven to be challenging. Given the difficulties of targeting RAS directly, drugs that impact the other components of pathways where mutant RAS operates may potentially be effective. However, the system-level features, including different localizations of RAS isoforms, competition between downstream effectors, and interlocking feedback and feed-forward loops, must be understood to fully grasp the opportunities and limitations of inhibiting specific targets. Mathematical modeling can help us discern the system-level impacts of these features in normal and cancer cells. New technologies enable the acquisition of experimental data that will facilitate development of realistic models of oncogenic RAS behavior. In light of the wealth of empirical data accumulated over decades of study and the advancement of experimental methods for gathering new data, modelers now have the opportunity to advance progress toward realization of targeted treatment for mutant RAS-driven cancers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Mathematical model of the metal mould surface temperature optimization

The article is focused on the problem of generating a uniform temperature field on the inner surface of shell metal moulds. Such moulds are used e.g. in the automotive industry for artificial leather production. To produce artificial leather with uniform surface structure and colour shade the temperature on the inner surface of the mould has to be as homogeneous as possible. The heating of the mould is realized by infrared heaters located above the outer mould surface. The conceived mathematical model allows us to optimize the locations of infrared heaters over the mould, so that approximately uniform heat radiation intensity is generated. A version of differential evolution algorithm programmed in Matlab development environment was created by the authors for the optimization process. For temperate calculations software system ANSYS was used. A practical example of optimization of heaters locations and calculation of the temperature of the mould is included at the end of the article.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

A Coupled Mathematical Model of the Intracellular Replication of Dengue Virus and the Host Cell Immune Response to Infection

Dengue virus (DV) is a positive-strand RNA virus of the Flavivirus genus. It is one of the most prevalent mosquito-borne viruses, infecting globally 390 million individuals per year. The clinical spectrum of DV infection ranges from an asymptomatic course to severe complications such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS), the latter because of severe plasma leakage. Given that the outcome of infection is likely determined by the kinetics of viral replication and the antiviral host cell immune response (HIR) it is of importance to understand the interaction between these two parameters. In this study, we use mathematical modeling to characterize and understand the complex interplay between intracellular DV replication and the host cells' defense mechanisms. We first measured viral RNA, viral protein, and virus particle production in Huh7 cells, which exhibit a notoriously weak intrinsic antiviral response. Based on these measurements, we developed a detailed intracellular DV replication model. We then measured replication in IFN competent A549 cells and used this data to couple the replication model with a model describing IFN activation and production of IFN stimulated genes (ISGs), as well as their interplay with DV replication. By comparing the cell line specific DV replication, we found that host factors involved in replication complex formation and virus particle production are crucial for replication efficiency. Regarding possible modes of action of the HIR, our model fits suggest that the HIR mainly affects DV RNA translation initiation, cytosolic DV RNA degradation, and naïve cell infection. We further analyzed the potential of direct acting antiviral drugs targeting different processes of the DV lifecycle in silico and found that targeting RNA synthesis and virus assembly and release are the most promising anti-DV drug targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Mathematical model and calculation algorithm of micro and meso levels of separation process of gaseous mixtures in molecular sieves

In the given article, the main thermodynamic aspects of the issue of modeling diffusion transfer in molecular sieves have been formulated. Dissipation function is used as a basic notion. The differential equation, connecting volume flow with the change of the concentration of catchable component has been derived. As a result, the expression for changing the concentration of the catchable component and the coefficient of membrane detecting has been received. As well, the system approach to describing the process of gases separation in ultra porous membranes has been realized and micro and meso-levels of mathematical modeling have been distinguished. The non-ideality of the shared system is primarily taken into consideration at the micro-level and the departure from the diffusion law of Fick has been taken into account. The calculation method of selectivity considering fractal structure of membranes has been developed at the meso level. The calculation algorithm and its software implementation have been suggested.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

Searching for fat tails in CRISPR-Cas systems: Data analysis and mathematical modeling

Understanding CRISPR-Cas systems—the adaptive defence mechanism that about half of bacterial species and most of archaea use to neutralise viral attacks—is important for explaining the biodiversity observed in the microbial world as well as for editing animal and plant genomes effectively. The CRISPR-Cas system learns from previous viral infections and integrates small pieces from phage genomes called spacers into the microbial genome. The resulting library of spacers collected in CRISPR arrays is then compared with the DNA of potential invaders. One of the most intriguing and least well understood questions about CRISPR-Cas systems is the distribution of spacers across the microbial population. Here, using empirical data, we show that the global distribution of spacer numbers in CRISPR arrays across multiple biomes worldwide typically exhibits scale-invariant power law behaviour, and the standard deviation is greater than the sample mean. We develop a mathematical model of spacer loss and acquisition dynamics which fits observed data from almost four thousand metagenomes well. In analogy to the classical ‘rich-get-richer’ mechanism of power law emergence, the rate of spacer acquisition is proportional to the CRISPR array size, which allows a small proportion of CRISPRs within the population to possess a significant number of spacers. Our study provides an alternative explanation for the rarity of all-resistant super microbes in nature and why proliferation of phages can be highly successful despite the effectiveness of CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

A mathematical model of extremely low frequency ocean induced electromagnetic noise

Magnetic Anomaly Detection (MAD) system uses the principle that ferromagnetic objects disturb the magnetic lines of force of the earth. These lines of force are able to pass through both water and air in similar manners. A MAD system, usually mounted on an aerial vehicle, is thus often employed to confirm the detection and accomplish localization of large ferromagnetic objects submerged in a sea-water environment. However, the total magnetic signal encountered by a MAD system includes contributions from a myriad of low to Extremely Low Frequency (ELF) sources. The goal of the MAD system is to detect small anomaly signals in the midst of these low-frequency interfering signals. Both the Range of Detection (R{sub d}) and the Probability of Detection (P{sub d}) are limited by the ratio of anomaly signal strength to the interfering magnetic noise. In this paper, we report a generic mathematical model to estimate the signal-to-noise ratio or SNR. Since time-variant electro-magnetic signals are affected by conduction losses due to sea-water conductivity and the presence of air-water interface, we employ the general formulation of dipole induced electromagnetic field propagation in stratified media [1]. As a first step we employ a volumetric distribution of isolated elementary magnetic dipoles, each having its own dipole strength and orientation, to estimate the magnetic noise observed by a MAD system. Numerical results are presented for a few realizations out of an ensemble of possible realizations of elementary dipole source distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Understanding the Connection between Nanoparticle Uptake and Cancer Treatment Efficacy using Mathematical Modeling

Nanoparticles have shown great promise in improving cancer treatment efficacy while reducing toxicity and treatment side effects. Predicting the treatment outcome for nanoparticle systems by measuring nanoparticle biodistribution has been challenging due to the commonly unmatched, heterogeneous distribution of nanoparticles relative to free drug distribution. We here present a proof-of-concept study that uses mathematical modeling together with experimentation to address this challenge. Individual mice with 4T1 breast cancer were treated with either nanoparticle-delivered or free doxorubicin, with results demonstrating improved cancer kill efficacy of doxorubicin loaded nanoparticles in comparison to free doxorubicin. We then developed a mathematical theory to render model predictions from measured nanoparticle biodistribution, as determined using graphite furnace atomic absorption. Model analysis finds that treatment efficacy increased exponentially with increased nanoparticle accumulation within the tumor, emphasizing the significance of developing new ways to optimize the delivery efficiency of nanoparticles to the tumor microenvironment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Estimating the population level impact of a gonococcal vaccine candidate: Predictions from a simple mathematical model

Neisseria gonorrhoeae cross-protection was suggested in a New Zealand meningitis B vaccine. We modeled the potential impact of similar vaccines on gonorrhea prevalence in heterosexuals in the United States. Here, our mathematical model incorporated infection, behavior, and vaccination dynamics. Approximate Bayesian Computation calibrated our model to US prevalence. Primary analyses assumed New Zealand vaccine characteristics: 30% efficacy and 2-year duration of protection. We estimated impact under two vaccine coverages (20%, 50%). Reduction in gonorrhea prevalence ranged from 4.8 to 39.4%, depending on vaccine coverage. Vaccine impact was correlated with both size of the highly sexually active subpopulation and sexual mixing between high and low activity subpopulations. A meningitis vaccine providing low efficacy cross-protection against gonorrhea acquisition and short duration of protection could result in a large reduction in gonorrhea prevalence in the United States. Potential dual protective effects can be considered when making vaccine recommendations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Mathematical Models for Optimization of Grid-Integrated Energy Storage Systems

Energy storage has been proven to yield positive effects on planning, operation and control of electric grids. It has become a crucial task to properly model the energy storage systems (ESS) under the framework of grid optimization on transmission and distribution networks including microgrids. This paper presents a review on mathematical models and test cases of ESSs used for grid optimization studies, where the network constraints of power systems are included. The existing ESS models are mainly classified into two categories – linear and nonlinear models. The two main categories are further divided into several subcategories respectively; such as mixed-integer linear and convex nonlinear sub-categories. Based on the review and discussions, this paper aims at providing suggestions for choosing proper ESS models for specific grid optimization studies considering the chosen power network model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

SU-F-E-20: A Mathematical Model of Linac Jaw Calibration Integrated with Collimator Walkout

Purpose: Accurate jaw calibration is possible, but it does not necessarily achieve good junctions because of collimator rotation walkout. We developed a mathematical model seeking to pick an origin for calibration that minimizes the collimator walkout effect. Methods: We use radioopaque markers aligned with crosshair on the EPID to determine the collimator walkout at collimator angles 0°, 90° and 270°. We can accurately calibrate jaws to any arbitrary origin near the radiation field centre. While the absolute position of an origin moves with the collimator walkout, its relative location to the crosshair is an invariant. We studied two approaches to select an optimal origin. One approach seeks to bring all three origin locations (0°–90°–270°) as close as possible by minimizing the perimeter of the triangle formed by these points. The other approach focuses on the gap for 0°–90° junctions. Results: Our perimeter cost function has two variables and non-linear behaviour. Generally, it does not have zero-perimeter-length solution which leads to perfect jaw matches. The zero solution can only be achieved, if the collimator rotates about a single fixed axis. In the second approach, we can always get perfect 0°–0° and 0°–90° junctions, because we ignore the 0°–270° situation. For our TrueBeams, both techniques for selecting an origin improved junction dose inhomogeneities to less than ±6%. Conclusion: Our model considers the general jaw matching with collimator rotations and proposes two potential solutions. One solution optimizes the junction gaps by considering all three collimator angles while the other only considers 0°–90°. The first solution will not give perfect matching, but can be clinically acceptable with minimized collimator walkout effect, while the second can have perfect junctions at the expense of the 0°–270° junctions. Different clinics might choose between these two methods basing on their clinical practices.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Experimental measurements and mathematical modeling of cold plate for aviation thermal management

Herein, this study, which has been motivated by the recent applications of the cold plate device in aviation thermal management, reports on physics-based mathematical models derived from the conservation laws of mass, momentum, and energy, and empiricism-based models. One of the objectives of the present work is to report on an elaborate and successful experimental work carried out on an additively manufactured device for the purpose of rigorously validating the numerical predictions. The excellent agreement between the numerical predictions and measured performance provides much needed confidence in the implementation, in the software package, of the offset-strip fin passage correlations, as well as in the software implementation of user-defined wavy fin correlations for aerospace heat exchangers and cold plates operating with ram air at Reynolds numbers below 8000. The contributions of this work can also be found in the development of a new and accurate thermal-hydraulic analysis procedure, referred to in this paper as plate-fin analogy. Results from this procedure are compared with those from thermal resistance network. The comparative study in this paper of the bulk and discrete enthalpy flux method is also new, as is the relative assessment of four off-set strip fin thermal-hydraulic models.

42 ENGINEERING↗
DOE OSTI2020

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗
DOE OSTI2016

Comment on “Mathematical modeling of unicellular microalgae and cyanobacteria metabolism for biofuel production” by Baroukh et al. [Curr Opin Biotechnol. 2015, 33:198–205]

This letter concerns an article recently published in Current Opinion in Biotechnology by Baroukh et al. entitled "Mathematical modeling of unicellular microalgae and cyanobacteria metabolism for biofuel production" (2015 Mar 26; 33:198-205). The issue we wish to bring to light is the authors’ claim that dynamic metabolic models including Hybrid Cybernetic Model (HCM) (Kim et al., 2008; Song et al., 2009) and Lumped HCM (L-HCM) (Song and Ramkrishna, 2010; 2011) are based on the balanced growth hypothesis so that they are unable to simulate the accumulation of intracellular metabolites. This is a misrepresentation of these models due to the following reasons.

balanced growth↗
DOE OSTI2017

RELAP-7 Progress Report: A Mathematical Model for 1-D Compressible, Single-Phase Flow Through a Branching Junction

In the literature, the abundance of pipe network junction models, as well as inclusion of dissipative losses between connected pipes with loss coefficients, has been treated using the incompressible flow assumption of constant density. This approach is fundamentally, physically wrong for compressible flow with density change. This report introduces a mathematical modeling approach for general junctions in piping network systems for which the transient flows are compressible and single-phase. The junction could be as simple as a 1-pipe input and 1-pipe output with differing pipe cross-sectional areas for which a dissipative loss is necessary, or it could include an active component, between an inlet pipe and an outlet pipe, such as a pump or turbine. In this report, discussion will be limited to the former. A more general branching junction connecting an arbitrary number of pipes with transient, 1-D compressible single-phase flows is also presented. These models will be developed in a manner consistent with the use of a general equation of state like, for example, the recent Spline-Based Table Look-up method [1] for incorporating the IAPWS-95 formulation [2] to give accurate and efficient calculations for properties for water and steam with RELAP-7 [3].

42 ENGINEERING↗