Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical analysis”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

SOMA: Observability, monitoring, and in situ analytics for exascale applications

With the rise of exascale systems and large, data-centric workflows, the need to observe and analyze high performance computing (HPC) applications during their execution is becoming increasingly important. HPC applications are typically not designed with online monitoring in mind, therefore, the observability challenge lies in being able to access and analyze interesting events with low overhead while seamlessly integrating such capabilities into existing and new applications. We explore how our service-based observation, monitoring, and analytics (SOMA) approach to collecting and aggregating both application-specific diagnostic data and performance data addresses these needs. Furthermore, we present our SOMA framework and demonstrate its viability with LULESH, a hydrodynamics proxy application. Then we focus on Astaroth, a multi-GPU library for stencil computations, highlighting the integration of the TAU and APEX performance tools and SOMA for application and performance data monitoring.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.

42 ENGINEERING↗
DOE OSTI2021

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Graphical Model for Fusing Diverse Microbiome Data

This paper develops a Bayesian graphical model for fusing disparate types of count data. The motivating application is the study of bacterial communities from diverse high-dimensional features, in this case, transcripts, collected from different treatments. In such datasets, there are no explicit correspondences between the communities and each corresponds to different factors, making data fusion challenging. We introduce a flexible multinomial-Gaussian generative model for jointly modeling such count data. This latent variable model jointly characterizes the observed data through a common multivariate Gaussian latent space that parameterizes the set of multinomial probabilities of the transcriptome counts. The covariance matrix of the latent variables induces a covariance matrix of co-dependencies between all the transcripts, effectively fusing multiple data sources. We present a computationally scalable variational Expectation-Maximization (EM) algorithm for inferring the latent variables and the parameters of the model. Here, the inferred latent variables provide a common dimensionality reduction for visualizing the data and the inferred parameters provide a predictive posterior distribution. In addition to simulation studies that demonstrate the variational EM procedure, we apply our model to a bacterial microbiome dataset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components

Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.

42 ENGINEERING↗
DOE OSTI2024

Case Study: Leveraging GenAI to Build AI-based Surrogates and Regressors for Modeling Radio Frequency Heating in Fusion Energy Science

This work presents a detailed case study on using Generative AI (GenAI) to develop AI surrogates for simulation models in fusion energy research. The scope includes the methodology, implementation, and results of using GenAI to assist in model development and optimization, comparing these results with previous manually developed models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gate Set Tomography

Gate set tomography (GST) is a protocol for detailed, predictive characterization of logic operations (gates) on quantum computing processors. Early versions of GST emerged around 2012-13, and since then it has been refined, demonstrated, and used in a large number of experiments. This paper presents the foundations of GST in comprehensive detail. The most important feature of GST, compared to older state and process tomography protocols, is that it is calibration-free . GST does not rely on pre-calibrated state preparations and measurements. Instead, it characterizes all the operations in a gate set simultaneously and self-consistently, relative to each other. Long sequence GST can estimate gates with very high precision and efficiency, achieving Heisenberg scaling in regimes of practical interest. In this paper, we cover GST's intellectual history, the techniques and experiments used to achieve its intended purpose, data analysis, gauge freedom and fixing, error bars, and the interpretation of gauge-fixed estimates of gate sets. Our focus is fundamental mathematical aspects of GST, rather than implementation details, but we touch on some of the foundational algorithmic tricks used in the pyGSTi implementation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Characterization and effects of impurities on carbonate quantification in heterogenous matrices

Mineral carbonation simulates a natural weathering phenomena by breaking down silicates and oxides to form Ca & Mg carbonates. Various mineralization methods have been demonstrated as a potential technique to improve the quality of slags and tailings through neutralization and stabilization of problematic species to yield a product better suited for use in concrete. This study aims to characterize, quantify and analyze carbonates in various carbonated products such as mineralized CaCO 3 , CO 2 mineralized Steel Slags and Mine Tailings. More than 10 samples were analyzed for carbonate measurement and verification from industrial and academic partners that pioneer commercial CO 2 mineralization technologies. The samples were characterized primarily by using X-ray diffraction (XRD), Thermogravimetric Analyses (TGA), and Scanning Electron Microscopy (SEM) to gain insights into CaCO 3 content. A baseline characterization of lab-grade CaCO 3 and MgCO 3 also revealed important considerations for CaCO 3 measurement using TGA alone. The experiments using synthetic lab-grade samples also revealed that the presence of MgCO 3 /MgO can accelerate the decomposition of CaCO 3 and thus can affect measurement parameters. Lab-grade CaCO 3 samples dosed into steel slag and mine tailing also showed significant deviation in their decomposition behavior. These insights are used to inform the development of a standardized protocol for the measurement and verification of carbonate-bearing products.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant At or Near Delayed Critical

