Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 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 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 OSTI2021

Mathematical model of a personalized neoantigen cancer vaccine and the human immune system

Cancer vaccines are an important component of the cancer immunotherapy toolkit enhancing immune response to malignant cells by activating CD4 + and CD8 + T cells. Multiple successful clinical applications of cancer vaccines have shown good safety and efficacy. Despite the notable progress, significant challenges remain in obtaining consistent immune responses across heterogeneous patient populations, as well as various cancers. We present a mechanistic mathematical model describing key interactions of a personalized neoantigen cancer vaccine with an individual patient’s immune system. Specifically, the model considers the vaccine concentration of tumor-specific antigen peptides and adjuvant, the patient’s major histocompatibility complexes I and II copy numbers, tumor size, T cells, and antigen presenting cells. We parametrized the model using patient-specific data from a clinical study in which individualized cancer vaccines were used to treat six melanoma patients. Model simulations predicted both immune responses, represented by T cell counts, to the vaccine as well as clinical outcome (determined as change of tumor size). This model, although complex, can be used to describe, simulate, and predict the behavior of the human immune system to a personalized cancer vaccine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Mathematical Model for Li-S Cell with Shuttling-Induced Capacity Loss Approximation

Lithium sulfur (Li-S) batteries have the potential to outperform the current lithium ion batteries and transform the technology of the future. However, dissolution, diffusion, and shuttling of the dissolved polysulfides result in parasitic reactions and substantial capacity loss. To provide a better understanding of the shuttling process, a 1D porous electrode mathematical model has been developed in this paper. An approximation method is used to account for the shuttling-induced capacity loss by adding an extra source/sink term in the material balance equations for the species involved in the parasitic reactions. Shuttling constants used in the source terms can be determined by fitting the model predictions to the experimental measurements. Here, the results showed that by including the approximation method, the model was able to predict the active material loss and the continuous decrease of volume fractions of Li 2 S on the cathode surface. The model sheds light on the capacity loss mechanism occurring inside the cell as a result of the shuttling of polysulfides.

25 ENERGY STORAGE↗
DOE OSTI2022

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗
DOE OSTI2020

Mathematical Modeling of Hydroxide-Exchange-Membrane Water Electrolyzer

Water electrolyzers can transform intermittent renewable energy like solar energy and wind energy into the chemical energy of hydrogen with zero greenhouse-gas emissions. The hydroxide-exchange membrane electrolyzer (HEME) combines the capability to produce pressurized hydrogen with the advantage of being able to use low or non-platinum group metal (PGM) electrocatalysts in the alkaline environment.1 Hydroxide salts, for example, KOH, are added to the HEME water feed on both anode and cathode to improve its performance. However, the specific mechanism of performance improvement still needs to be further understood. In addition, at high current densities, bubble evolution can result in mass-transport limitations, a less well studied phenomena. Mathematical modeling is ideal to explore these issues as it is cost and time efficient and can deconvolute the physics, processes, and observed phenomena and study the applied-voltage breakdown. In this work, we extend our previously developed 1D two-phase continuum model2 to study the varies processes in the HEME and provide insights on performance optimizations. First, the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) kinetics at different hydroxide concentrations have been studied by rotating disk electrodes (RDE) and implemented in the model. Then, the model is calibrated and validated against experimental HEME polarization curves for different KOH concentrations as a liquid electrolyte. The model clearly shows a performance increase with increasing KOH concentrations, which is consistent with the experimental results. The reduced ohmic resistance and increased electrochemical active surface area (ECSA) are the two main reasons for performance increase. The large amount of hydroxide in the liquid electrolyte not only helps to distribute the reactant hydroxide throughout the catalyst layer (CL), which reduces ohmic loss, but also enables reaction at the interface between the liquid electrolyte and electrocatalyst, which increases the ECSA. Applied-voltage breakdown demonstrates that the electrolyzer performance is dominated by anode kinetics and ohmic loss. A comparison with the DI water feed shows a more uniform current distribution in the anode CL when KOH is added, which indicates a higher utilization of the CL. Second, we present modeling on the effects of bubble coverage. As gas evolves, part of the ECSA is minimized due to bubble coverage. To account for this effect, an empirical relationship between the fractional bubble coverage and the current density is implemented in the model.3 The model shows this bubble coverage effect is more pronounced at large current densities with DI water feed. Acknowledgements This work was funded under the HydroGEN Consortium by the Energy Efficiency and Renewable Energy, Hydrogen and Fuel Cell Technologies Office, of the U. S. Department of Energy under contract number DE-AC02-05CH11231. References R. Abbasi, B. P. Setzler, S. Lin, J. Wang, Y. Zhao, H. Xu, B. Pivovar, B. Tian, X. Chen, G. Wu and Y. Yan, 31, 1805876 (2019). L. N. Stanislaw, M. R. Gerhardt and A. Z. Weber, ECS Transactions, 92, 767 (2019). H. Vogt and R. J. Balzer, Electrochimica Acta, 50, 2073 (2005).

Liu, Jiangjin↗
DOE OSTI2022

Minimum reflux calculation for multicomponent distillation in multi‐feed, multi‐product columns: Mathematical model

Abstract Multi‐feed, multi‐product distillation columns are ubiquitous in multicomponent distillation systems. The minimum reflux ratio of a distillation column is directly related to its energy consumption and capital cost. Thus, it is a key parameter for distillation systems design, operation, and comparison. In this series, we present the first accurate shortcut based algorithmic method to determine the minimum reflux condition for any general multi‐feed, multi‐product (MFMP) distillation column separating any ideal multicomponent mixture. The classic McCabe‐Thiele or Underwood method is a special case of this general approach. Compared with existing techniques, this method does not involve any rigorous tray‐by‐tray calculation, nor does it require guessing of key components. In this first part of the series, we present the mathematical model for a general MFMP column, derive constraints for feasible separation and minimum reflux condition, discuss their geometric interpretations, and present an illustrative example to demonstrate the effectiveness of our approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

Evaluation of the rheological and electrical percolation of high‐density polyethylene/carbon black composites using mathematical models

Abstract In this work, conductive polymer composites (CPCs) of bio‐based polyethylene (BioPe) containing different concentrations of carbon black (CB) were developed. By using oscillatory rheology analysis, a Newtonian plateau was observed in BioPe, and all BioPe/CB composites had a behavior of a pseudo‐solid and that composites with volume fractions ranging from 0.24 to 0.56 presented higher viscosity, storage, and loss modulus. This suggests the formation of a percolated network and by using the power‐law models, it was observed that the electrical percolation threshold was higher than the rheological percolation threshold. The electrical conductivity was measured using the four‐point probe method and a sigmoid model was used to predict the CPCs' electrical conductivity percolation threshold. The results indicated that the four‐point probe method presented satisfactory results according to the calculated standard deviations and voltage–current characteristics for each round of measurements considering the same ranging as used in rheology analysis. The analytical model used showed a coefficient of determination ( R 2 ) higher than 95%, allowing the prediction of the electrical conductivity of the CPC and the percolation threshold as a function of the volumetric fraction of the CB.

da SiIva, Moacy P.↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