Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICS - MATRICES”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Cooperative Education

Los Alamos National Laboratory (LANL) is a multidisciplinary national laboratory that conducts research and development in national security, engineering, materials science, computational modeling, and advanced manufacturing. The laboratory develops innovative technologies to address complex scientific and engineering challenges. This project focuses on the development and evaluation of high-performance absorbing structures through computational design, simulation, and engineering analysis. Absorbing structures are used in applications where damage mitigation, structural protection, and material efficiency are critical performance requirements. The increasing demand for lightweight, high-strength, and highly efficient structural systems has created a need for improved design methodologies capable of maximizing absorption while minimizing weight and material usage. The project utilizes advanced engineering software, including 3D CAD software and FEA, to generate and optimize structural concepts. Computational simulations are performed to evaluate structural behavior under loading conditions, while mathematical analyses are conducted using Python-based tools as well as established analytical equations from material and structural mechanics. The project benefits LANL by supporting the development of advanced design methodologies and improving the understanding of material and structural performance. During the internship term, a significant portion of the design development, simulation, and data analysis activities will be completed. Success of the project depends on collaboration among engineering mentors and technical staff members. Work will be conducted at Los Alamos National Laboratory using laboratory computing resources and engineering software.

42 ENGINEERING↗
DOE OSTI2026

From Ensemble Climate to Ensemble Impacts

Many climate-risk tools rely on ensemble mean projections or endpoint climate snapshots to characterize future hazards. Although convenient for communication, these representations remove the statistical, temporal, and physical information that real infrastructure systems respond to. Infrastructure degradation and failure arise from extremes, sequences, cumulative stress, compound hazards, and nonlinear fragility relationships, none of which survive ensemble averaging or temporal compression. Power-system failure statistics and cascading failure models further show that infrastructure risk is dominated by tail events and path-dependent dynamics rather than by mean conditions. This paper demonstrates why ensemble mean or endpoint-only climate representations are mathematically and physically inconsistent with engineering-grade risk analysis. We outline a model-resolved, time-series-based workflow that preserves extremes, variability, and sequencing by propagating each climate-model realization independently through hazard formation, exposure, fragility, and cascading failure mechanisms. Taking the ensemble of impacts—rather than the ensemble of climate—provides a defensible, physically coherent foundation for infrastructure resilience planning, regulatory compliance, and long-term investment decisions.

54 - ENVIRONMENTAL SCIENCES/GLOBAL CLIMATE CHANGE ↗
DOE OSTI2026

On Transforming the Specific Force Vector for an ECEF Inertial Navigation System

An inertial measurement unit accumulates the instantaneous specific force applied to the sensor over a small window of time. Importantly, this accumulation occurs in the sensor’s frame, yet the specific force must be mathematically integrated in a different navigation frame, and these two frames may be rotating with respect to each other. Consequently, a correction is needed to not only account for the frame difference but also for the time variation in the frame difference. Many references provide derivations of this correction [1], but the derivation is often abbreviated, and important assumptions are sometimes left unstated. This memo addresses this gap and provides a clear derivation of the specific force transform with all assumptions specifically stated. This derivation is specific to the Earth-centered, Earth-fixed navigation frame, but, hopefully, the derivation is sufficiently clear that readers can apply it to other navigation frames. Generally, notation follows the same format as described in both [1, 2].

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

High Purity Germanium Detector Normalization

The presentation summarizes a project completed by our group a couple years ago comparing two detectors and attempting to normalize the results using a mathematical equation and how those results impacted our current operations.

