Search NASA⌕ Search

SEARCH · Search NASA

Results for “Types Of Problems”

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 541 records · Page 30

A Computationally Efficient Algorithm for Computing Convex Hull Prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may execute side payments when a participant’s cost is not covered through energy sale settlements from locational marginal pricing schemes, or when a participant incurs lost opportunity costs to follow the dispatch signal. Convex hull pricing minimizes these and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. However, computing convex hull prices involves solving either a large-scale linear program - which in turn requires explicit descriptions of market participants’ convex hulls -or the Lagrangian dual of the corresponding non-convex scheduling problem. Here, we propose a computationally feasible and industrially scalable Benders decomposition approach to computing convex hull prices at least an order of magnitude faster than the current state-of-the-art while leveraging recent advances in convex hull formulations for thermal generating units.

61 RADIATION PROTECTION AND DOSIMETRY↗

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗

A Garnet-Type Solid-Electrolyte-Based Molten Lithium–Molybdenum–Iron(II) Chloride Battery with Advanced Reaction Mechanism

Solid-electrolyte-based molten-metal batteries have attracted considerable attention for grid-scale energy storage. Although ZEBRA batteries are considered one of the promising candidates, they still have the potential concern of metal particle growth and ion exchange with the β”-Al 2 O 3 electrolyte. Herein, a Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-electrolyte-based molten lithium–molybdenum–iron(II) chloride battery (denoted as Li–Mo–FeCl 2 ) operated at temperature of 250 °C, comprising a mixture of Fe and LiCl cathode materials, a Li anode, a garnet-type Li-ion ceramic electrolyte, and Mo additive, is designed to overcome these obstacles. Different from conventional battery reaction mechanisms, this battery revolutionarily synchronizes the reversible Fe–Mo alloying–dealloying reactions with the delithiation–lithiation processes, meaning that the porous Mo framework derived from Fe–Mo alloy simultaneously suppresses the growth of pure Fe particles. By adopting a Li anode and a Li-ion ceramic electrolyte, the corrosion problem between the cathode and the solid electrolyte is overcome. With similar battery cost ($12 kWh –1 ), the theoretical energy density of Li–Mo–FeCl 2 battery surpasses that of a Na–FeCl 2 ZEBRA battery over 25%, to 576 Wh kg –1 and 2216 Wh L –1 , respectively. Experimental results further prove this cell has excellent cycling performance (472 mAh g LiCl –1 after 300 cycles, 50 mg active material) and strong tolerance against the overcharge–overdischarge (3–1.6 V) and freezing–thawing (25–250 °C) incidents.

36 MATERIALS SCIENCE↗

Biocatalyst discovery and design for plastics deconstruction: A multi‐scale perspective

Plastic waste accumulation poses significant environmental challenges due to a lack of economical solutions for the molecular deconstruction of diverse synthetic polymers. Biological‐based degradation offers promise but is hindered by the crystallinity, hydrophobicity, and additive complexity of plastics, which restrict biocatalyst access and activity. To address these problems, we propose a multi‐scale framework that combines detailed materials characterization, optimization of plastic‐biomolecular interfacial interactions, and enhancement of biocatalytic kinetics to develop effective plastic‐deconstructing enzymes. This approach leverages principles from reaction kinetics, transport and interfacial phenomena, and enzyme engineering to systematically address barriers across diverse plastic types. Our framework aims to accelerate the discovery and optimization of biocatalysts capable of scalable, selective, and efficient deconstruction of plastic waste. These advances hold potential to enable sustainable biological recycling and upcycling pathways, contributing to global efforts in mitigating plastic pollution and promoting circular material economies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum algorithmic measurement

There has been recent promising experimental and theoretical evidence that quantum computational tools might enhance the precision and efficiency of physical experiments. However, a systematic treatment and comprehensive framework are missing. Here we initiate the systematic study of experimental quantum physics from the perspective of computational complexity. To this end, we define the framework of quantum algorithmic measurements (QUALMs), a hybrid of black box quantum algorithms and interactive protocols. We use the QUALM framework to study two important experimental problems in quantum many-body physics: determining whether a system’s Hamiltonian is time-independent or time-dependent, and determining the symmetry class of the dynamics of the system. We study abstractions of these problems and show for both cases that if the experimentalist can use her experimental samples coherently (in both space and time), a provable exponential speedup is achieved compared to the standard situation in which each experimental sample is accessed separately. Our work suggests that quantum computers can provide a new type of exponential advantage: exponential savings in resources in quantum experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancing the Survivability of Power Systems With Grid-Edge DERs Against DoS Attacks

