Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear systems”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Real Time Applications Using Linear State Estimation Technology (RTA/LSE)

Electric Power Group, LLC (EPG) was awarded DOE OE0000849 to design, develop, and demonstrate three real-time applications for monitoring power system stability. The applications use phasor measurements and processing by a linear state estimator (LSE) to improve system visibility and accuracy. These applications leverage the large and expanding phasor measurement systems being deployed by utilities throughout the world, and specifically in the United States. Phasor measurements provide a much more detailed view of the power systems than traditional SCADA systems, so can provide better visibility of events and an early warning of developing problems. These three applications are real-time contingency analysis (RTCA), area angle monitoring (AAM), and a voltage stability index (VSI). The project was fully completed as proposed. All three applications were successfully implemented and demonstrated. There were challenges the caused project delays, however. While the RTCA is based on well-known concepts and readily available tools, its deployment using phasor measurements meant it needed to be adapted to sometimes sparse measurements and limited coverage areas. EPG met these challenges with extensive modeling and testing. The AAM is based on unproven research, so EPG had to deal with a lot of unknowns and develop additional methods for practical implementation. These activities took more time than anticipated, so required extra time for development. In addition, there was difficulty getting adequate measurement locations for successful deployment and utility installation delays. Consequently, the project required a time extension to complete. These three real-time applications, RTCA, AAM, and VSI, use only phasor measurement input. Since they do not require the EMS or other measurement systems, they can serve as a backup and security check for the EMS. Phasor measurements are validated and extended by an LSE. The LSE will always produce an output even when some inputs are impaired, so the RT applications will operate even when the EMS state estimator does not converge and the applications it supports cannot function. Consequently, these RT applications can provide valuable additional reliability for control center operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A Relaxed PV Bus Model in Linear Power Flow

The emerging multimodal active energy resource units make it necessary to introduce constant voltage amplitude buses (PV bus) in the operation of modern distribution systems. To enable the analysis of PV bus using linear power flow models, this letter proposes a relaxed PV bus model for solving power flow related problems using an indirect modeling approach. The developed PV bus model can be applied to linear power flow models in both rectangular and polar coordinates. Besides, the approximation error of the relaxed PV bus model can also be controlled and quantified. Lastly, the proposed relaxed PV bus model is validated using different test systems. According to the numerical studies, the maximum error on the PV bus approximation in all case studies is only 0.61%, and the maximum error of the supporting power at the PV buses is 6.5%.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Quantum three-wave instability

For the three-wave interaction, the lowest-order nonlinear interaction in plasma dynamics, we describe how the quantum system, which is time-independent, nonchaotic, finite-dimensional, and Hermitian, can give rise to a linear instability corresponding to that in the classical system. Furthermore, we show that the instability is realized in the quantum regime as a cascade of the wave function in the space of occupation number states, and the unstable quantum system has a richer spectrum and a much longer recurrence time than the stable quantum system. The conditions for instability of the quantum three-wave interaction are described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗
DOE OSTI2023

Open-Loop Control Co-Design of Semisubmersible Floating Offshore Wind Turbines Using Linear Parameter-Varying Models

Here, this paper discusses a framework to design elements of the plant and control systems for floating offshore wind turbines in an integrated manner using linear parameter-varying models. Multiple linearized models derived from aero-elastic simulation software in different operating regions characterized by the incoming wind speed are combined to construct an approximate low-fidelity model of the system. The combined model is then used to generate open-loop, optimal control trajectories as part of a nested control co-design strategy that explores the system’s power production and stability using the platform pitch tilt as a proxy in the context of crucial plant and control design decisions. The radial distance between the central and outer columns and the diameter of the outer columns of the semisubmersible platform are the plant design variables. The platform stability and power production are studied for different plant design decisions. The effect of plant decisions on subsequent power production and stability response of the floating wind turbine is quantified in terms of the levelized cost of energy. The results show that the inner-loop constraints and the plant design decisions affect the turbine’s power and, subsequently, the cost of the system.

