-
Jacobi Set Simplification for Tracking Topological Features in Time-Varying Scalar Fields
Authors:
Dhruv Meduri,
Mohit Sharma,
Vijay Natarajan
Abstract:
The Jacobi set of a bivariate scalar field is the set of points where the gradients of the two constituent scalar fields align with each other. It captures the regions of topological changes in the bivariate field. The Jacobi set is a bivariate analog of critical points, and may correspond to features of interest. In the specific case of time-varying fields and when one of the scalar fields is tim…
▽ More
The Jacobi set of a bivariate scalar field is the set of points where the gradients of the two constituent scalar fields align with each other. It captures the regions of topological changes in the bivariate field. The Jacobi set is a bivariate analog of critical points, and may correspond to features of interest. In the specific case of time-varying fields and when one of the scalar fields is time, the Jacobi set corresponds to temporal tracks of critical points, and serves as a feature-tracking graph. The Jacobi set of a bivariate field or a time-varying scalar field is complex, resulting in cluttered visualizations that are difficult to analyze. This paper addresses the problem of Jacobi set simplification. Specifically, we use the time-varying scalar field scenario to introduce a method that computes a reduced Jacobi set. The method is based on a stability measure called robustness that was originally developed for vector fields and helps capture the structural stability of critical points. We also present a mathematical analysis for the method, and describe an implementation for 2D time-varying scalar fields. Applications to both synthetic and real-world datasets demonstrate the effectiveness of the method for tracking features.
△ Less
Submitted 7 June, 2024;
originally announced July 2024.
-
Model reference adaptive control for state and input constrained linear systems
Authors:
Sudipta Chattopadhyay,
Srikant Sukumar,
Vivek Natarajan
Abstract:
State and input constraints are ubiquitous in all engineering systems. In this article, we derive adaptive controllers for uncertain linear systems under pre-specified state and input constraints. Several modifications of the model reference adaptive control (MRAC) framework have been proposed to address input constraints in uncertain linear systems. Considering the infeasibility of arbitrary refe…
▽ More
State and input constraints are ubiquitous in all engineering systems. In this article, we derive adaptive controllers for uncertain linear systems under pre-specified state and input constraints. Several modifications of the model reference adaptive control (MRAC) framework have been proposed to address input constraints in uncertain linear systems. Considering the infeasibility of arbitrary reference trajectories, reference modification has been implemented in the case of input constraints in literature. The resulting conditions on the reference and input signals are difficult to verify online. Similar results on state and input constraints together have also been proposed, albeit resulting in more complex and unverifiable conditions on the control. The primary objective of this article is therefore to account for state and input constraints in uncertain linear systems by providing easily verifiable conditions on the control and reference. A combination of reference modification and barrier Lyapunov methods in adaptive control are employed to arrive at these results.
△ Less
Submitted 23 August, 2023;
originally announced August 2023.
-
Topological Feature Search in Time-Varying Multifield Data
Authors:
Tripti Agarwal,
Amit Chattopadhyay,
Vijay Natarajan
Abstract:
A wide range of data that appear in scientific experiments and simulations are multivariate or multifield in nature, consisting of multiple scalar fields. Topological feature search of such data aims to reveal important properties useful to the domain scientists. It has been shown in recent works that a single scalar field is insufficient to capture many important topological features in the data,…
▽ More
A wide range of data that appear in scientific experiments and simulations are multivariate or multifield in nature, consisting of multiple scalar fields. Topological feature search of such data aims to reveal important properties useful to the domain scientists. It has been shown in recent works that a single scalar field is insufficient to capture many important topological features in the data, instead one needs to consider topological relationships between multiple scalar fields. In the current paper, we propose a novel method of finding similarity between two multifield data by comparing their respective fiber component distributions. Given a time-varying multifield data, the method computes a metric plot for each pair of histograms at consecutive time stamps to understand the topological changes in the data over time. We validate the method using real and synthetic data. The effectiveness of the proposed method is shown by its ability to capture important topological features that are not always possible to detect using the individual component scalar fields.
△ Less
Submitted 2 November, 2019;
originally announced November 2019.
-
Minimal order controllers for output regulation of nonlinear systems
Authors:
Vivek Natarajan,
George Weiss
Abstract:
This paper is about the nonlinear local error feedback regulator problem. The plant is a nonlinear finite-dimensional system with a single control input and a single output and it is locally exponentially stable around the origin. The plant is driven, via a separate disturbance input, by a Lyapunov stable exosystem whose states are nonwandering. The reference signal that the plant output must trac…
▽ More
This paper is about the nonlinear local error feedback regulator problem. The plant is a nonlinear finite-dimensional system with a single control input and a single output and it is locally exponentially stable around the origin. The plant is driven, via a separate disturbance input, by a Lyapunov stable exosystem whose states are nonwandering. The reference signal that the plant output must track is a nonlinear function of the exosystem state. The local error feedback regulator problem is to design a dynamic feedback controller, with the tracking error as its input, such that (i) the closed-loop system of the plant and the controller is locally exponentially stable, and (ii) the tracking error tends to zero for all sufficiently small initial conditions of the plant, the controller and the exosystem. Under the assumption that the above regulator problem is solvable, we propose a nonlinear controller whose order is relatively small - typically equal to the order of the exosystem, and which solves the regulator problem. The emphasis is on the low order of the controller. The stability assumption on the plant (which can be relaxed to some extent) is crucial for making it possible to design a low order controller. We will show, under certain assumptions, that our proposed controller is of minimal order. Three examples are presented - the first illustrates our controller design procedure using an exosystem whose trajectories are periodic even though the state operator of the linearized exosystem contains a nontrivial Jordan block. The second example is more involved, and shows that sometimes a nontrivial immersion of the exosystem is needed in the design. The third example, based on output voltage regulation for a boost power converter, shows how the regulator equations may reduce to a first order PDE with no given boundary conditions, but which nevertheless has a locally unique solution.
△ Less
Submitted 5 January, 2019;
originally announced January 2019.
-
The global bifurcation picture for ground states in nonlinear Schrodinger equations
Authors:
Eduard Kirr,
Vivek Natarajan
Abstract:
In this paper, we propose a method of finding all coherent structures supported by a given nonlinear wave equation. It relies on enhancing the recent global bifurcation theory as developed by Dancer, Toland, Buffoni and others, by determining all the limit points of the coherent structure manifolds at the boundary of the Fredholm domain. Local bifurcation theory is then used to trace back these ma…
▽ More
In this paper, we propose a method of finding all coherent structures supported by a given nonlinear wave equation. It relies on enhancing the recent global bifurcation theory as developed by Dancer, Toland, Buffoni and others, by determining all the limit points of the coherent structure manifolds at the boundary of the Fredholm domain. Local bifurcation theory is then used to trace back these manifolds from their limit points into the interior of the Fredholm domain identifying the singularities along them. This way all coherent structure manifold are discovered except may be the ones which form loops and hence never reach the boundary. The method is then applied to the Schrodinger equation with a power nonlinearity for which all ground states are identified.
△ Less
Submitted 14 November, 2018;
originally announced November 2018.
-
Stability of the integral control of stable nonlinear systems
Authors:
George Weiss,
Vivek Natarajan
Abstract:
PI controllers are the most widespread type of controllers and there is an intuitive understanding that if their gains are sufficiently small and of the correct sign, then they always work. In this paper we try to give some rigorous backing to this claim, under specific assumptions. Let $\bf P$ be a nonlinear system described by $\dot x=f(x,u)$, $y=g(x)$, where the state trajectory $x$ takes value…
▽ More
PI controllers are the most widespread type of controllers and there is an intuitive understanding that if their gains are sufficiently small and of the correct sign, then they always work. In this paper we try to give some rigorous backing to this claim, under specific assumptions. Let $\bf P$ be a nonlinear system described by $\dot x=f(x,u)$, $y=g(x)$, where the state trajectory $x$ takes values in $R^n$, $u$ and $y$ are scalar and $f,g$ are of class $C^1$. We assume that there is a Lipschitz function $Ξ:[u_{min},u_{max}]\rightarrow R^n$ such that for every constant input $u_0\in[u_{min},u_{max}]$, $Ξ(u_0)$ is an exponentially stable equilibrium point of $\bf P$. We also assume that $G(u)=g(Ξ(u))$, which is the steady state input-output map of $\bf P$, is strictly increasing. Denoting $y_{min}=G(u_{min})$ and $y_{max}=G(u_{max})$, we assume that the reference value $r$ is in $(y_{min},y_{max})$. Our aim is that $y$ should track $r$, i.e., $y\rightarrow r$ as $t\rightarrow\infty$, while the input of $P$ is only allowed to be in $[u_{min},u_{max}]$. For this, we introduce a variation of the integrator, called the saturating integrator, and connect it in feedback with $\bf P$ in the standard way, with gain $k>0$. We show that for any small enough $k$, the closed-loop system is (locally) exponentially stable around an equilibrium point $(Xi(u_r),u_r)$, with a large region of attraction $X_T\subset R^n\times[u_{min},u_{max}]$. When the state $(x(t),u(t))$ of the closed-loop system converges to $(Ξ(u_r),u_r)$, then the tracking error $r-y$ tends to zero. The compact set $X_T$ can be made larger by choosing a larger parameter $T>0$, resulting in smaller $k$.
△ Less
Submitted 13 January, 2021; v1 submitted 16 October, 2016;
originally announced October 2016.
-
Almost global asymptotic stability of a grid-connected synchronous generator
Authors:
Vivek Natarajan,
George Weiss
Abstract:
We study the global asymptotic behavior of a grid-connected constant field current synchronous generator (SG). The grid is regarded as an "infinite bus", i.e. a three-phase AC voltage source. The generator does not include any controller other than the frequency droop loop. This means that the mechanical torque applied to this generator is an affine function of its angular velocity. The negative s…
▽ More
We study the global asymptotic behavior of a grid-connected constant field current synchronous generator (SG). The grid is regarded as an "infinite bus", i.e. a three-phase AC voltage source. The generator does not include any controller other than the frequency droop loop. This means that the mechanical torque applied to this generator is an affine function of its angular velocity. The negative slope of this function is the frequency droop constant. We derive sufficient conditions on the SG parameters under which there exist exactly two periodic state trajectories for the SG, one stable and another unstable, and for almost all initial states, the state trajectory of the SG converges to the stable periodic trajectory (all the angles are measured modulo $2π$). Along both periodic state trajectories, the angular velocity of the SG is equal to the grid frequency. Our sufficient conditions are easy to check computationally. An important tool in our analysis is an integro-differential equation called the {\em exact swing equation}, which resembles a forced pendulum equation and is equivalent to our fourth order model of the grid-connected SG. Apart from our objective of providing an analytical proof for a global asymptotic behavior observed in a classical dynamical system, a key motivation for this work is the development of synchronverters which are inverters that mimic the behavior of SGs. Understanding the global dynamics of SGs can guide the choice of synchronverter parameters and operation. As an application we find a set of stable nominal parameters for a 500 kW synchronverter.
△ Less
Submitted 16 October, 2016;
originally announced October 2016.
-
Felix: A Topology based Framework for Visual Exploration of Cosmic Filaments
Authors:
Nithin Shivshankar,
Pratyush Pranav,
Vijay Natarajan,
Rien van de Weygaert,
E G Patrick Bos,
Steven Rieder
Abstract:
The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which s…
▽ More
The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which span a relatively larger range. The large scale filamentary network thus forms an intricate part of the cosmic web.
In this paper, we describe Felix, a topology based framework for visual exploration of filaments in the cosmic web. The filamentary structure is represented by the ascending manifold geometry of the 2-saddles in the Morse-Smale complex of the density field. We generate a hierarchy of Morse-Smale complexes and query for filaments based on the density ranges at the end points of the filaments. The query is processed efficiently over the entire hierarchical Morse-Smale complex, allowing for interactive visualization.
We apply Felix to computer simulations based on the heuristic Voronoi kinematic model and the standard $Λ$CDM cosmology, and demonstrate its usefulness through two case studies. First, we extract cosmic filaments within and across cluster like regions in Voronoi kinematic simulation datasets. We demonstrate that we produce similar results to existing structure finders. Filaments that form the spine of the cosmic web, which exist in high density regions in the current epoch, are isolated using Felix. Also, filaments present in void-like regions are isolated and visualized. These filamentary structures are often over shadowed by higher density range filaments and are not easily characterizable and extractable using other filament extraction methodologies.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
