Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical modelling”

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 253 records · Page 14

Mathematical modeling of control subsystems for CELSS: Application to diet

The dynamic control of a Closed Ecological Life Support System (CELSS) in a closed space habitat is of critical importance. The development of a practical method of control is also a necessary step for the selection and design of realistic subsystems and processors for a CELSS. Diet is one of the dynamic factors that strongly influences, and is influenced, by the operational states of all major CELSS subsystems. The problems of design and maintenance of a stable diet must be obtained from well characterized expert subsystems. The general description of a mathematical model that forms the basis of an expert control program for a CELSS is described. The formulation is expressed in terms of a complete set of time dependent canonical variables. System representation is dynamic and includes time dependent storage buffers. The details of the algorithm are described. The steady state results of the application of the method for representative diets made from wheat, potato, and soybean are presented.

Waleh, Ahmad↗

A systems analysis of the erythropoietic responses to weightlessness. Volume 1: Mathematical model simulations of the erythropoietic responses to weightlessness

Theoretical responses to weightlessness are summarized. The studies include development and validation of a model of erythropoiesis regulation, analysis of the behavior of erythropoiesis under a variety of conditions, simulations of bed rest and space flight, and an evaluation of ground-based animal studies which were conducted as analogs of zero-g. A review of all relevant space flight findings and a set of testable hypotheses which attempt to explain how red cell mass decreases in space flight are presented. An additional document describes details of the mathematical model used in these studies.

Leonard, J. I.↗

