Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lyapunov function”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A Lyapunov Function Based Remedial Action Screening Tool Using Real-Time Data

This report summarizes the outcome of a research project that comprised the development of a Lyapunov function based remedial action screening tool using real-time data (L-RAS). The L-RAS is an advanced computational tool that is intended to assist system operators in making real-time redispatch decisions to preserve power grid stability. The tool relies on screening contingencies using a homotopy method based on Lyapunov functions to avoid, to the extent possible, the use of time domain simulations. This enables transient stability evaluation at real-time speed without the use of massively parallel computational resources. The project combined the following components. 1. Development of a methodology for contingency screening using a homotopy method based on Lyapunov functions and real-time data. 2. Development of a methodology for recommending remedial actions based on the screening results. 3. Development of a visualization and operator interaction interface. 4. Testing of screening tool, validation of control actions, and demonstration of project outcomes on a representative real system simulated on a Real-Time Digital Simulator (RTDS) cluster. The project was led by Michigan State University (MSU), where the theoretical models including homotopy-based screening, trajectory correction using real-time data, and remedial action were developed and implemented in the form of research-grade software. Los Alamos National Laboratory (LANL) contributed to the development of energy margin sensitivity dynamics, which constituted a part of the remedial action portfolio. Florida State University (FSU) and Southern California Edison (SCE) developed a model of the SCE system that was implemented on FSU's RTDS cluster to simulate real-time data that was streamed over the internet to MSU where the L-RAS tool was executed and remedial actions were communicated back to FSU to execute stabilizing controls on the simulated system. LCG Consulting developed the visualization and operator interaction interface, based on specifications provided by MSU. The project was performed from October 2012 to December 2016, at the end of which the L-RAS tool, as described above, was completed and demonstrated. The project resulted in the following innovations and contributions: (a) the L-RAS software prototype, tested on a simulated system, vetted by utility personnel, and potentially ready for wider testing and commercialization; (b) an RTDS-based test bed that can be used for future research in the field; (c) a suite of breakthrough theoretical contributions to the field of power system stability and control; and (d) a new tool for visualization of power system stability margins. While detailed descriptions of the development and implementation of the various project components have been provided in the quarterly reports, this final report provides an overview of the complete project, and is demonstrated using public domain test systems commonly used in the literature. The SCE system, and demonstrations thereon, are not included in this report due to Critical Energy Infrastructure Information (CEII) restrictions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Adaptive Path-Following Control for Ground Vehicles Using a Switching Non-Quadratic Lyapunov Function

The application of adaptive control techniques in the development of control systems for intelligent vehicles, especially for ground vehicle path-following controllers, has gained popularity due to their ability to handle large-scale parametric uncertainties. However, the use of a standard quadratic Lyapunov function in existing adaptive control-based path-following controllers can lead to poor transient performance, such as slow convergence and/or large overshoot. To address this limitation, this study proposes the use of a switching non-quadratic Lyapunov function to design a model reference adaptive path-following controller that aims to provide superior transient performance. The stability and signal convergence of the closed-loop system are demonstrated through a Lyapunov-like analysis. Through dSPACE ASM simulation, the effectiveness of the proposed controller is illustrated, which confirms improved tracking performance over a baseline solution.

Zhou, Xingyu↗
DOE OSTI2022

Provably-Stable Overload Ride-Through Control for Grid-Forming Inverters Using System-Wide Lyapunov Function Analysis

A key challenge associated with a grid-forming (GFM) inverter based resource (IBR) is its behavior during severe grid disturbances: since a GFM inverter regulates voltage in the fast timescale instead of current or power, it may experience a transient overload of current, power and/or energy during a severe grid disturbance. While many promising control strategies for overload ride-through have been proposed over the past two decades, transient stability of the system during and after the transition to an overload ride-through control mode remains difficult to guarantee. In this work, a novel overload ride-through control strategy is proposed for a system of grid-forming inverters that takes both self-protection and system-wide transient stability into account. A proposed system-level supervisory control uses slow communication to pre-emptively assign a set of local ride through control parameters to individual GFM IBR, including a current-limiting virtual reactance, that guarantees that synchronism is still preserved for any set of anticipated grid disturbances. At the core of the supervisory control lies a Lyapunov-function-based routine capable of establishing a strong, albeit conservative, transient stability guarantee for the system. Here, the proposed overload ride-through control strategy is validated via numerical integration of a reduced-order model, as well as through detailed electromagnetic transient (EMT) simulation.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Network analysis of memristive device circuits: dynamics, stability and correlations

