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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Lyapunov-Based Generalized Dc-Side Controller Design for PV-Connected Systems

The objective of this paper is to realize a universal dc-side controller for photovoltaic (PV) systems where the control is agnostic to the downstream converter configuration. To achieve this, the downstream power converter and its controls are manipulated into an effective power control loop that is then cast into a generalized multi-loop design framework. On the dc side, a nonlinear small-signal model of the PV input is realized exclusively in terms of PV datasheet parameters (i.e., open-circuit voltage, short-circuit current, and maximum power point). Finally, a linear controller is used to modulate the dc-side PV system with the generic downstream power controller. A Lyapunov candidate is proposed to analyze the stability of the interconnected system and provide a streamlined approach for the controller design. The proposed design is validated on a 1k VA experimental setup that interfaces a PV module to the grid.

generalized PV controller↗
DOE OSTI2022

A Lyapunov-Based Generalized Dc-Side Controller Design for PV-Connected Systems: Preprint

The objective of this paper is to realize a universal dc-side controller for photovoltaic (PV) systems where the control is agnostic to the downstream converter configuration. To achieve this, the downstream power converter and its controls are manipulated into an effective power control loop that is then cast into a generalized multi-loop design framework. On the dc side, a nonlinear small-signal model of the PV input is realized exclusively in terms of PV datasheet parameters (i.e., open-circuit voltage, short-circuit current, and maximum power point). Finally, a linear controller is used to modulate the dc-side PV system with the generic downstream power controller. A Lyapunov candidate is proposed to analyze the stability of the interconnected system and provide a streamlined approach for the controller design. The proposed design is validated on a 1-kVA experimental setup that interfaces a PV module to the grid.

dc-dc control↗
DOE OSTI2022

Two-Stage Optimization Framework for Detecting and Correcting Parameter Cyber-Attacks in Power System State Estimation

One major tool of Energy Management Systems for monitoring the status of the power grid is State Estimation. Since the results of state estimation are used within the energy management system, the security of the state estimation process is most important. The focus research in this area is on detecting False Data Injection attacks on measurements. While this is important, State Estimation also rely on database that are used to describe the relationship between measurements and systems' states. This paper presents a two-stage programming framework to detect and correct attacks in the parameters of the measurement model used by the state estimation process in the Energy Management System. In the first stage, an estimate of the line parameters ratios are obtained. In the second stage, the estimated ratios from stage I are used in a Bi-Level model for obtaining a final estimate of the measurements' model parameters. Hence, the presented framework does not only unify the detection and correction in a single optimization run, but also provide a monitoring scheme for the SE database that is typically considered static. In addition, in the two stages, linear programming framework is preserved. For validation, the IEEE 118 bus system is used for implementation. The results of this paper illustrate the effectiveness of the proposed model for detecting attacks in the database used in the state estimation process.

state estimation, two-stage optimization, cyber-ph↗
DOE OSTI2024

Dynamical instability in multi-orbiter systems with gas friction

ABSTRACT Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

Li, Jiaru (ORCID:0000000155507421)↗
DOE OSTI2026

Krylov winding and emergent coherence in operator growth dynamics

The operator wavefunction provides a fine-grained description of quantum chaos and of the irreversible growth of simple operators into increasingly complex ones. Remarkably, at finite temperature this wavefunction can acquire a phase that increases linearly with the operator’s size, a phenomenon called . Although size winding occurs naturally in a holographic setting, the emergence of a coherent phase in a scrambled operator remains mysterious from the standpoint of a thermalizing quantum many-body system. Here, in this article, we elucidate this phenomenon by introducing the related concept of , whereby the operator wavefunction acquires a phase which winds linearly with the Krylov index. We show that Krylov winding is a generic feature of quantum chaotic systems and is a direct consequence of the universal operator growth bound hypothesis. It gives rise to size winding under two additional conditions: (i) a low-rank mapping between the Krylov and size bases, which ensures phase alignment among operators of the same size, and (ii) the saturation of the "chaos-operator growth" bound 𝜆 𝐿 ≤ 2⁢𝛼 (with 𝜆 𝐿 the Lyapunov exponent and 𝛼 the growth rate), which ensures a linear phase dependence on size. For systems which do not saturate this bound, with ℎ = 𝜆 𝐿 /2⁢𝛼 < 1, the winding with Pauli size ℓ becomes superliner, behaving as ℓ 1/ℎ . We illustrate these results with two classes of microscopic models: the Sachdev-Ye-Kitaev (SYK) model and its variants, and a disordered 𝑘-local spin model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Optimal Power System Black start using Inverter-Based Generation

