Refine
Year of publication
Language
- English (1103)
- German (10)
- Multiple languages (1)
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)
We prove an existence and uniqueness theorem for weighted Delaunay triangulations (with non-intersecting site-circles) with prescribed combinatorial type and circle intersection angles. Such weighted Delaunay triangulations can also be interpreted as hyperbolic polyhedra with vertices beyond the infinite boundary. The proof is based on a variational principle. This extends similar work by Rivin on Delaunay triangulations and ideal polyhedra to weighted Delaunay triangulations and hyperideal polyhedra.
The slow processes of metastable stochastic dynamical systems are difficult to access by direct numerical simulation due the sampling problem. Here, we suggest an approach for modeling the slow parts of Markov processes by approximating the dominant eigenfunctions and eigenvalues of the propagator. To this end, a variational principle is derived that is based on the maximization of a Raleigh coefficient. It is shown that this Raleigh coefficient can be estimated from statistical observables that can be obtained from short distributed simulations starting from different parts of state space. The approach forms a basis for the development of adaptive and efficient computational algorithms for simulating and analyzing metastable Markov processes while avoiding the sampling problem. Since any stochastic process with finite memory can be transformed into a Markov process, the approach is applicable to a wide range of processes relevant for modeling complex real-world phenomena.
We analyze a rate-independent model for damage evolution in elastic bodies. The central quantities are a stored energy functional and a dissipation functional, which is assumed to be positively homogeneous of degree one. Since the energy is not simultaneously (strictly) convex in the damage variable and the displacements, solutions may have jumps as a function of time. The latter circumstance makes it necessary to recur to suitable notions of weak solution. However, the by-now classical concept of global energetic solution fails to describe accurately the behavior of the system at jumps.
Hence, we consider rate-independent damage models as limits of systems driven by viscous, rate-dependent dissipation. We use a technique for taking the vanishing viscosity limit, which is based on arc-length reparameterization. In this way, in the limit we obtain a novel formulation for the rate-independent damage model, which highlights the interplay of viscous and rate-independent effects in the jump regime, and provides a better description of the energetic behavior of the system at jumps.
We are concerned with a nonstandard phase field model of
Cahn-Hilliard type. The model, which was introduced by Podio-Guidugli (Ric. Mat. 2006), describes two-species phase segregation and consists of a system of two highly nonlinearly coupled PDEs. It has been recently investigated
by Colli, Gilardi, Podio-Guidugli, and Sprekels in a series of papers: see, in particular, SIAM J. Appl. Math. 2011,
and Boll. Unione Mat. Ital. 2012. In the latter contribution, the authors can treat the very general case in which the diffusivity coefficient of the parabolic PDE
is allowed to depend nonlinearly on both variables. In the same framework, this paper investigates the asymptotic limit of the solutions to the initial-boundary value problems as the diffusion coefficient sigma in the equation governing the evolution of the order parameter tends to zero. We prove that such a limit actually exists and solves the limit problem, which couples a nonlinear PDE of parabolic type with an ODE accounting for the phase dynamics. In the case of a constant diffusivity, we are able to show uniqueness and to improve the regularity of the solution.
It is known that for each combinatorial type of convex 3-dimensional
polyhedra, there is a representative with edges tangent to the unit sphere.
This representative is unique up to projective transformations that fix the unit
sphere. We show that there is a unique representative (up to congruence) with
edges tangent to the unit sphere such that the origin is the barycenter of the
points where the edges touch the sphere.
The subject of the present paper is a simplified model for a symmetric bistable system with memory
or delay, the reference model, which in the presence of noise exhibits a phenomenon similar to what
is known as stochastic resonance. The reference model is given by a one dimensional parametrized
stochastic differential equation with point delay, basic properties whereof we check.
With a view to capturing the effective dynamics and, in particular, the resonance-like behavior of
the reference model we construct a simplified or reduced model, the two state model, first in discrete
time, then in the limit of discrete time tending to continuous time. The main advantage of the
reduced model is that it enables us to explicitly calculate the distribution of residence times which
in turn can be used to characterize the phenomenon of noise-induced resonance.
Drawing on what has been proposed in the physics literature, we outline a heuristic method for
establishing the link between the two state model and the reference model. The resonance characteristics
developed for the reduced model can thus be applied to the original model.
We develop a model for the dynamic evolution of default-free and defaultable interest rates in a LIBOR framework. Utilizing the class of affine processes, this model produces positive LIBOR rates and spreads, while the dynamics are analytically tractable under defaultable forward measures. This leads to explicit formulas for CDS spreads, while semi-analytical formulas are derived for other credit derivatives. Finally, we give an application to counterparty risk.
In this paper we give an overview of model
order reduction techniques for coupled
systems. We consider linear time-invariant
control systems that are coupled through
input-output relations and discuss model
reduction of such systems using moment
matching and balanced truncation.
Structure-preserving approaches to model
order reduction of coupled systems are also
presented. Numerical examples are given.
Supercomputers can simulate complex molecular systems. However, there is a very large gap between the fastest oscillations of covalent bonds of a molecule and the time-scale of the dominant processes. In order to extract the dominant time-scales and to identify the dominant processes, a clustering of information is needed. This thesis shows that only the subspace-based Robust Perron Cluster Analysis (PCCA+) can solve this problem correctly by the construction of a Markov State Model. PCCA+ allows for time-extrapolation in molecular kinetics. This thesis shows the difference between molecular dynamics and molecular kinetics. Only in the molecular kinetics framework a definition of transition rates is possible. In this context, the existence of an infinitesimal generator of the dynamical processes is discussed. If the existence is assumed, the Theorem of Gauß can be applied in order to compute transition rates efficiently. Molecular dynamics, however, is not able to provide a suitable statistical basis for the determination of the transition pattern.
We present structure preserving algorithms for the numerical com-
putation of structured staircase forms of skew-symmetric/symmetric
matrix pencils along with the Kronecker indices of the associated skew-
symmetric/symmetric Kronecker-like canonical form. These methods
allow deflation of the singular structure and deflation of infinite eigenvalues with index greater than one. Two algorithms are proposed: one
for general skew-symmetric/symmetric pencils and one for pencils in
0
which the skew-symmetric matrix is a direct sum of 0 and J = −I I .
0
We show how to use the structured staircase form to solve boundary
value problems arising in control applications and present numerical
examples.
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.
We propose a robust and efficient numerical discretization scheme for the infinitesimal generator of a diffusion process based on a finite volume approximation. The resulting discrete-space operator can be interpreted as a jump process on the mesh whose invariant measure is precisely the cell approximation of the Boltzmann distribution of the original process. Moreover the resulting jump process preserves the detailed balance property of the original stochastic process.
A Structure-Preserving Method for Generalized Algebraic RiccatiEquations Based on Pencil Arithmetic
(2004)
This paper describes a numerical method for extracting the stable
right deflating subspace of a matrix pencil Z Y using
a spectral projection method. It has several advantages compared
to other spectral projection methods like the sign function
method. In particular it avoids the rounding error induced
loss of accuracy associated with matrix inversions. The new algorithm
is particularly well adapted to solving continuous-time
algebraic Riccati equations. In numerical examples, it solves
Riccati equations to high accuracy.
In this paper we propose a Libor model with a high-dimensional specially structured system of
driving CIR volatility processes. A stable calibration procedure which takes into account
a given local correlation structure is presented. The calibration algorithm is FFT based, so fast and easy
to implement.
Trajectory- or mesh-based methods for analyzing the dynamical behavior of large
molecules tend to be impractical due to the curse of dimensionality - their computational cost increases
exponentially with the size of the molecule. We propose a method to break the curse by a
novel square root approximation of transition rates, Monte Carlo quadrature and a discretization
approach based on solving linear programs. With randomly sampled points on the molecular energy
landscape and randomly generated discretizations of the molecular conguration space as our initial
data, we construct a matrix describing the transition rates between adjacent discretization regions.
This transition rate matrix yields a Markov State Model of the molecular dynamics. We use Perron
cluster analysis and coarse-graining techniques in order to identify metastable sets in conguration
space and approximate the transition rates between the metastable sets. Application of our method
to a simple energy landscape on a two-dimensional conguration space provides proof of concept and
an example for which we compare the performance of dierent discretizations. We show that the
computational cost of our method grows only polynomially with the size of the molecule. However,
nding discretizations of higher-dimensional conguration spaces in which metastable sets can be
identied remains a challenge.
Laplace transforms which admit a holomorphic extension to some sector strictly
containing the right half plane and exhibiting a potential behavior are considered. A spectral order,
parallelizable method for their numerical inversion is proposed. The method takes into account the
available information about the errors arising in the evaluations. Several numerical illustrations are
provided.
We study the expansion of the eigenfunctions of Schrödinger operators with smooth confinement potentials in Hermite functions; confinement potentials are potentials that become unbounded at infinity. The key result is that such eigenfunctions and all their derivatives decay more rapidly than any exponential function under some mild growth
conditions to the potential and its derivatives. Their expansion in Hermite functions converges therefore very fast, super-algebraically.
Short term climate events such as the sea surface temperature anomaly known as El Nino are financial risk sources leading to incomplete markets. To make such risk tradable, we use a market model in which a climate index provides an extra investment opinion. Given one possible market price of risk each agent can maximize the exponential utility from three sources of income: capital market, additional security, and individual risk exposure. Under an equilibrium condition the market price of risk is uniquely determined by a backward stochastic differential equation. We translate these stochastic equations into semi-linear partial differential equations for the simulation of which numerical schemes are available. We choose two simple models for sea surface temperature, and with ENSO risk exposed fisher and farmer and a nonh-exposed bank three toy agents. By simulating their optimal investment into the climat index we obtain first insight into the dynamics of the market.
Bovine fertility is the subject of extensive research in animal sciences,
especially because fertility of dairy cows has declined during the last
decades. The regulation of estrus is controlled by the complex interplay
of various organs and hormones. Mathematical modeling of the bovine
estrous cycle could help in understanding the dynamics of this complex
biological system. In this paper we present a mechanistic mathematical
model of the bovine estrous cycle that includes the processes of follicle
and corpus luteum development and the key hormones that interact to
control these processes. The model generates successive estrous cycles of
21 days, with three waves of follicle growth per cycle. The model contains
12 differential equations and 54 parameters. Focus in this paper is on
development of the model, but also some simulation results are presented,
showing that a set of equations and parameters is obtained that describes
the system consistent with empirical knowledge. Even though the majority
of the mechanisms that are included in the model are based on relations
that in literature have only been described qualitatively (i.e. stimulation
and inhibition), the output of the model is surprisingly well in line with
empirical data. This model of the bovine estrous cycle could be used
as a basis for more elaborate models with the ability to study effects of
external manipulations and genetic differences.
We present a linear time approximation algorithm with a performance ratio of 1/2 for finding a maximum weight matching in an arbitrary graph. Such a result is already known and is due to Preis [STACS'99, Lecture Notes in Comput. Sci., Vol. 1563, 1999, pp. 259–269]. Our algorithm uses a new approach which is much simpler than the one given by Preis and needs no amortized analysis for its running time.