(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↗

Development of mathematical models of environmental physiology

Selected articles concerned with mathematical or simulation models of human thermoregulation are presented. The articles presented include: (1) development and use of simulation models in medicine, (2) model of cardio-vascular adjustments during exercise, (3) effective temperature scale based on simple model of human physiological regulatory response, (4) behavioral approach to thermoregulatory set point during exercise, and (5) importance of skin temperature in sweat regulation.

Stolwijk, J. A. J.↗

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↗

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↗

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↗

Fully Automatic Guidance and Control for Rotorcraft Nap-of-the-earth Flight Following Planned Profiles. Volume 2: Mathematical Model

Developing a single-pilot, all-weather nap-of-the-earth (NOE) capability requires fully automatic NOE (ANOE) navigation and flight control. Innovative guidance and control concepts are investigated in a four-fold research effort that: (1) organizes the on-board computer-based storage and real-time updating of NOE terrain profiles and obstacles in course-oriented coordinates indexed to the mission flight plan; (2) defines a class of automatic anticipative pursuit guidance algorithms and necessary data preview requirements to follow the vertical, lateral, and longitudinal guidance commands dictated by the updated flight profiles; (3) automates a decision-making process for unexpected obstacle avoidance; and (4) provides several rapid response maneuvers. Acquired knowledge from the sensed environment is correlated with the forehand knowledge of the recorded environment (terrain, cultural features, threats, and targets), which is then used to determine an appropriate evasive maneuver if a nonconformity of the sensed and recorded environments is observed. This four-fold research effort was evaluated in both fixed-base and moving-base real-time piloted simulations; thereby, providing a practical demonstration for evaluating pilot acceptance of the automated concepts, supervisory override, manual operation, and re-engagement of the automatic system. Volume one describes the major components of the guidance and control laws as well as the results of the piloted simulations. Volume two describes the complete mathematical model of the fully automatic guidance system for rotorcraft NOE flight following planned flight profiles.

Clement, Warren F.↗

Building Mathematical Models Of Solid Objects

Solid Modeling Program (SMP) version 2.0 provides capability to model complex solid objects mathematically through aggregation of geometric primitives (parts). System provides designer with basic set of primitive parts and capability to define new primitives. Six primitives included in present version: boxes, cones, spheres, paraboloids, tori, and trusses. Written in VAX/VMS FORTRAN 77.

Randall, Donald P.↗

Time and frequency-domain identification and verification of BO-105 dynamic models

Mathematical models for the dynamics of the DLR BO 105 helicopter are extracted from flight test data using two different approaches: frequency-domain and time-domain identification. Both approaches are reviewed. Results from an extensive data consistency analysis are given. Identifications for 6 degrees of freedom (DOF) rigid body models are presented and compared in detail. The extracted models compare favorably and their prediction capability is demonstrated in verification results. Approaches to extend the 6 DOF models are addressed and first results are presented. System identification is broadly defined as the deduction of system characteristics from measured data. It provides the only possibility to extract both non-parametric (e.g., frequency responses) and parametric (e.g., state space matrices) aircraft models from flight test data and therefore gives a reliable characterization of the dynamics of the actually existing aircraft. Main applications of system identification are seen in areas where higher accuracies of the mathematical models are required: Simulation validation, control system design (in particular model-following control system design for in-flight simulation), and handling qualities.

Kaletka, Juergen↗

Time and frequency-domain identification and verification of BO-105 dynamic models

Mathematical models for the dynamics of the DLR BO 105 helicopter are extracted from flight test data using two different approaches: frequency-domain and time-domain identification. Both approaches are reviewed. Results from an extensive data consistency analysis are given. Identifications for 6 degrees of freedom (DOF) rigid body models are presented and compared in detail. The extracted models compare favorably and their prediction capability is demonstrated in verification results. Approaches to extend the 6 DOF models are addressed and first results are presented. System identification is broadly defined as the deduction of system characteristics from measured data. It provides the only possibility to extract both non-parametric (e.g., frequency responses) and parametric (e.g., state space matrices) aircraft models from flight test data and therefore gives a reliable characterization of the dynamics of the actually existing aircraft. Main applications of system identification are seen in areas where higher accuracies of the mathematical models are required: Simulation validation, control system design (in particular model-following control system design for in-flight simulation), and handling qualities.

Kaletka, Juergen↗

A mathematical model for the doubly fed wound rotor generator

A mathematical analysis of a doubly-fed wound rotor machine used as a constant frequency generator is presented. The purpose of this analysis is to derive a consistent set of circuit equations which produce constant stator frequency and constant stator voltage. Starting with instantaneous circuit equations, the necessary rotor voltages and currents are derived. The model, thus obtained, is assumed to be valid, since the resulting relationships between mechanical power and active volt-amperes agrees with the results of others. In addition, the model allows for a new interpretation of the power flow in the doubly-fed generator.

Brady, F. J.↗

Mathematical Model Development and Experimental Verification of Electromagnetic Force Actuators for use in Space Applications

This technical memorandum shows the development of the mathematics needed to model several types of electromagnetic force actuators that may be of interest for use in space applications. The actuators described are capable of providing attractive, repulsive, and longitudinal forces between the actuator and a variety of metal surfaces. Potential applications for these propellant free actuators include spacecraft docking, in-space fabrication, robotic spacecraft inspection and servicing, deflection of metallic debris, and interaction with iron core asteroids (e.g., landing, repelling, and relocating).

Robert C. Youngquist↗

Thermal mathematical modeling and system simulation of Space Shuttle less subsystem

Applications, validation tests, and upgrades of the two- and three-dimensional system level thermal mathematical system simulation models (TMSSM) used for thermal protection system (TPS) analyses are described. The TMSSM were developed as an aid to predicting the performance requirements and configurations of the Shuttle wing leading edge (WLE) and nose cone (NC) TPS tiles. The WLE and its structure were subjected to acoustic, thermal/vacuum, and air loads tests to simulate launch, on-orbit, and re-entry behavior. STS-1, -2 and -5 flight data led to recalibration of on-board instruments and raised estimates of the thermal shock at the NC and WLE. Baseline heating data are now available for the design of future TPS.

Chao, D. C.↗

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↗

A Mathematical Model of Oxygen Transport in Skeletal Muscle During Hindlimb Unloading

During hindlimb unloading (HU) dramatic fluid shifts occur within minutes of the suspension, leading to a less precise matching of blood flow to O2 demands of skeletal muscle. Vascular resistance directs blood away from certain muscles, such as the soleus (SOL). The muscle volume gradually reduces in these muscles so that eventually the relative blood flow returns to normal. It is generally believed that muscle volume change is not due to O2 depletion, but a consequence of disuse. However, the volume of the unloaded rat muscle declines over the course of weeks, whereas the redistribution of blood flow occurs immediately. Using a Krogh Cylinder Model, the distribution of O2 was predicted in two skeletal muscles: SOL and gastrocnemius (GAS). Effects of the muscle blood flow, volume, capillary density, and O2 uptake, are included to calculate the pO2 at rest and after 10 min and 15 days of unloading. The model predicts that 32 percent of the SOL muscle tissue has a pO2 1.25 mm Hg within 10 min, whereas the GAS maintains normal O2 levels, and that equilibrium is reached only as the SOL muscle cells degenerate. The results provide evidence that there is an inadequate O2 supply to the mitochondria in the SOL muscle after 10 min HU.

Mathematical models↗