Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Mathematical 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 19 records

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

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

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

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)

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

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry

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

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

System Identification of a DC-DC Buck Converter Based on Two-Channel Relay Method

This work presents a system identification approach based on a two-channel relay applied to a DC-DC buck converter operating in closed-loop with a PI controller to regulate its output voltage. The implemented two-relay method allows system identification while the system is running online. Moreover, the algorithm works in a closed-loop configuration while introducing a minimal perturbation in the converter's operation. The presented algorithm was used to identify the frequency response of the converter with no-prior knowledge about the system structure. Furthermore, the identified transfer function parameters has a good agreement with the mathematical model of the converter. Simulation results of the comparison between the switching model, mathematical model and the identified model are provided to validate the theoretical analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

System Identification for Integrated Aircraft Development and Flight Testing [l'Identification Des Systemes Pour le Developpement Integre des Aeronefs et les Essais en Vol]

Over the last decades flight vehicles such as aircraft and helicopters entering service and requiring increased operational effectiveness have with few exceptions experienced prolonged flight test development to achieve full certification. In many cases the original requirements had later to be reduced to enable release to service. The impact on the customer, and manufacturer has been considerable leading to increased costs and or reduced operational capabilities. These costly experiences are largely a result of the flight vehicle not behaving as modelled and designed. The evaluation of flight test data can be used as a tool for validating windtunnel results and mathematical models describing the flight dynamical behaviour. In this sense the uncertainty of important aerodynamic stability and control parameters can be reduced and the confidence of aircraft mathematical models improved. An additional important factor comes from the implementation of active control systems offering the promise of significantly increased flight vehicle performance and operational capability. This approach extends the traditional trade-offs between aerodynamics, structures and propulsion systems to include full- time, full-authority fly-by-wire/light systems. It is imperative that the aerodynamic stability and control parameters of such integrated flight and propulsion control systems have to turn out inflight as predicted, since inherent stability margins will be lower and the flight control system must correct these deficiencies to provide flight critical redundancy and safety. With the methodology of system identification from flight tests it is possible to sense the control inputs and the flight vehicle reactions Such as accelerations, rates and attitudes. The mathematical model, e.g. the model structure and parameters, has to be determined from the relationship of the measured control inputs and the system's responses. The aim of this symposium was to review the present state of the art of flight vehicle system and parameter identification techniques, and to provide a critical appraisal of current methods developed and applied to flight test data in a number of NATO nations. Particular emphasis was placed on practical aspects and lessons learned in order to generate information useful to the flight test community in industry and government agencies. The technical papers share invaluable experience and emphasize the advances of flight vehicle system identification over the last years to the point where confidence and robustness level is now reasonably high. The symposium covered overviews of identification methodologies, flight test techniques, recent aircraft and helicopter application programs, and a session of short papers covering up-to-the-minute flight test results. A final discussion included prepared comments from experts and concluded with key issues learned in the application of system identification and future research needs. The essential benefits to NATO nations can be condensed as follows: More accurate mathematical models for high bandwidth flight control systems, Improved assessment and evaluation of flying qualities, High fidelity mathematical models for flight vehicle development and mission training simulators, and generally, Reduced flight test time and costs.

Advisory Group for Aerospace Research and Developm

Implementation and evaluation of multi-dual mode counter-current chromatography in the CUP Modeler software

Counter-current chromatography (CCC) is a separation technique that utilizes immiscible solvent pairs as stationary and mobile phases, which imparts numerous benefits compared to solid-liquid chromatography including the ability to treat either the more-dense or less-dense solvent layer as the mobile phase. Multi-dual mode (MDM) is a CCC elution mode capable of improving the separation of closely eluting compounds by alternating upper- and lower-layer solvent flows in opposing directions within the same separation. While some effort has been made to model MDM, implementation of these models in experimental design has yet to be widely adopted. Accordingly, we further developed our previously published cell utilized partitioning (CUP) model to include MDM predictions with CCC and packaged the full suite of CUP modeling capabilities into a user-friendly, open-source tool called the CUP Modeler. The mathematical model for MDM CCC was derived and validated with experimental separation of ethyl guaiacol (EG) and ethyl phenol (EP), two compounds that co-elute in our previously demonstrated reductive catalytic fractionation (RCF) lignin monomer isolation method. The developed MDM model provided insights into the effect of multiple operating parameters - including stationary phase retention, flow rate, column efficiency, feed concentration ratio, selectivity factor, and solute distribution ratios - on the separation yields, productivity, and purities. Our model agreed with prevailing understanding of MDM but also revealed new insights including that the ideal distribution ratios for co-eluting solutes to be separated by MDM is between 1.1 and 1.5, with the lower value ideally close to 1.25. Overall, this work provides fundamental insights for MDM process design and enables broader adoption of general liquid-liquid chromatography with a new, open-source user-friendly interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Viral Dynamic Models During COVID‐19: Are We Ready for the Next Pandemic?

Mathematical models have been used for about 30 years to improve our understanding of virus-host interaction, in particular during chronic infections. During the COVID-19 pandemic, these models have been used to provide insights into the natural history of acute SARS-CoV-2 infection, optimize antiviral treatment strategies, understand factors associated with transmission, and optimize surveillance systems. The impact of modeling has been accelerated by the availability of unprecedented multidimensional immune data from animal and human systems, which enhanced partnerships between experimentalists and theorists and led to exciting new modeling and statistical developments. In this mini review, we examine the lessons learned from the COVID-19 pandemic and discuss the main insights provided by mathematical models of viral dynamics at the different stages of the outbreak. Although we focus on respiratory infection, we also consider the new areas for development in anticipation of future acute infections from new or reemerging pathogens.

59 BASIC BIOLOGICAL SCIENCES

Towards modeling phage therapy

Patients infected with life-threatening multi-drug resistant (MDR) bacteria have been treated with cocktails of bacteriophages. This is a complicated form of personalized medicine as the phages given to a patient have to be selected beforehand on the basis of their lytic capacity of the infecting bacteria. Because bacteria rapidly become resistant, the evolution of resistance to a diverse cocktail of phages is a complicated dynamical process, during which competing bacterial strains replace one another by accumulating several resistance mechanisms, each of which may involve a fitness cost. As a consequence, it is typically not known why a particular phage therapy succeeded or failed, and how one can optimize the composition of the cocktails to maximize the rate of success. To improve upon this, we extend an existing in vivo -calibrated mouse model into a novel mathematical model for the human situation, and include multiple phages infecting multiple bacterial strains, differing in their resistance to each of the phages. We adjust several parameter estimates of the bacterial model to the human situation, and use the model to describe a successful case of phage therapy involving several cocktails, each containing several phages. In the model, treatment success crucially depended on pretreatment resistance levels, and on the diversity and the timing of the cocktails. Once an appropriate cocktail is found, it is less important to further optimize the infection rates of the phages. Resistant bacterial strains expand rapidly when sensitive strains decline, and the higher the infectivity of the phages, the faster resistant strains expand. Because resistance evolves rapidly, it is best to provide a diverse set of phages right from the start of therapy, i.e., to hit hard and early, and create a high genetic barrier to bacterial resistance.

59 BASIC BIOLOGICAL SCIENCES

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE

Mechanical, Electrochemical & Thermal Models - Training (CRADA CRD-19-00798 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Renewable Energy Laboratory (NLR) will host Dr. Jaeyoung Lim from HMC for a period of one year to jointly develop mathematical models for battery cells and modules subject to mechanical crush. NLR will assist with the development of mathematical models that Dr. Lim will incorporate into his research effort on new concepts of mobility with electric vehicles.

33 ADVANCED PROPULSION SYSTEMS

Adopting Code Verification Methodology Based on Model Form

Code verification is an essential part of credibility analysis for computational models. It assesses whether the mathematical model is implemented correctly into the code and whether the numerical methods behave consistently, and is done before solution verification and validation. Robust guidance for code verification exists in the literature. However, there is no known, concise guide for selecting the approach based on the model form that also presents an overview of the common elements. This document was written to address this gap as an accessible reference for beginning a code-verification effort.

97 MATHEMATICS AND COMPUTING

Adapting Code Verification Methodology to Model Form

Code verification is an essential part of credibility analysis for computational models. It assesses whether the mathematical model is implemented correctly into the code and whether the numerical methods behave consistently, and is done before solution verification and validation. Robust guidance for code verification exists in the literature. However, there is no known, concise guide for selecting the approach based on the model form that also presents an overview of the common elements. This document was written to address this gap as an accessible reference for beginning a code-verification effort.

97 MATHEMATICS AND COMPUTING