Refine
Year of publication
- 2010 (70) (remove)
Language
- English (70)
Keywords
- dissipativity (3)
- optimal control (3)
- backward error (2)
- behavior approach (2)
- even pencil (2)
- even polynomial (2)
- finite elements (2)
- linear quadratic optimal control (2)
- problem (2)
- structured eigenvalue problem (2)
This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study.
We propose a generalization of the Structured Doubling Algorithm (SDA) to compute invariant subspaces
of structured matrix pencils
that arise in the context of solving linear quadratic optimal control problems.
The new algorithm is
designed to attain better accuracy when the classical Riccati equation approach for the solution of the optimal control problem is not well suited because
the stable and unstable invariant subspaces are not well separated (due to eigenvalues near or on the imaginary
axis) or in the case when the Riccati solution does not exist at all. We analyze the convergence
of the method and compare the new method with the classical SDA algorithm as well as some recent structured QR-methods.
Recently, the format of TT tensors
\cite{hackbuschHT,osele1,tyrtosele2,tyrtosele3} has turned out to be
a promising new format for the approximation of solutions of high
dimensional problems. In this paper, we prove some new results for
the TT representation of a tensor $U \in \R^{n_1\times \ldots\times
n_d}$ and for the manifold of tensors of TT-rank $\underline{r}$.\As a first result, we prove that the TT (or compression) ranks $r_i$
of a tensor $U$ are unique and equal to the respective separation
ranks of $U$ if the components of the TT decomposition are required to
fulfil a certain maximal rank condition. We then show that the set
$\mathcal{T}$ of TT tensors of fixed rank $\underline{r}$ forms an embedded
manifold in $\R^{n^d}$, therefore preserving the essential theoretical
properties of the Tucker format, but often showing an improved scaling
behaviour. Extending a similar approach for matrices \cite{conte_lub},
we introduce certain gauge conditions to obtain a unique
representation of the tangent space $\cT_U\mathcal{T}$ of $\mathcal{T}$
and deduce a
local parametrization of the TT manifold. The parametrisation of
$\cT_{U}\mathcal{T}$ is often crucial for an algorithmic treatment of
high-dimensional time-dependent PDEs and minimisation problems
\cite{lubuch_blau}. We conclude with remarks on those applications and
present some numerical examples.
We consider the numerical solution of projected Lyapunov equations using Krylov subspace iterative methods. Such equations
play a fundamental role in balanced truncation model reduction of descriptor systems. We present generalizations of
the extended block and global Arnoldi methods to projected Lyapunov equations and compare these methods with the alternating direction implicit method with respect to performance on different examples.
A deflation strategy is also proposed to overcome possible breakdown in the
recurrence.
The PSurface Library
(2010)
We describe psurface, a C++ library that allows to store and access piecewise linear mappings between simplicial surfaces in $\R^2$ and $\R^3$. These mappings are stored in a graph data structure and can be constructed explicitly, by projection, or by surface simplification. Piecewise linear maps can be used, e.g., to construct boundary
approximations for finite element grids, and grid intersections for domain decomposition methods. In computer graphics the mappings allow to build level-of-detail representations as well as texture- and bump maps. We document the data structures and algorithms used and show how \psurface is used in the numerical analysis framework Dune
and the visualization software Amira.
We describe the appearance and stability of spatio-temporal periodic
patterns (rotating waves) in unidirectional rings of coupled oscillators
with delayed couplings. We show how delays in the coupling lead
to a splitting of each rotating wave into several new ones. The appearance
of rotating waves is mediated by Hopf bifurcations of the symmetric
equilibrium.
We also conclude that the coupling delays can be effectively
replaced by increasing the number of oscillators in the chain.
The phenomena are shown for Stuart-Landau
oscillators as well as for coupled FitzHugh-Nagumo systems interacting
via excitatory chemical synapses.
Classical surface parameterization algorithms often place singularities
in order to enhance the quality of the resulting parameter map. Unfortunately, singularities of positive integral index (as the north pole of a sphere) were not handled since they cannot be described with piecewise linear parameter functions on a triangle mesh. Preprocessing is needed to adapt the mesh connectivity. We present an extension to the QuadCover parameterization algorithm [KNP07], which allows to handle those singularities. A singularity of positive integral index can be resolved using bilinear parameter functions on quadrilateral elements. This generalization
of piecewise linear functions for quadrilaterals enriches the space of parameterizations. The resulting parameter map can be visualized by textures using a rendering system which supports quadrilateral elements, or it can be used for remeshing into a pure quad mesh.
We study the class of linear differential-algebraic m-input m-output systems
which have a transfer function with proper inverse.
A sufficient condition for the transfer function to have proper inverse
it that the system has 'strict and non-positive relative degree'.
We present two main results:
First, a so called 'zero dynamics form' is derived: this form is - within the class of system equivalence - a simple "almost normal" form of the DAE; it is a counterpart to the well-known Byrnes-Isidori form for
ODE systems with strictly proper transfer function.
The 'zero dynamics form' is exploited to characterize structural properties such as
asymptotically stable zero dynamics,
minimum phase, and high-gain stabilizability.
The zero dynamics are characterized by (A,E,B)-invariant subspaces.
Secondly, it is shown that the 'funnel controller' (that is a static nonlinear output error feedback) achieves, for all DAE systems with asymptotically stable zero dynamics and transfer function with proper inverse, tracking of a reference signal by the output signal within a pre-specified funnel. This funnel determines the transient behaviour.
We study discrete curvatures computed from nets of curvature lines on a given smooth surface and prove their uniform convergence to smooth principal curvatures. We provide explicit error bounds, with constants depending only on properties of the smooth limit surface and the shape regularity of the discrete net.
We present a novel algorithm for automatic parameterization of tube-like surfaces of arbitrary genus such as the surfaces of knots, trees, blood vessels, neurons, or any tubular graph with a globally consistent stripe texture. We use the principal curvature frame field of the underlying tube-like surface to guide the creation of a global, topologically consistent stripe parameterization of the surface. Our algorithm extends the QuadCover algorithm and is based, first, on the use of so-called projective vector fields instead of frame fields, and second, on different types of branch points. That does not only simplify the mathematical theory, but also reduces computation time by the decomposition of the underlying stiffness matrices.