Search NASA⌕ Search

SEARCH · Search NASA

Results for “alternating minimization”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Matrix Completion Using Alternating Minimization for Distribution System State Estimation

This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Matrix Completion Using Alternating Minimization for Distribution System State Estimation: Preprint

This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A realistic theory of E6 unification through novel intermediate symmetries

Abstract We propose a non-supersymmetric E 6 GUT with the scalar sector consisting of650⨁351′ ⨁27. Making use of the first representation for the initial symmetry breaking to an intermediate stage, and the latter two representations for second-stage breaking to the Standard Model and a realistic Yukawa sector, this theory represents the minimal E 6 GUT that proceeds through one of the intermediate stages that are novel compared to SU(5) or SO(10) GUT: trinification SU(3) C × SU(3) L × SU(3) R , SU(6) × SU(2) and flipped SO(10) × U(1). We analyze these possibilities under the choice of vacuum that preserves a ℤ 2 “spinorial parity”, which disentangles the chiral and vector-like fermions of E 6 and provides a dark matter candidate in the form of a (scalar) inert doublet. Three cases are shown to consistently unify under the extended survival hypothesis (with minimal fine-tuning): trinification symmetry SU(3) C × SU(3) L × SU(3) R with eitherLRorCRparity, and SU(6) CR × SU(2) L . Although the successful cases give a large range for proton lifetime estimates, all of them include regions consistent with current experimental bounds and within reach of forthcoming experiments. The scenario investigated in this paper essentially represents the unique (potentially) viable choice in the class of E 6 GUTs proceeding through a novel-symmetry intermediate stage, since non-minimal alternatives seem to be intrinsically non-perturbative.

Physics↗
DOE OSTI2023

Plug-and-Play Methods for Integrating Physical and Learned Models in Computational Imaging: Theory, algorithms, and applications

Plug-and-play (PnP) priors constitute one of the most widely used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data fidelity term to promote data consistency and imposing a learned regularizer in the form of an image denoiser. Recent highly successful applications of PnP algorithms include biomicroscopy, computerized tomography (CT), magnetic resonance imaging (MRI), and joint ptychotomography. This article presents a unified and principled review of PnP by tracing its roots, describing its major variations, summarizing main results, and discussing applications in computational imaging. Additionally, we also point the way toward further developments by discussing recent results on equilibrium equations that formulate the problem associated with PnP algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A note on minimizing time assurance tests for repairable systems

We consider assurance testing for repairable systems when supplementary information is available in addition to the data collected in the assurance test. Here, the supplementary information is incorporated using a Bayesian inferential framework. Here we consider assurance testing for a homogeneous Poisson process. In this note we consider an alternative criterion that minimizes the test time while ensuring that the requirements on the producer's and consumer's risks are met. We illustrate the use of this alternative criterion with an example.

42 ENGINEERING↗
DOE OSTI2026

The Kinetic Consequences of Water on Catalytic Methane Pyrolysis

Hydrogen production from biomass and natural gas has emerged as a prominent research area in response to the growing demand for energy from alternative sources that minimize CO 2 emissions. In this study, we investigate the impact of water, which is present in and generated from biomass-derived streams, on carbon nanotube (CNT) growth and hydrogen production during methane decomposition using Ni–Mo/MgO as a catalyst. We reveal here that the role of water on CNT growth is highly complex; its effect depends on the stage of growth at which the water is incorporated. When water is introduced at the beginning of methane decomposition ( t = 0 h), methane conversion rates are negatively impacted. We hypothesize that water inhibits the significant phase changes the Ni–Mo/MgO catalyst undergoes during catalyst carburization. In contrast, the incorporation of a small percentage of water after a stabilization period ( t = 3 h) results in methane conversion rate enhancements that scale with the introduced water partial pressure as water selectively reacts with amorphous carbon deposits that lead to catalyst deactivation, thus prolonging the lifetime of some of the most active sites. Moreover, water incorporation after stabilization significantly reduces the apparent activation energy. Density Functional Theory (DFT) calculations reveal that water preferentially interacts with carbon fragments on the catalyst surface to remove carbon deposits with a barrier lower than that required for methane activation, further supporting its role in cleaning active sites on the catalyst surface. Characterization of the resulting carbon nanotubes reveals the formation of more graphitic materials produced in the presence of water, highlighting the impact of water on nanotube properties. These results provide clarity toward the many ways in which water, or cofeeding of biomass-derived materials, may impact catalytic methane pyrolysis rates.

carbon nanotubes↗
DOE OSTI2026

Technoeconomic and Life Cycle Analysis of an Integrated Fermentation and Microbial Electrochemical Process for Volatile Fatty Acid Production from Food Waste

Techno-economic analysis (TEA) and life cycle assessment (LCA) were conducted for an integrated system designed for the production of volatile fatty acid (VFA) from food waste. The TEA estimated a production cost of $\$$3.12/kg VFA, and the LCA predicted negative greenhouse gas (GHG) emissions of -0.4 kg CO 2 e/kg VFA, driven primarily by diverting organic waste from landfills and avoiding methane emissions while producing valuable chemical products. Hotspot analysis showed arrested methanogenesis (AM) fermentation as the largest contributor to costs (37%) and environmental burden (47%), driven by high sodium hydroxide (NaOH) consumption. Distillation and microbial electrosynthesis (MES) units were the next-largest environmental contributors (28% and 18%). Major cost drivers also included residuals management (biosolids and wastewater) and the equipment and operating costs for AM, MES, and sonication pretreatment units. Although the new integrated system is environmentally benign, its costs and environmental impacts can be further reduced by integrating alternative energy sources, minimizing chemical and energy inputs through process optimization, and improving efficiency. In conclusion, this work highlighted the viability of waste-derived VFA production and provided a clear, data-driven strategy to accelerate the commercialization of waste valorization technology.

