Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical 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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reducing Conservativeness of Polytopic Linear-Parameter-Varying Robust Vehicle Sideslip Angle Observer Through Minimum-Area Convex Quadrilateral Design

The polytopic linear-parameter-varying (LPV) method has become increasingly popular for designing intelligent control and estimation systems for ground vehicles, particularly for the vehicle sideslip angle observer. In vehicle lateral dynamics, the vehicle longitudinal velocity-induced nonlinearities are conventionally outer-approximated with polytopes such as rectangles or triangles to obtain a polytopic LPV system. Yet, such polytopic approximations tend to be conservative and may lead to inadequate observer performance. To address this issue, a minimum-area convex quadrilateral construction is proposed in this paper to reduce design conservatism. The suggested design is demonstrated through the synthesis of an LPV H ∞ robust vehicle sideslip angle observer. Furthermore, a dSPACE-ASM simulation study is conducted to demonstrate the effectiveness and the advantage of the proposed polytopic construction over a baseline approach.

Zhou, Xingyu↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through various methods, two of them being thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine thermal quadrant powers. The N-16 system utilizes a beta detector that outputs Nitrogen activation levels to calculate lobe power through an algorithm called Lobe Power Calculation and Indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. These correct indications are crucial for maintaining safety limits because operators rely on this information for decision making. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have increased the accuracy of power readings by making the system more dynamic. The new proposed coefficient and multiplier method implements a cycle specific and physics-based model intended for changing coefficients during operation. This calibration experiment focused on power splitting, outer shim and neck shim, as well as fuel burning into the reactor digital acquisition system (RDAS) weighting factors. Results demonstrated that the physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Quantum time dynamics employing the Yang-Baxter equation for circuit compression

Quantum time dynamics (QTD) is considered a promising problem for quantum supremacy on near-term quantum computers. However, QTD quantum circuits grow with increasing time simulations. This study focuses on simulating the time dynamics of one-dimensional (1D) integrable spin chains with nearest-neighbor interactions. We have proved the existence of a reflection symmetry in the quantum circuit employed for simulating the time evolution of certain classes of 1D Heisenberg model Hamiltonians by virtue of the quantum Yang-Baxter equation, and how this symmetry can be exploited to compress and produce a shallow quantum circuit. With this compression scheme, the depth of the quantum circuit becomes independent of step size and only depends on the number of spins. We show that the depth of the compressed circuit is rigorously a linear function of the system size for the studied Heisenberg model Hamiltonians in the present work. As a consequence, the number of CNOT gates in the compressed circuit only scales quadratically with the system size, which allows for the simulations of time dynamics of very large 1D spin chains. We derive the compressed circuit representations for different special cases of the Heisenberg Hamiltonian. We compare and demonstrate the effectiveness of this approach by performing simulations on quantum computers.

1-dimensional spin chains↗
DOE OSTI2022

Optimal Sizing and Operation of a Hybrid Clean Energy Center

In this paper, we consider the deployment of battery energy storage systems (BESSs) paired with photovoltaics (PV) to form a hybrid clean energy center (HCEC), where the BESS is equipped for two main purposes: i) reducing the volatility of PV output and ii) shifting energy to meet an export limit. We develop a linear dynamical model to characterize the output of an HCEC and express the volatility and export power as functions of the battery's dispatch operations. We also propose an optimal sizing method via second-order cone programming (SOCP) to maximize the net benefits of the HCEC. Numerical studies are performed to illustrate and validate the proposed method.

Battery energy storage systems, photovoltaics, sec↗
DOE OSTI2023

Synthetically non-Hermitian nonlinear wave-like behavior in a topological mechanical metamaterial

