Refine
Year of publication
Language
- English (1103) (remove)
Keywords
- optimal control (27)
- stability (14)
- integer programming (11)
- Stochastic programming (9)
- finite elements (9)
- mixed integer programming (9)
- Hamiltonian matrix (8)
- finite element method (8)
- model reduction (8)
- state constraints (8)
This paper discusses the stability of quasi-static paths for a continuous elastic-plastic system with hardening in a one-dimensional (bar) domain. Mathematical formulations, as well as existence and uniqueness results for dynamic and quasi-static problems involving elastic-plastic systems with linear kinematic hardening are recalled in the paper. The concept of stability of quasi-static paths used here is essentially a continuity property of the system dynamic solutions relatively to the quasi-static ones, when (as in Lyapunov stability) the size of initial perturbations is decreased and the rate of application of the forces (which plays the role of the small parameter in singular perturbation problems) is also decreased to zero. The stability of the quasi-static paths of these elastic-plastic systems is the main result proved in the paper.
Mathematical results on existence for viscoelastodynamic problems with unilateral constraints
(2011)
We study a damped wave equation and the evolution of a Kelvin-Voigt material, both problems have unilateral boundary conditions. Under appropriate regularity assumptions on the initial data, both problems possess a weak solution which is obtained as the limit of a sequence of penalized problems; the functional properties of all the traces are precisely identified through Fourier analysis, and this enables us to infer the existence of a strong solution.
Energetic solutions to rate-independent processes are usually constructed via time-incremental minimization problems. In this work we show that all energetic solutions can be approximated by incremental problems if we allow approximate minimizers, where the error in minimization has to be of the order of the time step. Moreover, we study sequences of problems where the energy functionals have a Gamma limit.
A model for the evolution of damage that allows for complete disintegration is addressed. Small strains and a linear response function are assumed. The ``flow rule'' for the damage parameter is rate-independent. The stored energy involves the gradient of the damage variable, which determines an internal length-scale. Quasi-static fully rate-independent evolution is considered as well as rate-dependent evolution including viscous/inertial effects. Illustrative 2-dimensional computer simulations are presented, too.
Quasistatic small-strain plasticity in the limit of small hardening and its numerical approximation
(2011)
The quasistatic rate-independent evolution of the Prager-Ziegler-type model of linearized plasticity with hardening is shown to converge to the rate-independent evolution of the Prandtl-Reuss elastic/perfectly plastic model. Based on the concept of energetic solutions we study the convergence of the solutions in the limit for hardening coefficients converging to 0 by using the abstract method of Gamma-convergence for rate-independent systems. An unconditionally convergent numerical scheme is devised and 2D and 3D numerical experiments are presented. A two-sided energy inequality is a posteriori verified to document experimental convergence rates.
Deflated and augmented Krylov subspace methods: Basic Facts and a Breakdown-free deflated MINRES
(2011)
In this paper we consider deflation and augmentation techniques for accelerating
the convergence of Krylov subspace methods for the solution of nonsingular linear
algebraic systems. The two techniques are conceptually different from
preconditioning. Deflation "removes" certain parts from the operator, while
augmentation adds a subspace to the Krylov subspace. Both approaches have been
used in a variety of methods and settings. For Krylov subspace methods that
satisfy a (Petrov-) Galerkin condition we show that augmentation can in general
be achieved implicitly by projecting the residuals appropriately and correcting
the approximate solutions in a final step. In this context, we analyze known
methods to deflate CG, GMRes and MinRes. Our analysis reveals that the recently
proposed RMinRes method can break down. We show how such breakdowns can be
avoided by choosing a special initial guess, and we derive a breakdown-free
deflated MinRes method. In numerical experiments we study the properties of
different variants of MinRes analyzed in this paper.
We introduce a~numerical method for the numerical solution of the Lur'e matrix equations that arise, for instance, in linear-quadratic infinite time horizon optimal control. The method is based on the characterization of the solutions in terms of deflating subspaces of a suitable even matrix pencil. Via a Cayley transformation, the problem is transformed to the discrete-time case. This leaves us with a symplectic problem with several Jordan blocks of eigenvalue 1 and even size, which arise from the remaining eigenvalues at infinity of the original problem. For the solution of this modified problem, we use the {\em structure-preserving doubling algorithm} (SDA), an iterative scheme for the solution of dense continuous- and discrete-time algebraic Riccati equations. Unlike other iterative schemes, this algorithm converges also when the pencil has eigenvalues on the unit circle, as is the case in our problem. Implementation issues such as the choice of the parameter $\gamma$ in the Cayley transform are discussed. The numerical examples presented confirm the effectiveness of this method.
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 investigate the regularity of the weak solution to elliptic transmission problems that involve two layered anisotropic materials separated by a boundary intersecting interface. Under a compatibility condition for the angle of contact of the two surfaces and the boundary data, we prove the existence of square-integrable second derivatives, and the global Lipschitz continuity of the solution. We show that the second weak derivatives remain integrable to a certain power less than two if the compatibility condition is violated.
While it is well-known that the standard integral operator K of (stationary) diffuse-gray radiation, as it occurs in the radiosity equation, is compact if the domain of radiative interaction is sufficiently regular, we show noncompactness of the operator if the domain is polyhedral. We also show that a stationary operator is never compact when reinterpreted in a transient setting. Moreover, we provide new proofs, which do not use the compactness of K, for 1 being a simple eigenvalue of K for connected enclosures, and for I-(1-e)K being invertible, provided the emissivity e does not vanish identically.
In this paper we will consider elliptic boundary value problems with
oscillatory diffusion coefficient, say A. We will derive regularity
estimates in Sobolev norms which are weighted by certain derivatives of A.
The constants in the regularity estimates then turn out to be independent of
the variations in A.
These regularity results will be employed for the derivation of error
estimates for hp-finite element discretizations which are explicit with
respect to the local variations of the diffusion coefficient.
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.