Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICS - BOUNDARY VALUES”

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

Power Management Solution for Growing Loads at Airport Rental Facilities

Growing electrical loads across airport facilities could expand beyond local grid capacities. Coordinating flexible loads with generation and storage assets can reduce coincident peaks, mitigate costly upgrades, and reduce electricity costs. This report outlines how The Athena Team has developed a suite of tools, resources, and solutions to leverage load flexibility with energy assets to harness this value potential.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Multi-channel, multi-template event reconstruction for SuperCDMS data using machine learning

SuperCDMS SNOLAB uses kilogram-scale germanium and silicon detectors to search for dark matter. Each detector has Transition Edge Sensors (TESs) patterned on the top and bottom faces of a large crystal substrate, with the TESs electrically grouped into six phonon readout channels per face. Noise correlations are expected among a detector's readout channels, in part because the channels and their readout electronics are located in close proximity to one another. Moreover, owing to the large size of the detectors, energy deposits can produce vastly different phonon propagation patterns depending on their location in the substrate, resulting in a strong position dependence in the readout-channel pulse shapes. Both of these effects can degrade the energy resolution and consequently diminish the dark matter search sensitivity of the experiment if not accounted for properly. We present a new algorithm for pulse reconstruction, mathematically formulated to take into account correlated noise and pulse shape variations. This new algorithm fits N readout channels with a superposition of M pulse templates simultaneously - hence termed the N$\times$M filter. We describe a method to derive the pulse templates using principal component analysis (PCA) and to extract energy and position information using a gradient boosted decision tree (GBDT). We show that these new N$\times$M and GBDT analysis tools can reduce the impact from correlated noise sources while improving the reconstructed energy resolution for simulated mono-energetic events by more than a factor of three and for the 71Ge K-shell electron-capture peak recoils measured in a previous version of SuperCDMS called CDMSlite to $<$ 50 eV from the previously published value of $\sim$100 eV. These results lay the groundwork for position reconstruction in SuperCDMS with the N$\times$M outputs.

Albakry, M. F. [British Columbia U.; TRIUMF]↗
DOE OSTI2026

Validation of an Erythema-Weighted UV Model Using Broadband Solar Irradiance Measurements From Eleven U.S. Sites: Preprint

Erythema-weighted UV solar irradiance (UV-E) has a potential impact on human health if the recommended maximum exposure times are exceeded. In spite of this, it is not measured at most sites that measure Global Horizontal Irradiance (GHI). However, since UV-E is highly correlated with GHI and total ozone content it can be estimated from this information with sufficient accuracy to assisst in public health recommendations. The Power Model (PM) provides a simple method to estimate the erythema-weighted UV irradiance (UV-E) from measured GHI, total ozone column and air mass. In this work, the performance of the PM method is assessed using high-quality data from 11 sites in the continental U.S. (part of SURFRAD and SOLRAD networks) and total ozone estimates publicly available from the MERRA-2 re-analysis database. A three year period (2021- 2023) at 1-minute frequency is considered. The results (for time aggregations of 5 and 60 minutes) show a high Pearson's correlations (> 0.99), consistently positive mean bias deviations (below 13%) and dispersions in the 8-19% range at all sites. Relative values are expressed in terms of the corresponding measurement mean. The performance indicators remain consistent across time resolutions (5 or 60 minutes), suggesting that the model's performance is robust and not significantly affected by short-term variability (which is captured by GHI). Spatial patterns reveal higher biases and RMSD in northern and eastern locations. These values represent a significant improvement over widely used satellite-based global UV-E estimates and open the possibility of using the PM with satellite-based GHI estimates for operational UV-E mapping over the contiguous U.S. territory.

14 SOLAR ENERGY↗
DOE OSTI2025

Selection of Global Climate Model Data for Downscaling With Generative Machine Learning and Use in the Power Planning for Alignment of Climate and Energy Systems Project