Power system survivability, defined as the ability of a system to maintain steady-state functionality under varying operational conditions, reflects its resilience against disturbances. While existing research primarily focuses on physical-layer disturbances, the increasing prevalence of grid-edge DERs, which are primarily used for integrating renewable energy, has significantly expanded the cyber attack surface. As a result, operational disruptions caused by cyber threats are posing significant challenges to system survivability and cannot be overlooked. To fill this gap, we redefine system survivability to incorporate the cyber layer’s status and propose a Distributionally Robust Optimization (DRO) approach to enhance power system survivability against potential cyber-physical threats. In this paper, we first analyze the operational guidelines of systems with a high penetration of DERs under various cyber network conditions and redefine survivability in this context. Next, we focus on the most common cyber threat, Denial-of-Service (DoS) attacks, and develop a corresponding attack model. This model allows for the creation of a kernel-based ambiguity set that captures attack uncertainties using historical data. Finally, we transform the proposed DRO model as a tractable optimization problem, with its solution providing an optimal cyber redundancy plan to enhance system survivability in DoS attack scenarios. Simulation results on the IEEE 13-node and 123-node test feeders demonstrate the effectiveness of our proposed model in improving system survivability. This model can also be expanded to include other types of common attacks and serve as a comprehensive planning tool to improve overall cyber physical survival of the system.

cybersecurity↗

Minimum Quench Energy of Nb 3 Sn Wires With High Specific Heat Tape

A major problem of state-of-the-art Nb3Snaccelerator magnets is their long training due to thermo-mechanical perturbations. Increasing the specific heat, Cp, of the Rutherford cable would reduce and/or eliminate training by limiting the coils temperature rise. This paper studies feasibility of increasing the Cp of Rutherford-type cables by using thin composite Cu/Gd 2 O 3 and Cu/Gd 2 O 2 S tapes produced by HyperTech Research, Inc. The tape can be either wrapped around the cable, placed on the cable wide faces under the insulation, and/or inserted as a core. Wire samples outfitted with these high-Cp ribbons, or tapes, were prepared and tested at FNAL for their Minimum Quench Energy (MQE). At 90%Ic and 15 T, the average gain of MQE of the Nb 3 Sn wire soldered to the Cu/Gd 2 O 2 S 55 μm thick ribbon was 2.5, and further increased at larger transport currents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Switching Device-Cognizant Sequential Distribution System Restoration

This paper presents an optimization framework for sequential reconfiguration using an assortment of switching devices and repair process in distribution system restoration. Compared to existing studies, this paper considers types, capabilities and operational limits of different switching devices, making it applicable in practice. We develop a novel multi-phase method to find the optimal sequential operation of various switching devices and repair faulted areas. We consider circuit breakers, reclosers, sectionalizers, load breaker switches, and fuses. The switching operation problem is decomposed into two mixed-integer linear programming (MILP) subproblems. The first subproblem determines the optimal network topology and estimates the number of steps to reach that topology, while the second subproblem generates a sequence of switching operations to coordinate the switches. For repairing the faults, we design an MILP model that dispatches repair crews to clear faults and replace melted fuses. After clearing a fault, we update the topology of the network by generating a new sequence of switching operations, and the process continues until all faults are cleared. To improve the computational efficiency, a network reduction algorithm is developed to group line sections, such that only switchable sections are present in the reduced network. The proposed method is validated on the IEEE 123-bus and 8500-bus systems.

distribution system↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗

Exploration of Options for Ultra-Fast Mass Spectroscopy for Pulsed Power Applications

Mega-ampere class pulsed power machines drive intense currents into small volumes to study high energy and density environments. Power lost during these events is a difficult and paramount problem to solve. For example, facilities such as Sandia National Laboratories’ Z machine experience meaningful power loss, which can be linked to non-linear ohmic heating at high currents (i.e., 26 MA on Z) leading to thermal desorption of contaminants and subsequent shunt plasma formation. Characterizing and understanding this type of thermal desorption is key to design optimizations necessary to minimize current loss, which will be even more important for next generation pulsed power. This type of characterization requires the ability to identify and determine concentration of analytes with nanosecond resolution given the pulse width of Z is on the order of 100 ns. This report summarizes progress on a small exploratory project focused on investigating options to meet this challenge using mass spectrometry. The main focus of these efforts utilized an Energy and Velocity Analyzer for Distributions of Electric Rockets intending to determine how quickly transient data could be resolved. This probe combines an electrostatic analyzer with a Wien velocity filter (ExB) to obtain ion energy and velocity distributions. Primary results from this exploratory project indicate significant additional work is needed to demonstrate a nanosecond time scale mass spectrometer for this application and also highlight that alternative detection methods such as laser-based diagnostics should be considered to meet the need for ultra-fast detection.