Power system black start readiness is part of the system planning. Utility planners perform periodic studies to assess if their power system is capable of total restoration following a black out. Hydro and diesel generators are the most commonly used black start capable resources by power utilities. However, with increasing penetration of solar generation, inverter-based resources can be considered to provide black start capability. Since the solar inverters can be located at multiple locations throughout the power system, and in view of their unique characteristics, an optimal real-time capable plan is helpful for system operators for faster black start. Black start optimization is a multi stage mixed-integer non-linear optimization which is extremely hard to solve. In this paper, we propose an optimal black start methodology that is easier to solve and scalable in real-time. We demonstrate the proposed methodology on two test systems and illustrate how inverter-based resources can contribute and improve power system restoration.

power system restoration, blackstart, Inverter-bas↗
DOE OSTI2022

Time dependent supervisory control update with FARM using rolling window

This report describes improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution to the dispatch problem involves economically optimal dispatches that satisfy limits on production variables and their rates of variation (explicit constraints) as well as process variables tied to the service life of equipment (implicit constraints). FARM serves to validate or confirm that a HERON solution for explicit constraints also satisfies the implicit constraints. FARM-alpha was released by Argonne National Laboratory in January 2021 followed by FARM-Beta in January 2022 with the latter providing increased flexibility for the user. In this report, FARM-Gamma, the latest version of the code, is described. The major improvement is the implementation of a system identification algorithm based on the Dynamic Mode Decomposition with Control (DMDc) coupled with a “Rolling Window” scheme that allows obtaining linear time-varying state-space models. This feature equips FARM with the most accurate approximation of system dynamics, and it relieves the user from the burden of performing an exhaustive off-line characterization of the dynamics. FARM-Gamma capabilities are assessed by solving the power dispatch problem for a representative IES unit. The simulation times corresponding to the different releases are estimated and compared. These values capture the increasing computational burden of the successively higher fidelity state-space models adopted by FARM-Alpha, FARM-Beta and FARM-Gamma. The code implementation provides significant flexibility, i.e., the user can always select the most suitable version of FARM according to the problem to be solved and the available computational resources. It is anticipated that FARM will play a role in addressing several future IES applications. We outline how it can support the coordinated management and safe operation of a nuclear plant coupled to industrial processes to produce hydrogen and synfuels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Tunable valleytronics with symmetry-retaining high polarization degree in SnS x Se 1–x model system

SnS has recently been shown to possess unique valleytronic capability with a large polarization degree, where non-degenerate valleys can be accessed using linearly polarized light, bestowed upon by the unique anisotropy and wavefunction symmetry. It is thus of utmost importance to demonstrate the extension of such effects for the IV–VI system in general, thereby elucidating the generality and tunability of such valleytronics. In this study, we show the highly tunable valleytronics via gradual compositional control of the tin(II) sulfo-selenide (SnS x Se 1–x ) alloy system with excellent retainment of symmetry-determined selection rules. We show the presence of both ΓY and ΓX valleys in all alloy compositions via selectivity in absorption and emission of linearly polarized light by optical reflection (R)/transmission (T) and photoluminescence measurements and tuned the bandgaps of the valleys within a range of 1.28 eV–1.05 eV and 1.48 eV–1.24 eV, respectively. This simultaneous tuning of non-degenerate valleys agrees well with theoretical calculations. We then fitted the bandgap values in compositional space, obtaining bowing parameters as a useful database. We further demonstrated the feasibility of using IV–VI valleytronics systems in general by elucidating the retainment of strong polarization degrees of as high as 91% across all compositions. The generalization of such purely symmetry-dependent valleytronics also opens up opportunities for the discovery of more multi-functional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Geometric Interpretation of a Non-Linear Extension of Quantum Mechanics

