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 presents three different adaptive algorithms for eigenvalue problems
associated with non-selfadjoint partial differential
operators. The basis for the developed algorithms is a homotopy method.
The homotopy method starts from a well-understood selfadjoint problem,
for which well-established adaptive methods are available.
Apart from the adaptive grid refinement, the progress of the homotopy as
well as the solution of the iterative method are adapted to balance the contributions
of the different error sources.
The first algorithm balances the homotopy, discretization and approximation errors with respect
to a fixed step-size $\tau$ in the homotopy.
The second algorithm combines the adaptive step-size control for the homotopy
with an adaptation in space that ensures an error below a fixed tolerance $\varepsilon$.
The third algorithm allows the complete adaptivity in space,
homotopy step-size as well as the iterative algebraic eigenvalue solver.
All three algorithms are compared in numerical examples.
Given a system of closed convex domains (inclusions) in n-dimensional Euclidean space, new computational meshes are introduced which partition the convex hull of the inclusion set into simple geometric objects. These partitions generalize the concept of Delaunay triangulations by interpreting the inclusions as generalized vertices while the remaining elements of the partition serve as connections between generalized vertices and therefore assume the classical role of edges, faces, etc. The proposed partitions are derived in two different ways: by exploiting duality with respect to certain generalized Voronoi partitions and by generalizing the well known Delaunay (empty circumcircle) criterion.
Generalized Delaunay partitions are of practical importance for the modeling of particle- and fiber-reinforced composite materials
since they enable an efficient conforming resolution of the highly complicated component geometries. The number of elements in the partitions is proportional to the number of inclusions which is minimal.
Functional Magnetic Resonance Imaging inherently involves noisy measurements and a severe multiple test
problem. Smoothing is usually used to reduce the effective number of multiple
comparisons and to locally integrate the signal and hence increase the
signal-to-noise ratio. Here, we provide a new structural adaptive segmentation
algorithm (AS)
that naturally combines the signal detection with noise reduction in one procedure.
Moreover, the new method
is closely related to a recently proposed structural adaptive smoothing
algorithm and preserves shape and spatial extent of activation areas without
blurring the borders.
MINLP Solver Software
(2010)
Motivated by an obstacle problem for a membrane
subject to cohesion forces, constrained minimization problems involving
a non-convex and non-differentiable objective functional
representing the total potential energy are considered. The associated
first order optimality system leads to a hemi-variational inequality,
which can also be interpreted as a special complementarity
problem in function space. Besides an analytical investigation
of first-order optimality, a primal-dual active set solver is introduced.
It is associated to a limit case of a semi-smooth Newton
method for a regularized version of the underlying problem class.
For the numerical algorithms studied in this paper, global as well as
local convergence properties are derived and verified numerically.
A duality based semismooth Newton framework for solving variational inequalities of the second kind
(2010)
In an appropriate function space setting, semismooth Newton methods are proposed
for iteratively computing the solution of a rather general class of variational inequalities (VIs) of the
second kind. The Newton scheme is based on the Fenchel dual of the original VI problem which
is regularized if necessary. In the latter case, consistency of the regularization with respect to the
original problem is studied. The application of the general framework to specific model problems
including Bingham flows, simplified friction, or total variation regularization in mathematical imaging
is described in detail. Finally, numerical experiments are presented in order to verify the theoretical
results.
In this paper a general form of the infinite-horizon linear quadratic control problem is considered. We will discuss quadratic cost functionals which involve not only the state and input-variables but also derivatives of the state and input-variables of arbitrary order under constraints given by linear systems of higher order. We will examine two results that relate the linear quadratic control problem to an optimality system, which is given through a para-Hermitian matrix polynomial. The results can be applied to general rectangular descriptor systems (see Subsection 6.1) to obtain results which so far were only known for quadratic descriptor systems. Also we will see that the notion of dissipativity (when introduced in the proper way) is equivalent to the solvability of the linear quadratic control problem.
The behavior approach and the problem of dissipativity have both been introduced and studied extensively by Willems et al. However, a computationally feasible method to check dissipativity is missing. Current methods will mostly rely on symbolic representations of rational functions. We will discuss a new characterization for linear systems in behavior form that allows to check dissipativity via the solution of a para-Hermitian, polynomial eigenvalue problem. Thus, we can employ standard methods of cubic complexity.
Under market frictions like illiquidity or transaction costs, contingent claims
can incorporate some inevitable intrinsic risk that cannot be completely hedged
away but remains with the holder. In general, they cannot be synthesized by
dynamical trading in liquid assets and hence not be priced by no-arbitrage arguments alone. Still, an agent can determine a valuation with respect to her
preferences towards risk. The utility indifference value for a variation in the
quantity of illiquid assets held by the agent is defined as the compensating variation
of wealth, under which her maximal expected utility remains unchanged.
Arrow Debreu Prices
(2010)
Arrow Debreu prices are the prices of ‘atomic’ time and state contingent
claims which deliver one unit of a specific consumption good if a specific uncertain
state realizes at a specific future date. For instance, claims on the good
‘ice cream tomorrow’ are split into different commodities depending whether the
weather will be good or bad, so that good-weather and bad-weather ice cream
tomorrow can be traded separately. Such claims were introduced by K.J. Arrow
and G. Debreu in their work on general equilibrium theory under uncertainty,
to allow agents to exchange state and time contingent claims on goods. Thereby
the general equilibrium problem with uncertainly can be reduced to a conventional
one without uncertainty. In finite state financial models, Arrow-Debreu
securities delivering one unit of the numeraire good can be viewed as natural
atomic building blocks for all other state-time contingent financial claims; their
prices determine a unique arbitrage-free price system.
We present and compare two different approaches to conditional
risk measures. One approach draws from vector space based convex analysis
and presents risk measures as functions on L^p spaces while the other approach
utilizes module based convex analysis where conditional risk measures are defined on L^p type modules. Both approaches utilize general duality theory for
vector valued convex functions in contrast to the current literature in which
we fi nd ad hoc dual representations. By presenting several applications such
as monotone and sub(cash) invariant hulls with corresponding examples we
illustrate that module based convex analysis is well suited to the concept of
conditional risk measures.
In this paper, we propose and investigate numerical methods based on QR factorization for computing all or some Lyapunov or Sacker-Sell spectral intervals for
linear differential-algebraic equations.
Furthermore, a perturbation and error analysis for these methods is presented. We
investigate how errors in the data and in the numerical integration affect the
accuracy of the approximate spectral intervals. Although we need to integrate
numerically some differential-algebraic systems on usually very long
time-intervals, under certain assumptions, it is shown that the error of the
computed spectral intervals can be controlled by the local error of numerical
integration and the error in solving the algebraic constraint.
Some numerical examples are presented to illustrate the theoretical results.
We estimate potential energy savings in IP-over-WDM networks achieved by switching off router line cards in low-demand hours. We compare three approaches to react on dynamics in the IP traffic over time, FUFL,
DUFL and DUDL. They provide different levels of freedom in adjusting the routing of lightpaths in the WDM layer and the routing of demands in the IP layer. Using MILP models based on realistic network topologies and node architectures as well as realistic demands, power, and cost values, we show that already a simple monitoring of the lightpath utilization in order to deactivate empty line cards (FUFL) brings substantial
benefits. The most significant savings, however, are achieved by rerouting traffic in the IP layer (DUFL), which allows emptying and deactivating lightpaths together with the corresponding line cards. A
sophisticated reoptimization of the virtual topologies and the routing in the optical domain for every demand scenario (DUDL) yields nearly no additional profits in the considered networks.
A new implicitly-restarted Krylov subspace method
for real symmetric/skew-symmetric generalized eigenvalue problems
is presented. The new method improves and generalizes the SHIRA method
to the case where the skew symmetric matrix is singular.
It computes a few eigenvalues and eigenvectors of the matrix pencil
close to a given target point. Several applications from control theory are
presented and the properties of the new method are illustrated by benchmark
examples.
In this paper, we investigate the interconversion processes of the major flame retardant - 1,2,5,6,9,10-hexabromocyclododecane (HBCD) - by the means of statistical thermodynamics based on classical force-fields. Three ideas will be presented. First, the application of classical hybrid Monte-Carlo simulations for quantum mechanical processes will be justified. Second, the problem of insufficient convergence properties of hybrid Monte-Carlo methods for the generation of low temperature canonical ensembles will be solved by an interpolation approach. Furthermore, it will be shown how free energy differences can be used for a rate matrix computation. The results of our numerical simulations will be compared to experimental results.
Keywords: Markov process, molecular dynamics, rate matrix
Given a general mixed integer program (MIP), we automatically detect block structures in the constraint matrix together with the coupling by capacity constraints arising from multi-commodity flow formulations. We
identify the underlying graph and generate cutting planes based on cuts in the detected network. Our implementation adds a separator to the branch-and-cut libraries of SCIP and CPLEX. We make use of the complemented mixed integer rounding framework (cMIR) but provide a special purpose aggregation heuristic that exploits the network structure. Our separation scheme speeds-up the computation for a large set of MIPs coming from network design problems
by a factor of two on average.
We study the perturbation theory of structured matrices under
structured rank one perturbations, and then focus on several classes of complex
matrices. Generic Jordan structures of perturbed matrices are identified.
It is shown that the perturbation
behavior of the Jordan structures
is substantially different from the corresponding theory for unstructured generic
rank one perturbations.
An one-dimensional Kohn-Sham system for spin particles is considered which effectively describes semiconductor nanostructures and which is investigated at zero temperature. We prove the existence of solutions and derive a priori estimates. For this purpose we find estimates for eigenvalues of the Schrödinger operator with effective Kohn-Sham potential and obtain $W^1,2$-bounds of the associated particle density operator. Afterwards, compactness and continuity results allow to apply Schauder's fixed point theorem. In case of vanishing exchange-correlation potential uniqueness is shown by monotonicity arguments. Finally, we investigate the behavior of the system if the temperature approaches zero.
We regard drift-diffusion equations for semiconductor devices in Lebesgue spaces. To that end we reformulate the (generalized) van Roosbroeck system as an evolution equation for the potentials to the driving forces of the currents of electrons and holes. This evolution equation falls into a class of quasi-linear parabolic systems which allow unique, local in time solution in certain Lebesgue spaces. In particular, it turns out that the divergence of the electron and hole current is an integrable function. Hence, Gauss' theorem applies, and gives the foundation for space discretization of the equations by means of finite volume schemes. Moreover, the strong differentiability of the electron and hole density in time is constitutive for the implicit time discretization scheme. Finite volume discretization of space, and implicit time discretization are accepted custom in engineering and scientific computing. ---This investigation puts special emphasis on non-smooth spatial domains, mixed boundary conditions, and heterogeneous material compositions, as required in electronic device simulation.