Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution operators”

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 181 records · Page 10

A Scalable Method for Decarbonizing Modular Building Solutions: Preprint

The decarbonization movement emphasizes the shift in focus from energy efficiency to directly reducing global-warming impact. Blokable, LLC, a vertically integrated modular builder with an all-electric portfolio, worked with NREL on a roadmap to decarbonize its high-performance building product at a relative cost advantage by utilizing the learning curves of mass production. Previous decarbonization literature focused on either i) lifecycle assessments, or ii) efficiency measures. These decarbonization exercises were bespoke to individual building projects and did not consider positive feedback loops of builder experience or process repetition. Vertically integrated, prefab builders possess the unique ability to leverage learning and repetition to decarbonize their design-build-operate process. This collaboration between Blokable and NREL resulted in a decarbonization strategy utilizing the company's scaling and production efficiencies, on-site renewable energy and storage, and the projected evolution of building components over time based on trends and emerging legislation. The method developed here encompasses a growing business model, lifecycle carbon assessment, and projected changes in product and grid emissions over time due to existing trends and emerging legislation. This methodology incorporates learning-curve efficiencies gleaned from scaled manufacturing, as well as open-source tools integration for energy and carbon accounting. The output projects and compares cost and carbon savings per modular unit as production increases to 10,000 dwelling units annually over 15 years. The resulting roadmap illustrates a path to roughly 60% carbon savings and beyond-net-zero-energy performance at no incremental cost by 2030. The methodology can be mapped to other integrated or productized builders for methodical decarbonization.

affordable housing↗

A Scalable Method for Decarbonizing Modular Building Solutions

The decarbonization movement emphasizes the shift in focus from energy efficiency to directly reducing global-warming impact. Blokable, LLC, a vertically integrated modular builder with an all-electric portfolio, worked with NREL on a roadmap to decarbonize its high-performance building product at a relative cost advantage by utilizing the learning curves of mass production. Previous decarbonization literature focused on either i) lifecycle assessments, or ii) efficiency measures. These decarbonization exercises were bespoke to individual building projects and did not consider positive feedback loops of builder experience or process repetition. Vertically integrated, prefab builders possess the unique ability to leverage learning and repetition to decarbonize their design-build-operate process. This collaboration between Blokable and NREL resulted in a decarbonization strategy utilizing the company's scaling and production efficiencies, on-site renewable energy and storage, and the projected evolution of building components over time based on trends and emerging legislation. The method developed here encompasses a growing business model, lifecycle carbon assessment, and projected changes in product and grid emissions over time due to existing trends and emerging legislation. This methodology incorporates learning-curve efficiencies gleaned from scaled manufacturing, as well as open-source tools integration for energy and carbon accounting. The output projects and compares cost and carbon savings per modular unit as production increases to 10,000 dwelling units annually over 15 years. The resulting roadmap illustrates a path to roughly 60% carbon savings and beyond-net-zero-energy performance at no incremental cost by 2030. The methodology can be mapped to other integrated or productized builders for methodical decarbonization.

affordable housing↗

Hybrid quantum-classical approach for coupled-cluster Green's function theory

The three key elements of a quantum simulation are state preparation, time evolution, and measurement. While the complexity scaling of time evolution and measurements are well known, many state preparation methods are strongly system-dependent and require prior knowledge of the system's eigenvalue spectrum. Here, we report on a quantum-classical implementation of the coupled-cluster Green's function (CCGF) method, which replaces explicit ground state preparation with the task of applying unitary operators to a simple product state. While our approach is broadly applicable to many models, we demonstrate it here for the Anderson impurity model (AIM). The method requires a number of T gates that grows as O ( N 5 ) per time step to calculate the impurity Green's function in the time domain, where N is the total number of energy levels in the AIM. Since the number of T gates is analogous to the computational time complexity of a classical simulation, we achieve an order of magnitude improvement over a classical CCGF calculation of the same order, which requires O ( N 6 ) computational resources per time step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiobjective optimization of nuclear microreactor reactivity control system operation with swarm and evolutionary algorithms

To improve the marketability of novel microreactor designs, there is a need for automated and optimal control of these reactors. This paper presents a methodology for performing multiobjective optimization of control drum operation for a microreactor under normal and off-nominal conditions. Here, two different case studies are used where the control drum configuration is optimized for the reactor to be critical with some desired power distribution that would satisfy peaking limits. A surrogate model for power distribution is developed based on a feedforward neural network. The process for determining weights for scalarization of the multiobjective optimization problem is also detailed. Six optimization algorithms: evolutionary strategies, differential evolution, grey wolf optimization, Harris hawks optimization, moth flame optimization and particle swarm optimization, are all applied to these cases and the results analyzed. Although all these algorithms will demonstrate optima-seeking behavior, for real-time control it is necessary to identify the best algorithm to efficiently provide reasonable optima without operator interference. The moth flame optimization algorithm was found to perform particularly well on both cases. Overall, it was found that the algorithms capable of supplying the best optima were also the most consistent. Finally, the found optima were verified with the original model used to train surrogates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evolution of parton distribution functions in the short-distance factorization scheme