17 WIND ENERGY↗
DOE OSTI2021

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

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

Protonic nickelate device networks for spatiotemporal neuromorphic computing

Computation in biological neural circuits arises from the interplay of nonlinear temporal responses and spatially distributed dynamic network interactions. Replicating this richness in hardware has remained challenging, as most neuromorphic devices emulate only isolated neuron- or synapse-like functions. Here we introduce an integrated neuromorphic computing platform in which both nonlinear spatiotemporal processing and programmable memory are realized within a single perovskite nickelate material system. By engineering symmetric and asymmetric hydrogenated NdNiO 3 junction devices on the same wafer, we combine ultrafast, proton-mediated transient dynamics with stable multilevel resistance states. Networks of symmetric NdNiO 3 junctions exhibit emergent spatial interactions mediated by proton redistribution, while each node simultaneously provides short-term temporal memory, enabling nanosecond-scale operation with an energy cost of ~0.2 nJ per input. When interfaced with asymmetric output units serving as reconfigurable long-term weights, these networks allow both feature transformation and linear classification in the same material system. Leveraging these emergent interactions, the platform enables real-time pattern recognition and achieves high accuracy in spoken digit classification and early seizure detection, outperforming temporal-only or uncoupled architectures. These results position protonic nickelates as a compact, energy-efficient, CMOS-compatible platform that integrates processing and memory for scalable intelligent hardware.

Electrical and electronic engineering↗
DOE OSTI2021

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Universal framework for simultaneous tomography of quantum states and SPAM noise

We present a general denoising algorithm for performing simultaneous tomography of quantum states and measurement noise. This algorithm allows us to fully characterize state preparation and measurement (SPAM) errors present in any quantum system. Our method is based on the analysis of the properties of the linear operator space induced by unitary operations. Given any quantum system with a noisy measurement apparatus, our method can output the quantum state and the noise matrix of the detector up to a single gauge degree of freedom. We show that this gauge freedom is unavoidable in the general case, but this degeneracy can be generally broken using prior knowledge on the state or noise properties, thus fixing the gauge for several types of state-noise combinations with no assumptions about noise strength. Such combinations include pure quantum states with arbitrarily correlated errors, and arbitrary states with block independent errors. This framework can further use available prior information about the setting to systematically reduce the number of observations and measurements required for state and noise detection. Our method effectively generalizes existing approaches to the problem, and includes as special cases common settings considered in the literature requiring an uncorrelated or invertible noise matrix, or specific probe states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Systems Analysis Approach to Polyethylene Terephthalate and Olefin Plastics Supply Chains in the Circular Economy: A Review of Data Sets and Models

The environmental and economic impacts of implementing a circular economy in plastic waste supply chains are not well understood. The proposed systems analysis framework assesses environmental, social, and economic impacts of plastic waste supply chains in a circular economy. The first objective of this article is to identify datasets, models and knowledge gaps associated with waste plastic supply chain processes, mainly in the U.S. Our literature review indicated that the best datasets exist for virgin plastic resin production, mechanical recycling, landfilling, and incineration, with materials recovery facility being intermediate, and with chemical recycling the lowest. The second objective of this article is to develop an illustrative application of the framework by conducting a preliminary systems analysis of PET bottles with closed-loop recycling. Here, the preliminary systems analysis of PET bottles utilized a linear programming optimization method. Our optimization model indicated that both chemical and mechanical recycling processes are needed to achieve a true circular economy of PET bottles with the least greenhouse gas emissions, specifically reductions of 24% when compared with the linear economy. Good quality and standardized life cycle assessment and techno-economic analysis studies are needed to better understand the environmental, economic, and social impacts of advanced sorting and chemical recycling technologies.

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

Strange Metal in Magic-Angle Graphene with near Planckian Dissipation