We recently introduced a particular non-linear generalization of quantum mechanics that has the property that it is exactly solvable in terms of the eigenvalues and eigenfunctions of the Hamiltonian of the usual linear quantum mechanics problem. In this paper, we suggest that the two components of the wave function represent the system described by the Hamiltonian H in two different asymptotic regions of spacetime and we show that the non-linear terms can be viewed as giving rise to gravitational effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Applying Quantum Computing to Simulate Power System Dynamics

Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. This paper demonstrates the potential use of quantum computing algorithms to model the power system dynamics. Leveraging a symbolic programming framework, we equivalently convert the power system dynamics’ differential algebraic equations (DAEs) into ordinary differential equations (ODEs), where the data of the state vector can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion, the quantum state tensor, and Hamiltonian simulation, whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can simulate the dynamics of the power system with high accuracy, whereas its complexity is polynomial in the logarithm of the system dimension. Our work also illustrates the use of scientific machine learning tools for implementing scientific computing concepts, e.g., Taylor expansion, DAEs/ODEs transform, and quantum computing solver, in the field of power engineering.

Tran, Huynh↗
DOE OSTI2023

Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction

Understanding dynamics in complex systems is challenging because there are many degrees of freedom, and those that are most important for describing events of interest are often not obvious. The leading eigenfunctions of the transition operator are useful for visualization, and they can provide an efficient basis for computing statistics, such as the likelihood and average time of events (predictions). Here, we develop inexact iterative linear algebra methods for computing these eigenfunctions (spectral estimation) and making predictions from a dataset of short trajectories sampled at finite intervals. We demonstrate the methods on a low-dimensional model that facilitates visualization and a high-dimensional model of a biomolecular system. Implications for the prediction problem in reinforcement learning are discussed.

Chemistry↗
DOE OSTI2023

Design Considerations of a Coordinative Demand Charge Mitigation Strategy

This paper presents a coordinative demand charge mitigation (DCM) strategy for reducing electricity consumption during system peak periods. Available DCM resources include batteries, diesel generators, controllable loads, and conservation voltage reduction. All resources are directly controlled by load serving entities. A mixed integer linear programming based energy management algorithm is developed to optimally coordinate of DCM resources considering the load payback effect. To better capture system peak periods, two different kinds of load forecast are used: the day-ahead load forecast and the peak-hour probability forecast. Five DCM strategies are compared for reconciling the discrepancy between the two forecasting results. The DCM strategies are tested using actual utility data. Simulation results show that the proposed algorithm can effectively mitigate the demand charge while preventing the system peak from being shifted to the payback hours. We also identify the diminishing return effect, which can help load serving entities optimize the size of their DCM resources.

Hu, Rongxing↗
DOE OSTI2024

Role of Bottlebrush Additives on the Structure of Block Copolymers in the Bulk and Thin Films

Blending block copolymers (BCP) with additives is a useful approach for controlling BCP morphology and properties. In athermal systems, blends of BCPs with polymer additives having very high molecular (M n ) mass generally result in macrophase separation. Bottlebrush polymers, which consist of a linear backbone and grafted side chains, present an interesting alternative where the overall M n of the system can be very large but the low M n side chains may drive miscibility with the BCP. Here, in this study, a bottlebrush with a polynorbornene backbone and polystyrene (PS) side chains is blended with PS-b-poly(methyl methacrylate) (PS-b-PMMA) of varying M n , and the resulting morphologies are examined in both the bulk and thin films. Two different M n of PS-b-PMMA were used in the bulk study, and the analysis of small-angle X-ray scattering data shows that the blends were miscible and lamellar at all concentrations. This deviates from reference series of both low and high M n linear polymer additives, which either showed morphological transitions from lamellae to cylinders (low M n ) or were immiscible at all mass fractions studied (high M n ). The relative molecular mass of the side chain (N SC ) and the corresponding component in the BCP (N A ) dictate the distribution of the bottlebrush throughout the BCP, analogous to BCP/linear blends or grafted nanoparticles in a homopolymer matrix. The studies on thin films show a thickness dependence for bottlebrush mass fractions at or above 0.17, a behavior which may be driven by conformational changes of the bottlebrush upon confinement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of a new control rod drive mechanism design for the ISU AGN-201M reactor

