SEARCH · Engineering Papers
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.
NREL researcher Andrey Bernstein discusses asynchronous and distributed tracking of time-varying fixed points as part of the Innovative Optimization and Control Methods for Highly Distributed Autonomous Systems workshop. NREL hosted the workshop to exchange ideas on state-of-the-art control and optimization strategies as well as get feedback on work being done at NREL.
In this paper, we consider an optimal coordination problem for distributed energy resources (DERs) including distributed generators and energy storage devices. We first propose an algorithm based on the push-sum and gradient method to optimally coordinate storage devices and distributed generators in a distributed manner. In the proposed algorithm, each DER only maintains a set of variables and updates them through information exchange with a few neighboring DERs over a time-varying directed communication network. We show that the proposed distributed algorithm with appropriately chosen step-sizes solve the optimal DER coordination problem if the time-varying directed communication network is uniformly jointly strongly connected, which is a mild condition on the connectivity of communication topologies. The proposed distributed algorithm is illustrated and validated by numerical simulations.
We study synchronization of coupled linear systems over networks with weak connectivity and time-varying delays. We focus on the case that the internal dynamics are time-varying but non-expansive. Both uniformly connected and infinitely connected communication topologies are considered. A new concept of P-synchronization is introduced and we first show that global asymptotic P-synchronization can be achieved over directed networks with uniform joint connectivity and arbitrarily bounded delays. We then study the case of the infinitely jointly connected communication topology. In particular, for the undirected communication topologies, it turns out that the existence of a uniform time interval for the communication topology is not necessary and P-synchronization can be achieved when the time varying delays are arbitrarily bounded. Simulations are given to validate the theoretical results.
Abstract not provided.
Light fields carrying orbital angular momentum (OAM) provide powerful capabilities for applications in optical communications, microscopy, quantum optics, and microparticle manipulation. We introduce a property of light beams, manifested as a temporal OAM variation along a pulse: the self-torque of light. Although self-torque is found in diverse physical systems (i.e., electrodynamics and general relativity), it was not realized that light could possess such a property. We demonstrate that extreme-ultraviolet self-torqued beams arise in high-harmonic generation driven by time-delayed pulses with different OAM. We monitor the self-torque of extreme-ultraviolet beams through their azimuthal frequency chirp. Lastly, this class of dynamic-OAM beams provides the ability for controlling magnetic, topological, and quantum excitations and for manipulating molecules and nanostructures on their natural time and length scales.
In the analysis of observational data with mathematical models, it is often desirable to let one or more model parameters vary with time to account for changing environmental conditions or to offer more flexibility. A standard approach for dynamic models is to allow for the number of parameters to grow as a function of time, leading to a high-dimensional inference problem for long data records. We propose using Hamiltonian Monte Carlo, a gradient-based Markov chain Monte Carlo (MCMC) method, to perform Bayesian inference for time-varying parameters (TVP) in hydrology models. As derivatives of model error functions with regard to parameters are not available in closed form, we implemented the GR4J rainfall-runoff models in Theano/PyMC3, and Jax/NumPyro, allowing for the application of automatic differentiation techniques. A simulation experiment assessing the viability of these methods for recovering underlying temporal variation in parameters indicates that a discrete Gaussian random walk-based prior appears to be best suited across a range of scenarios while the continuous-time Gaussian process compares relatively poorly. In an analysis of 20 years of daily real-world streamflow records from the Model Parameter Estimation Experiment (MOPEX) hydrology data set with increasing data sparsity and at multiple temporal resolutions, we find that TVP inference with gradient-based MCMC is a flexible and reliable approach for analyzing hydrology models with dynamic parameter sets.
Simple models of therapy for viral diseases such as hepatitis C virus (HCV) or human immunodeficiency virus assume that, once therapy is started, the drug has a constant effectiveness. More realistic models have assumed either that the drug effectiveness depends on the drug concentration or that the effectiveness varies over time. Here a previously introduced varying-effectiveness (VE) model is studied mathematically in the context of HCV infection. We show that while the model is linear, it has no closed-form solution due to the time-varying nature of the effectiveness. We then show that the model can be transformed into a Bessel equation and derive an analytic solution in terms of modified Bessel functions, which are defined as infinite series, with time-varying arguments. Fitting the solution to data from HCV infected patients under therapy has yielded values for the parameters in the model. We show that for biologically realistic parameters, the predicted viral decay on therapy is generally biphasic and resembles that predicted by constant-effectiveness (CE) models. We introduce a general method for determining the time at which the transition between decay phases occurs based on calculating the point of maximum curvature of the viral decay curve. For the parameter regimes of interest, we also find approximate solutions for the VE model and establish the asymptotic behavior of the system. We show that the rate of second phase decay is determined by the death rate of infected cells multiplied by the maximum effectiveness of therapy, whereas the rate of first phase decline depends on multiple parameters including the rate of increase of drug effectiveness with time.
Interactive time-varying volume visualization is challenging due to its complex spatiotemporal features and sheer size of the dataset. Recent works transform the original discrete time-varying volumetric data into continuous Implicit Neural Representations (INR) to address the issues of compression, rendering, and super-resolution in both spatial and temporal domains. However, training the INR takes a long time to converge, especially when handling large-scale time-varying volumetric datasets. In this work, we proposed F-Hash, a novel feature-based multi-resolution Tesseract encoding architecture to greatly enhance the convergence speed compared with existing input encoding methods for modeling time-varying volumetric data. The proposed design incorporates multi-level collision-free hash functions that map dynamic 4D multi-resolution embedding grids without bucket waste, achieving high encoding capacity with compact encoding parameters. Our encoding method is agnostic to time-varying feature detection methods, making it a unified encoding solution for feature tracking and evolution visualization. Experiments show the F-Hash achieves state-of-the-art convergence speed in training various time-varying volumetric datasets for diverse features. We also proposed an adaptive ray marching algorithm to optimize the sample streaming for faster rendering of the time-varying neural representation.
In this paper, the optimal distributed linear time varying H2 control problem for cone causal spatially invariant linear time varying (SILTV) systems is considered. This class of systems was initially defined in the work of Voulgaris et al. for cone causal spatially invariant (SI) linear time invariant (LTI) systems. First, the optimal time varying H2 problem is posed using a version of the Youla parameterization. This allows to transform the problem into an affine form which results in a convex but infinite dimensional distributed optimization. Next, the problem is solved by computing an operator projection on a class of time varying Youla parameters with a cone causal distributed structure. The cone causal spatially invariant systems are viewed as input-output multiplicative operators with a mixed causal time varying structure with respect to time, and spatially invariant with respect to the space dimension.
Quantifying the time-varying value of energy efficiency is necessary to properly account for all of its benefits and costs and to identify and implement efficiency resources that contribute to a low-cost, reliable electric system. Historically, most quantification of the benefits of efficiency has focused largely on the economic value of annual energy reduction. Due to the lack of statistically representative metered end-use load shape data in Michigan (i.e., the hourly or seasonal timing of electricity savings), the ability to confidently characterize the time-varying value of energy efficiency savings in the state, especially for weather-sensitive measures such as central air conditioning, is limited. Still, electric utilities in Michigan can take advantage of opportunities to incorporate the time-varying value of efficiency into their planning. For example, end-use load research and hourly valuation of efficiency savings can be used for a variety of electricity planning functions, including load forecasting, demand-side management and evaluation, capacity planning, long-term resource planning, renewable energy integration, assessing potential grid modernization investments, establishing rates and pricing, and customer service (KEMA 2012). In addition, accurately calculating the time-varying value of efficiency may help energy efficiency program administrators prioritize existing offerings, set incentive or rebate levels that reflect the full value of efficiency, and design new programs.
Simulating quantum channels is a fundamental primitive in quantum computing, since quantum channels define general (trace-preserving) quantum operations. An arbitrary quantum channel cannot be exactly simulated using a finite-dimensional programmable quantum processor, making it important to develop optimal approximate simulation techniques. In this paper, we study the challenging setting in which the channel to be simulated varies adversarially with time. We propose the use of matrix exponentiated gradient descent (MEGD), an online convex optimization method, and analytically show that it achieves a sublinear regret in time. Through experiments, we validate the main results for time-varying dephasing channels using a programmable generalized teleportation processor.
An experiment on the propagation of electromagnetic (EM) signals in continuous time-varying plasma is described. The time-varying characteristics of plasma are considered to cause a parasitic modulation in both amplitude and phase, and the strength of this modulation, which carries the information of the electron density profile, is closely related to the plasma frequency and the incident wave frequency. Through theoretical analysis, we give an explanation and mechanism of the interaction between the continuous time-varying plasma and EM waves, which is verified by a comparative analysis with experiments performed under the same conditions. The effects of this modulation on the EM signals in the plasma sheath cannot be ignored.
This paper focuses on a time-varying constrained nonconvex optimization problem, and considers the synthesis and analysis of online regularized primal-dual gradient methods to track a Karush-Kuhn-Tucker (KKT) trajectory. The proposed regularized primal-dual gradient method is implemented in a running fashion, in the sense that the underlying optimization problem changes during the execution of the algorithms. In order to study its performance, we first derive its continuous-time limit as a system of differential inclusions. We then study sufficient conditions for tracking a KKT trajectory, and also derive asymptotic bounds for the tracking error (as a function of the time-variability of a KKT trajectory). Further, we provide a set of sufficient conditions for the KKT trajectories not to bifurcate or merge, and also investigate the optimal choice of the parameters of the algorithm. Illustrative numerical results for a time-varying nonconvex problem are provided.
Electric energy efficiency resources save energy and may reduce peak demand. Historically, quantification of energy efficiency benefits has largely focused on the economic value of energy savings during the first year and lifetime of the installed measures. Due in part to the lack of publicly available research on end-use load shapes (i.e., the hourly or seasonal timing of electricity savings) and energy savings shapes, consideration of the impact of energy efficiency on peak demand reduction (i.e., capacity savings) has been more limited. End-use load research and the hourly valuation of efficiency savings are used for a variety of electricity planning functions, including load forecasting, demand-side management and evaluation, capacity and demand response planning, long-term resource planning, renewable energy integration, assessing potential grid modernization investments, establishing rates and pricing, and customer service. This study reviews existing literature on the time-varying value of energy efficiency savings, provides examples in four geographically diverse locations of how consideration of the time-varying value of efficiency savings impacts the calculation of power system benefits, and identifies future research needs to enhance the consideration of the time-varying value of energy efficiency in cost-effectiveness screening analysis. Findings from this study include: -The time-varying value of individual energy efficiency measures varies across the locations studied because of the physical and operational characteristics of the individual utility system (e.g., summer or winter peaking, load factor, reserve margin) as well as the time periods during which savings from measures occur. -Across the four locations studied, some of the largest capacity benefits from energy efficiency are derived from the deferral of transmission and distribution system infrastructure upgrades. However, the deferred cost of such upgrades also exhibited the greatest range in value of all the components of avoided costs across the locations studied. -Of the five energy efficiency measures studied, those targeting residential air conditioning in summer-peaking electric systems have the most significant added value when the total time-varying value is considered. -The increased use of rooftop solar systems, storage, and demand response, and the addition of electric vehicles and other major new electricity-consuming end uses are anticipated to significantly alter the load shape of many utility systems in the future. Data used to estimate the impact of energy efficiency measures on electric system peak demands will need to be updated periodically to accurately reflect the value of savings as system load shapes change. -Publicly available components of electric system costs avoided through energy efficiency are not uniform across states and utilities. Inclusion or exclusion of these components and differences in their value affect estimates of the time-varying value of energy efficiency. -Publicly available data on end-use load and energy savings shapes are limited, are concentrated regionally, and should be expanded.
A fault isolation, estimation, and fault-tolerant control (FTC) scheme for nonlinear time-varying delay stochastic distribution control systems was presented in this paper. The Takagi-Sugeno fuzzy model was adopted to approach the nonlinear dynamics of time-varying delay systems. According to the output equivalence principle and Laplace transformation, an augmented state vector was given to solve the time-varying delay problem. When multiple actuator faults and interference occur simultaneously, fault detection, isolation and fault estimation was designed to obtained the fault information. To decouple faults and obtain the value and location information of the fault, the system was separated into two parts through the designed multiple conversion matrices, in which one subsystem was only affected by one actuator fault. This has simplified the design of fault isolation and estimation. A adaptive observer for fault estimation was given. Then, fault information such as the time, location, and size was determined. The observer gain matrices were calculated using linear matrix inequality (LMI). When a fault was detected and diagnosed, a FTC algorithm was devised using the proportional-integral control scheme to compensate the fault as much as possible. It has been shown that even if multiple faults actuator occurred simultaneously, the FTC controller still ensured the output probability density function of the system traced the desired probability density function when a fault occurred. Finally, the expected results were obtained through the simulation example, which confirmed the effectiveness of the method.
Wireless networks are characterized by nodes mobility, which makes the propagation environment time-varying and subject to fading. As a consequence, the statistical characteristics of the received signal vary continuously, giving rise to a Doppler power spectral density (DPSD) that varies from one observation instant to the next. This paper is concerned with dynamical modeling of time-varying wireless fading channels, their estimation and parameter identification, and optimal power control from received signal measurement data. The wireless channel is characterized using a stochastic state-space form and derived by approximating the time-varying DPSD of the channel. The expected maximization and Kalman filter are employed to recursively identify and estimate the channel parameters and states, respectively, from online received signal strength measured data. Moreover, we investigate a centralized optimal power control algorithm based on predictable strategies and employing the estimated channel parameters and states. The proposed models together with the estimation and power control algorithms are tested using experimental measurement data and the results are presented.