Recent experiments on magic-angle twisted bilayer graphene have discovered correlated insulating behavior and superconductivity at a fractional filling of an isolated narrow band. Here we show that magic-angle bilayer graphene exhibits another hallmark of strongly correlated systems—a broad regime of $\textit{T{$-linear resistivity above a small density-dependent crossover temperature—for a range of fillings near the correlated insulator. This behavior is reminiscent of similar behavior in other strongly correlated systems, often denoted “strange metals,” such as cuprates, iron pnictides, ruthenates, and cobaltates, where the observations are at odds with expectations in a weakly interacting Fermi liquid. We also extract a transport “scattering rate,” which satisfies a near Planckian form that is universally related to the ratio of $(k_{\text{B}}T/ℏ)$. In this work, our results establish magic-angle bilayer graphene as a highly tunable platform to investigate strange metal behavior, which could shed light on this mysterious ubiquitous phase of correlated matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

COZMIC. III. Cosmological Zoom-in Simulations of Self-interacting Dark Matter with Suppressed Initial Conditions

We present eight cosmological dark matter (DM)-only zoom-in simulations of a Milky Way–like system that include suppression of the linear matter power spectrum P(k), and/or velocity-dependent DM self-interactions, as the third installment of the COZMIC suite. We consider a model featuring a massive dark photon that mediates DM self-interactions and decays into massless dark fermions. The dark photon and dark fermions suppress linear matter perturbations, resulting in dark acoustic oscillations in P(k), which ultimately affect dwarf galaxy scales. The model also features a velocity-dependent elastic self-interaction between DM particles (SIDM), with a cross section that can alleviate small-scale structure anomalies. For the first time, our simulations test the impact of P(k) suppression on gravothermal evolution in an SIDM scenario that leads to core collapse in (sub)halos with present-day virial masses below ≈10 9 M ⊙ . In simulations with P(k) suppression and self-interactions, the lack of low-mass (sub)halos and the delayed growth of structure reduce the fraction of core-collapsed systems relative to SIDM simulations without P(k) suppression. In particular, P(k) suppression that saturates current warm DM constraints almost entirely erases core collapse in isolated halos. Models with less extreme P(k) suppression produce core collapse in ≈20% of subhalos and ≈5% of isolated halos above 10 8 M ⊙ , and also increase the abundance of extremely low-concentration isolated low-mass halos relative to SIDM. These results reveal a complex interplay between early and late-Universe DM physics, revealing new discovery scenarios in the context of upcoming small-scale structure measurements.

dark matter↗
DOE OSTI2022

Rapid Load Transition for Integrated Solid Oxide Fuel Cell - Gas Turbine (SOFC-GT) Energy Systems: A Demonstration of the Potential for Grid Response

Rapid load transition is an essential requirement for integrated energy systems to maintain grid resilience as more renewable resources are added to the grid. Integrated solid oxide fuel cell - gas turbine (SOFC-GT) systems can provide high efficiency and low carbon emissions over a broad range of turndown. These hybrids also have the potential to enable rapid grid response. The challenge has been to demonstrate effective control strategies to manage load transitions. In the present study, a load transition of ~50% was achieved in 10 s using a novel but simple strategy. Power demand on the SOFC and the GT were ramped down concurrently. During this transition, the SOFC anode fuel was manipulated to maintain SOFC fuel utilization while the cathode inlet air flow and temperature were also manipulated to thermally protect the SOFC. This study was conducted using the Hybrid Performance (Hyper) facility at the National Energy Technology Laboratory (NETL) in a co-simulation environment with the Idaho National Laboratory (INL)'s grid-simulation. The load ramping strategy was tested using a hardware-based cyber-physical simulation methodology. The results demonstrate a high-fidelity representation of SOFC-GT hybrid dynamics and validation of the control strategy. Thermal and electrochemical transients indicated that the SOFC was well protected during rapid load turndown without violating operability constraints. This demonstration revealed the non-linear nature of tightly coupled SOFC-GT system components, especially the non-linear response of SOFC cathode air flow and inlet temperature controls. These results highlight the needs and challenges in developing adaptive automatic controls for autonomous rapid load transitions. This work demonstrates that SOFC-GT hybrids are a viable option to provide the fast-ramping characteristics essential to accommodate high levels of variable renewable power while maintaining grid resilience, reliability, and environmental performance. The results also demonstrate the utility of co-simulation in advancing the tightly-coupled integrated energy systems needed to meet goals for zero-carbon power generation.

DIRECT ENERGY CONVERSION,POWER TRANSMISSION AND DI↗
DOE OSTI2021

Reservoir Computing as a Tool for Climate Predictability Studies

Reduced-order dynamical models play a central role in developing our understanding of predictability of climate irrespective of whether we are dealing with the actual climate system or surrogate climate models. In this context, the linear inverse modeling (LIM) approach, by capturing a few essential interactions between dynamical components of the full system, has proven valuable in providing insights into predictability of the full system. We demonstrate that reservoir computing (RC), a form of learning suitable for systems with chaotic dynamics, provides an alternative nonlinear approach that improves on the predictive skill of the LIM approach. We do this in the example setting of predicting sea surface temperature in the North Atlantic in the preindustrial control simulation of a popular earth system model, the Community Earth System Model so that we can compare the performance of the new RC-based approach with the traditional LIM approach both when learning data are plentiful and when such data are more limited. The improved predictive skill of the RC approach over a wide range of conditions—larger number of retained EOF coefficients, extending well into the limited data regime, etc.—suggests that this machine-learning technique may have a use in climate predictability studies. While the possibility of developing a climate emulator—the ability to continue the evolution of the system on the attractor long after failing to be able to track the reference trajectory—is demonstrated in the Lorenz-63 system, it is suggested that further development of the RC approach may permit such uses of the new approach in more realistic predictability studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
DOE OSTI2021

Tuning successive linear programming to solve AC optimal power flow problem for large networks

Successive linear programming (SLP) is a practical approach for solving large-scale nonlinear optimization problems. Alternating current optimal power flow (ACOPF) is no exception, particularly the large size of real-world networks. However, in order to achieve tractability, it is essential to tune the SLP algorithm presented in the literature. This paper presents a modified SLP algorithm to solve the ACOPF problem, specified by the U.S. Department of Energy’s (DOE) Grid Optimization (GO) Competition Challenge 1, within strict time limits. The algorithm first finds a near-optimal solution for the relaxed problem (i.e., Stage 1). Then, it finds a feasible solution in the proximity of the near-optimal solution (i.e., Stage 2 and Stage 3). The numerical experiments on test cases ranging from 500-bus to 30,000-bus systems show that the algorithm is tractable. Here the results show that our proposed algorithm is tractable and can solve more than 80% of test cases faster than the well-known Interior Point Method while significantly reduce the number of iterations required to solve ACOPF. The number of iterations is considered an important factor in the examination of tractability which can drastically reduce the computational time required within each iteration.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Dynamics and energy coupling of gas puff Z-pinches on a fast linear transformer driver

Gas puff Z-pinch experiments with annular Ar and Ne gas shells have been conducted on the Compact Experimental System for Z-pinch and Ablation Research (CESZAR) linear transformer driver (LTD) with 500 kA current and 160 ns rise time. Here, we present results from the first systematic gas puff Z-pinch experiments using a fast (≤ 200 ns) LTD as a driver, in which we show that 7% of the stored energy in the capacitors is coupled to plasma kinetic energy as estimated via self-emission and laser schlieren images. 0D and 1D simulations—which do not allow instability growth and thus reach greater maximum average velocities—using initial conditions inferred from experimental implosion trajectories predict coupling in excess of 10% of the stored energy. Additionally, the Ar and Ne implosions were comparably massed and thus achieved similar maximum kinetic energies, though the Ne pinches were more stable and the x-ray pulses were longer and produced higher yield: 2–5 ns and 0.21–0.52 J (0.15–0.37 J/cm) of Ar K-shell and 12–25 ns and 2.2–3.9 J (1.6–2.6 J/cm) of Ne K-shell, respectively. The difference in stability is most likely attributed to variations in initial conditions such as density distribution and gas breakdown initiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