Topological mechanical metamaterials have enabled new ways to control stress and deformation propagation. Exemplified by Maxwell lattices, they have been studied extensively using a linearized formalism. Herein, we study a two-dimensional topological Maxwell lattice by exploring its large deformation quasi-static response using geometric numerical simulations and experiments. We observe spatial nonlinear wave-like phenomena such as harmonic generation, localized domain switching, amplification-enhanced frequency conversion, and solitary waves. We further map our linearized, homogenized system to a non-Hermitian, nonreciprocal, one-dimensional wave equation, revealing an equivalence between the deformation fields of two-dimensional topological Maxwell lattices and nonlinear dynamical phenomena in one-dimensional active systems. Our study opens a regime for topological mechanical metamaterials and expands their application potential in areas including adaptive and smart materials and mechanical logic, wherein concepts from nonlinear dynamics may be used to create intricate, tailored spatial deformation and stress fields greatly transcending conventional elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Networked Microgrid Topology Reconfiguration to Promote Fairness in Proactive Load Shedding

Increasing occurrences of natural disasters and grid emergency events consistently challenge the safe and reliable operations of power systems. During such emergency situations, system operators may proactively shed load to mitigate risks. However, uncoordinated implementation of load shedding may disrupt electricity supply and even lead to cascading failures. Meanwhile, it is crucial to address potential biases affecting different customers when executing load shedding. This paper addresses the dynamic topology reconfiguration problem for networked microgrids with distributed energy resources under emergency conditions. Specifically, we propose a novel rolling-horizon optimization model that integrates fairness-aware constraints into the networked microgrid topology reconfiguration. Unlike existing approaches that focus solely on efficiency or apply fairness considerations in static settings, our method explicitly incorporates temporal fairness constraints to restrict repeated or excessive load curtailment for load blocks. Moreover, the fairness-aware constraints are specifically developed for the context of dynamic networked microgrid topology reconfiguration, and are designed to be convex or amenable to linear reformulations, which offers a more tractable alternative to traditional models with non-convex formulations. Numerical studies on a modified IEEE 13-bus system and a larger-sized SMART-DS networked microgrid system demonstrate the performance of the proposed algorithm towards more fairness-aware networked microgrid topology reconfiguration decision-making.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

How the dynamical properties of globular clusters impact their $γ$-ray and X-ray emission

The X-ray and $γ$-ray emission of globular clusters (GCs) is attributed to their large fraction of compact binary systems, especially those with millisecond pulsars (MSPs). We analyse a population of 124 Galactic GCs to investigate how their dynamical properties affect the formation and evolution of compact binary systems and how this can be translated into the clusters’ observed X-ray and γ-ray emission. We use mainly Chandra X-ray Observatory and Fermi Large Area Telescope observations to achieve our goals and start by detecting 39 GCs in γ rays, seven of which are not listed in previous Fermi -LAT catalogs. Additionally, we find that the total number of X-ray sources within a GC and its γ-ray luminosity are linearly correlated with the stellar encounter rate, indicating that compact binary systems are mainly formed via close stellar encounters. We also find an unexpected rise in the number of X-ray sources for GCs with low rates of stellar encounters, suggesting that there is a dynamical threshold where the formation of X-ray sources is dominated by stellar encounters. Furthermore, we use the Heggie-Hills law to find that subsequent stellar encounters in these compact binaries will, on average, make the binaries even harder, with basically no possibility of binary ionization. Finally, we find that all GCs are point-like sources in γ rays, indicating that the MSPs are concentrated in the clusters’ cores, likely due to dynamical friction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

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

Multi-spatio-temporal scales PIV in a turbulent buoyant jet discharging in a linearly stratified environment

Time-resolved particle image velocimetry is implemented with a camera array and several laser sheets; this results in a multi-spatio-temporal scale system that is modular and flexible. The setup is optimized to capture global flow features, while locally resolving in space and time near the Kolmogorov scale. The apparatus relies extensively on machine vision CMOS cameras; they are inexpensive and stream to computer hard drive with virtually continuous time-resolved records (up to one hour for the current system). This leads to statistically converged data and also helps in minimizing down time between experimental runs. Details of the implementation and design of experiment are reported. The system instruments a vertical buoyant jet discharging in a linearly stratified environment. Refractive index matched aqueous solutions enable precise optical deployment. The density difference is 3% and the fluids have similar dynamic viscosity. The jet Reynolds number is 2.00 x 10 4 and is above the mixing transition. Such flows are typically challenging to instrument and few velocity field data are available to date. The velocity statistics are reported as well as first insights gained from the campaign.