Lattice QCD offers the possibility of computing parton distributions from first principles, although not in the usual $\overline{MS}$ factorization scheme. Calculations are therefore matched to $\overline{MS}$ using a perturbative procedure which is the source of significant uncertainty within the currently accessible kinematics. We present the possibility of computing the z 2 evolution of non-singlet pseudo-parton distribution functions within the short factorization scheme in a numerically improvable way. The goal is to have tools to evolve a calculation to a scale where perturbative uncertainties are less pronounced. We compare a numerical extraction of the evolution operator from lattice data to the computation of z 2 dependence in perturbation theory. Finally, we discuss how this numerical work may be extended to address the two-scale problem that arises when the Ioffe time range must be made large to extend the reach of the calculation of the pseudo-PDF to smaller values of the momentum fraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗

Modeling gas release from a Bingham plastic slurry and deconvoluting measured data

Physics-based models for predicting the off-gassing characteristics of agitated slurries are presented. This approach decomposes the system into two separate but coupled slurry and headspace models. The physics driving bubble transport through the slurry and gas mixing within the headspace are discussed. An analytical expression for predicting the time evolution of the headspace concentration using first principles theory is also presented. Predictions from the numerical models, as well as expectations from the analytical solution, both agree with measured off gassing data for two different experimental operating conditions. After benchmarking the numerical model predictions against experimental data, the model was used to make predictions for gas release from a full-scale vessel and to perform sensitivity analyses to examine the sensitivity of gas release to parameters such as yield stress, consistency index, slurry density, bubble size, bubble concentration, and impeller speed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗

Control-Oriented Modeling of Cycle-to-Cycle Combustion Variability at the Misfire Limit in SI Engines

A control-oriented model is presented that can capture the prior-cycle correlation of combustion cycles during conditions with high levels of exhaust gas recirculation (EGR). Combustion events are modeled in discrete time and the dynamic evolution is captured by the residual air, fuel, and inert gas trapped in the combustion chamber. The mathematical formulation of the model is presented together with the calibration procedure to emulate a particular engine operating condition. A cycle-to-cycle system identification methodology is described which allows regressing model parameters from experimental data. Simulations are presented and compared to real engine measurements to show the modeling potential for analysis and control of combustion events.

Maldonado Puente, Bryan↗

Modes and Mechanisms of Pacific Decadal-Scale Variability

The modes of Pacific decadal-scale variability (PDV), traditionally defined as statistical patterns of variance, reflect to first order the ocean's integration (i.e., reddening) of atmospheric forcing that arises from both a shift and a change in strength of the climatological (time-mean) atmospheric circulation. While these patterns concisely describe PDV, they do not distinguish among the key dynamical processes driving the evolution of PDV anomalies, including atmospheric and ocean teleconnections and coupled feedbacks with similar spatial structures that operate on different timescales. In this review, we synthesize past analysis using an empirical dynamical model constructed from monthly ocean surface anomalies drawn from several reanalysis products, showing that the PDV modes of variance result from two fundamental low-frequency dynamical eigenmodes: the North Pacific–central Pacific (NP-CP) and Kuroshio–Oyashio Extension (KOE) modes. Both eigenmodes highlight how two-way tropical–extratropical teleconnection dynamics are the primary mechanisms energizing and synchronizing the basin-scale footprint of PDV. While the NP-CP mode captures interannual- to decadal-scale variability, the KOE mode is linked to the basin-scale expression of PDV on decadal to multidecadal timescales, including contributions from the South Pacific.

54 ENVIRONMENTAL SCIENCES↗

GSAS Tools

SAND2023-06684O GSAS Tools is a web application that manages user access to modeling and simulation tools and promotional material. This software, which is a spiking neural network (SNN) simulator, represents an SNN as a graph of stochastic differential equations and simulates the time-evolution of these equations. It has the capability of reading inputs from file and writing outputs to file, and generally supports experimentation, analysis, and algorithm development using SNNs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Noel, Todd↗

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗

Data-driven models of nonautonomous systems

Nonautonomous dynamical systems are characterized by time-dependent inputs, which complicates the discovery of predictive models describing the spatiotemporal evolution of the state variables of quantities of interest from their temporal snapshots. When dynamic mode decomposition (DMD) is used to infer a linear model, this difficulty manifests itself in the need to approximate the time-dependent Koopman operators. Our approach is to approximate the original nonautonomous system with a modified system derived via a local parameterization of the time-dependent inputs. The modified system comprises a sequence of local parametric systems, which are subsequently approximated by a parametric surrogate model using the DRIPS (dimension reduction and interpolation in parameter space) framework. The offline step of DRIPS relies on DMD to build a linear surrogate model, endowed with reduced-order bases for the observables mapped from training data. The online step interpolates on suitable manifolds to construct a sequence of iterative parametric surrogate models; the target/test parameter points on these manifolds are specified by a local parameterization of the test time-dependent inputs. Here, we use numerical experimentation to demonstrate the robustness of our method and compare its performance with that of deep neural networks.

