Two Problems in Applied Stochastic Processes: Scaling Limits of the Multiple Range in 1D Random Walk and Uncertainty Quantification in Satellite Imagery
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “Applied Mathematics”
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.
Analytic bounds on late-time axion-scalar cosmologies
The cosmological dynamics of multiple scalar/pseudoscalar fields are difficult to solve, especially when the field-space metric is curved. This presents a challenge in determining whether a given model can support cosmic acceleration, without solving for the on-shell solution. In this work, we present bounds on late-time FLRW-cosmologies in classes of theories that involve arbitrary numbers of scalar and pseudoscalar fields coupled both kinetically (leading to a curved field space metric) and through scalar potentials. Such bounds are proven analytically, independently of initial conditions, with no approximation in the field equations and without referring to explicit solutions. Besides their broad applications to cosmological model building, our bounds can be applied to studying asymptotic cosmologies of certain classes of string compactifications.
79 ASTRONOMY AND ASTROPHYSICS↗
Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools
This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity
This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.
43 PARTICLE ACCELERATORS↗
Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems
Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.
97 MATHEMATICS AND COMPUTING↗
A retrospective review of von Neumann’s analysis of hidden variables in quantum mechanics
This article reviews the history of J. von Neumann’s analysis of hidden variables in quantum mechanics and the subsequent analysis by others. In his book The Mathematical Foundations of Quantum Mechanics , published in 1932, von Neumann performed an analysis of the consequences of introducing hidden parameters (hidden variables) into quantum mechanics. He arrived at two principal conclusions: first, hidden variables cannot be incorporated into the existing theory of quantum mechanics without major modifications, and second, if they did exist, the theory would have already failed in situations where it has been successfully applied. This analysis has been taken as an “incorrect proof” against the existence of hidden variables, possibly due to a mistranslation of the German word prufen . von Neumann’s so-called proof isn’t even wrong as such a proof does not exist, but it is an examination of the limitations imposed by internal consistency of the Hilbert space formulation of the theory. One of the earliest attempts to eliminate uncertainty, by D. Bohm, requires a major modification of quantum mechanics (observables are not represented by Hermitian operators), which supports von Neumann’s first principal conclusion. However, testing the Bohm theory requires constructing a physically impossible initial state. As such, the theory has no experimental consequences, so W. Pauli referred to it as an “uncashable check”. As there are no observable consequences, the Bohm theory is possibly a counterexample to von Neumann’s second conclusion that hidden variables in particular would have already led to a failure of the theory.
density matrix↗
An Advanced Machine Learning and Artificial Intelligence System for Demonstrating Radiation Regulatory Compliance in DOE Accelerator Facilities
In this Phase II proposal, Applied Research LLC (ARLLC), Thomas Jefferson National Accelerator Facility (Jefferson Lab), and Old Dominion University (ODU) propose the combination of domain knowledge (beam characteristics, fixed structural shielding, earthen burden (the soil and foliage added to the dome of the experimental halls as additional shielding), etc.), machine learning (ML) and/or artificial intelligence (AI) to correlate a variety of multi-modal onsite signals and the radiation fields seen in accessible areas of the accelerator site and the site boundary. The ML/AI will consider the complex influence of environmental parameters affecting the radon contribution of the measurements, focusing on actual data obtained from Jefferson Lab. In Phase I, the coded beam and location data were fed into a deep learning model to predict doses at several designated locations in Jefferson Lab’s facility. Moreover, a dense radiation map was generated using only a sparse collection of the samples in a facility. In Phase II, we will develop a software prototype containing a radiation prediction algorithm, dense radiation map algorithms, and background noise prediction algorithms, with actual data used to evaluate the prototype. This work will provide a framework for evaluation of radiation measurement results around the site based on learned responses. In addition, the proposed approach allows more granular mapping of radiation levels. Better understanding and communication of these levels is related to the overall approach in keeping doses to personnel ALARA.
43 PARTICLE ACCELERATORS↗
Interesting Solar Energy System Simulations in 1990 and 2025
This presentation uses a simulation by the author in 1991 and newer developements in 2025 to illustrate strategies to address problems that arise when steady state assumption is applied in time series simulation: 1) high resolution time series data; 2) distirbution functiion; 3) machine learning. The presentation does not report new findings (previously published material is cited).
14 SOLAR ENERGY↗
Aggregated DER_A Model Parameterization via Online Moving Horizon Estimation
Here, this paper introduces a methodology for parameterizing the DER_A model using a novel smooth mathematical representation, simplifying the process and preserving accuracy in modeling inverter-based generator (IBG). The methodology employs an online parameterization process that can operate in real-time. The model parameterization process is structured into five sequential steps, each targeting a specific aspect of the DER_A model through moving horizon estimation. This approach adapts to systems with varying voltage and frequency support requirements by selectively applying each step. Simulation results on systems with both known and unknown parameters validate the methodology’s effectiveness. The online moving horizon estimation technique accurately captures the dynamics of the overall system and ensures that the parameterized DER_A model closely mirrors the real system’s voltage, current, and power dynamics. The findings highlight the potential of this methodology to substantially improve and simplify the dynamic modeling of power systems, paving the way for more reliable and robust IBG and grid integration.
42 ENGINEERING↗
Half-closed discontinuous Galerkin discretisations
Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.
97 MATHEMATICS AND COMPUTING↗
An overview of HR-EBSD techniques for mapping local stress and dislocations in crystalline materials at sub-micron resolution
High resolution electron backscatter diffraction (HR-EBSD) is a technique used to map elastic strain, crystallographic orientation and dislocation density in a scanning electron microscope. Here, this review covers the background and mathematics of this technique, contextualizing it within the broader landscape of EBSD techniques and other materials characterization methods. Several case studies are presented showing the application of HR-EBSD to the study of plasticity in metals, failure analysis in microelectronics and defect quantification in thin films. This is intended to be a comprehensive resource for researchers developing this technique as well as an introduction to those wishing to apply it.
Ruggles, Timothy J. [Sandia National Laboratories ↗
Developing Data-Driven Synthetic Infrastructure Models for Resilience Analysis
Research on infrastructure resilience has produced promising methods to simulate and optimize complex networks to improve performance. However, restrictions on sharing infrastructure models and the steep cost of developing and maintaining infrastructure models presents a roadblock to adoption. To overcome this limitation, this research focuses on methods to create data-driven infrastructure models that will help improve infrastructure resilience and security. The analysis couples incomplete utility data, geospatial data, machine learning, and synthetic network generation methods to rapidly develop and update infrastructure models. The methods are validated using realistic utility models and site-specific data, with a focus on Puerto Rico due to its unique infrastructure challenges and available data. This research highlights promising opportunities for the use of synthetic network generation and machine learning to create infrastructure models when very little data is available. Results demonstrate that hybrid methods, which combine sparse utility data with synthetic models, can enhance model accuracy, and machine learning can predict model attributes using training data from other models. However, the complexity of infrastructure systems means that even minor changes in network connectivity can significantly impact simulation results. Resilience analysis using synthetic infrastructure models shows that while some system behaviors are preserved, the magnitude of disruptions may not be accurately represented, indicating the need for more research and validation before using synthetic models for critical infrastructure investment decisions. The framework outlined in this report represents a significant advance to infrastructure model development and could be applied to additional domains and sites. Future research will continue to streamline and validate methods to help reduce roadblocks to resilience analysis.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Erratum to The Scaling and Units of the Elastic Response Term for Rayleigh Waves that is Output by Computer Programs in Seismology (CPS)
This erratum corrects two notational errors and disambiguates a statement in the original article. These corrections apply to equations that follow Equation 9 in the text of the original.
58 GEOSCIENCES↗
Bridging the time scale in exascale computing of chemical systems (Final Technical Report)
This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.
08 HYDROGEN↗
Toward a sustainable circular economy of multilayer plastic films: Life cycle and techno-economic assessment with a focus on end-of-life treatment and multiple recovery cycles
This study presents a life cycle assessment (LCA) and techno-economic analysis (TEA) of end-of-life technologies for treating polyethylene–polyamide barrier film waste, focusing on quality degradation across recovery cycles. Novel treatment methods are experimentally validated, while others are drawn from literature and industry consultations. A displacement approach, assuming no quality loss, is first applied. Results show that solvent-based recycling via the solvent-targeted recovery and precipitation (STRAP) process outperforms alternatives across environmental indicators, reducing global warming potential (GWP) by 40% compared to landfilling. Incineration performs worst in most categories, particularly eutrophication (80% higher than landfilling), due to nitrogen emissions. Experimentally validated downcycling (pelletizing) proves more economically viable. The assumption of infinite recoverability is overly optimistic. To address this, we propose a mathematical framework accounting for a finite number of recovery cycles. This refined model shows reduced GWP and cost savings for solvent recovery, making its benefits less pronounced than initially estimated. Sensitivity and uncertainty analyses reveal strong dependence on recovered material quality and solvent recovery efficiency, underscoring the need for optimized process design. Finally, hotspot analysis identifies greenhouse gas emissions from the polyamide supply chain as the dominant GWP contributor. In conclusion, the results underscore potential trade-offs across pathways and show that solvent-based recovery’s sustainability depends heavily on process conditions.
36 MATERIALS SCIENCE↗
Towards verifiable cancer digital twins: tissue level modeling protocol for precision medicine
Cancer exhibits substantial heterogeneity, manifesting as distinct morphological and molecular variations across tumors, which frequently undermines the efficacy of conventional oncological treatments. Developments in multiomics and sequencing technologies have paved the way for unraveling this heterogeneity. Nevertheless, the complexity of the data gathered from these methods cannot be fully interpreted through multimodal data analysis alone. Mathematical modeling plays a crucial role in delineating the underlying mechanisms to explain sources of heterogeneity using patient-specific data. Intra-tumoral diversity necessitates the development of precision oncology therapies utilizing multiphysics, multiscale mathematical models for cancer. This review discusses recent advancements in computational methodologies for precision oncology, highlighting the potential of cancer digital twins to enhance patient-specific decision-making in clinical settings. We review computational efforts in building patient-informed cellular and tissue-level models for cancer and propose a computational framework that utilizes agent-based modeling as an effective conduit to integrate cancer systems models that encode signaling at the cellular scale with digital twin models that predict tissue-level response in a tumor microenvironment customized to patient information. Furthermore, we discuss machine learning approaches to building surrogates for these complex mathematical models. These surrogates can potentially be used to conduct sensitivity analysis, verification, validation, and uncertainty quantification, which is especially important for tumor studies due to their dynamic nature.
60 APPLIED LIFE SCIENCES↗
Using PyBioNetFit to leverage qualitative and quantitative data in biological model parameterization and uncertainty quantification
Data generated in studies of cellular regulatory systems are often qualitative. For example, measurements of signaling readouts in the presence and absence of mutations may reveal a rank ordering of responses across conditions but not the precise extents of mutation-induced differences. Qualitative data are often ignored by mathematical modelers or are considered in an ad hoc manner, as in the study of Kocieniewski and Lipniacki (2013) [Phys Biol 10: 035006], which was focused on the roles of MEK isoforms in ERK activation. In this earlier study, model parameter values were tuned manually to obtain consistency with a combination of qualitative and quantitative data. This approach is not reproducible, nor does it provide insights into parametric or prediction uncertainties. Here, starting from the same data and the same ordinary differential equation (ODE) model structure, we generate formalized statements of qualitative observations, making these observations more reusable, and we improve the model parameterization procedure by applying a systematic and automated approach enabled by the software package PyBioNetFit. We also demonstrate uncertainty quantification (UQ), which was absent in the original study. Our results show that PyBioNetFit enables qualitative data to be leveraged, together with quantitative data, in parameterization of systems biology models and facilitates UQ. These capabilities are important for reliable estimation of model parameters and model analyses in studies of cellular regulatory systems and reproducibility.
59 BASIC BIOLOGICAL SCIENCES↗
AIR Framework for Physics-Inspired Artificial Intelligence in High Energy Physics (Final Report)
The FAIR4HEP project was a collaboration between Argonne National Laboratory, University of Illinois Urbana-Champaign, Massachusetts Institute of Technology, University of Minnesota, and the University of California San Diego funded by the US Department of Energy, Office of Science, Office of Advanced Scientific Research (ASCR) from 2020 to 2023. The primary focus of the FAIR4HEP project was to advance our understanding of the relationship between data and artificial intelligence (AI) models by exploring relationships among them through the development of findable, accessible, interoperable, and reusable (FAIR) frameworks. Using high-energy physics (HEP) as the science driver, this project developed a FAIR framework to advance our understanding of AI, provide new insights to apply AI techniques and provide an environment where novel approaches to AI can be explored. This final report summarizes the accomplishments of the MIT group.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