42 ENGINEERING↗

PyApprox: A software package for sensitivity analysis, Bayesian inference, optimal experimental design, and multi-fidelity uncertainty quantification and surrogate modeling

PyApprox is a Python-based one-stop-shop for probabilistic analysis of numerical models such as those used in the earth, environmental and engineering sciences. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent a wide range of methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of new or existing algorithms for a wide range of model analyses. Here, this paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in groundwater.

54 ENVIRONMENTAL SCIENCES↗

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

The high resolution Rydberg absorption spectrum of 2-butyne C 4 H 6 recorded previously at the SOLEIL synchrotron facility has been interpreted using multichannel quantum defect theory (MQDT). The calculations are based on the continuum scattering calculations of Xu et al., J. Chem. Phys. 136, 154303 (2012) and of Jacovella et al., J. Phys. Chem. A 119, 12339 (2015) pertaining to the dipole-allowed excited state symmetries in absorption from the ground state. In contrast to the traditional approach of calculating low-lying electronic states first and then attempting to extend the calculations to ever higher energy, here the analysis proceeds through the extension of these detailed calculations of the electronic continuum scattering down into the discrete region of the spectrum. The continuum reaction matrices and dipole transition moments are adapted to the discrete Rydberg region via the use of an energy-modified formulation of MQDT theory and associated energy dependences of the quantum defects. The analysis reproduces more than 40 Rydberg states from n ≈ 10 down to the 3 d and 4 s levels with an rms error of better than 20 cm −1 . Further, these belong to five Rydberg series with three different molecular symmetries. While the approach predicts many additional series, most of these are calculated and observed to carry only little oscillator strength. The analysis shows that the Rydberg spectrum is dominated by the excitation of an e″ symmetry electron of fδ and gπ type, in line with what previous studies of the above-threshold shape resonance of 2-butyne have shown. The present study is intended to serve as an example showing how first principles continuum calculations may be useful for the interpretation of highly bound discrete states in a range that poses problems for the standard ab initio techniques. The quantitative treatment of the dipole absorption cross sections is deferred to a future paper.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Options for Achieving Cost Reduction in Advanced Reactors through Open Architecture

A key contributor to high capital costs and schedule overruns for new nuclear power plants is lack of standardization, driven by site-specific customization and construction of multiple designs by competing vendors rather than commitment to a single standardized program. While advanced reactor vendors typically individually target repeat construction of standardized units, the many competing designs could exacerbate the problem. “Open Architecture”, the open specification of requirements and interfaces for structures, systems and components (SSCs), has been proposed as a means of promoting standardization, by facilitating existing non-nuclear suppliers to enter the industry and/or allowing SSCs to be configured for more than one reactor within the same technology type. Contracting mechanisms that facilitate information sharing and alignment of incentives between stakeholders may complement such an approach. A preliminary scheme is presented for selection of SSCs for which such strategies could be adopted, based on a vendor make/buy decision model and stakeholder interviews. SSCs are categorized according to number of suppliers and their contribution to the reactor’s competitive edge. SSCs with many potential suppliers and a high contribution to competitive edge may be attractive for widening the supply chain via open specification of system requirements and interfaces, e.g., SSCs in the power island. SSCs with few suppliers and low contribution to competitive edge may be potential avenues for common system specification between vendors, e.g., some of the auxiliary SSCs. Potential cost reductions from such strategies will depend upon the size of the build program and the reactor type.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Data-Driven Framework for Power System Event Type Identification via Safe Semi-Supervised Techniques

Herein this paper investigates the use of phasor measurement unit (PMU) data with deep learning techniques to construct real-time event identification models for transmission networks. Increasing penetration of distributed energy resources represents a great opportunity to achieve decarbonization, as well as challenges in systematic situational awareness. When high-resolution PMU data and sufficient manually recorded event labels are available, the power event identification problem is defined as a statistical classification problem that can be solved by numerous cutting-edge classifiers. However, in real grids, collecting tremendous high-quality event labels is quite expensive. Utilities frequently have a large number of event records without in-depth details (i.e., unlabeled events). To bridge this gap, we propose a novel semi-supervised learning-based method to improve the performance of event classifiers trained with a limited number of labeled events by exploiting the information from massive unlabeled events. In other words, compared to existing data-driven methods, our method requires only a small portion of labeled data to achieve a similar level of accuracy. Meanwhile, this work discusses and addresses the performance degradation caused by class distribution mismatch between the training set and the real applications. Based on the proposed safe learning mechanism, our model does not directly use all unlabeled events during model training, but selectively uses them through a comprehensive evaluation procedure. Numerical studies on a sizable PMU dataset have been used to validate the performance of the proposed method.