The derivation of a non-invasive prompt neutron decay constant sensitivity coefficient is provided in this work. The computation of the sensitivity coefficient derived in this work does not require modification of Monte Carlo source code and is based on capabilities available in Monte Carlo N-Particle R© Code Version 6.2. The prompt neutron decay constant sensitivity coefficients are calculated for 44-group and 252-group energy structures for specific nuclide-reaction pairs in the Jezebel benchmark experiment. The nuclide-reaction pairs investigated in this work include Pu- 239(n,f), Pu-240(n,f), and Pu-241(n,f). Physical explanations of the sensitivity profiles exhibited by the 252-group energy structure are investigated for the prompt neutron multiplication factor, mean neutron lifetime, and prompt neutron decay constant. The prompt neutron decay constant sensitivity coefficients calculated for the 44-group and 252-group energy structure of Pu-239(n,f) are compared. Lastly, the 44-group energy structure sensitivity coefficients calculated are used for nuclear-data induced uncertainty quantification of the neutron multiplication factor. This work shows that a reduction in the nuclear data-induced uncertainty of the neutron multiplication factor is possible for all nuclide-reaction pairs investigated when prompt neutron decay constant sensitivity coefficients are utilized. This is important for new critical experiment design optimization studies of measurement configurations. This work provides a basis for more detailed sensitivity analysis and uncertainty quantification of nuclide, reaction, and energy-specific cross section data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Beyond PCA: Additional Dimension Reduction Techniques to Consider in the Development of Climate Fingerprints

Abstract Dimension reduction techniques are an essential part of the climate analyst’s toolkit. Due to the enormous scale of climate data, dimension reduction methods are used to identify major patterns of variability within climate dynamics, to create compelling and informative visualizations, and to quantify major named modes such as El Niño–Southern Oscillation. Principal components analysis (PCA), also known as the method of empirical orthogonal functions (EOFs), is the most commonly used form of dimension reduction, characterized by a remarkable confluence of attractive mathematical, statistical, and computational properties. Despite its ubiquity, PCA suffers from several difficulties relevant to climate science: high computational burden with large datasets, decreased statistical accuracy in high dimensions, and difficulties comparing across multiple datasets. In this paper, we introduce several variants of PCA that are likely to be of use in climate sciences and address these problems. Specifically, we introduce non-negative , sparse , and tensor PCA and demonstrate how each approach provides superior pattern recognition in climate data. We also discuss approaches to comparing PCA-family results within and across datasets in a domain-relevant manner. We demonstrate these approaches through an analysis of several runs of the E3SM climate model from 1991 to 1995, focusing on the simulated response to the Mt. Pinatubo eruption; our findings are consistent with a recently identified stratospheric warming fingerprint associated with this type of stratospheric aerosol injection.

Weylandt, Michael↗
DOE OSTI2023

Addressing Data Center Cooling Needs through the Use of Subsurface Thermal Energy Storage Systems

This study aims to evaluate the feasibility of addressing the cooling needs for information technology (IT) equipment in data centers by using reservoir thermal energy storage (RTES) to provide reliable and sustainable low-temperature fluid. This project focuses on the technical viability of operating such a system in Houston, Texas, which is a representative location for crypto mining data centers. An analysis has been performed to investigate the technical feasibility with climate data for Houston. Results show that data centers on the scale of 30 MW, and operating at a temperature of 27 degrees C, can be reliably cooled by a combined RTES and dry coolers setup. A techno-economic analysis will be performed and energy/water saving benefits will be quantified in the future. In addition, the study will be extended to other representative data center locations with different climates and geographical locations.

cooling↗
DOE OSTI2025

Transforming Energy Through Computational Excellence: A View From NREL

At the National Renewable Energy Laboratory (NREL)—a U.S. Department of Energy laboratory—computational science, high-performance computing, applied mathematics, advanced computer science, visualization, and data play a pivotal role in advancing energy abundance, affordability, security, and reliability. From fundamental scientifc discovery to systems engineering and analysis, NREL researchers tackle market-relevant challenges to develop solutions for an independent energy system that is reliable, resilient and secure. Collaborative partnerships with industry, government, and academia ensure that our research remains cutting edge, impactful, applicable, and aligned with real-world energy needs. This special issue of Computing in Science & Engineering highlights exemplary NREL projects where computational tools and methodologies drive discovery and accelerate innovation in scalable and integrated energy systems. The featured articles explore the role of computational modeling, high-performance computing, generative AI, and adaptive computing in advancing independent energy solutions, optimizing sustainability research, and enhancing decision-making for energy solutions using a broad mix of energy technologies. Here, these contributions demonstrate how NREL’s computational research bridges the gap between theoretical advancements and practical implementation, emphasizing interdisciplinary collaboration and a commitment to innovation, with a focus on translating computational excellence into real-world impact, thus accelerate progress toward national energy goals. By showcasing cutting-edge research at the intersection of computational science and energy systems, this issue aims to inspire and inform researchers, practitioners, and policymakers dedicated to shaping a more reliable energy future.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Solar Thermal Energy Planner (STEP 1): A New Decision Support Tool for Solar Industrial Process Heat Applications

