Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated two new models into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model and a Euler-Euler model. The Euler-Euler model offers higher fidelity, while the mixture drift-flux model provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model's accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase models implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, both models are applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of both models for capturing the two-phase flow dynamics critical to MSR operations.

42 - ENGINEERING↗
DOE OSTI2024

FORESTR: Finding, Organizing, Representing, Explaining, Summarizing, and Thinning Random forests

Random forests have become popular models used for data driven predictions. As a result, random forests are currently used or being considered for high-consequence mission applications in national security, such as the prediction of yield from optical signals and malware detection. While random forests may provide accurate predictions, the complexity of the algorithm causes a lack of interpretability. Random forests are an ensemble of regression or decision trees. Individual regression and decision trees are interpretable, but ensembles are inherently difficult to interpret due to the compilation of many models. We aim to increase the interpretability of random forests by finding patterns in the ensemble of trees that can be used to “thin” (or remove) trees. As a starting point, in this report, we develop a new distance metric for quantifying the similarity between trees based on their topologies (i.e., shapes). We base the metric on a novel distance metric for graphs that is a proper mathematical distance, is invariant to transformations, has registration between graphs, and computes topological evolutions between graphs. We use the tree distance metric to compute tree statistics such as a “mean tree” and to identify clusters of trees. We apply the developed methodology to a toy dataset and a mission relevant product inspection dataset to demonstrate how the metric can provide insight into random forests. Furthermore, we discuss the limitations of the approach and ideas for future research into how the metric could be used as a thinning tool to develop less complex models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Resilient Consensus State Observer for Nonlinear Systems and Against Attacks

Herein, we investigate the design of a resilient consensus observer that yields correct and continuous state estimates of a nonlinear control system whose state observations are subject to attacks. The observations are assumed to be obtained from multiple sensors and/or communication channels employed to provide redundancy. It is shown that, in an attack-free setting, individual globally convergent observers can be designed for each sensor/channel using the corresponding state-dependent Jacobian system. To counter attacks, a resilient consensus observer is then designed by merging the observers estimates to produce a smooth consensus estimate. Resilient convergence of the consensus estimate to the true state is ensured under the condition that the sparsity of the attacks being less than a half of the number of redundant channels. As an illustrative example, the proposed scheme is successfully applied to the chaotic circuit synchronization problem in which their synchronization has been attacked.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling suggests SARS-CoV-2 rebound after nirmatrelvir-ritonavir treatment is driven by target cell preservation coupled with incomplete viral clearance

In a subset of SARS-CoV-2-infected individuals treated with the antiviral nirmatrelvir-ritonavir, the virus rebounds following treatment. The mechanisms driving this rebound are not well understood. We used a mathematical model to describe the longitudinal viral load dynamics of 51 individuals treated with nirmatrelvir-ritonavir, 20 of whom rebounded. Target cell preservation, either by a robust innate immune response or initiation of N-R near the time of symptom onset, coupled with incomplete viral clearance, appears to be the main factor leading to viral rebound. Moreover, the occurrence of viral rebound is likely influenced by the time of treatment initiation relative to the progression of the infection, with earlier treatments leading to a higher chance of rebound. A comparison with an untreated cohort suggests that early treatments with nirmatrelvir-ritonavir may be associated with a delay in the onset of an adaptive immune response. Nevertheless, our model demonstrates that extending the course of nirmatrelvir-ritonavir treatment to a 10-day regimen may greatly diminish the chance of rebound in people with mild-to-moderate COVID-19 and who are at high risk of progression to severe disease. Altogether, our results suggest that in some individuals, a standard 5-day course of nirmatrelvir-ritonavir starting around the time of symptom onset may not completely eliminate the virus. Thus, after treatment ends, the virus can rebound if an effective adaptive immune response has not fully developed. These findings on the role of target cell preservation and incomplete viral clearance also offer a possible explanation for viral rebounds following other antiviral treatments for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Accidental symmetries, Hilbert series, and friends

Accidental symmetries in effective field theories can be established by computing and comparing Hilbert series. This invites us to study them with the tools of invariant theory. Applying this technology, we spotlight three classes of accidental symmetries that hold to all orders for non-derivative interactions. They are broken by derivative interactions and become ordinary finite-order accidental symmetries. To systematically understand the origin and the patterns of accidental symmetries, we introduce a novel mathematical construct — a (non-transitive) binary relation between subgroups that we call friendship. Equipped with this, we derive new criteria for all-order accidental symmetries in terms of friends, and criteria for finite-order accidental symmetries in terms of friends ma non troppo. They allow us to verify and identify accidental symmetries more efficiently without computing the Hilbert series. We demonstrate the success of our new criteria by applying them to a variety of sample accidental symmetries, including the custodial symmetry in the Higgs sector of the Standard Model effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Entanglement Throughput Measurements and Capacity Estimates for Aerial-Inground Fiber

