An official website of the United States government
The ability of a nonlinear system to recover from a large disturbance to a desired stable equilibrium point depends on system parameter values, which are often uncertain and time varying. A particular disturbance acting for a finite time can be modeled as an implicit map that takes a parameter value to its corresponding post disturbance initial condition in state space. The system recovers when the post-disturbance initial condition lies inside the region of attraction of the stable equilibrium point. Critical parameter values are defined to be parameter values whose corresponding post-disturbance initial condition lies on the boundary of the region of attraction. Computing such values is important in numerous applications because they represent the boundary between desirable and undesirable system behavior. Many realistic system models involve controller clipping limits and other forms of switching. Furthermore, these hybrid dynamics are closely linked to the ability of a system to recover from disturbances. The paper develops theory which underpins a novel algorithm for numerically computing critical parameter values for nonlinear systems with clipping limits and switching. For an almost generic class of vector fields with event-selected discontinuities, it is shown that the boundary of the region of attraction is equal to a union of the stable manifolds of the equilibria and periodic orbits it contains, and that this decomposition persists and the boundary varies continuously under small changes in parameter.
more »
« less
Determining Parameter Space Margins for Fault Recovery By Inverse Sensitivity Minimization
Power system model parameter values are becoming increasingly uncertain and time-varying. Therefore, it is important to determine the margin in parameter space between a given set of parameter values for which the system will recover from a particular fault, and the nearest parameter values for which it will not recover from that fault. This work presents an efficient method for computing parameter space recovery margins by exploiting the property that the trajectory becomes infinitely sensitive to small changes in parameter value along the operating point’s region of attraction boundary. Consequently, along this boundary the inverse sensitivity of the trajectory approaches zero. The method proceeds by varying parameter values so as to minimize the inverse sensitivity of the system trajectory. Recent results provide theoretical justification for the approach. The efficacy of the method is demonstrated using a modified IEEE 39-bus New England power system test case.
more »
« less
- Award ID(s):
- 1810144
- PAR ID:
- 10126163
- Date Published:
- Journal Name:
- 21st Power Systems Computation Conference
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Geodetic observations of post-seismic deformation due to afterslip and viscoelastic relaxation can be used to infer fault and lithosphere rheologies by combining the observations with mechanical models of post-seismic processes. However, estimating the spatial distributions of rheological parameters remains challenging because it requires solving a nonlinear inverse problem with a high-dimensional parameter space and potentially computationally expensive forward model. Here we introduce an inversion method to estimate spatially varying fault and lithospheric rheological parameters in a mechanical model of post-seismic deformation using geodetic time series. The forward model combines afterslip and viscoelastic relaxation governed by a velocity-strengthening frictional rheology and a power-law Burgers rheology, respectively, and incorporates the mechanical coupling between coseismic slip, afterslip and viscoelastic relaxation. The inversion method estimates spatially varying fault frictional parameters, viscoelastic constitutive parameters and coseismic stress change. We formulate the inverse problem in a Bayesian framework to quantify the uncertainties of the estimated parameters. To solve this problem with reasonable computational costs, we develop an algorithm to estimate the mean and covariance matrix of the posterior probability distribution based on an ensemble Kalman filter. We validate our method through numerical tests using a 2-D forward model and synthetic post-seismic GNSS time-series. The test results suggest that our method can estimate the spatially varying rheological parameters and their uncertainties reasonably well with tolerable computational costs. Our method can also recover spatially and temporally varying afterslip, viscous strain and effective viscosities and can distinguish the contributions of afterslip and viscoelastic relaxation to observed post-seismic deformation.more » « less
-
Abstract In this paper, we consider the inverse scattering problem associated with an anisotropic medium with a conductive boundary. We will assume that the corresponding far–field pattern is known/measured and we consider two inverse problems. First, we show that the far–field data uniquely determines the boundary coefficient. Next, since it is known that anisotropic coefficients are not uniquely determined by this data we will develop a qualitative method to recover the scatterer. To this end, we study the so–called monotonicity method applied to this inverse shape problem. This method has recently been applied to some inverse scattering problems but this is the first time it has been applied to an anisotropic scatterer. This method allows one to recover the scatterer by considering the eigenvalues of an operator associated with the far–field operator. We present some simple numerical reconstructions to illustrate our theory in two dimensions. For our reconstructions, we need to compute the adjoint of the Herglotz wave function as an operator mapping intoH1of a small ball.more » « less
-
Abstract Cloud and convective parameterizations strongly influence uncertainties in equilibrium climate sensitivity. We provide a proof‐of‐concept study to constrain these parameterizations in a perturbed parameter ensemble of the atmosphere‐only version of the Goddard Institute for Space Studies Model E2.1 simulations by evaluating model biases in the present‐day runs using multiple satellite climatologies and by comparing simulated δ18O of precipitation (δ18Op), known to be sensitive to parameterization schemes, with a global database of speleothem δ18O records covering the Last Glacial Maximum (LGM), mid‐Holocene (MH) and pre‐industrial (PI) periods. Relative to modern interannual variability, paleoclimate simulations show greater sensitivity to parameter changes, allowing for an evaluation of model uncertainties over a broader range of climate forcing and the identification of parts of the world that are parameter sensitive. Certain simulations reproduced absolute δ18Opvalues across all time periods, along with LGM and MH δ18Opanomalies relative to the PI, better than the default parameterization. No single set of parameterizations worked well in all climate states, likely due to the non‐stationarity of cloud feedbacks under varying boundary conditions. Future work that involves varying multiple parameter sets simultaneously with coupled ocean feedbacks will likely provide improved constraints on cloud and convective parameterizations.more » « less
-
This paper investigates a control strategy in which the state of a dynamical system is driven slowly along a trajectory of stable equilibria. This trajectory is a continuum set of points in the state space, each one representing a stable equilibrium of the system under some constant control input. Along the continuous trajectory of such constant control inputs, a slowly varying control is then applied to the system, aimed to create a stable quasistatic equilibrium that slowly moves along the trajectory of equilibria. As a stable equilibrium attracts the state of system within its vicinity, by moving the equilibrium slowly along the trajectory of equilibria, the state of system travels near this trajectory alongside the equilibrium. Despite the disadvantage of being slow, this control strategy is attractive for certain applications, as it can be implemented based only on partial knowledge of the system dynamics. This feature is in particular important for the complex systems for which detailed dynamical models are not available.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- The DOI is not currently available.