97 MATHEMATICS AND COMPUTING↗

Quantum Tensor-Product Decomposition from Choi-State Tomography

The Schmidt decomposition is the go-to tool for measuring bipartite entanglement of pure quantum states. Similarly, it is possible to study the entangling features of a quantum operation using its operator-Schmidt or tensor-product decomposition. While quantum technological implementations of the former are thoroughly studied, entangling properties on the operator level are harder to extract in the quantum computational framework because of the exponential nature of sample complexity. Here, we present an algorithm for unbalanced partitions into a small subsystem and a large one (the environment) to compute the tensor-product decomposition of a unitary the effect of which on the small subsystem is captured in classical memory, while the effect on the environment is accessible as a quantum resource. This quantum algorithm may be used to make predictions about operator nonlocality and effective open quantum dynamics on a subsystem, as well as for finding low-rank approximations and low-depth compilations of quantum circuit unitaries. We demonstrate the method and its applications on a time-evolution unitary of an isotropic Heisenberg model in two dimensions. Published by the American Physical Society 2024

Mansuroglu, Refik (ORCID:000000017352513X)↗

ATR-SEIRAS Method to Measure Interfacial pH during Electrocatalytic Nitrate Reduction on Cu

This study reports the accuracy and applications of an attenuated total reflectance–surface-enhanced infrared absorption spectroscopy (ATR–SEIRAS) technique to indirectly measure the interfacial pH of the electrolyte within 10 nm of the electrocatalyst surface. This technique can be used in situ to study aqueous electrochemical reactions with a calibration range from pH 1–13, time resolution down to 4 s, and an average 95% confidence interval of 14% that varies depending on the pH region (acidic, neutral, or basic). The method is applied here to electrochemical nitrate reduction at a copper cathode to demonstrate its capabilities, but is broadly applicable to any aqueous electrochemical reaction (such as hydrogen evolution, carbon dioxide reduction, or oxygen evolution) and the electrocatalyst may be any SEIRAS-active thin film (e.g., silver, gold, or copper). The time-resolved results show a dramatic increase in the interfacial pH from pH 2–7 in the first minute of operation during both constant current and pulsed current experiments where the bulk pH is unchanged. Attempts to control the pH polarization at the surface by altering the electrochemical operating conditions—lowering the current or increasing the pulse frequency—showed no significant change, demonstrating the challenge of controlling the interfacial pH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A reinforcement learning approach to long-horizon operations, health, and maintenance supervisory control of advanced energy systems

In this work, we develop a Reinforcement Learning (RL) approach to the supervisory control problem for advanced energy systems, such as novel nuclear reactors and other demand-driven, mission-critical, and component-health-sensitive energy plants. The inclusive problem landscape considered captures the stochastic confluence of plant performance, component health evolution, power demand from the grid, diverse maintenance actions, and operator-defined goals and constraints, all considered over meaningfully long-enough reasoning horizons. Key aspects of the proposed approach are a receding horizon control-inspired technique dictating time- or event-triggered supervisory policy (re-)constructions, as well as additional capability-enabling contributions such as timescale compression, to handle long reasoning horizons and uncertainty in parts of the problem, and practical yet demonstrably-effective handling of hybrid action spaces with continuous and discrete decision variables. The resulting algorithm consists of a simulation-based RL agent constructing stochastic supervisory control policies over nontrivial action spaces and for long horizons, applying the learned policy to the system for a much shorter interval, and perpetually repeating, to construct the next long-horizon policy. That next policy will only be applied, again, for a short interval, yet originally far-in-time events move progressively closer, their associated uncertainty decreases, and new events and aspects enter the reasoning horizon. The proposed methodology bridges fundamental receding horizon concepts with the unequivocally stronger and more scalable reasoning of contemporary RL. Numerical examples using Soft Actor–Critic Deep RL illustrate the operation and efficacy of the proposed technique for a power plant tasked with health-aware load following missions in a dynamic electricity market landscape.

97 MATHEMATICS AND COMPUTING↗

Quantum Fokker-Planck modeling of degenerate electrons

In this work, an implicit and conservative numerical scheme is proposed for the isotropic quantum Fokker-Planck equation describing the evolution of degenerate electrons subject to elastic collisions with other electrons and ions. The electron-ion and electron-electron collision operators are discretized using a discontinuous Galerkin method, and the electron energy distribution is updated by an implicit time integration method. The numerical scheme is designed to satisfy all conservation laws exactly. Numerical tests and comparisons with other modeling approaches are shown to demonstrate the accuracy and conservation properties of the proposed method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