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.