Drumm, Natalie Daphne [Sandia National Laboratorie↗
DOE OSTI2026

Reconstruction of 6D Phase Space

A ‘Phase Space’ is a mathematical representation of all possible physical states of matter of a physical volume or system. Particle beams at Fermilab exist in a six-dimensional phase space: three position dimensions, and three momentum dimensions. To ‘reconstruct’ that space means to take measurement data from virtual detectors along those beams and map out where particles are and how they’re moving, using computational methods. In this experiment, we simulated a particle beam (Fig 1) and used the position distribution of each particle at each detector to map out the magnitude of the beam’s initial momentum.

Shirlee, Jermaine [DuPage Coll.]↗
DOE OSTI2026

SIMBER: the Simula Berkeley Education and Research Collaboration (CRADA Final Report)

The SIMBER project is centered around advancing the state-of-the-art in the science of modelling of the human heart and leveraging the collaborations between leading research groups at the University of California Berkeley (UC Berkeley), the Lawrence Berkeley National Laboratory (Berkeley Lab), and Simula Research Laboratory (Simula). In particular, the following areas of expertise are shared and expanded through this project: “heart-on-chip” experimental systems from UC Berkeley, mathematical modelling of the heart from Simula, and supercomputing software from Simula and Berkeley Lab. This intersection of expertise is already enabling the development of novel tools and knowledge that produce more accurate models of the human heart in both health and disease, which in turn are leading to drug screening technologies that make cardiac drug development faster, cheaper and more humane.

Li, Xiaoye Sherry [Lawrence Berkeley National Labo↗
DOE OSTI2025

What exactly does Bekenstein bound?

The Bekenstein bound posits a maximum entropy for matter with finite energy confined to a spatial region. It is often interpreted as a fundamental limit on the information that can be stored by physical objects. In this work, we test this interpretation by asking whether the Bekenstein bound imposes constraints on a channel's communication capacity, a context in which information can be given a mathematically rigorous and operationally meaningful definition. We study specifically the Unruh channel that describes a stationary Alice exciting different species of free scalar fields to send information to an accelerating Bob, who is confined to a Rindler wedge and exposed to the noise of Unruh radiation. We show that the classical and quantum capacities of the Unruh channel obey the Bekenstein bound that pertains to the decoder Bob. In contrast, even at high temperatures, the Unruh channel can transmit a significant number of zero-bits , which are quantum communication resources that can be used for quantum identification and many other primitive protocols. Therefore, unlike classical bits and qubits, zero-bits and their associated information processing capability are generally not constrained by the Bekenstein bound. However, we further show that when both the encoder and the decoder are restricted, the Bekenstein bound does constrain the channel capacities, including the zero-bit capacity.

Hayden, Patrick [Stanford University, CA (United S↗
DOE OSTI2024

Is the exquisite specificity of lymphocytes generated by thymic selection or due to evolution?

We have previously argued that the antigen receptors of T and B lymphocytes evolved to be sufficiently specific to avoid massive deletion of clonotypes by negative selection. Their optimal ‘specificity’ level, i.e., probability of binding any particular epitope, was shown to be inversely related to the number of self-antigens that the cells have to be tolerant to. Experiments have demonstrated that T lymphocytes also become more specific during negative selection in the thymus, because cells expressing the most crossreactive receptors have the highest likelihood of binding a self-antigen, and hence to be tolerized (i.e., deleted, anergized, or diverted into a regulatory T cell phenotype). Thus, there are two —not mutually exclusive— explanations for the exquisite specificity of T cells, one involving evolution and the other thymic selection. To better understand the impact of both, we extend a previously developed mathematical model by allowing for T cells with very different binding probabilities in the pre-selection repertoire. We confirm that negative selection tends to tolerize the most crossreactive clonotypes. As a result, the average level of specificity in the functional post-selection repertoire depends on the number of self-antigens, even if there is no evolutionary optimization of binding probabilities. However, the evolutionary optimal range of binding probabilities in the pre-selection repertoire also depends on the number of self-antigens. Species with more self antigens need more specific pre-selection repertoires to avoid excessive loss of T cells during thymic selection, and hence mount protective immune responses. We conclude that both evolution and negative selection are responsible for the high level of specificity of lymphocytes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Towards interpretable Cryo-EM: disentangling latent spaces of molecular conformations

Molecules are essential building blocks of life and their different conformations (i.e., shapes) crucially determine the functional role that they play in living organisms. Cryogenic Electron Microscopy (cryo-EM) allows for acquisition of large image datasets of individual molecules. Recent advances in computational cryo-EM have made it possible to learn latent variable models of conformation landscapes. However, interpreting these latent spaces remains a challenge as their individual dimensions are often arbitrary. The key message of our work is that this interpretation challenge can be viewed as an Independent Component Analysis (ICA) problem where we seek models that have the property of identifiability. That means, they have an essentially unique solution, representing a conformational latent space that separates the different degrees of freedom a molecule is equipped with in nature. Thus, we aim to advance the computational field of cryo-EM beyond visualizations as we connect it with the theoretical framework of (nonlinear) ICA and discuss the need for identifiable models, improved metrics, and benchmarks. Moving forward, we propose future directions for enhancing the disentanglement of latent spaces in cryo-EM, refining evaluation metrics and exploring techniques that leverage physics-based decoders of biomolecular systems. Moreover, we discuss how future technological developments in time-resolved single particle imaging may enable the application of nonlinear ICA models that can discover the true conformation changes of molecules in nature. The pursuit of interpretable conformational latent spaces will empower researchers to unravel complex biological processes and facilitate targeted interventions. This has significant implications for drug discovery and structural biology more broadly. More generally, latent variable models are deployed widely across many scientific disciplines. Thus, the argument we present in this work has much broader applications in AI for science if we want to move from impressive nonlinear neural network models to mathematically grounded methods that can help us learn something new about nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Variance Preserving Spectral Subsampling

Generating statistically faithful short-duration gamma-ray spectra from a single long measurement is essential in nuclear safeguards, supporting tasks such as algorithm development and machine-learning applications, especially when list-mode data are unavailable. Existing subsampling methods often distort the statistical characteristics of genuine short-duration measurements, leading to biased or unreliable analytical outcomes and thereby undermining downstream tasks. In this work, we compare five subsampling approaches using a benchmark set of 156 genuine replicate spectra collected with a high-purity germanium detector. We evaluate each method with respect to run-to-run variance, channel-to-channel variance, and preservation of total counts (losslessness). Across a wide range of subsampling ratios, only binomial subsampling without replacement consistently reproduces the statistical properties of genuine short-duration spectra, maintaining proper dispersion even in sparse spectral regions and perfectly preserving total counts. These results provide a mathematically principled and practically validated framework for generating synthetically shortened spectra when true short-duration measurements are unavailable.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Application Potential of a Dew-Point Cooling Tower in Selected Energy Intensive Applications in Temperate Climate

In the article, the application potential of the dew-point cooling tower (DPCT) in selected energy-intensive applications in temperate climates was analyzed and discussed. The applications selected for analysis are power generation with natural gas turbines and chilled water air conditioning systems. The study is based on a mathematical model derived from a modified ε-NTU model. The model was validated against experimental results and showed satisfactory agreement with the experimental data. DPCT was compared with a typical cooling tower limited by the wet-bulb temperature (wet-bulb cooling tower, WBCT). The simulation results showed that DPCT is able to provide significant energy savings in energy-intensive applications; therefore, its application potential in temperate climates can be considered justified. In the case of gas turbines, DPCT was able to generate 2 to 10 percentage points more capacity than operating on outdoor air and 1.8 to 5 percentage points more than operating with WBCT. In the case of air conditioning systems, the system equipped with DPCT achieved EERs (energy efficiency ratios) higher by 1 to 7.2 compared to dry cooling and by 0.3 to 5.1 compared to systems equipped with WBCT. The annual energy savings obtained by the system with DPCT were 14.7 MWh compared to WBCT and 30 MWh compared to dry cooling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Modeling of Induction Motors and Variable Speed Drives for Multi-Domain System Simulations Using Modelica and the OpenIPSL Library

This paper introduces an innovative method for characterizing, implementing, and validating both three-phase and single-phase induction motor models, accompanied by a variable speed drive model. The primary goal is to investigate interactions between the electrical power grid and other dynamic domains (e.g., thermofluidic) that impact motor/load drive behavior. Our approach involves establishing a mechanical interface based on a physically meaningful equation linking motor torque/speed to the electrical model in the phasor domain. This allows seamless integration of diverse domain subsystems into a unified multi-domain model using Modelica v4.0.0 and the OpenIPSL library v3.0.1, overcoming co-simulation limitations. The proposed model, which requires only one Modelica-compliant tool for simulation, introduces additional dynamics through the mechanical interface, enabling explicit simulation of load disturbances based on constitutive physics. This deepens our understanding of dynamic interactions between the electrical power domain and other subsystems connected through the motor. We detail the modeled components using mathematical equations and textual descriptions, emphasizing the Modelica modeling approach. Simulation examples validate the implementation, demonstrating the multi-domain modeling capabilities of the newly developed components.

Computer Science↗
DOE OSTI2024

An Overview of Electric Vehicle Load Modeling Strategies for Grid Integration Studies

The adoption of electric vehicles (EVs) has emerged as a solution to reduce greenhouse gas emissions in the transportation sector, which has motivated the implementation of public policies to promote their use in several countries. However, the high adoption of EVs poses challenges for the electricity sector, as it would imply an increase in energy demand and possible impacts on the power quality (PQ) of the power grid. Therefore, it is important to conduct EV integration studies in the power grid to determine the amount that can be incorporated without causing problems and identify the areas of the power sector that will require reinforcements. Accurate EV load patterns are required for this type of study that, through mathematical modeling, reflect both the dynamic behavior and the factors that influence the decision to recharge EVs. This article aims to present an overview of EVs, examine the different factors considered in the literature for modeling EV load patterns, and review modeling methods. EV load modeling methods are classified into deterministic, statistical, and machine learning. The article shows that each modeling method has its advantages, disadvantages, and data requirements, ranging from simple load modeling to more accurate models requiring large datasets.

Computer Science↗
DOE OSTI2024

Trust-Based Detection and Mitigation of Cyber Attacks in Distributed Cooperative Control of Islanded AC Microgrids

In this study, we address the challenge of detecting and mitigating cyber attacks in the distributed cooperative control of islanded AC microgrids, with a particular focus on detecting False Data Injection Attacks (FDIAs), a significant threat to the Smart Grid (SG). The SG integrates traditional power systems with communication networks, creating a complex system with numerous vulnerable links, making it a prime target for cyber attacks. These attacks can lead to the disclosure of private data, control network failures, and even blackouts. Unlike machine learning-based approaches that require extensive datasets and mathematical models dependent on accurate system modeling, our method is free from such dependencies. To enhance the microgrid’s resilience against these threats, we propose a resilient control algorithm by introducing a novel trustworthiness parameter into the traditional cooperative control algorithm. Our method evaluates the trustworthiness of distributed energy resources (DERs) based on their voltage measurements and exchanged information, using Kullback-Leibler (KL) divergence to dynamically adjust control actions. We validated our approach through simulations on both the IEEE-34 bus feeder system with eight DERs and a larger microgrid with twenty-two DERs. The results demonstrated a detection accuracy of around 100%, with millisecond range mitigation time, ensuring rapid system recovery. Additionally, our method improved system stability by up to almost 100% under attack scenarios, showcasing its effectiveness in promptly detecting attacks and maintaining system resilience. These findings highlight the potential of our approach to enhance the security and stability of microgrid systems in the face of cyber threats.

Computer Science↗
DOE OSTI2024

Convex Relaxations of Maximal Load Delivery for Multi-Contingency Analysis of Joint Electric Power and Natural Gas Transmission Networks

Recent increases in gas-fired power generation have engendered increased interdependencies between natural gas and power transmission systems. These interdependencies have amplified existing vulnerabilities in gas and power grids, where disruptions can require the curtailment of load in one or both systems. Although typically operated independently, coordination of these systems during severe disruptions can allow for targeted delivery to lifeline services, including gas delivery for residential heating and power delivery for critical facilities. To address the challenge of estimating maximum joint network capacities under such disruptions, we consider the task of determining feasible steady-state operating points for severely damaged systems while ensuring the maximal delivery of gas and power loads simultaneously, represented mathematically as the nonconvex joint Maximal Load Delivery (MLD) problem. To increase its tractability, we present a mixed-integer convex relaxation of the MLD problem. Then, to demonstrate the relaxation’s effectiveness in determining bounds on network capacities, exact and relaxed MLD formulations are compared across various multi-contingency scenarios on nine joint networks ranging in size from 25 to 1191 nodes. The relaxation-based methodology is observed to accurately and efficiently estimate the impacts of severe joint network disruptions, often converging to the relaxed MLD problem’s globally optimal solution within ten seconds.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

A Review of the Multiple-Readout Concept and Its Application in an Integrally Active Calorimeter

A comprehensive multi-jet physics program is anticipated for experiments at future colliders. Key physics processes necessitate detectors that can distinguish signals from W and Z bosons and the Higgs boson. Typical examples include channels with or pairs and processes involving new physics in those cases where neutral particles must be disentangled from charged ones due to the presence of W or Z bosons in their final states. Such a physics program demands calorimetric energy resolution at or beyond the limits of traditional calorimetric techniques. Multiple-readout calorimetry, which aims to reduce fluctuations in energy measurements of hadronic showers, is a promising approach. The first part of this article reviews dual- and triple-readout calorimetry within a mathematical framework describing the underlying compensating mechanism. The second part proposes a potential implementation using an integrally active and total absorption detector. This model serves as the basis for several Monte Carlo studies, illustrating how the response of a multiple-readout calorimeter depends on construction parameters. Among the layouts considered, one configuration operating in triple-readout mode shows the potential to achieve an energy resolution approaching .

47 OTHER INSTRUMENTATION↗