The range of results from climate models and scenarios is important to the understanding of uncertainty in power planning analysis. A U.S. Department of Energy-funded analytic project called Power Planning for Alignment of Climate and Energy Systems is developing data and analytic methods to reflect the effects of climate change on key variables for power system planning, as part of the Grid Modernization Lab Consortium. This project will select and prepare global climate model results for use in power system planning models. A related report (Evaluation of Global Climate Models for Use in Energy Analysis) assesses the performance of various global climate models from the Coupled Model Intercomparison Project Phase 6 data archive for their historical skill with respect to energy system performance and for their future projections under multiple climate change scenarios. Building from that report, we describe the selection of a climate scenario (Shared Socioeconomic Pathway [SSP] 2-4.5) and five climate models: TaiESM1, EC-Earth3-CC, GFDL-CM4, EC-Earth3-Veg, and MPI-ESM1-2-HR. We describe the model selection criteria, which were based on the quality of the match between model results under historical conditions and on the representation of the range of future values for several variables. These results will be downscaled via an open-source generative machine learning method called Super-Resolution for Renewable Energy Resource Data with Climate Change Impacts.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2027

Nodal capacity expansion planning with flexible large-scale load siting

We propose explicitly incorporating large-scale load siting into a stochastic nodal power system capacity expansion planning model that concurrently co-optimizes generation, transmission, and storage expansion. The potential operational flexibility of some of these large loads is also taken into account by considering them as consisting of a set of tranches with different reliability requirements, which are modeled as a constraint on expected served energy across operational scenarios. We implement our model as a two-stage stochastic mixed-integer optimization problem with cross-scenario expectation constraints. To overcome the challenge of scalability, we build upon existing work to implement this model on a high performance computing platform and exploit scenario parallelization using an augmented Progressive Hedging Algorithm. The algorithm is implemented using the bounding features of mpisppy, which have shown to provide satisfactory provable optimality gaps despite the absence of theoretical guarantees of convergence. We test our approach and assess the value of this proactive planning framework on total system cost and reliability metrics using realistic testcases geographically assigned to San Diego and South Carolina, with datacenter and direct air capture facilities as large loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

DiffESM: Conditional Emulation of Temperature and Precipitation in Earth System Models With 3D Diffusion Models

Earth system models (ESMs) are essential for understanding the interaction between human activities and the Earth's climate. However, the computational demands of ESMs often limit the number of simulations that can be run, hindering the robust analysis of risks associated with extreme weather events. While low-cost climate emulators have emerged as an alternative to emulate ESMs and enable rapid analysis of future climate, many of these emulators only provide output on at most a monthly frequency. This temporal resolution is insufficient for analyzing events that require daily characterization, such as heat waves or heavy precipitation. We propose using diffusion models, a class of generative deep learning models, to effectively downscale ESM output from a monthly to a daily frequency. Trained on a handful of ESM realizations, reflecting a wide range of radiative forcings, our DiffESM model takes monthly mean precipitation or temperature as input, and is capable of producing daily values with statistical characteristics close to ESM output. Combined with a low-cost emulator providing monthly means, this approach requires only a small fraction of the computational resources needed to run a large ensemble. We evaluate model behavior using a number of extreme metrics, showing that DiffESM closely matches the spatio-temporal behavior of the ESM output it emulates in terms of the frequency and spatial characteristics of phenomena such as heat waves, dry spells, or rainfall intensity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An Uncertainty-Informed and High-Fidelity Performance Forecasting Framework for Heliostat Fields

Concentrating Solar Thermal (CST) tower systems employ heliostat fields to direct solar energy to a central receiver, which then transfers the heat either directly to a thermal process (e.g., steam production) or to a thermal energy storage system for future use. Heliostat fields compose a significant proportion of the project costs of a CST tower system and the performance of the heliostats determines a plant's productivity at a given location. While CST characterization tools such as SolarPILOT and System Advisor Model (SAM) include a large collection of inputs that influence the performance of a CST tower system, many are uncertain prior to the development of the project and may have a significant impact on the overall energy delivery and profitability of a project; moreover, the fidelity of these models under default conditions may be insufficient to determine the value of component improvements such as those under development in the Heliostat Consortium. This work introduces a Monte Carlo simulation framework that incorporates uncertainty in key performance parameters to generate confidence intervals and percentile estimates for a CST solar field's energy delivery.

