SEARCH · Search NASA

Results for “mathematical foundation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. 1: Theoretical foundations

The theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions (3-D) is discussed. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 solar radii where the flow is supersonic) is approximately described by the nonlinear, single fluid, polytropic (magneto-) hydrodynamic equations. Efficient numerical techniques for solving this complex system of coupled, hyperbolic partial differential equations are outlined. The formulation is inviscid and nonmagnetic, but methods allow for the potential inclusion of both features with only modest modifications. One simple, highly idealized, hydrodynamic model stream is examined to illustrate the fundamental processes involved in the 3-D dynamics of stream evolution. Spatial variations in the rotational stream interaction mechanism were found to produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction.

Pizzo, V. J.↗

DOE OSTI2025

Basic Research Needs for Inverse Methods for Complex Systems under Uncertainty

Inverse problems, which aim to infer unknown properties of a system using experimental and observational data, are central to addressing many of the U.S. Department of Energy’s (DOE) most critical scientific and engineering challenges. Accurate, computationally efficient, and data-efficient solutions to inverse problems are essential for advancing DOE mission-critical science drivers, including analyzing data from large-scale experimental facilities, optimizing fusion reactor performance, accelerating materials discovery, enhancing geophysical imaging, improving wildfire predictions, and enabling autonomous systems and digital twins. However, these problems are becoming increasingly complex, often involving nonlinear, highdimensional, and interconnected systems and models that span multiple physics and scales, while relying on data with varying quantity, quality, and information content. Compounding these challenges is the uncertainty inherent in DOE-relevant systems, where errors in inputs, noise in data, incompleteness of data, and discrepancies between models and reality constrain the accuracy and precision of solutions. At the same time, the convergence of recent scientific computing trends—scientific machine learning, artificial intelligence, and computing advances such as exascale computing—is creating unprecedented opportunities for tackling these challenges. The cross-cutting nature of inverse problems, combined with their growing complexity and rapidly evolving data and algorithmic demands, strongly motivates the formulation of a prioritized research agenda to maximize their capabilities and impact. In response to this need, DOE’s Advanced Scientific Computing Research (ASCR) program in the Office of Science convened the Workshop on Basic Research Needs for Inverse Problems for Complex Systems Under Uncertainty in June 2025. This workshop brought together experts across disciplines to identify grand challenges and major opportunities in the field. Through collaborative discussions, the workshop defined transformative research directions aimed at addressing the mathematical, statistical, and computational challenges posed by inverse problems under uncertainty. As a result of these efforts, four priority research directions (PRDs) were identified to guide future research and development in this area. These PRDs, summarized below, represent a roadmap for advancing the foundational science and mathematics of inverse problems, enabling robust, scalable, and uncertainty-aware solutions that are critical for DOE applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches, Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models (GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models (CRMs) and large-eddy simulation (LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations: multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions; discrete vs “continuous” representation of hydrometeor types; turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation; theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology; and approaches for developing bulk microphysics parameterizations. Also presented are the spectral bin scheme and particle-based scheme (especially, super-droplet method) for representing explicit microphysics. Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations. Particle-resolved direct numerical simulation (DNS) models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds. Outstanding challenges and future research directions are explored as well.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

buhito

buhito is a Python library for graph analysis and machine learning. Graphs can represent networks with objects as nodes and their relationships as edges. buhito focuses on graphlet methods that study graphs through enumerating their component subgraphs to enable interpretable and fast models of complex systems. The package provides tools for different algorithmic designs for computing, analyzing, and applying graphlets to research problems such as machine learning, data compression, and anomaly detection in graph-structured data. A central feature is performing decomposition data analysis on graphs for machine learning models. Implemented in Python and built upon open-source scientific libraries such as NetworkX, NumPy, and SciPy, buhito provides high-performance methods for researchers exploring the mathematical and computational foundations of graphlet analysis applicable to systems of different sizes.

Pimonova, Yulia↗

NASA NTRS1976

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗

NASA NTRS1990

AMOEBA clustering revisited

A description of the clustering, classification, and image display program AMOEBA is presented. Using a difficult high resolution aircraft-acquired MSS image, the steps the program takes in forming clusters are traced. A number of new features are described here for the first time. Usage of the program is discussed. The theoretical foundation (the underlying mathematical model) is briefly presented. The program can handle images of any size and dimensionality.

Bryant, Jack↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

DOE OSTI2024

Classifying photonic topology using the spectral localizer and numerical K -theory

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Mathematics Integrated with Computer Science through Scratch: Curricular Modules for the Middle Grades

A project funded by the National Science Foundation, Computer Science Integrated with Mathematics in Middle Schools (CSIMMS) (DRL-1640039), brought middle-school mathematics teachers together with university computer science (CS) faculty and STEM education faculty to design, develop, and test curriculum modules in which CS is integrated into instruction for middle-school general-mathematics courses. The project developed integrated math/CS curriculum modules (for grades 6, 7, and 8), complete with student tasks and teacher materials to guide classroom implementation. Across the project, sixteen teachers from four middle schools representing different school contexts (both urban and suburban, serving students from a range of demographic and socioeconomic backgrounds) participated as members of the design team and trial testers of the math/CS modules. All modules underwent multiple years of testing and refinement, with extensive input from the teachers who co-designed and implemented them. Approximately 50-100 middle-school students participated in the classes of those who taught each year. The modules in this volume feature a range of models for integrating mathematics and computer science at the 6th and 7th grade levels. All modules foreground the teaching of grade-level mathematics content, with computer science functioning to motivate and/or reinforce these ideas.

Andrews larson, Christine↗

DOE OSTI2024

Foundation Models for the Electric Power Grid

Foundation models (FMs) currently dominate news headlines. They employ advanced deep learning architectures to extract structural information autonomously from vast datasets through self-supervision. The resulting rich representations of complex systems and dynamics can be applied to many downstream applications. Therefore, advances in FMs can find uses in electric power grids, challenged by the energy transition and climate change. This paper calls for the development of FMs for electric grids. We highlight their strengths and weaknesses amidst the challenges of a changing grid. It is argued that FMs learning from diverse grid data and topologies, which we call grid foundation models (GridFMs), could unlock transformative capabilities, pioneering a new approach to leveraging AI to redefine how we manage complexity and uncertainty in the electric grid. Finally, we discuss a practical implementation pathway and road map of a GridFM-v0, a first GridFM for power flow applications based on graph neural networks, and explore how various downstream use cases will benefit from this model and future GridFMs.

AI-based power flow simulation↗

NASA NTRS1984

An operator calculus for surface and volume modeling

The mathematical techniques which form the foundation for most of the surface and volume modeling techniques used in practice are briefly described. An outline of what may be termed an operator calculus for the approximation and interpolation of functions of more than one independent variable is presented. By considering the linear operators associated with bivariate and multivariate interpolation/approximation schemes, it is shown how they can be compounded by operator multiplication and Boolean addition to obtain a distributive lattice of approximation operators. It is then demonstrated via specific examples how this operator calculus leads to practical techniques for sculptured surface and volume modeling.

Gordon, W. J.↗

NASA NTRS1984

The Interaction of Global Biochemical Cycles

The global biosphere in an exceedingly complex system. To gain an understanding of its structure and dynamic features, it is necessary not only to increase the knowledge about the detailed processes but also to develop models of how global interactions take place. Attempts to analyze the detailed physical, chemical and biological processes in this context need to be guided by an advancement of understanding of the latter. It is necessary to develop a strategy of data gathering that severs both these purposes simultaneously. The following papers deal with critical aspects in the global cycles of carbon, nitrogen, phosphorus and sulfur in details as well as the cycle of water and the flow of energy in the Earth's environment. The objective is to set partly the foundation for the development of mathematical models that allow exploration of the coupled dynamics of the global cycles of carbon, nitrogen, phosphorus, sulfur, as well as energy and water flux.

Moore, B., III↗

NASA NTRS1999

A Methodology for Investigating Adaptive Postural Control

Our research on postural control and human-environment interactions provides an appropriate scientific foundation for understanding the skill of mass handling by astronauts in weightless conditions (e.g., extravehicular activity or EVA). We conducted an investigation of such skills in NASA's principal mass-handling simulator, the Precision Air-Bearing Floor, at the Johnson Space Center. We have studied skilled movement-body within a multidisciplinary context that draws on concepts and methods from biological and behavioral sciences (e.g., psychology, kinesiology and neurophysiology) as well as bioengineering. Our multidisciplinary research has led to the development of measures, for manual interactions between individuals and the substantial environment, that plausibly are observable by human sensory systems. We consider these methods to be the most important general contribution of our EVA investigation. We describe our perspective as control theoretic because it draws more on fundamental concepts about control systems in engineering than it does on working constructs from the subdisciplines of biomechanics and motor control in the bio-behavioral sciences. At the same time, we have attempted to identify the theoretical underpinnings of control-systems engineering that are most relevant to control by human beings. We believe that these underpinnings are implicit in the assumptions that cut across diverse methods in control-systems engineering, especially the various methods associated with "nonlinear control", "fuzzy control," and "adaptive control" in engineering. Our methods are based on these theoretical foundations rather than on the mathematical formalisms that are associated with particular methods in control-systems engineering. The most important aspects of the human-environment interaction in our investigation of mass handling are the functional consequences that body configuration and stability have for the pick up of information or the achievement of overt goals. It follows that an essential characteristic of postural behavior is the effective maintenance of the orientation and stability of the sensory and motor "platforms" (e.g., head or shoulders) over variations in the human, the environment and the task. This general skill suggests that individuals should be sensitive to the functional consequences of body configuration and stability. In other words, individuals should perceive the relation between configuration, stability, and performance so that they can adaptively control their interaction with the surroundings. Human-environment interactions constitute robust systems in that individuals can maintain the stability of such interactions over uncertainty about and variations in the dynamics of the interaction. Robust interactions allow individuals to adopt orientations and configurations that are not optimal with respect to purely energetic criteria. Individuals can tolerate variation in postural states, and such variation can serve an important function in adaptive systems. Postural variability generates stimulation which is "textured" by the dynamics of the human-environment system. The texture or structure in stimulation provides information about variation in dynamics, and such information can be sufficient to guide adaption in control strategies. Our method were designed to measure informative patterns of movement variability.

McDonald, P. V.↗

NASA NTRS1994

Integration of the Response Surface Methodology with the Compromise Decision Support Problem in Developing a General Robust Design Procedure

In this paper we introduce a comprehensive and rigorous robust design procedure to overcome some limitations of the current approaches. A comprehensive approach is general enough to model the two major types of robust design applications, namely, robust design associated with the minimization of the deviation of performance caused by the deviation of noise factors (uncontrollable parameters), and robust design due to the minimization of the deviation of performance caused by the deviation of control factors (design variables). We achieve mathematical rigor by using, as a foundation, principles from the design of experiments and optimization. Specifically, we integrate the Response Surface Method (RSM) with the compromise Decision Support Problem (DSP). Our approach is especially useful for design problems where there are no closed-form solutions and system performance is computationally expensive to evaluate. The design of a solar powered irrigation system is used as an example. Our focus in this paper is on illustrating our approach rather than on the results per se.

Chen, Wei↗

NASA NTRS2011

Kennedy Educate to Innovate (KETI) Microgravity Powerpoint Presentation

The purpose of this presentation is to define and explain microgravity and show how microgravity can help students learn about the phenomena of the world. The presentation is designed to provide teachers of science, technology, engineering, and mathematics at many levels with a foundation in microgravity science and applications.

Source record↗

NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

DOE OSTI2021

A probabilistic graphical model foundation for enabling predictive digital twins at scale

A unifying mathematical formulation is needed to move from one-off digital twins built through custom implementations to robust digital twin implementations at scale. This work proposes a probabilistic graphical model as a formal mathematical representation of a digital twin and its associated physical asset. We create an abstraction of the asset–twin system as a set of coupled dynamical systems, evolving over time through their respective state spaces and interacting via observed data and control inputs. The formal definition of this coupled system as a probabilistic graphical model enables us to draw upon well-established theory and methods from Bayesian statistics, dynamical systems and control theory. The declarative and general nature of the proposed digital twin model make it rigorous yet flexible, enabling its application at scale in a diverse range of application areas. Here, we demonstrate how the model is instantiated to enable a structural digital twin of an unmanned aerial vehicle (UAV). The digital twin is calibrated using experimental data from a physical UAV asset. Its use in dynamic decision-making is then illustrated in a synthetic example where the UAV undergoes an in-flight damage event and the digital twin is dynamically updated using sensor data. The graphical model foundation ensures that the digital twin calibration and updating process is principled, unified and able to scale to an entire fleet of digital twins.

42 ENGINEERING↗

DOE OSTI2022

Dynamics and growth rate implications of ribosomes and mRNAs interaction in E. coli

Understanding how cells grow and adapt under various nutrient conditions is pivotal in the study of biological stoichiometry. Recent studies provide empirical evidence that cells use multiple strategies to maintain an optimal protein production rate under different nutrient conditions. Mathematical models can provide a solid theoretical foundation that can explain experimental observations and generate testable hypotheses to further our understanding of the growth process. In this study, we generalize a modeling framework that centers on the translation process and study its asymptotic behaviors to validate algebraic manipulations involving the steady states. Using experimental results on the growth of E. coli under C-, N-, and P-limited environments, we simulate the expected quantitative measurements to show the feasibility of using the model to explain empirical evidence. Our results support the findings that cells employ multiple strategies to maintain a similar protein production rate across different nutrient limitations. Moreover, we find that the previous study underestimates the significance of certain biological rates, such as the binding rate of ribosomes to mRNA and the transition rate between different ribosomal stages. Furthermore, our simulation shows that the strategies used by cells under C- and P-limitations result in a faster overall growth dynamics than under N-limitation. In conclusion, the general modeling framework provides a valuable platform to study cell growth under different nutrient supply conditions, which also allows straightforward extensions to the coupling of transcription, translation, and energetics to deepen our understanding of the growth process.

59 BASIC BIOLOGICAL SCIENCES↗