The Aerojet General Nucleonics (AGN) model 201-Modified, known as the AGN-201M reactor, plays an essential role in the educational and research activities at Idaho State University (ISU). The ISU AGN-201M's original Control Rod Drive Mechanism (CRDM) has been in operation for more than fifty years with no large-scale redesigns. The CRDM is required to eject the fuel rods within one second during a SCRAM event (also known as a 'reactor trip') and adjust the control rods' insertion speed, and keeps the rod insertion sequence correct. The existing control rod drive mechanisms meet these criteria but experience a few concerns due to aging. Concerns include complex maintenance and costly repairs for old electromechanical components, rod position feedback errors, and the impediment of the plate during a SCRAM due to binding of the lead screws of the existing mechanism. During a binding event, the drive mechanism becomes locked, preventing the control rod's magnetic plate from moving in or out under the reactor's normal and emergency operating conditions. Although, the binding has no effect on the ability of the rod to exit the core during the SCRAM. To counteract the concerns and issues with the current CRDM, a new design has been proposed using newer components and a simplified design. The new design utilizes more advanced electric and mechanical components that are commercially available. The new CRDM system is divided into four main aspects: (1) control rod movement design (motor, lead screw, guide rods), (2) control rod ejection (springs, electromagnet), (3) control rod position and feedback (position transducer, microswitches), and (4) material selection and structural analysis. The new design aims to reduce the overall complexity and probability of failure to improve the reactor's overall reliability. With proper material selection and improved structural design, the new drives are lighter with little to no change in structural integrity. The new control rod drive mechanism eliminates binding scenarios by using a single lead screw and implementing additional guide rods. An advanced linear position sensor and microswitches replace the existing and aging synchro system for accurate rod position feedback resulting in better reactivity control. The new design meets the reactor's operational limits by having an average reactivity insertion of 0.065% Δk/k per second, which corresponds to a total control rod insertion time of 19.23 s, while the control rod's ejection time remains less than one second during a SCRAM event. The new design ensures the reactor's long-term viability for educational and research activities by increasing the reliability and safety of operation for years to come.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Parallel-in-Time Integration for Nonlinear Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to diffusion-dominated PDEs, with some recent success in scalar linear hyperbolic PDEs. Given the ubiquity of numerical simulations of hyperbolic PDEs throughout the sciences, in particular, nonlinear hyperbolic systems, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems beyond simple scalar and linear cases, which is the main topic of this subcontract. The main focus of the work was to further develop and perfect coarse-grid operators for the Multigrid Reduction-inTime (MGRIT) method applied to hyperbolic PDEs that were recently proposed in PhD thesis, based on a modified semi-Lagrangian approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

ToPolyAgent: AI agents for coarse-grained bead-spring topological polymer simulations

We introduce ToPolyAgent, a multi-agent AI framework for performing coarse-grained molecular dynamics (MD) simulations of topological polymers through natural language instructions. By integrating large language models (LLMs) with domain-specific computational tools, ToPolyAgent supports both interactive and autonomous simulation workflows across diverse polymer architectures, including linear, ring, brush, and star polymers, as well as dendrimers. The system consists of four LLM-powered agents: a Config Agent for generating initial polymer–solvent configurations, a Simulation Agent for executing LAMMPS-based MD simulations and conformational analyses, a Report Agent for compiling markdown reports, and a Workflow Agent for streamlined autonomous operations. Interactive mode incorporates user feedback loops for iterative refinements, while autonomous mode enables end-to-end task execution from detailed prompts. We demonstrate ToPolyAgent's versatility through case studies involving diverse polymer architectures under varying solvent conditions, thermostats, and simulation lengths. Furthermore, we highlight its potential as a research assistant by directing it to investigate the effect of interaction parameters on the linear polymer conformation, and the influence of grafting density on the persistence length of the brush polymer. By coupling natural language interfaces with rigorous simulation tools, ToPolyAgent lowers barriers to complex computational workflows and advances AI-driven materials discovery in polymer science. It lays the foundation for autonomous and extensible multi-agent scientific research ecosystems.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2022

Probing ultrafast excited-state dynamics using EUV-IR six-wave-mixing emission spectroscopy

Non-linear processes such as four-wave-mixing have become instrumental in attosecond EUV spectroscopy. Using EUV high harmonics in conjunction with collinear near-infrared and mid-infrared fields, we extended this approach to high-order-mixing between three colors. Specifically, we find that atomic resonances in neon exhibit a significant cross section for six-wave mixing. The MIR probe frequency tunability in our multicolor scheme is employed to access several optically dark resonances and probe the quantum beat of a coherent electronic wavepacket using background-free EUV emission as a diagnostic. This technique can be easily extended to other atomic and molecular systems, and opens the door to multi-dimensional non-linear spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