Solar thermal technologies are a promising technology to supply low-cost thermal energy to industrial processes, but there are often significant barriers to entry to industrial owners considering these technologies for their energy demands. To overcome this barrier and convey economic value to customers, NREL and Sandia National Laboratories developed Solar Thermal Energy Planner (STEP 1), a new web-based decision support tool for solar industrial process heat systems. At SolarPACES 2024, the STEP 1 tool was still under development; progress, methodologies, and a preliminary case study was presented. With the STEP 1 tool launch in May 2025, in this work, the initial version of the full public tool will be presented with demonstrations of its capabilities using a few case studies. First, the user's process heat needs such as location, process media (e.g., steam, air), process temperature, land availability, electricity and fuel costs, among other parameters. STEP 1 features a mapping interface that allows users to draw land and roof boundaries. The process media and temperature inform technology selection criteria modules that determine the appropriate solar thermal collection technologies, as well as congruent heat transfer media (e.g., hot water, oil, salt). Once the solar thermal technology selected, its nominal thermal production for the given site is characterized using NREL's System Advisor Model (SAM). Then, a modified version of NREL's REopt optimal sizing and dispatch optimization tool determines cost-optimal sizing. Within minutes, the user receives the results of the technoeconomics analysis, including the size and performance of the cost-optimal solar-plus-storage system. The cost of the system is compared to business-as-usual (e.g., an existing, standalone natural gas boiler). Users can download key results to store for sensitivity analyses. Examples of flat plate collector, parabolic trough, and molten salt tower applications with and without PV hybridization for different industrial facility types are presented in this work. The STEP 1 tool aims to reduce barriers to the adoption of solar heating solutions stemming from a lack of familiarity and technical background with solar system design options and costs among industry stakeholders.

14 SOLAR ENERGY↗
DOE OSTI2023

Decomposition of multi-particle azimuthal correlations in Q-cumulant analysis*

Abstract The method of Q -cumulants is a powerful tool for studying the fine details of azimuthal anisotropies in high energy nuclear collisions. This paper presents a new method, based on mathematical induction, to evaluate the analytical form of high-order Q -cumulants. The capability of this method is demonstrated via a toy model that uses the elliptic power distribution to simulate the anisotropic emission of particles, quantified in terms of Fourier flow harmonics . The method can help in studying the large amount of event statistics that can be collected in the future and allow measurements of the very high central moments of the distribution. This can, in turn, facilitate progress in understanding the initial geometry, the input to the hydrodynamic calculations of medium expansion in high energy nuclear collisions, and the constraints on it.

Physics↗
DOE OSTI2021

Creation of the VADER Code in SCALE [Abstract]

The VADER (Validation Analysis Data Evaluation Resource) is a new module in SCALE 6.3 that has been derived from the legacy USLSTATS program. VADER is a tool that allows the determination of bias and bias uncertainty for criticality safety computational methods. The older USLSTATS program, written in Java, existed outside of SCALE and provided tools to calculate only the confidence band with administrative margin (sometimes called USL-1) and the single-sided uniform width closed interval (USL-2). For normality testing it only offered a crude chi- squared normality test that had no user-configurable options and presented a simple pass/no-pass functionality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Validation of non-negative matrix factorization for rapid assessment of large sets of atomic pair distribution function data

The use of the non-negative matrix factorization (NMF) technique is validated for automatically extracting physically relevant components from atomic pair distribution function (PDF) data from time-series data such as in situ experiments. The use of two matrix-factorization techniques, principal component analysis and NMF, on PDF data is compared in the context of a chemical synthesis reaction taking place in a synchrotron beam, applying the approach to synthetic data where the correct composition is known and on measured PDFs from previously published experimental data. The NMF approach yields mathematical components that are very close to the PDFs of the chemical components of the system and a time evolution of the weights that closely follows the ground truth. Lastly, it is discussed how this would appear in a streaming context if the analysis were being carried out at the beamline as the experiment progressed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interactive Supercomputing With Jupyter

Rich user interfaces like Jupyter have the potential to make interacting with a supercomputer easier and more productive, consequently attracting new kinds of users and helping to expand the application of supercomputing to new science domains. For the scientist-user, the ideal rich user interface delivers a familiar, responsive, introspective, modular, and customizable platform upon which to build, run, capture, document, re-run, and share analysis workflows. From the provider or system administrator perspective, such a platform would also be easy to configure, deploy securely, update, customize, and support. Jupyter checks most if not all of these boxes. But from the perspective of leadership computing organizations that provide supercomputing power to users, such a platform should also make the unique features of a supercomputer center more accessible to users and more composable with high performance computing (HPC) workflows. Project Jupyter’s core design philosophy of extensibility, abstraction, and agnostic deployment, has allowed HPC centers like NERSC to bring in advanced supercomputing capabilities that can extend the interactive notebook environment. This has enabled a rich scientific discovery platform, particularly for experimental facility data analysis and machine learning problems.

97 MATHEMATICS AND COMPUTING↗