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 the effect of generic but structured low rank perturbations on the Jordan structure and sign
characteristic of matrices that have structure in an indefinite inner product space.
The paper is a follow-up of earlier papers in which the effect of rank one perturbations
was considered. Several results that are in contrast to the case of unstructured low rank
perturbations of general matrices are presented here.
In large-area Organic Light-Emitting Diodes (OLEDs) spatially inhomogeneous luminance at high power due to inhomogeneous current flow and electrothermal feedback can be observed. To describe these self-heating effects in organic semiconductors we present a stationary thermistor model based on the heat equation for the temperature coupled to a p-Laplace-type equation for the electrostatic potential with mixed boundary conditions. The p-Laplacian describes the non-Ohmic electrical behavior of the organic material. Moreover, an Arrhenius-like temperature dependency of the electrical conductivity is considered.
We introduce a finite-volume scheme for the system and discuss its relation to recent network models for OLEDs. In two spatial dimensions we derive a priori estimates for the temperature and the electrostatic potential and prove the existence of a weak solution by Schauder's fixed point theorem.
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.
Several classes of optimal control of electromagnetic fields are considered. Special emphasis is
laid on a non-standard $H$-based formulation of the equations of electromagnetism in multiply connected conductors. By this technique, the Maxwell equations can be solved with reduced computational complexity. While the magnetic field $H$ in the conductor is obtained from an elliptic equation
with the $\curl \sigma^{-1} \curl$ operator, an elliptic equation with the $\div \mu \nabla$ operator is set up for a potential $\psi$ in the isolator.
Both equations are coupled by appropriate interface conditions. In all problems, the
electrical current is controlled in the conducting domain. Several types of control functions are discussed. In particular, the problem of sparse optimal control is investigated in a package of electrical wires. For all problems, the associated sensitivity
analysis is performed.
A survey of methods from numerical linear algebra for linear constant coefficient differential-algebraic equations (DAEs) and descriptor control systems is presented. We discuss numerical methods to check the solvability properties of DAEs as well as index reduction and regularization techniques. For descriptor systems we discuss controllability and observability properties and how these can be checked numerically. These methods are based on staircase forms and derivative arrays, transformed with real orthogonal transformations that are discussed in detail. Then we use the reformulated problems in several control applications for differential-algebraic equations ranging from regular and singular linear-quadratic optimal and robust control to dissipativity checking. We discuss these applications and give a systematic overview over the theory and the numerical solution methods. In particular, we show that all these applications can be treated with a common approach that is based on the computation of eigenvalues and deflating subspaces of even matrix pencils. The unified approach allows to generalize and improve several techniques that are currently in use in systems and control.
The efficient and reliable computation of guided modes in photonic crystal wave-guides is of great importance for designing optical devices. Transparent boundary conditions based on Dirichlet-to-Neumann operators allow for an exact computation of well-confined modes and modes close to the band edge in the sense that no modelling error is introduced. The well-known super-cell method, on the other hand, introduces a modelling error which may become prohibitively large for guided modes that are not well-confined. The Dirichlet-to-Neumann transparent boundary conditions are, however, not applicable for all frequencies as they are not uniquely defined and their computation is unstable for a countable set of frequencies that correspond to so called Dirichlet eigenvalues. In this work we describe how to overcome this theoretical difficulty introducing Robin-to-Robin transparent boundary conditions whose construction do not exhibit those forbidden frequencies. They seem, hence, well suited for an exact and reliable computation of guided modes in photonic crystal wave-guides.
Regularization and Numerical Solution of the Inverse Scattering Problem using Shearlet Frames
(2014)
Regularization techniques for the numerical solution of nonlinear inverse scattering
problems in two space dimensions are discussed. Assuming that the boundary of a scatterer is its most prominent feature, we exploit as model the class of cartoon-like functions.
Since functions in this class are asymptotically optimally sparsely approximated by shearlet frames, we consider shearlets as a means for the regularization in a Tikhonov method.
We examine both directly the nonlinear problem and a linearized problem obtained by
the Born approximation technique. As problem classes we study the acoustic inverse
scattering problem and the electromagnetic inverse scattering problem. We show that
this approach introduces a sparse regularization for the nonlinear setting and we present
a result describing the behavior of the local regularity of a scatterer under linearization,
which shows that the linearization does not affect the sparsity of the problem. The analytical results are illustrated by numerical examples for the acoustic inverse scattering problem that highlight the effectiveness of this approach.
An eigenvalue perturbation theory under rank-one perturbations is developed for classes
of real matrices that are symmetric with respect to a non-degenerate bilinear form,
or Hamiltonian with respect to a non-degenerate skew-symmetric form.
In contrast to the case of complex matrices, the sign characteristic is a crucial feature
of matrices in these classes. The behavior of the sign characteristic under generic
rank-one perturbations is analyzed in each of these two classes of matrices.
Partial results are presented, but some questions remain open. Applications
include boundedness and robust boundedness for solutions of structured systems
of linear differential equations with respect to general perturbations as well
as with respect to structured rank perturbations of the coefficients.
Using a standard first-order optimality condition for nonsmooth optimization problems, a general framework for a descent method is developed. This setting is applied to a class of mathematical programs with equilibrium constraints in function space from which a new algorithm is derived. Global convergence of the algorithm is demonstrated in function space and the results are then illustrated by numerical experiments.
We derive formulae for the second-order subdifferential of polyhedral norms. These formulae are fully explicit in terms of initial data. In a first step we rely on the explicit formula for the coderivative of normal cone mapping to polyhedra. Though being explicit, this formula is quite involved and difficult to apply. Therefore, we derive simple formulae for the 1-norm and - making use of a recently obtained formula for the second-order subdifferential of the maximum function - for the maximum norm.
In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides. Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone. We show the derivation and computation of the group velocity, the group velocity dispersion, and any higher derivative of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows
for the identification of crossings and anti-crossings of dispersion curves, respectively. This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.
Structured eigenvalue backward errors of matrix pencils and polynomials with palindromic structures
(2014)
We derive formulas for the backward error of an approximate eigenvalue of a *-palindromic
matrix polynomial with respect to *-palindromic perturbations. Such formulas are also obtained
for complex T-palindromic pencils and quadratic
polynomials. When the T-palindromic polynomial is real, then we derive the backward error
of a real number considered as an approximate eigenvalue of the matrix polynomial with
respect to real T-palindromic perturbations.
In all cases the corresponding minimal structure preserving perturbations are obtained as well.
The results are illustrated by numerical experiments. These show that there is
significant difference between the backward errors with respect to structure
preserving and arbitrary perturbations in many cases.
A mechanical equilibrium problem for a material consisting of two components with dierent densities is considered. Due to the heterogeneous material densities, the
outer shape of the underlying workpiece can be changed by shifting the interface between
the subdomains. In this paper, the problem is modeled as a shape design problem for optimally compensating unwanted workpiece changes. The associated control variable is
the interface. Regularity results for transmission problems are employed for a rigorous
derivation of suitable first-order optimality conditions based on the speed method. The paper concludes with several numerical results based on a spline approximation of the
interface.
Multi-valued network models can be described by their topology and a set of parameters capturing the effects of the regulators for each component. Dynamics can then be derived and represented as state transition systems.
Different network models may lead to the same transition system, meaning dynamics analysis of a representative model covers a larger class of models. While rather clear in the Boolean case, the properties contributing to this effect become more involved for multi-valued models. We analyse these properties and present a mathematical description of the resulting model equivalence classes.
In this paper a new theorem is formulated
which allows a rigorous proof of the shape differentiability without the usage of the material derivative; the domain expression is automatically obtained and the boundary expression is easy to derive.
Furthermore, the theorem is applied to a cost function which depends on a quasi-linear transmission
problem. Using a Gagliardo penalization the existence of optimal shapes is established.
A coupling of discrete and continuous optimization to solve kinodynamic motion planning problems
(2014)
This paper studies the relationship between the material derivative method, the shape derivative method, the min-max formulation of Correa and Seeger, and the Lagrange method introduced by Cea. A theorem is formulated
which allows a rigorous proof of the shape differentiability without the usage of material derivative;
the domain expression is automatically obtained and the boundary expression is easy to derive.
Furthermore, the theorem is applied to a cost function which depends on a quasi-linear transmission
problem. Using a Gagliardo penalization the existence of optimal shapes is established.
An algorithm to detect collisions between robots moving along given trajectories is presented. The method is a combination of the adaptive
dynamic collision checking developed by Schwarzer et al. and Lin and Canny's algorithm, which computes efficiently the distance between two polyhedra. The resulting algorithm is part of a global model that computes the optimal task assignment, sequencing and kinodynamic
motion planning in a robotic work-cell.
Two optimal control problems for instationary magnetization
processes are considered in 3D spatial domains that
include electrically conducting and nonconducting regions. The magnetic
fields are generated by induction coils. In the first model, the induction coil
is considered as part of the conducting region and the electrical current is taken
as control. In the second, the coil is viewed as part of the nonconducting region and the
electrical voltage is the control. Here, an integro-differential equation accounts
for the magnetic induction law that couples the given
electrical voltage with the induced electrical current in the
induction coil.
We derive first-order necessary
optimality condition for the optimal controls of both problems. Based on them,
numerical methods of gradient type are applied. Moreover, we report on the application
of model reduction by POD that lead to tremendous savings. Numerical tests are
presented for academic 3D geometries but also for a real-world application.
We study the fundamental problem of scheduling bidirectional traffic across machines arranged on a path. The main feature of the problem is that jobs traveling in the same direction can be scheduled in quick succession on a machine, while jobs in the other direction have to wait for an additional transit time. We show that this tradeoff makes the problem significantly harder than the related flow shop problem, by showing that it is NP-hard even for jobs with identical processing and transit times. We give polynomial algorithms for a single machine and any constant number of machines. In contrast, we show the problem to be NP-hard on a single machine and with identical processing and transit times if some pairs of jobs in different directions are allowed to run on the machine concurrently. We generalize a PTAS of Afrati et al. [1999] for one direction and a single machine to the bidirectional case on any constant number of machines.
The analysis of adaptive finite element methods in practice immediately leads to eigenvalue clusters which requires the simultaneous marking in adaptive finite element methods. A first analysis for multiple eigenvalues of the recent work [Dai, He, Zhou, arXiv Preprint 1210.1846v2] introduces an adaptive method whose marking strategy is based on the element-wise sum of local error estimator contributions for multiple eigenvalues. This paper proves optimality of a practical adaptive algorithm for eigenvalue clusters for the eigenvalues of the Laplace operator in terms of nonlinear approximation classes. All estimates are explicit in the initial mesh-size, the eigenvalues and the cluster width to clarify the dependence of the involved constants.