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 1,189 records · Page 66

Topology of the $O(3)$ non-linear sigma model under the gradient flow

The $O(3)$ non-linear sigma model (NLSM) is a prototypical field theory for QCD and ferromagnetism, featuring topological qualities. Though the topological susceptibility should vanish in physical theories, lattice simulations of the NLSM find that the topological susceptibility diverges in the continuum limit. We introduce a $\theta$-term and analyze the topological charge as a function of $\theta$ under the gradient flow. Our results show that divergence persists in the presence of the flow, even at non-zero $\theta$.

Thomas, Stuart↗

Using a Grid-Forming Inverter to Stabilize a Low-Inertia Power System - Maui Hawaiian Island

As power systems around the world integrate greater amounts of wind and solar photovoltaic power, periods of very high instantaneous power shares of inverters, the primary interfacing technology for these generation sources, are complicating system stability and control. The contemporary, primary mode of inverter operation, grid-following, which explicitly assumes the presence of a local, stable voltage waveform, yields operational inadequacy at high instantaneous power shares potentially leading to instability due to the low-inertia conditions, as well as the correlated reduced voltage forming capacity on the respective system. Parallel connected grid-forming inverters, which directly regulate the local voltage, are a solution that is expected to bolster system stability and mitigate the shortcomings of the grid-following technology. In this paper, grid-following and two types of grid-forming inverter control, the traditional linear droop and the recently introduced nonlinear exponential droop (Droop-e), are simulated on a low-inertia, H = 0.48s, high inverter-based resource scenario, 97%, with a validated electromagentic transient domain model of the Hawaiian island of Maui power system. The benefit of a single grid-forming device over its grid-following counterpart is significant, both in terms frequency deviation and voltage stability. Further, the superiority of the Droop-e and the associated secondary power sharing control over both grid-following and linear droop grid-forming technologies is displayed, with improved nadir and rate of change of frequency over the linear droop control.

droop control↗

Toward witnessing molecular exciton entanglement from spectroscopy

Entanglement is a defining feature of quantum mechanics that can be a resource in engineered and natural systems, but measuring entanglement in experiment remains elusive especially for large chemical systems. Most practical approaches require determining and measuring a suitable entanglement witness which provides some level of information about the entanglement structure of the probed state. A fundamental quantity of quantum metrology is the quantum Fisher information (QFI), which is a rigorous witness of multipartite entanglement that can be evaluated from linear response functions for certain states. Here, in this paper, we explore measuring the QFI of molecular exciton states of the first-excitation subspace from spectroscopy. In particular, we utilize the fact that the linear response of a pure state subject to a weak electric field over all possible driving frequencies encodes the variance of the collective dipole moment in the probed state, which is a valid measure for QFI. The systems that are investigated include the molecular dimer, N-site linear aggregate with nearest-neighbor coupling, and N-site circular aggregate, all modeled as a collection of interacting qubits. Our theoretical analysis shows that the variance of the collective dipole moment in the brightest dipole-allowed eigenstate is the maximum QFI. The optical response of a thermally equilibrated state in the first-excitation subspace is also a valid QFI. Theoretical predictions of the measured QFI for realistic linear dye aggregates as a function of temperature and energetic disorder due to static variations of the host matrix show that two- to three-partite entanglement is realizable. This paper lays some groundwork and inspires measurement of multipartite entanglement of molecular excitons with ultrafast pump-probe experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum tensor network algorithms for evaluation of spectral functions on quantum computers

We investigate quantum algorithms derived from tensor networks to simulate the static and dynamic properties of quantum many-body systems. Using a sequentially prepared quantum circuit representation of a matrix product state (MPS) that we call a quantum tensor network (QTN), we demonstrate algorithms to prepare ground and excited states on a quantum computer and apply them to molecular nanomagnets (MNMs) as a paradigmatic example. In this setting, we develop two approaches for extracting the spectral correlation functions measured in neutron-scattering experiments: (a) a generalization of the SWAP test for computing wave function overlaps and, (b) a generalization of the notion of matrix product operators to the QTN setting which generates a linear combination of unitaries. The latter method is discussed in detail for translationally invariant spin-half systems, where it is shown to reduce the qubit resource requirements compared with the SWAP method and may be generalized to other systems. We demonstrate the versatility of our approaches by simulating spin-1/2 and spin-3/2 MNMs, with the latter being an experimentally relevant model of a Cr$^{3+}_{8}$ ring. Here, our approach has qubit requirements that are independent of the number of constituents of the many-body system and scale only logarithmically with the bond dimension of the MPS representation, making them appealing for implementation on near-term quantum hardware with mid-circuit measurement and reset.

Neutron scattering↗

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗

Development of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an Integrated Energy System (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. A new RAVEN (Risk Analysis Virtual ENvironment) plugin for grid and capacity optimization (HERON) has been developed for optimizing the production variables of the IES given the demand profile. Currently, only the limits that affect the production variables and their corresponding time rates of change are considered (explicit constraints). However, other variables are additionally subject to constraints, but the associated limits are not accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response effects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. The Reference Governor (RG) algorithm is traditionally used for enforcing state and control constraints by modifying the set-point trajectories supplied to the feedback regulators. In our application, the RG is coupled within an iterative loop with the HERON-power dispatcher to generate optimal trajectories that ensure the operational constraints are met. A data-driven procedure to derive a representation of the dynamics of the controlled system was developed. First, the variables that represented the state of the system are selected (PCA-based approach), and then state-space representation matrices are derived from the collected measurements (DMDc algorithm). A preliminary version of the developed workflow based on Linear Time Invariant matrices was assessed by adopting a two-unit test case. More sophisticated versions of this workflow foreseeing the on-line derivation of system matrices will be deployed in FY 2021. Finally, a “plug-and-play” library of controllers and state observers was developed in Dymola. Some aspects of the current configuration of the IES unit components, e.g., the encapsulation of the control schemes into dedicated blocks, are consistent with the “plug-and-play” philosophy. Other features, e.g., the system buses collecting the input and the output variables, are not. For this reason, once listed and described the limits of the current configuration, necessary modifications to the plant model interface are presented. As a test-case, the interfaces of the SES model in the RAVEN-based HYBRID framework were reworked accordingly, and two different control schemes were applied to the same plant model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Microgrid Assisted Design for Remote Areas

In this work, we present a three-stage multiobjective mixed-integer linear programming (MILP) for the optimal expansion planning and operation of isolated multienergy microgrids in remote areas. By selecting the optimal distributed generators (DGs) and energy storage systems (ESSs) mix selection, siting, sizing, and scheduling in the remote microgrid, the proposed model is targeted to minimize the annualized total cost of microgrids while enhancing the performance of the system, i.e., minimizing the voltage deviations and line power loss. To represent the electricity and heat flow between generation resources and various electrical, heating, and cooling loads in the isolated microgrid, linearized power flow, and heat flow constraints are employed in the proposed optimization model. The available capacity of DGs and ESSs are modeled as discrete constants instead of continuous variables for practical purpose. Numerical simulation results on a remote microgrid consisting of DGs, ESSs, and various loads validate the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Influence of Communication on the Complexity of Connected Lighting Systems

Experience in the Next Generation Lighting Systems (NGLS) ”Living Lab” at the Parsons School of Design, has clearly shown the impacts of communication on the configuration complexity of the connected lighting systems installed there. Since 2017, NGLS has invited manufacturers to submit connected lighting systems for assessment against a set of specifications. To date, 14 systems of wirelessly connected LED luminaires and controls have been installed in classrooms with existing lighting. In nine of the rooms, new LED linear pendants or recessed troffers replaced existing luminaires; in the other rooms, LED “kits” were retrofitted into the existing troffers. Eleven of the installations featured sensors integrated in the luminaires or retrofit kits; the others used ceiling mounted sensors. All systems provided wall control devices for manual operation. For half of the installations, the complete systems were furnished by single companies; for the other half, luminaires and controls were furnished by separate companies. A total of 18 different companies provided products, with four of these located outside the United States including, Korea, Poland, Germany, and the Netherlands. In addition, two different electrical contractors in New York City provided multiple installation teams totaling 10 individuals, none of whom were familiar with any of the systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Optimal Sizing and Dispatch of Solar Power with Storage

Designers of utility-scale solar plants with storage, seeking to maximize some aspect of plant performance, face multiple challenges. In many geographic locations, there is significant penetration of photovoltaic generation, which depresses energy prices during the hours of solar availability. An energy storage system affords the opportunity to dispatch during higher-priced time periods, but complicates plant design and dispatch decisions. Solar resource variability compounds these challenges, because determining optimal system sizes requires simultaneously considering how the plant will be operated under the imposed market and weather conditions. We develop an approach to analyze the economic performance of hybrid and single-technology solar power plants, which incorporates optimal dispatch, and considers the expected electricity market and weather conditions. We utilize the System Advisor Model software package to simulate the operation of multiple renewable generation and energy storage technologies, in conjunction with hourly-fidelity generation decisions determined by a revenue-maximizing, mixed-integer linear program. We show that, under our assumed market and weather conditions, the lifetime benefit-to-cost ratio can be improved by 6 to 19 percent, relative to a baseline design without optimizing, and that a concentrating solar power with thermal energy storage design produces significantly more energy per year, but is less profitable under our cost assumptions.

black-box optimization↗

Aqueous titanium redox flow batteries—State-of-the-art and future potential

Market-driven deployment of inexpensive (but intermittent) renewable energy sources, such as wind and solar, in the electric power grid necessitates grid-stabilization through energy storage systems Redox flow batteries (RFBs), with their rated power and energy decoupled (resulting in a sub-linear scaling of cost), are an inexpensive solution for the efficient electrochemical storage of large amounts of electrical energy. Titanium-based RFBs, first developed by NASA in the 1970s, are an interesting albeit less examined chemistry and are the focus of the present review. Ti, constituting 0.6% of the Earth’s crust and an ingredient in inexpensive white paints, is amongst the few elements (V and Mn being some others) which exhibit multiple soluble oxidation states in aqueous electrolytes. Further, the very high (approaching 10 M) solubility of Ti in low pH solutions suggests the possibility of developing exceptionally high energy density aqueous Redox Flow Batteries systems. With these advantages in mind, we present the state-of the-art in Ti-RFBs with a focus on Ti/Mn, Ti/Fe and Ti/Ce couples and systems that use Ti as an additive (such as Ti/V/Mn). The inherent advantages of inexpensive Ti actives and relatively high energy density is contrasted with potential side-reactions resulting in reduced energy efficiency. Technological pathways are presented with a view to overcoming critical bottlenecks and a vision is presented for the future development of Ti-RFBs.

25 ENERGY STORAGE↗

Block encoding bosons by signal processing

Block Encoding (BE) is a crucial subroutine in many modern quantum algorithms, including those with near-optimal scaling for simulating quantum many-body systems, which often rely on Quantum Signal Processing (QSP). Currently, the primary methods for constructing BEs are the Linear Combination of Unitaries (LCU) and the sparse oracle approach. In this work, we demonstrate that QSP-based techniques, such as Quantum Singular Value Transformation (QSVT) and Quantum Eigenvalue Transformation for Unitary Matrices (QETU), can themselves be efficiently utilized for BE implementation. Specifically, we present several examples of using QSVT and QETU algorithms, along with their combinations, to block encode Hamiltonians for lattice bosons, an essential ingredient in simulations of high-energy physics. We also introduce a straightforward approach to BE based on the exact implementation of Linear Operators Via Exponentiation and LCU (LOVE-LCU). We find that, while using QSVT for BE results in the best asymptotic gate count scaling with the number of qubits per site, LOVE-LCU outperforms all other methods for operators acting on up to qubits, highlighting the importance of concrete circuit constructions over mere comparisons of asymptotic scalings. Using LOVE-LCU to implement the BE, we simulate the time evolution of single-site and two-site systems in the lattice theory using the Generalized QSP algorithm and compare the gate counts to those required for Trotter simulation.

Kane, Christopher F↗

Linear-frequency conversion with time-varying metasurfaces

Frequency conversion is a hallmark of nonlinearity. The spectral manifestations, emergent within a system, can typically be attributed to a marked nonlinearity within the material properties, complex geometric configurations, and/or the unique functional form of interactions taking place in the constitutive subsystems. These phenomena, irrespective of their origins, have been harnessed and exploited in applications ranging from the generation of entangled photons, a cornerstone in quantum technologies, to nanomechanical frequency mixing, advancing subsurface scanning probe microscopy. Here, we propose a frequency conversion mechanism based on time-varying metasurfaces, an emerging frontier in metamaterial research. We show how temporal properties of metasurfaces can effectively emulate a nonlinear medium, thereby facilitating frequency conversion. The proposed material configuration has the potential not only to advance integrated photonics and quantum optics, but also to create opportunities in quantum sensing, quantum materials, and crucially quantum communications. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stable electrospray signal on a microfabricated glass chip with three-dimensional open edge and tiered depth geometries

This paper presents the microfabrication and performance of a three-dimensional electrospray ionization (ESI) emitter tip made from glass, which achieves stable current signals important for chemical analysis. Our novel microfabrication process and custom-built signal conditioning hardware provides the advantage of providing accurate features and steady signals. The fabrication process relies on standard microfabrication techniques (i.e., deposition, photolithography, and wet etching). This fabrication method involves the novel application of two layers of positive and negative photoresists in addition to Parafilm® wax tape. Open edge and tiered depth details were successfully created from a multilayer planar mask. This is a benefit for integrated miniaturized and microfluidic systems that often require micro features for their functionality but relatively large millimeter size features for their physical periphery. We demonstrate the fundamental performance of electrospray with this microfluidic chip. The emitter tip was fixed on a linear axis stage with high resolution (10 μm) to finely control the tip distance from a metal counter electrode plate. Further, a custom printed circuit board system was built to safely control four voltages applied to the microchip ports from a single high voltage power supply. To readily form the electrospray, non-aqueous solvents were used for their low viscosity and a constant voltage of +2.7 kV was applied to the sheath electrospray microchannel. The liquid being sprayed was 80/20 (v/v) methanol/acetonitrile with 0.1% acetic acid in the sheath microchannel and with ammonium acetate (10–40 mM) in its remaining microchannels. The electrospray signal was measured in response to varying the distance (1.4 to 1.6mm) between the electrospray emitter tip and a metal counter electrode plate in addition to the varying concentration of the background electrolyte, ammonium acetate. Stable and repeatable electrospray signal showed linear relationships with emitter tip distance and concentration (r 2 ≥ 0.95).

42 ENGINEERING↗

Conformational and oligomeric states of SPOP from small-angle X-ray scattering and molecular dynamics simulations

Speckle-type POZ protein (SPOP) is a substrate adaptor in the ubiquitin proteasome system, and plays important roles in cell-cycle control, development, and cancer pathogenesis. SPOP forms linear higher-order oligomers following an isodesmic self-association model. Oligomerization is essential for SPOP’s multivalent interactions with substrates, which facilitate phase separation and localization to biomolecular condensates. Structural characterization of SPOP in its oligomeric state and in solution is, however, challenging due to the inherent conformational and compositional heterogeneity of the oligomeric species. Here, we develop an approach to simultaneously and self-consistently characterize the conformational ensemble and the distribution of oligomeric states of SPOP by combining small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. We build initial conformational ensembles of SPOP oligomers using coarse-grained molecular dynamics simulations, and use a Bayesian/maximum entropy approach to refine the ensembles, along with the distribution of oligomeric states, against a concentration series of SAXS experiments. Our results suggest that SPOP oligomers behave as rigid, helical structures in solution, and that a flexible linker region allows SPOP’s substrate-binding domains to extend away from the core of the oligomers. Additionally, our results are in good agreement with previous characterization of the isodesmic self-association of SPOP. In the future, the approach presented here can be extended to other systems to simultaneously characterize structural heterogeneity and self-assembly.

59 BASIC BIOLOGICAL SCIENCES↗

Antibacterial ADP-ribosyl cyclase toxins inhibit bacterial growth by rapidly depleting NAD(P) +

In metazoans, enzymes belonging to the bifunctional ADP-ribosyl cyclase/cyclic ADP-ribose (cADPr) hydrolase family regulate diverse cellular processes by synthesizing and hydrolyzing the intracellular second messenger cADPr, derived from the electron carrier NAD+. However, bacterial enzymes belonging to this family have not been characterized. Here, we identify a bacterial ADP-ribosyl cyclase that is associated with the type VII secretion system and functions as an antibacterial toxin. This enzyme, which we name Tac1, inhibits bacterial growth by rapidly hydrolyzing NAD+ and NADP+. We determine the X-ray crystal structure of Tac1 to a resolution of 1.4 Å, which reveals that this protein adopts the core catalytic fold of metazoan ADP-ribosyl cyclase enzymes such as CD38. Using a combination of biochemical and mutagenesis approaches, we identify catalytic residues within the active site of Tac1, which are responsible for the formation of a cADPr catalytic intermediate and subsequent hydrolysis of this intermediate into linear ADP-ribose. A bioinformatic analysis reveals that Tac1 is the founding member of a widespread family of bacterial ADP-ribosyl cyclase enzymes, many of which are associated with interbacterial conflict systems. We also identify enzymes in this family that are not associated with biological conflict systems and demonstrate that they produce cADPr as their major product rather than linear ADP-ribose, suggesting that these enzymes serve a biological function distinct from interbacterial antagonism. Together, these findings demonstrate that ADP-ribosyl cyclase/cADPr hydrolase enzymes function as toxins in diverse bacterial conflict systems and suggest that cADPr may play a previously overlooked role in bacterial physiology.

Colautti, Jake↗

Neural network-based single material beam-hardening correction for X-ray CT in Additive Manufacturing

Beam-hardening (BH) artifacts are ubiquitous in X-ray CT scans of dense metal additively manufactured (AM) parts. While linearization approaches are useful for correcting beam-hardened data from single material objects, they either require a calibration scan or detailed system and material composition information. In this paper, we introduce a neural network-based, material-agnostic method to correct beam-hardening artifacts. We train a neural network to map the acquired beam-hardened projection values and the corresponding estimated thickness of the part based on an initial segmentation to beam-hardening related parameters, which can be used to compute the coefficients of a linearizing correction polynomial. A key strength of our approach is that, once the network is trained, it can be used for correcting beam hardening from a variety of materials without any calibration scans or detailed system and material composition information. Furthermore, our method is robust to errors in the estimated thickness due to the typical challenge of obtaining an accurate initial segmentation from reconstructions impacted by BH artifacts. We demonstrate the utility of our method to obtain high-quality CT reconstructions from a collection of AM parts -- suppressing cupping and streaking artifacts

Rahman, Obaid↗

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