The throughput of entangled qubit pairs per second (eqps) is a basic performance metric of quantum networks that provide the entanglement distribution capability. Over fiber connections, it is measured using specialized instruments, including photonic entanglement sources and single-photon detectors. Extensive theory has been developed to estimate the capacity of a generic quantum channel, which in turn is applied to estimate the maximum achievable eqps over a fiber connection. However, there is a gap in relating these two performance metrics, in part due to their disparate nature, namely, mathematical formulae of the channel capacity and specialized instrumentation for eqps measurements. We describe eqps measurements collected over testbed connections composed of aerial and inground fiber of lengths up to 45 km. We estimate the normalized capacity using the transmissivity parameter derived from single photon detector measurements. These estimates are then converted to capacity bounds on the connection's eqps using the source eqps rate derived from co-located detector measurements. The results indicate a consistency between eqps measurements and their capacity estimates, and provide insights into relating the parameters of analytic capacity estimates to physical measurements.

Rao, Nageswara [ORNL] (ORCID:0000000234085941)↗
DOE OSTI2025

Reliable and Efficient Machine Learning (Final Technical Report)

Modern scientific experiments generate massive amounts of data at a pace much faster than humans can manually analyze. While machine learning has revolutionized commercial data analysis (such as recommending movies or recognizing faces), applying these tools to complex scientific discovery is challenging because scientific answers must be precise, interpretable, and adhere to physical laws. The research under this project aims to develop new mathematical tools and computer algorithms specifically designed for scientific applications. Major progress has been made in automatically cleaning and deconstructing messy experimental data, analyzing the visual information of physical phenomena, determining the underlying physical variables, and providing rig orous mathematical analysis of interesting algorithms and concepts widely used in machine learning. This project addressed the critical gap between our ability to generate massive scientific data and our ability to extract interpretable information from it. We established mathematical foundations for Scientific Machine Learning (SciML) aimed at effective data analytics and automated discovery. Our work focused on three core objectives: (1) developing reliable feature extraction methods for dynamic high-dimensional data, (2) establishing mathematical foundations for discovering dynamics via neural networks, and (3) creating rigorous optimization techniques for these models. Key outcomes come from two fronts. On the practical side, they include the development of algorithms that significantly enhance the extraction of signals from field data, as well as the capability to handle situations that exhibit smooth variations or physical stretching due to temperature changes. They also include the creation of an automated framework for discovering fundamental state variables from raw experimental data, demonstrating the ability to identify intrinsic physical dimensions without prior knowledge of the governing laws. On the theoretical front, the research results in theoretical advances in Optimal Transport, a widely used notion in SciML, specifically regarding functions with fixed-size nodal sets, provide sharp bounds relevant to uncertainty quantification. Meanwhile, the outcomes also include the establishment of convergence theories for nonlocal gradient descent methods, enabling robust optimization with noisy data in high-dimensional settings commonly encountered in scientific modeling. The project also helps creating opportunities to train the next generation of researchers, equipping them with the necessary technical skills for today’s workplace and preparing them for future advances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modeling Analysis of Ball-Milling Process for Battery-Electrode Synthesis

The mechanical alloying process is a promising method for synthesizing electrode materials for batteries owing to its benefits such as the ability to produce nanostructured, high-performing electrode alloys, no adverse effects on the solid electrolyte for solid-state batteries, stable production of thick electrodes, simple processing steps, and low processing costs. It is gaining intensive attention in the battery industry as one of the best methods to replace the conventional wet-slurry-solvent method, and its application is rapidly increasing these days. However, the operation is currently conducted purely based on trial-and-error methods without fully utilizing the features of its functions. Here, this may be attributed to a lack of understanding of the effect of operating parameters on the alloying process and final products. Surprisingly, there is a scarcity of the literature conducting fundamental research to comprehend the underlying physics of the entire mechanical alloying process, resulting in a significant knowledge gap. To address this knowledge gap, extensive research was conducted. The existing literature on mechanical alloying was reviewed to comprehend the current state of understanding and to discuss the direction for future research. Mathematical expressions were developed to create physics-based models capable of capturing the entire mechanical alloying process, including milling kinetics and defect-enhanced phase evolution. These methods were then applied to investigate the impact of operating parameters such as milling frequency, initial mole ratio of the alloyed materials, density of grinding balls, and energy required for the powders to become amorphous (i.e., the amorphization energy threshold). This research aimed not only to comprehend the direct effects of these operating parameters but also to unveil the physics underlying the ball-milling process. The results of our study can serve as crucial information for the battery industry in designing or operating the ball-milling process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

An Introduction to Variational Calculus

Variational Calculus is an advanced topic in mathematics. This document exists to compile the basics of Variational Calculus in a readable and digestible manner for application by engineers. This document is a living document and will be updated as time moves forward. It is the hope that by reading this document, the basics of Variational Calculus can be readily applied to any system/function of interest. This does not serve as a replacement for textbooks or other learned sources but should act as a companion piece to better serve the reader in times of confusion and distress.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A tri-level optimization model for interdependent infrastructure network resilience against compound hazard events