42 ENGINEERING↗

PQ axiverse

We show that the strong CP problem is solved in a large class of compactifications of string theory. The Peccei-Quinn mechanism solves the strong CP problem if the CP-breaking effects of the ultraviolet completion of gravity and of QCD are small compared to the CP-preserving axion potential generated by low-energy QCD instantons. We characterize both classes of effects. To understand quantum gravitational effects, we consider an ensemble of flux compactifications of type IIB string theory on orientifolds of Calabi-Yau hypersurfaces in the geometric regime, taking a simple model of QCD on D7-branes. We show that the D-brane instanton contribution to the neutron electric dipole moment falls exponentially in N 4 , with N the number of axions. In particular, this contribution is negligible in all models in our ensemble with N > 17. We interpret this result as a consequence of large N effects in the geometry that create hierarchies in instanton actions and also suppress the ultraviolet cutoff. We also compute the CP breaking due to high-energy instantons in QCD. In the absence of vectorlike pairs, we find contributions to the neutron electric dipole moment that are not excluded, but that could be accessible to future experiments if the scale of supersymmetry breaking is sufficiently low. The existence of vectorlike pairs can lead to a larger dipole moment. Finally, we show that a significant fraction of models are allowed by standard cosmological and astrophysical constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Generalized fiducial inference on differentiable manifolds

We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. Here, this new approach uses generalized fiducial inference (GFI) to obtain a posterior-like distribution on the manifold, without needing to know local parameterizations that map to the constrained space from an unconstrained Euclidean space. Using mathematical tools from Riemannian geometry, we construct a constrained generalized fiducial distribution (CGFD). A Bernstein-von Mises-type result for the CGFD, which provides intuition for how the desirable asymptotic qualities of the unconstrained generalized fiducial distribution are inherited by the CGFD, is provided. To illustrate the practical use of the CGFD, we provide a proof-of-concept example in the context of a linear logspline density estimation problem, and demonstrate that CGFD-based confidence sets exhibit desirable coverage properties via simulation. As an application, we fit a CGFD to COVID-19 case count data from North Carolina, USA.

97 MATHEMATICS AND COMPUTING↗

Fallowed agricultural lands dominate anthropogenic dust sources in California

Air pollution remains a major problem in many parts of California, significantly impacting public health and regional climate. However, the contribution of anthropogenic dust from agricultural sources, among major pollutants in California’s semi-arid Central Valley, remains largely unclear. Here, we used the Cropland Data Layer from the U.S. Department of Agriculture to identify crop types and land use/cover and leveraged satellite-derived estimates of major dust events between 2008 and 2022 over California. We identified fallowed land—an unplanted agricultural land parcel—as a key anthropogenic dust source in California. Specifically, we find that the Central Valley accounts for about 77% of total fallowed land areas in California, where they are associated with about 88% of major anthropogenic dust events. We also find that the geographic coverage of these fallowed lands expanded between 2008 and 2022 with associated increasing anthropogenic dust activities. Additionally, these anthropogenic dust activities are sensitive to the drought severity over the fallowed lands, with potential cumulative effects on downstream dust burden during prolonged multi-year drought conditions. Overall, our results have important implications for public health, including increased risk for Valley fever and for regional climates, such as increases in extreme precipitation and snowmelt over the Sierra Nevada.

Climate sciences↗

Computational design of thermoelectric alloys through optimization of transport and dopability

Alloying is a common technique to optimize the functional properties of materials for thermoelectrics, photovoltaics, energy storage etc. Designing thermoelectric (TE) alloys is especially challenging because it is a multi-property optimization problem, where the properties that contribute to high TE performance are interdependent. In this work, we develop a computational framework that combines first-principles calculations with alloy and point defect modeling to identify alloy compositions that optimize the electronic, thermal, and defect properties. We apply this framework to design n-type Ba 2(1–x) Sr 2x CdP 2 Zintl thermoelectric alloys. Our predictions of the crystallographic properties such as lattice parameters and site disorder are validated with experiments. To optimize the conduction band electronic structure, we perform band unfolding to sketch the effective band structures of alloys and find a range of compositions that facilitate band convergence and minimize alloy scattering of electrons. Here, we assess the n-type dopability of the alloys by extending the standard approach for computing point defect energetics in ordered structures. Through the application of this framework, we identify an optimal alloy composition range with the desired electronic and thermal transport properties, and n-type dopability. Such a computational framework can also be used to design alloys for other functional applications beyond TE.

36 MATERIALS SCIENCE↗