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 review some known results for POD model reduction applied to ODEs. Then, these
results are generalized to several types of DAEs. We provide algorithms for the
model reduction and error bounds for the reduced order models. Some limits of
the approach are pointed out and alternative methods for reduced order subspace approximation
are presented. The POD approach is tested and evaluated for a medium sized DAE example
from multibody dynamics.
We consider a semilinear parabolic equation subject to a nonlinear dynamical boundary condition that is related to the so-calles Wentzell boundary condition. First, we prove the existence and uniqueness of global solutions as well as the existence of a global attractor. Then we derive a suitable Lojasiewicz-Simon-type inequality to show the convergence of global solutions to single steady states as time tends to infinity under the assumption that the nonlinear terms $f$, $g$ are real analytic. Moreover, we provide an estimate for the convergence rate.
A nonlocal quasilinear multi-phase system with nonconstant specific heat and heat conductivity
(2012)
In this paper, we prove the existence
and global boundedness from above for a solution to an
integrodifferential model for nonisothermal multi-phase
transitions under nonhomogeneous third type boundary conditions.
The system couples a quasilinear internal energy balance
ruling the evolution of the absolute temperature with a vectorial
integro-differential inclusion governing the vectorial
phase-parameter dynamics. The specific heat and the heat
conductivity k are allowed to depend both on the order parameter
$\chi$ and on the absolute temperature $\teta$ of the system, and
the convex component of the free energy may or may not be
singular. Uniqueness and continuous data dependence are
also proved under additional assumptions.
We prove the existence, uniqueness, thermodynamic consistency,
global boundedness from both above and below, and continuous data
dependence for a strong solution to an
integrodifferential model for nonisothermal phase transitions
under nonhomogeneous mixed boundary conditions.
The specific heat is allowed to depend on the order parameter,
and the convex component of the free energy may or may not
be singular.
For a nice Markov process such as Brownian motion on a bounded domain, we introduce a non-linear potential operator defined in terms of running
suprema, and we prove a non-linear Riesz representation of a given function as
the sum of a harmonic function and a non-linear potential. The proof involves
a family of optimal stopping problems in analogy to the general construction
of Bank and El Karoui [3], but here the analysis is carried out in terms of
probabilistic potential theory.
The numerical simulation of the band structure of three-dimensional dispersive metallic photonic crystals with face-centered cubic lattices leads to large-scale nonlinear eigenvalue problems, which are very challenging due to a high dimensional subspace associated with the eigenvalue zero and the fact that the desired eigenvalues (with smallest real part) cluster near the zero eigenvalues. For
the solution of the eigenvalue problem, a Newton-type iterative method is proposed and the nullspace-free method is applied to exclude the zero eigenvalues from the associated generalized eigenvalue problem. To find the successive eigenvalue/eigenvector pairs, we propose a new non-equivalence deflation method to transform converged eigenvalues to infinity, while all other eigenvalues remain unchanged. The deflated problem is then solved by the same Newton-type method, which uses a hybrid method that combines the Jacobi-Davidson, the shift-invert residual Arnoldi and nonlinear Arnoldi methods to compute the clustered eigenvalues. Numerical results illustrate that the method is robust even for the case of computing many eigenvalues in very large problems.
In this paper we propose a new finite element realization of the Perfectly Matched
Layer method (PML-method). Our approach allows to deal with a wide class of
polygonal domains and with certain types of inhomogeneous exterior domains.
Among the covered inhomogeneities are open waveguide structures playing an essential
role in integrated optics. We give a detailed insight into implementation
aspects. Numerical examples show exponential convergence behavior to the exact
solution with the thickness of the PML sponge layer.
Modelling incompressible ideal fluids as a finite collection of
vortex filaments is important in physics (super-fluidity, models for the
onset of turbulence) as well as for numerical algorithms used in computer
graphics for the real time simulation of smoke. Here we introduce
a time-discrete evolution equation for arbitrary closed polygons in 3-
space that is a discretisation of the localised induction approximation of
filament motion. This discretisation shares with its continuum limit the
property that it is a completely integrable system. We apply this polygon
evolution to a significant improvement of the numerical algorithms
used in Computer Graphics.
A generalization of the method of Chu, Liu and Mehrmann
for the computation of the Hamiltonian real Schur form is presented.
The new method avoids some of the difficulties that may arise when
a Hamiltonian matrix has tightly clustered groups of eigenvalues.
A detailed analysis of the method is presented and several numerical examples demonstrate the superior behavior of the method.
O&D revenue management (RM)
– either leg-based or PNR-based – has become
a standard in the airline industry. In this paper,
we present a new approach to O&D RM which
does not make any assumptions on demand distributions
or on the correlations of the booking
process. Protection levels are determined for all
origin destination itineraries, fare classes, points
of sale and data collection points (DCPs). This
approach to the seat inventory problem is modelled
as a multistage stochastic program, where
its stages correspond to the DCPs of the booking
horizon. The stochastic passenger demand
process is approximated by a scenario tree generated
from historical data by a recursive scenario
reduction procedure. The stochastic program
represents a specially structured large scale
LP that may be solved by standard LP software
(e.g. CPLEX). Preliminary numerical experience
is reported.
Flux coupling analysis is a method to identify blocked and coupled reactions
in a metabolic network at steady state. We present a new approach to
flux coupling analysis, which uses a minimum set of generators of the steady state
flux cone. Our method does not require to reconfigure the network by splitting
reversible reactions into forward and backward reactions.
By distinguishing different types of reactions (irreversible, pseudo-irreversible,
fully reversible), we show that reaction coupling relationships can only hold
between certain reaction types. Based on this mathematical analysis, we propose
a new algorithm for flux coupling analysis.
We propose a new approach to competitive analysis by introducing the novel concept of online approximation schemes. Such scheme algorithmically constructs an online algorithm with a competitive ratio arbitrarily close to the best possible competitive ratio for any online algorithm. We study the problem of scheduling jobs online to minimize the weighted sum of completion times on parallel, related, and unrelated machines, and we derive both deterministic and randomized algorithms which are almost best possible among all online algorithms of the respective settings. Our method relies on an abstract characterization of online algorithms combined with various simplifications and transformations. We also contribute algorithmic means to compute the actual value of the best possible competitive ratio up to an arbitrary accuracy. This strongly contrasts all previous manually obtained competitiveness results for algorithms and, most importantly, it reduces the search for the optimal competitive ratio to a question that a computer can answer. We believe that our method can also be applied to many other problems and yields a completely new and interesting view on online algorithms.
Dominance constraints are logical descriptions of trees. Efficient algorithms for the subclass of normal dominance constraints were recently proposed. We present a new and simpler graph algorithm solving these constraints more efficiently, in quadratic time per solved form. It also applies to weakly normal dominance constraints as needed for an application to computational linguistics. Subquadratic running time can be achieved employing decremental graph biconnectivity algorithms.
A novel Finite Element Method (FEM) for the computational simulation in particle reinforced composite materials with many inclusions is presented. It is based on an adapted mesh which consists of triangles and parametric quadrilaterals in 2D. The number of elements and, hence, the number of degrees of freedom are proportional to the number of inclusions. The error of the method is independent of the distance of the neighboring inclusions. While being related to network methods, the approach can tackle more general settings. We present an efficient residual a posteriori error estimator which enables to compute reliable upper and lower error bounds. Several numerical examples illustrate the performance of the method and the error estimator. Moreover, it is demonstrated that the assumption of a lattice structure of inclusions can easily lead to incorrect predictions about material properties.
We consider the linear-quadratic optimal control problem for a controlled differential
algebraic equation (DAE). Under minimal assumptions on the DAE concerning index and regularity
it will be possible to prove that the sufficient optimality condition given in papers of G. Kurina
and R. März is also a necessary condition. This condition includes the solution of an appropriate
boundary value problem and, in the special case of an explicit ordinary differential equation, the
condition is equal to the well known necessary and sufficient condition of the classical linear-quadratic
optimal control problem.
The purpose of the paper is to apply monotone multigrid methods
to static and dynamic biomechanical contact problems.
In space, a finite element method involving a mortar
discretization of the contact conditions is used.
In time, a new contact--stabilized Newmark scheme is presented.
Numerical experiments for a two body Hertzian contact problem
and a biomechanical knee problem are reported.
In this paper we discuss the numerical
solution of projected generalized Lyapunov
equations using the matrix sign function
method. Such equations arise in stability
analysis and control problems for
descriptor systems including model reduction
based on balanced truncation. It is known
that the matrix sign function method applied
to a matrix pencil $\lambda E-A$ converges
if and only if $\lambda E-A$ is of index at
most two. The convergence is quadratic if
$E$ is nonsingular, and it is linear,
otherwise. We will propose a modification
of the matrix sign function method that
converges quadratically for pencils of
arbitrary index. Numerical examples will be
presented to demonstrate the properties of
the modified method.