385 research outputs found
Nonlinear cross Gramians and gradient systems
We study the notion of cross Gramians for non-linear gradient systems, using the characterization in terms of prolongation and gradient extension associated to the system. The cross Gramian is given for the variational system associated to the original nonlinear gradient system. We obtain linearization results that precisely correspond to the notion of a cross Gramian for symmetric linear systems. Furthermore, first steps towards relations with the singular value functions of the nonlinear Hankel operator are studied and yield promising results.
emgr - The Empirical Gramian Framework
System Gramian matrices are a well-known encoding for properties of
input-output systems such as controllability, observability or minimality.
These so-called system Gramians were developed in linear system theory for
applications such as model order reduction of control systems. Empirical
Gramian are an extension to the system Gramians for parametric and nonlinear
systems as well as a data-driven method of computation. The empirical Gramian
framework - emgr - implements the empirical Gramians in a uniform and
configurable manner, with applications such as Gramian-based (nonlinear) model
reduction, decentralized control, sensitivity analysis, parameter
identification and combined state and parameter reduction
Model Reduction for Complex Hyperbolic Networks
We recently introduced the joint gramian for combined state and parameter
reduction [C. Himpe and M. Ohlberger. Cross-Gramian Based Combined State and
Parameter Reduction for Large-Scale Control Systems. arXiv:1302.0634, 2013],
which is applied in this work to reduce a parametrized linear time-varying
control system modeling a hyperbolic network. The reduction encompasses the
dimension of nodes and parameters of the underlying control system. Networks
with a hyperbolic structure have many applications as models for large-scale
systems. A prominent example is the brain, for which a network structure of the
various regions is often assumed to model propagation of information. Networks
with many nodes, and parametrized, uncertain or even unknown connectivity
require many and individually computationally costly simulations. The presented
model order reduction enables vast simulations of surrogate networks exhibiting
almost the same dynamics with a small error compared to full order model.Comment: preprin
Reduced order models for control of fluids using the Eigensystem Realization Algorithm
In feedback flow control, one of the challenges is to develop mathematical
models that describe the fluid physics relevant to the task at hand, while
neglecting irrelevant details of the flow in order to remain computationally
tractable. A number of techniques are presently used to develop such
reduced-order models, such as proper orthogonal decomposition (POD), and
approximate snapshot-based balanced truncation, also known as balanced POD.
Each method has its strengths and weaknesses: for instance, POD models can
behave unpredictably and perform poorly, but they can be computed directly from
experimental data; approximate balanced truncation often produces vastly
superior models to POD, but requires data from adjoint simulations, and thus
cannot be applied to experimental data.
In this paper, we show that using the Eigensystem Realization Algorithm (ERA)
\citep{JuPa-85}, one can theoretically obtain exactly the same reduced order
models as by balanced POD. Moreover, the models can be obtained directly from
experimental data, without the use of adjoint information. The algorithm can
also substantially improve computational efficiency when forming reduced-order
models from simulation data. If adjoint information is available, then balanced
POD has some advantages over ERA: for instance, it produces modes that are
useful for multiple purposes, and the method has been generalized to unstable
systems. We also present a modified ERA procedure that produces modes without
adjoint information, but for this procedure, the resulting models are not
balanced, and do not perform as well in examples. We present a detailed
comparison of the methods, and illustrate them on an example of the flow past
an inclined flat plate at a low Reynolds number.Comment: 22 pages, 7 figure
Feedback control of unstable steady states of flow past a flat plate using reduced-order estimators
We present an estimator-based control design procedure for flow control,
using reduced-order models of the governing equations, linearized about a
possibly unstable steady state. The reduced models are obtained using an
approximate balanced truncation method that retains the most controllable and
observable modes of the system. The original method is valid only for stable
linear systems, and we present an extension to unstable linear systems. The
dynamics on the unstable subspace are represented by projecting the original
equations onto the global unstable eigenmodes, assumed to be small in number. A
snapshot-based algorithm is developed, using approximate balanced truncation,
for obtaining a reduced-order model of the dynamics on the stable subspace. The
proposed algorithm is used to study feedback control of 2-D flow over a flat
plate at a low Reynolds number and at large angles of attack, where the natural
flow is vortex shedding, though there also exists an unstable steady state. For
control design, we derive reduced-order models valid in the neighborhood of
this unstable steady state. The actuation is modeled as a localized body force
near the leading edge of the flat plate, and the sensors are two velocity
measurements in the near-wake of the plate. A reduced-order Kalman filter is
developed based on these models and is shown to accurately reconstruct the flow
field from the sensor measurements, and the resulting estimator-based control
is shown to stabilize the unstable steady state. For small perturbations of the
steady state, the model accurately predicts the response of the full
simulation. Furthermore, the resulting controller is even able to suppress the
stable periodic vortex shedding, where the nonlinear effects are strong, thus
implying a large domain of attraction of the stabilized steady state.Comment: 36 pages, 17 figure
Model reduction of controlled Fokker--Planck and Liouville-von Neumann equations
Model reduction methods for bilinear control systems are compared by means of
practical examples of Liouville-von Neumann and Fokker--Planck type. Methods
based on balancing generalized system Gramians and on minimizing an H2-type
cost functional are considered. The focus is on the numerical implementation
and a thorough comparison of the methods. Structure and stability preservation
are investigated, and the competitiveness of the approaches is shown for
practically relevant, large-scale examples
- …