Abstract Networks with memristive devices are a potential basis for the next generation of computing devices. They are also an important model system for basic science, from modeling nanoscale conductivity to providing insight into the information-processing of neurons. The resistance in a memristive device depends on the history of the applied bias and thus displays a type of memory. The interplay of this memory with the dynamic properties of the network can give rise to new behavior, offering many fascinating theoretical challenges. But methods to analyze general memristive circuits are not well described in the literature. In this paper we develop a general circuit analysis for networks that combine memristive devices alongside resistors, capacitors and inductors and under various types of control. We derive equations of motion for the memory parameters of these circuits and describe the conditions for which a network should display properties characteristic of a resonator system. For the case of a purely memresistive network, we derive Lyapunov functions, which can be used to study the stability of the network dynamics. Surprisingly, analysis of the Lyapunov functions show that these circuits do not always have a stable equilibrium in the case of nonlinear resistance and window functions. The Lyapunov function allows us to study circuit invariances, wherein different circuits give rise to similar equations of motion, which manifest through a gauge freedom and node permutations. Finally, we identify the relation between the graph Laplacian and the operators governing the dynamics of memristor networks operators, and we use these tools to study the correlations between distant memristive devices through the effective resistance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Deriving GENERIC from a Generalized Fluctuation Symmetry

Much of the structure of macroscopic evolution equations for relaxation to equilibrium can be derived from symmetries in the dynamical fluctuations around the most typical trajectory. For example, detailed balance as expressed in terms of the Lagrangian for the path-space action leads to gradient zero-cost flow. We expose a new such fluctuation symmetry that implies GENERIC, an extension of gradient flow where a Hamiltonian part is added to the dissipative term in such a way as to retain the free energy as Lyapunov function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Explicit Weierstrass traveling wave solutions and bifurcation analysis for dissipative Zakharov–Kuznetsov modified equal width equation

We consider the dissipative Zakharov–Kuznetsov modified equal width (ZK-MEW) equation and discuss the effect of the dissipation on the existence and nature of traveling wave solutions of the equation. We use Lyapunov function and dynamical system theory in order to show that when viscosity term is added to the ZK-MEW equation, yet there exists still bounded traveling wave solutions in certain regions. Subsequently, we obtain general solution of the ZK-MEW equation in the presence and absence of viscosity in terms of Weirstrass ℘ functions and Jacobi elliptic functions. We also derive a new form of kink-type solution by exploring a factorization method based on functional transformation and the Abel’s first order nonlinear equation. Finally, we use the phase plane analysis and examine the stability of the viscous waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Global dynamics of a vector-borne disease model with infection ages and general incidence rates

A vector-borne disease model with general incidence rates is proposed and investigated in this paper, where both vector and host are stratied by infection ages in the form of a hyperbolic system of partial differential equations coupled with ordinary differential equations.The existence, uniqueness, nonnegativeness, and boundedness of solution of the model are studied for biologically reasonable purpose. Furthermore, a global threshold dynamics of the system is established by constructing suitable Lyapunov functionals, which is determined by the basic reproduction number R{sub 0}: the infection-free equilibrium is globally asymptotically stable when R{sub 0} < 1 while the endemic equilibrium is globally asymptotically stable when R{sub 0} > 1.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Cluster synchronization induced by one-node clusters in networks with asymmetric negative couplings

This paper deals with the problem of cluster synchronization in networks with asymmetric negative couplings. By decomposing the coupling matrix into three matrices, and employing Lyapunov function method, sufficient conditions are derived for cluster synchronization. The conditions show that the couplings of multi-node clusters from one-node clusters have beneficial effects on cluster synchronization. Based on the effects of the one-node clusters, an effective and universal control scheme is put forward for the first time. The obtained results may help us better understand the relation between cluster synchronization and cluster structures of the networks. The validity of the control scheme is confirmed through two numerical simulations, in a network with no cluster structure and in a scale-free network.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Dynamics of the stochastic Lorenz chaotic system with long memory effects

Little seems to be known about the ergodic dynamics of stochastic systems with fractional noise. This paper is devoted to discern such long time dynamics through the stochastic Lorenz chaotic system (SLCS) with long memory effects. By a truncation technique, the SLCS is proved to generate a continuous stochastic dynamical system Λ. Based on the Krylov-Bogoliubov criterion, the required Lyapunov function is further established to ensure the existence of the invariant measure of Λ. Meanwhile, the uniqueness of the invariant measure of Λ is proved by examining the strong Feller property, together with an irreducibility argument. Therefore, the SLCS has exactly one adapted stationary solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A large scale system approach for stability analysis of multi-unit SMR plants - 248