14 SOLAR ENERGY↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantum Computing Market Report

The quantum computing market is entering a rapid growth phase. Growth was originally driven by private investment interest in high value applications such as drug discovery, materials science, and finance. Currently the main market drivers include technical advances, continued private investment, and rising public funding. The quantum computing market is poised for explosive growth — with its market size estimated at $1.42 billion USD in 2024, and an expected compound annual growth rate of 20.5% from 2025 to 2030 [4]. Superconducting qubits are the largest qubit type by revenue. However, alternative approaches such as trapped ion, neutral atom, spin based, and photonic qubits are also attracting investment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anticipating Technical Expertise and Capability Evolution in Research Communities Using Dynamic Graph Transformers

The ability to anticipate global technical expertise and capability evolution trends is essential for national and global security, especially in safety-critical domains such as nuclear nonproliferation (NN) and rapidly emerging fields like artificial intelligence (AI). Here, in this work, we extend traditional statistical relational learning approaches (e.g., link prediction in collaboration networks) and formulate a problem of anticipating technical expertise and capability evolution using dynamic heterogeneous graph representations. We develop novel capabilities to forecast collaboration patterns, authorship behavior, and technical capability evolution at different granularities (e.g., scientist and institution levels) in two distinct research fields. We implement a dynamic graph transformer (DGT) neural architecture, which pushes the state-of-the-art graph neural network models by: 1) forecasting heterogeneous (rather than homogeneous) nodes and edges; and 2) relying on both discrete- and continuous-time inputs. We demonstrate that our DGT models predict collaboration, partnership, and expertise patterns with 0.26, 0.73, and 0.53 mean reciprocal rank values for AI and 0.48, 0.93, and 0.22 for NN domains. DGT model performance exceeds the best-performing static graph baseline models by 30%–80% across AI and NN domains. Our findings demonstrate that DGT models boost inductive task performance when previously unseen nodes appear in the test data for the domains with emerging collaboration patterns (e.g., AI). Specifically, models accurately predict which established scientists will collaborate with early career scientists and vice versa in the AI domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Commercial Building Prototypes Based on ANSI/ASHRAE/IES Standard 90.1-2019 Appendix G PRM: Technical Support Document

The two paths for documenting compliance with ANSI/ASHRAE/IES Standard 90.1 are the prescriptive path and the performance path. Beyond code programs and rating systems (for example, USGBC-LEED ) are primarily known to use a third path – the Appendix G Performance Rating Method. An update in the 2016 edition of Standard 90.1 approved the Appendix G Performance Rating System for code compliance, extending its application and allowing for greater consistency of modeling rules for code and beyond code building energy modeling. The Appendix G PRM provides rules for the development of whole building energy models of baseline and proposed models for calculating the “performance cost index target” value using the simulated energy results of the baseline and proposed models and the building performance factors published in Table 4.2.1.1 of the Standard. This report documents (1) the methodology used for development of the baseline and proposed energy models of the Pacific Northwest National Laboratory and U.S. Department of Energy commercial building prototypes using the Appendix G Performance Rating Method; and (2) the building performance factors that were calculated using those models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On MCNP Stochastic Volume Estimation Normalization

This paper resolves a perennial point of confusion regarding the source-weighting normalization factor recommended in the MCNP manual (𝜋⁢𝑟 2 ) to stochastically estimate the volume of a region within an enclosing inward-directed spherical surface source with radius 𝑟. The normalization factor arises from the relationship between a sphere’ s mean chord length, its volume, and the values estimated by MCNP track-length tallies. A brief derivation is given that relates these quantities and results in the stated normalization. The correctness of this factor is demonstrated by estimating the volume of a variety of convex and nonconvex volumes. A heuristic demonstration of how biasing the inward-directed source reduces the statistical uncertainty of the stochastic volume estimate is also given, but a rigorous analysis of this improvement is left as future work.

42 ENGINEERING↗