42 ENGINEERING↗
DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗
DOE OSTI2025

Gapless dispersive continuum in a modulated quantum kagome antiferromagnet

The pursuit of quantum spin liquid (QSL) states in condensed matter physics has drawn attention to kagome antiferromagnets (AFM) where a two-dimensional corner-sharing network of triangles frustrates conventional magnetic orders. While quantum kagome AFMs based on Cu 2+ (3d 9 , s = 1/2) ions have been extensively studied, there is so far little work beyond copper-based systems. Here we present our bulk magnetization, specific heat and neutron scattering studies on single crystals of a new titanium fluorides Cs 8 RbK 3 Ti 12 F 48 where Ti 3+ (3d 1 , s = 1/2) ions form a modulated quantum kagome antiferromagnet that does not order magnetically down to 1.5 K. Our comprehensive map of the dynamic response function $S$($Q$,$\hslash$$\omega$) acquired at 1.5 K where the heat capacity is T-linear reveals a dispersive continuum emanating from soft lines that extend along (100). The data indicate fractionalized spinon-like excitations with quasi-one-dimensional dispersion within a quasi-two-dimensional spin system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Resilient Distribution System Restoration with Equitable Load Shedding

A methodology is proposed for the improvement of electric distribution system resilience to high-impact, low-probability catastrophic events. An approach for dynamic network reconfiguration and coordination of distributed energy resources is introduced to assist in restoration efforts. The problem is formulated as a mixed-integer linear program that minimizes generation costs, the cost of lost load, and costs associated with equitable load shedding, while respecting operational limits of generation, loads, and the network. Constraints are imposed on binary switching variables to ensure equitable load shedding in emergency situations. Numerical validation of the proposed approach is conducted on an example distribution feeder, and case studies are performed to analyze the impact of various parameters in the optimization problem formulation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Parameter Reduction of Composite Load Model Using Active Subspace Method

Over the past decades, the increasing penetration of distributed energy resources (DERs) has dramatically changed the power load composition in the distribution networks. The traditional static and dynamic load models can hardly capture the dynamic behavior of modern loads especially for fault-induced delayed voltage recovery (FIDVR) events. Thus, a more comprehensive composite load model with combination of static load, different types of induction motors, single-phase A/C motor, electronic load and DERs has been proposed by Western Electricity Coordinating Council (WECC). However, due to the large number of parameters and model complexity, the WECC composite load model (WECC CMLD) raises new challenges to power system studies. To overcome these challenges, in this paper, a cutting-edge parameter reduction (PR) approach for WECC CMLD based on active subspace method (ASM) is proposed. Firstly, the WECC CMLD is parameterized in a discrete-time manner for the application of the proposed method. Then, parameter sensitivities are calculated by discovering the active subspace, which is a lower-dimensional linear subspace of the parameter space of WECC CMLD in which the dynamic response is most sensitive. The interdependency among parameters can be taken into consideration by our approach. Finally, the numerical experiments validate the effectiveness and advantages of the proposed approach for WECC CMLD model.

active subspace↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonlinear Optics Measurements in IOTA

Nonlinear integrable optics is a recently proposed accelerator lattice design approach which allows to generate an amplitude dependent tune shift which is needed in high brightness accelerators to mitigate fast coherent instabilities. Whereas usually octupoles are used to achieve this task, this concept allows doing so without exciting any resonances, in turn preventing any particle loss. The concept is based around a special magnet design, together with specific constraints on the optics of the accelerator. To study such a system, the Integrable Optics Test Accelerator (IOTA) was recently constructed and commissioned at Fermilab. For the assessment of the performance of this concept, good knowledge of the optics and the (non-)linear dynamics without the special magnet is of key importance. As such, measurements were conducted in the IOTA ring, using the captured turn-by-turn data by the beam position monitors after excitation to infer quantities such as amplitude detuning and resonance driving terms. In this note, first results of these measurements are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