In this paper we present a framework based on large scale system theory for analyzing the stability of coupled advanced high temperature molten salt small modular reactors. It has been suggested that small modular reactors be coupled through a common molten salt energy repository, called a salt vault. Using vector Lyapunov functions the connective stability of multiple reactors coupled through a salt vault can be guaranteed. In this theory, scalar Lyapunov functions are assigned to portions of the composite system state space. Functions are constructed for each small modular reactor (SMR), each Brayton cycle, and the salt vault. These functions are then aggregated, becoming components of a vector Lyapunov function. If this vector function possesses special properties, the composite system is stable. Bounds on the connective stability of a multi-reactor plant can be found in terms of a free parameter, in this case the mass of molten salt in the salt vault. Examining bounds on the connective stability allows one to use control theory as a tool for making performance and economic design decisions related to the size of the salt vault. This technique can produce a constraint that is used in a design optimization problem for finding the smallest-possible salt vault that ensures stability. It can also be used as an indicator of interconnection and its effect on system behavior. This, in turn, could be used to guide multi-unit control system design to ensure that highly coupled systems behave favorably. The theory for multi-unit control and operation is discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Asymptotic Behavior of Memristive Circuits

The interest in memristors has risen due to their possible application both as memory units and as computational devices in combination with CMOS. This is in part due to their nonlinear dynamics, and a strong dependence on the circuit topology. We provide evidence that also purely memristive circuits can be employed for computational purposes. In the present paper we show that a polynomial Lyapunov function in the memory parameters exists for the case of DC controlled memristors. Such a Lyapunov function can be asymptotically approximated with binary variables, and mapped to quadratic combinatorial optimization problems. This also shows a direct parallel between memristive circuits and the Hopfield-Little model. In the case of Erdos-Renyi random circuits, we show numerically that the distribution of the matrix elements of the projectors can be roughly approximated with a Gaussian distribution, and that it scales with the inverse square root of the number of elements. This provides an approximated but direct connection with the physics of disordered system and, in particular, of mean field spin glasses. Using this and the fact that the interaction is controlled by a projector operator on the loop space of the circuit. We estimate the number of stationary points of the approximate Lyapunov function and provide a scaling formula as an upper bound in terms of the circuit topology only.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Synchronization of Lienard-Type Oscillators in Uniform Electrical Networks

This paper presents a condition for global asymptotic synchronization of Lienard-type nonlinear oscillators in uniform LTI electrical networks with series R-L circuits modeling interconnections. By uniform electrical networks, we mean that the per-unit-length impedances are identical for the interconnecting lines. We derive conditions for global asymptotic synchronization for a particular feedback architecture where the derivative of the oscillator output current supplements the innate current feedback induced by simply interconnecting the oscillator to the network. Our proof leverages a coordinate transformation to a set of differential coordinates that emphasizes signal differences and the particular form of feedback permits the formulation of a quadratic Lyapunov function for this class of networks. This approach is particularly interesting since synchronization conditions are difficult to obtain by means of quadratic Lyapunov functions when only current feedback is used and for networks composed of series R-L circuits. Our synchronization condition depends on the algebraic connectivity of the underlying network, and reiterates the conventional wisdom from Lyapunov- and passivity-based arguments that strong coupling is required to ensure synchronization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Global properties and bifurcation analysis of an HIV-1 infection model with two target cells

In this paper, the authors consider a four-compartmental HIV epidemiological model, which describes the interaction between HIV virus and two target cells, CD4 T cells and macrophages in vivo. It is proved that the bilinear incidence can cause the backward bifurcation, where a locally asymptotically stable disease-free equilibrium co-exists with a locally asymptotically stable endemic equilibrium when the basic reproduction number (R{sub 0}) is less than unity. It is shown that a sequence of Hopf bifurcations occur at the endemic equilibrium by choosing one parameter of the model as the bifurcation parameter. Meanwhile, the global asymptotic stabilities of the equilibria are established by constructing suitable Lyapunov functions under some conditions. Furthermore, the authors develop an extended model by incorporating with the intracellular delays and derive global asymptotic stability of the delayed model by constructing Lyapunov functions. Some numerical simulations for justifying the theoretical analysis results are also given.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

A mean-field model of memristive circuit interaction

