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 deals with error estimates for space-time discretizations in the context of evolutionary variational inequalities of rate-independent type. After introducing a general abstract evolution problem, we address a fully-discrete approximation and provide a priori error estimates. The application of the abstract theory to a semilinear case is detailed. In particular, we provide explicit space-time convergence rates for the isothermal Souza-Auricchio model for shape-memory alloys.
This paper presents a combined adaptive finite element method with an iterative algebraic eigenvalue solver for the Laplace eigenvalue problem of quasi-optimal computational complexity. The analysis is based on a direct approach for eigenvalue problems and allows the use of higher order conforming finite element spaces with fixed polynomial degree k>0. The optimal adaptive finite element eigenvalue solver (AFEMES) involves a proper termination criterion for the algebraic eigenvalue solver and does not need any coarsening. Numerical evidence illustrates the optimal computational complexity.
On probabilistic constraints induced by rectangular sets and multivariate normal distributions
(2009)
In this paper, we consider optimization problems under probabilistic constraints which are defined by two-sided
inequalities for the underlying normally distributed random vector. As a main step
for an algorithmic solution of such problems, we derive a derivative formula for (normal) probabilities
of rectangles as functions of their lower or upper bounds. This formula allows to reduce the calculus
of such derivatives to the calculus of (normal) probabilities of rectangles themselves thus generalizing a
similar well-known statement for multivariate normal distribution functions. As an application, we consider
a problem from water reservoir management. One of the outcomes of the problem solution is that the
(still frequently encountered) use of simple individual probabilistic can completely fail. In contrast, the
(more difficult) use of joint probabilistic constraints which heavily depends on the derivative formula mentioned
before yields very reasonable and robust solutions over the whole time horizon considered.
Alternating matrix polynomials, that is, polynomials whose coefficients
alternate between symmetric and skew-symmetric matrices,
generalize the notions of even and odd scalar polynomials.
We investigate the Smith forms of alternating matrix polynomials,
showing that each invariant factor is an even or odd scalar polynomial.
Necessary and sufficient conditions
are derived for a given Smith form to be that of an alternating matrix polynomial.
These conditions allow a characterization of the possible Jordan structures
of alternating matrix polynomials,
and also lead to necessary and sufficient conditions
for the existence of structure-preserving strong linearizations.
Most of the results are applicable to singular as well
as regular matrix polynomials.
The paper proposes goal-oriented error estimation and mesh refinement
for optimal control problems with elliptic PDE constraints using the value
of the reduced cost functional as quantity of interest. Error representation,
hierarchical error estimators, and greedy-style error indicators are derived and
compared to their counterparts when using the all-at-once cost functional as
quantity of interest. Finally, the efficiency of the error estimator and generated
meshes are demonstrated on numerical examples.
We introduce geodesic finite elements as a new way to discretize
the nonlinear configuration space of a geometrically exact Cosserat rod.
These geodesic finite elements naturally generalize standard one-dimensional
finite elements to spaces of functions with values in a Riemannian manifold.
For the special orthogonal group, our approach reproduces the
interpolation formulas of [Crisfield/Jelenic:1999].
Geodesic finite elements are
conforming and lead to objective and path-independent problem formulations.
We introduce geodesic finite elements for general Riemannian manifolds,
discuss the relationship between geodesic finite elements and
coefficient vectors, and estimate the interpolation error.
Then we use them to find static equilibria of hyperelastic Cosserat rods.
Using the Riemannian trust-region algorithm of [Absil/Mahony/Sepulchre:2008]
we show numerically that the discretization error depends optimally on
the mesh size.
he generalized Langevin equation is useful for modeling a wide
range of physical processes. Unfortunately its parameters,
especially the memory function, are difficult to determine for
nontrivial processes. In this paper, relations between a
time-discrete generalized Langevin model and discrete multivariate
autoregressive (AR) or autoregressive moving average models (ARMA)
are established. This allows a wide range of discrete linear
methods known from time series analysis to be applied. In
particular, the determination of the memory function {\it via} the
order of the respective AR or ARMA model is addressed. The method
is illustrated on a one-dimensional test system and subsequently
applied to the molecular dynamics of a biomolecule which exhibits
an interesting relationship between the solvent method used, the
molecular conformation and the depth of the memory.
e propose an algorithm for the fast and efficient simulation of polymers represented by
chains of hard spheres. The particles are linked by holonomic bond constraints.
While the motion of the polymers is free (i.e., no collisions occur) the equations
of motion can be easily integrated using a collocation-based partitioned Gauss-Runge-Kutta method.
The method is reversible, symplectic and preserves energy. Moreover the numerical scheme allows the integration using much longer time steps than any explicit integrator such as the popular Verlet method. If polymers collide the point of impact can be determined to arbitrary precision by simple nested intervals. Once the collision point is known the impulsive contribution can be computed analytically. We illustrate our approach by means of a suitable numerical example.
We study Balanced Truncation for stochastic differential equations. In doing so, we adopt ideas from large deviations theory and discuss notions of controllability and obervability for dissipative Hamiltonian systems with degenerate noise term, also known as Langevin equations. For partially-observed Langevin equations, we illustrate model reduction by balanced truncation with an example from molecular dynamics and discuss aspects of structure-preservation.
or stable linear input-output systems, the method of balanced truncation (B.C.~Moore, {\em IEEE Trans. Auto. Contr.} {\bf AC-26}, 17--32, 1981) consists in finding a coordinate transformation such that modes which are least sensitive to the external input (controllability) also give the least output (observability) and therefore can be neglected. A drawback is that projecting the original equations of motion onto the subspace of interest typically fails to preserve the problem's physical structure, e.g., if the original equations are of second-order form or Hamiltonian. For Hamiltonian systems, a natural way of restricting a system to a subspace is by means of constraints, and we show, employing singular perturbation arguments, that balanced truncation can be done in a structure-preserving fashion. The thus obtained reduced Hamiltonian system preserves stability and passivity and satisfies the usual balanced truncation error bound.
We present a formal procedure for structure-preserving model reduction of linear second-order control problems that appear in a variety of physical contexts, e.g., vibromechanical systems or electrical circuit design. Typical balanced truncation methods that project onto the subspace of the largest Hankel singular values fail to preserve the problem's physical structure and may suffer from lack of stability. In this paper we adopt the framework of generalized Hamiltonian systems that covers the class of relevant problems and that allows for a generalization of balanced truncation to second-order problems.
It turns out that the Hamiltonian structure, stability and passivity are preserved if the truncation is done by imposing a holonomic constraint on the system rather than standard Galerkin projection.
The Steiner connectivity problem is a generalization of
the Steiner tree problem. It consists in finding a minimum cost set of
simple paths to connect a subset of nodes in an undirected graph.
We show that important polyhedral and algorithmic results on the
Steiner tree problem carry over to the Steiner connectivity problem,
namely, the Steiner cut and the Steiner partition inequalities, as
well as the associated polynomial time separation algorithms, can be
generalized. Similar to the Steiner tree case, a certain directed
formulation, which is stronger than the natural undirected one,
plays a central role.
For many fundamental cooperative cost sharing games, especially when costs are supermodular, it is known that Moulin mechanisms inevitably suffer from poor budget balance factors. Mehta, Roughgarden, and Sundararajan recently introduced acyclic mechanisms, which achieve a slightly weaker notion of group-strategyproofness, but leave more flexibility to improve upon the approximation guarantees with respect to budget balance and social cost.
In this paper, we provide a very simple but powerful method for turning any rho-approximation algorithm for a combinatorial optimization problem into a rho-budget balanced acyclic mechanism. Hence, we show that there is no gap between the best possible approximation guarantees of full-knowledge approximation algorithms and weakly group-strategyproof cost sharing mechanisms.
The applicability of our method is demonstrated by deriving mechanisms for scheduling and network design problems which beat the best possible budget balance factors of Moulin mechanisms. By elaborating our framework, we provide means to construct weakly group-strategyproof mechanisms with approximate social cost. The mechanisms we develop for completion time scheduling problems perform surprisingly well by achieving the first constant budget balance and social cost factors.
When managing energy or weather related risk often only imperfect hedging instruments are available. In the first part we illustrate problems arising with imperfect hedging by studying a toy model. We consider an airline’s problem with covering income risk due to fluctuating kerosene prices by investing into futures written on heating oil with closely correlated price dynamics. In the second part we outline recent results on exponential utility based cross hedging concepts. They highlight in a generalization of the Black-Scholes delta hedge formula to incomplete markets. Its derivation is based on a purely stochastic approach of utility maximization. It interprets stochastic control problems in the BSDE language, and profits from the power of the stochastic calculus of variations.
We consider backward stochastic differential equations with drivers of quadratic growth (qgBSDE). We prove several statements concerning path regularity and stochastic smoothness of the solution processes of the qgBSDE, in particular we prove an extension of Zhang's path regularity theorem to the quadratic growth setting. We give explicit convergence rates for the difference between the solution of a qgBSDE and its truncation, filling an important gap in numerics for qgBSDE. We give an alternative proof of second order Malliavin differentiability for BSDE with drivers that are Lipschitz continuous (and differentiable), and then derive the same result for qgBSDE.
We show that the spectrum of linear delay differential equations with
large delay splits into two different parts. One part, called the
strong spectrum, converges to isolated points when the delay parameter
tends to infinity. The other part, called the pseudocontinuous spectrum,
accumulates near criticality and converges after rescaling to a set
of spectral curves, called the asymptotic continuous spectrum. We
show that the spectral curves and strong spectral points provide a
complete description of the spectrum for sufficiently large delay
and can be comparatively easily calculated by approximating expressions.
Local existence, uniqueness and smooth dependence for nonsmooth quasilinear parabolic problems
(2009)
We prove local existence, uniqueness, Hölder regularity in space and time, and smooth dependence in Hölder spaces for a general class of quasilinear parabolic initial boundary value problems with nonsmooth data.
As a result the gap between low smoothness of the data, which is typical for
many applications, and high smoothness of the solutions, which is necessary
for the applicability of differential calculus to abstract formulations of the
initial boundary value problems, has been closed. The theory works for any
space dimension, and the nonlinearities are allowed to be nonlocal and to
have any growth. The main tools are new maximal regularity results [19, 20]
in Sobolev–Morrey spaces for linear parabolic initial boundary value problems
with nonsmooth data, linearization techniques and the Implicit Function Theorem.
In this paper we study a certain cardinality constrained packing integer program which is motivated by the problem of dimensioning a cut in a two-layer network. We prove NP-hardness and consider the facial structure of the corresponding polytope. We provide a complete description for the smallest nontrivial case and develop two general classes of facet-defining inequalities. This approach extends the
notion of the well known cutset inequalities to two network layers.
In this paper, we present a model-based optimization approach for the design of multi-layer networks. The proposed framework is based on a series of increasingly abstract models – from a general technical system model to a problem specific mathematical model – which are used in a planning cycle to optimize the multi-layer networks. In a case study we show how central design questions for an IP-over-WDM network architecture can be answered
using this approach. Based on reference networks from the German research project EIBONE, we investigate the influence of various planning parameters on the total design cost. This includes a comparison of point-to-point vs. transparent optical layer architectures, different traffic distributions, and the use of PoS vs. Ethernet interfaces.
We study superhedging of contingent claims with physical delivery in a discrete-time market model with convex transaction costs. Our model extends Kabanov's currency market model by allowing for nonlinear illiquidity effects. We show that an appropriate generalization of Schachermayer's robust no arbitrage condition implies that the set of claims hedgeable with zero cost is closed in probability. Combined with classical techniques of convex analysis, the closedness yields a dual characterization of premium processes that are sufficient to superhedge a given claim process. We also extend the fundamental theorem of asset pricing for general conical models.