Resilient operation of interdependent infrastructures against compound hazard events is essential for maintaining societal well-being. To address consequence assessment challenges in this problem space, we propose a novel policy-guided tri-level optimization model applied to a proof-of-concept case study with fuel distribution and transportation networks – encompassing one realistic network; one fictitious, yet realistic network; as well as networks drawn from three synthetic distributions. Mathematically, our approach takes the form of a defender-attacker-defender (DAD) model—a multi-agent tri-level optimization, comprised of a defender, attacker, and an operator acting in sequence. Here, in this study, our notional operator may choose proxy actions to operate an interdependent system comprised of fuel terminals and gas stations (functioning as supplies) and a transportation network with traffic flow (functioning as demand) to minimize unmet demand at gas stations. A notional attacker aims to hypothetically disrupt normal operations by reducing supply at the supply terminals, and the notional defender aims to identify best proxy defense policy options which include hardening supply terminals or allowing alternative distribution methods such as trucking reserve supplies. We solve our DAD formulation at a metropolitan scale and present practical defense policy insights against hypothetical compound hazards. We demonstrate the generalizability of our framework by presenting results for a realistic network; a fictitious, yet realistic network; as well as for three networks drawn from synthetic distributions. Additionally, we demonstrate the scalability of the framework by investigating runtime performance as a function of the network size. Steps for future research are also discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗
DOE OSTI2025

“When I talk about it, my eyes light up!” Impacts of a national laboratory internship on community college student success

Participation in technical/research internships may improve undergraduate graduation rates and persistence in science, technology, engineering, and mathematics (STEM), yet little is known about the benefits of these activities a) for community college students, b) when hosted by national laboratories, and c) beyond the first few years after the internship. We applied Social Cognitive Career Theory (SCCT) to investigate alumni perspectives about how CCI at Lawrence Berkeley National Laboratory (LBNL) impacted their academic/career activities. We learned that alumni had low confidence and expectations of success in STEM as community college students. Participation in CCI increased their professional networks, expectations of success, and STEM skills, identity, and self-efficacy/confidence. Hispanic/Latinx alumni recalled the positive impact of mentors who prioritized personal connections, and women valued "warm" social environments. We propose several additions to the SCCT model, to better reflect the supports and barriers to STEM persistence for community college students.

careers↗
DOE OSTI2024

Distributed Augmentation, Hypersweeps, and Branch Decomposition of Contour Trees for Scientific Exploration

Contour trees describe the topology of level sets in scalar fields and are widely used in topological data analysis and visualization. A main challenge of utilizing contour trees for large-scale scientific data is their computation at scale using highperformance computing. To address this challenge, recent work has introduced distributed hierarchical contour trees for distributed computation and storage of contour trees. However, effective use of these distributed structures in analysis and visualization requires subsequent computation of geometric properties and branch decomposition to support contour extraction and exploration. In this work, we introduce distributed algorithms for augmentation, hypersweeps, and branch decomposition that enable parallel computation of geometric properties, and support the use of distributed contour trees as query structures for scientific exploration. Finally, we evaluate the parallel performance of these algorithms and apply them to identify and extract important contours for scientific visualization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗
DOE OSTI2025

Auroral and Non‐Auroral H 3 + Ion Winds at Uranus With Keck‐NIRSPEC and IRTF‐iSHELL

Abstract To date, no investigation has documented ionospheric flows at Uranus. Previous investigations of Jupiter and Saturn have demonstrated that mapping ion winds can be used to understand ionospheric currents and how these connect to magnetosphere‐ionosphere coupling. We present a study of Uranus's near infrared emissions (NIR) using data from the Keck II Telescope's Near InfraRed SPECtrograph (NIRSPEC) and the InfraRed Telescope Facility's iSHELL spectrograph. H 3 + emission lines were used to derive dawn‐to‐dusk intensity, ionospheric temperatures and ion densities to identify auroral emissions, with their Doppler shifts used to measure ion velocities. We confirm the presence of the southern NIR aurora in 2016, driven by elevated H 3 + column densities up to 6.0 × 10 16 m −2 . While no auroral emissions were detected in 2014, we find a 14%–20% super rotation across the planet's disk in 2014 and a 7%–18% super rotation in 2016.

Thomas, Emma M. [Department of Mathematics Physics↗
DOE OSTI2026

BoBa

BoBa is a C++ software library for working with large matrices, tensors, and tensor decompositions. The library provides tools for dense matrix and tensor operations, tensor decompositions, and tensor decomposition methods that support modern CPU and GPU architectures. It includes portable abstractions for linear algebra, tensor algebra, and multidimensional computation. BoBa is intended for scientific computing applications that involve large multidimensional data sets or high dimensional mathematical models. Its capabilities support tasks such as data compression, linear algebra, efficient numerical computation, and the development of scalable algorithms for heterogeneous hardware. Tutorials, tests, and example applications are included to help users learn and apply the library.

Yao, Jin [Lawrence Livermore National Laboratory (↗