We construct an exactly solvable circuit of interacting memristors and study its dynamics and fixed points. This simple circuit model interpolates between decoupled circuits of isolated memristors, and memristors in series, for which exact fixed points can be obtained. Here, we introduce a Lyapunov functional that is found to be minimized along the non-equilibrium dynamics and which resembles a long-range Ising Hamiltonian with non-linear self-interactions. We use the Lyapunov functional as a Hamiltonian to calculate, in the mean-field theory approximation, the average asymptotic behavior of the circuit given a random initialization, yielding exact predictions for the case of decay to the lower resistance state, and reasonable predictions for the case of a decay to the higher resistance state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data-driven modeling and control of dynamical systems using Koopman and Perron-Frobenius operators

This dissertation studies the data-driven modeling and control problem of nonlinear systems by exploiting the linear operator theoretic framework involving Koopman and Perro-Frobenius operator. A systematic linear-operator based controller design procedure has been established, which can be used to solve a variety of nonlinear control problems, including feedback stabilization using control Lyapunov functions, optimal quadratic regulation using Koopman eigenfunctions and convex optimization formulation of optimal control problem using P-F and Koopman operator approximation. As the core of data-driven modeling, we first propose a new algorithm for the finite-dimensional approximation of the linear transfer Koopman and Perron-Frobenius operator from time-series data. We argue that the existing approach for the finite-dimensional approximation of these transfer operators such as Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD) do not capture two important properties of these operators, namely positivity and Markov property. The algorithm we propose preserves these two properties. We call the proposed algorithm as naturally structured DMD (NSDMD) since it retains the inherent properties of these operators. Naturally structured DMD algorithm leads to a better approximation of the steady-state dynamics of the system regarding computing Koopman and Perron- Frobenius operator eigenfunctions and eigenvalues. However, preserving positivity property is critical for capturing the real transient dynamics of the system. This positivity property of the transfer operators and it's finite-dimensional approximation play an important role for controller and estimator design of nonlinear systems. To solve the feedback stabilization problem for nonlinear control systems, we tried to take advantage of the Koopman operator framework. The Koopman operator approach provides a linear representation for a nonlinear dynamical system and a bilinear representation for a nonlinear control system. The problem of feedback stabilization of a nonlinear control system is then transformed to the stabilization of a bilinear control system. We propose a control Lyapunov function (CLF)-based approach for the design of stabilizing feedback controllers for the bilinear system. The search for finding a CLF for the bilinear control system is formulated as a convex optimization problem. This leads to a schematic procedure for designing CLF-based stabilizing feedback controllers for the bilinear system and hence the original nonlinear system. Another advantage of the proposed controller design approach outlined in this dissertation is that it does not require explicit knowledge of system dynamics. In particular, the bilinear representation of a nonlinear control system in the Koopman eigenfunction space can be obtained from time-series data. Next, we study the optimal quadratic regulation problem for nonlinear systems. The linear operator theoretic framework involving the Koopman operator is used to lift the dynamics of nonlinear control system to an infinite-dimensional bilinear system. The optimal quadratic regulation problem for nonlinear system is formulated in terms of the finite-dimensional approximation of the bilinear system. A convex optimization-based approach is proposed for solving the quadratic regulator problem for bilinear system. We applied a variety of examples and compared the simulation results between our framework and conventional LQR control using linearized model. For more general optimal control problems, we provide a density-function based convex formulation for the optimal control problem of the nonlinear system. The convex formulation relies on the duality result in the stability theory of a dynamical system involving density function and Perron-Frobenius operator. The optimal control problem is formulated as an infinite-dimensional convex optimization program. The finite-dimensional approximation of the optimization problem relies on the recent advances made in the data-driven computation of the Koopman operator, which is dual to the Perron-Frobenius operator. Simulation results are presented to demonstrate the application of the developed framework.

Huang, Bowen↗
DOE OSTI2021

A Convex Data-Driven Approach for Nonlinear Control Synthesis

We consider a class of nonlinear control synthesis problems where the underlying mathematical models are not explicitly known. We propose a data-driven approach to stabilize the systems when only sample trajectories of the dynamics are accessible. Our method is built on the density-function-based stability certificate that is the dual to the Lyapunov function for dynamic systems. Unlike Lyapunov-based methods, density functions lead to a convex formulation for a joint search of the control strategy and the stability certificate. This type of convex problem can be solved efficiently using the machinery of the sum of squares (SOS). For the data-driven part, we exploit the fact that the duality results in the stability theory can be understood through the lens of Perron–Frobenius and Koopman operators. This allows us to use data-driven methods to approximate these operators and combine them with the SOS techniques to establish a convex formulation of control synthesis. The efficacy of the proposed approach is demonstrated through several examples.

97 MATHEMATICS AND COMPUTING↗