Carboxylic Acid Production↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

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

Diamond channel-cut crystals for high-heat-load beam-multiplexing narrow-band X-ray monochromators

Next-generation high-brilliance X-ray photon sources call for new X-ray optics. Here we demonstrate the possibility of using monolithic diamond channel-cut crystals as high-heat-load beam-multiplexing narrow-band mechanically stable X-ray monochromators with high-power X-ray beams at cutting-edge high-repetition-rate X-ray free-electron laser (XFEL) facilities. The diamond channel-cut crystals fabricated and characterized in these studies are designed as two-bounce Bragg reflection monochromators directing 14.4 or 12.4 keV X-rays within a 15 meV bandwidth to 57 Fe or 45 Sc nuclear resonant scattering experiments, respectively. The crystal design allows out-of-band X-rays transmitted with minimal losses to alternative simultaneous experiments. Only ≲2% of the incident ∼100 W X-ray beam is absorbed in the 50 µm-thick first diamond crystal reflector, ensuring that the monochromator crystal is highly stable. Other X-ray optics applications of diamond channel-cut crystals are anticipated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Metallurgical and mechanical characterization of ATP-5 aluminum alloys for minimizing machining distortion

ATP-5 aluminum alloy was considered as an alternative to 6061-T651 aluminum alloy to minimize machining distortion. ATP-5 is a proprietary cast alloy that is compositionally similar to 5083 aluminum and is purported to have excellent machinability, stability, and corrosion response. Dimensional stability tests, mechanical testing, chemical analysis, microstructural analysis, and fractography were completed to understand the metallurgy of the ATP-5 alloy and assess its potential as a 6061-T651 alloy alternative. Ultimately, the ATP-5 was found to retain dimensional stability issues, albeit to a lesser extent, relative to the 6061-T651 alloy. Combined with the residual cast structure observed microstructurally and lot-to-lot mechanical property variation, the alloy was deemed not suitable for a structural application. Alternative uses of ATP-5 include tooling, molds, vacuum chucks, fixtures, and jigs, where the material can shine in non-load bearing applications with a need for dimensional control and machinability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Alternant Hydrocarbon Diradicals as Optically Addressable Molecular Qubits

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical–radical interactions in the ground state, generating π-systems with high diradical characters. We call them m-dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m-diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. Furthermore, the molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m-diradicals as molecular color centers.

Chemistry↗
DOE OSTI2020

Virtual Metering for Monitoring Building Energy Consumption

The United States Department of Energy (DOE) has standard metering requirements of commercial buildings for optimizing energy performance. The guiding principles are to continuously track and optimize energy performance and install building-level meters for electricity, natural gas, and steam. Some buildings at Los Alamos National Laboratory (LANL) have physical submeters monitoring their energy consumption, but these meters have proven to be unreliable. And, in most cases, replacing them has proven to be a slow process. Installing new submeters also requires a temporary lockout of the circuit on which they are being installed. Many buildings at LANL contain laboratories with ongoing experiments or data centers, which makes an equipment power outage nearly impossible to plan. This inability to plan power outages results in long-term submeter failures. Although most submeters are eventually replaced, failures lead to missing consumption data for some unpredictable, extended time. A building automation system (BAS) is a system that provides control and monitoring on a building to maintain the operational performance of the building and occupancy comfort. Many buildings at LANL currently have a BAS, and all new renovations and installs will include installing a BAS if one does not already exist. The intended purpose for a BAS is primarily to monitor the health and efficiency of a building; however, it is also possible to calculate equipment power and energy consumption using BAS information. This project aims to use virtual meters to monitor building energy consumption as a cost-effective and minimally labor-intensive alternative to installing physical submeters. The fault detection and diagnostics tool, SkySpark, provides a centralized database for all the data from the various BAS that are active at LANL. This data includes the information that is needed to create virtual meters for heating, ventilating, and air conditioning (HVAC) systems in most buildings, including heating and cooling loads. 9 This report begins with a detailed summary of the project, including the reasoning, procedure, and results. The specific processes of creating the various virtual meters are then identified. Then the limitations are discussed. And, lastly, the results and future potential are presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

A research program to measure the lifetime of spin polarized fuel

The use of spin polarized fuel could increase the deuterium-tritium (D-T) fusion cross section by a factor of 1.5 and, owing to alpha heating, increase the fusion power by an even larger factor. Issues associated with the use of polarized fuel in a reactor are identified. Theoretically, nuclei remain polarized in a hot fusion plasma. The similarity between the Lorentz force law and the Bloch equations suggests polarization can be preserved despite the rich electromagnetic spectrum present in a magnetic fusion device. The most important depolarization mechanisms can be tested in existing devices. The use of polarized deuterium and 3 He in an experiment avoids the complexities of handling tritium, while encompassing the same nuclear reaction spin-physics, making it a useful proxy to study issues associated with full D-T implementation. 3 He fuel with 65% polarization can be prepared by permeating optically-pumped 3 He into a shell pellet. Dynamically polarized 7 Li-D pellets can achieve 70% vector polarization for the deuterium. Cryogenically-frozen pellets can be injected into fusion facilities by special injectors that minimize depolarizing field gradients. Alternatively, polarized nuclei could be injected as a neutral beam. Once injected, the lifetime of the polarized fuel is monitored through measurements of escaping charged fusion products. Multiple experimental scenarios to measure the polarization lifetime in the DIII-D tokamak and other magnetic-confinement facilities are discussed, followed by outstanding issues that warrant further study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