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)
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.
For regular matrix pencils the distance in norm to the nearest singular pencil
under low rank perturbation is studied. Characterizations of this distance are derived via the Weyl function of the perturbation. Special attention is paid to the Hermitian pencil case.
Estimates for the distance of a given pencil to the set of singular pencils are obtained.
A Branch-Cut-And-Price Algorithm for a Combined Buy-at-Bulk Network Design Facility Location Problem
(2014)
In the combined buy-at-bulk network design facility location
problem we are given an undirected network with a set of potential facilities and a set of clients with demands. In addition, we are also provided a set of cable types with each cable type having a capacity and cost per unit length. The cost to capacity ratio decreases from large to small
cables following economies of scale. A network planner is expected to determine a set of facilities to open and install cables along the edges of the network in order to route the demands of the clients to some open facility. The capacities of the installed cable must be able to support the demands of the clients routed along that edge.The objective is to minimize the cost that is paid for opening facilities and installing cables. We model the problem as an integer program and propose a branch-cut-and-price algorithm for solving it. We study the effect of two family of valid inequalities that naturally emerge from the model. We present the results of our implementation that were tested on a set of large real world instances.
Asymptotic behavior is often of particular interest when analyzing asynchronous Boolean networks representing biological systems such as signal transduction or gene regulatory networks. Methods based on a generalization of the steady state notion, the so-called symbolic steady states, can be exploited to investigate attractor properties as well as for model reduction techniques conserving attractors. In this paper, we propose a novel optimization-based method for computing all maximal symbolic steady states and motivate their use. In particular, we add a new result yielding a lower bound for the number of cyclic attractors and illustrate the methods with a short study of a MAPK pathway model.
Constraint-based modeling of genome-scale metabolic network reconstructions has become a widely used approach in computational biology. Flux coupling analysis is a constraint-based method that analyses the impact of single reaction knockouts on other reactions in the network.We present an extension of flux coupling analysis for double and multiple gene or reaction knockouts, and develop corresponding algorithms for an in silico simulation. To evaluate our method, we perform a full single and double knockout analysis on a selection of genome-scale metabolic network reconstructions and compare the results.
Flux coupling analysis (FCA) has become a useful tool for aiding metabolic reconstructions and guiding genetic manipulations. Originally, it was introduced for constraint-based models of metabolic networks that are based on the steady-state assumption. Recently, we have shown that the steady-state assumption can be replaced by a much weaker lattice-theoretic property related to the supports of metabolic fluxes. In this paper, we further extend our approach and delevelop an efficient algorithm for general qualitative flux coupling analysis (QFCA). We illustrate our method by thermodynamic flux coupling analysis (tFCA), which allows studying steady-state metabolic models with loop-law thermodynamic constraints. These models do not satisfy the lattice-theoretic properties required in our previous work. For a selection of genome-scale metabolic network reconstructions, we discuss both theoretically and practically, how thermodynamic constraints strengthen the coupling results that can be obtained with classical FCA.
Transitive reductions and minimal equivalent subgraphs have proven to be a powerful concept to simplify networks and to measure their redundancy. Here we consider a generalization of the minimal equivalent subgraph problem where a set of arcs is already given. For two digraphs D = (V,A), D′ = (V,A′) with A′ ⊆ A, we ask for the minimal set of edges of D that have to be added to D′ such that the transitive closure of D equals the transitive closure of D′.
We present a method to compute such an extension and show that if D is transitively closed, this problem can be solved in polynomial time.
The long standing problem is discussed of how to deflate the part associated with the eigenvalue infinity in a structured matrix pencil using structure preserving unitary transformations. We derive such a deflation procedure and apply this new technique to symmetric, Hermitian or alternating pencils and in a modified form to (anti)-palindromic pencils. We present a detailed error and perturbation analysis of this and other deflation procedures and demonstrate the properties of the new algorithm with several numerical examples.
We discuss Möbius transformations for general matrix polynomials over arbitrary
elds, analyzing their in
uence on regularity, rank, determinant, constructs such as compound
matrices, and on structural features including sparsity and symmetry. Results on
the preservation of spectral information contained in elementary divisors, partial multiplicity
sequences, invariant pairs, and minimal indices are presented. The eect on canonical
forms such as Smith forms and local Smith forms, on relationships of strict equivalence
and spectral equivalence, and on the property of being a linearization or quadratication
are investigated. We show that many important transformations are special instances
of Möbius transformations, and analyze a Möbius connection between alternating and
palindromic matrix polynomials. Finally, the use of Möbius transformations in solving
polynomial inverse eigenproblems is illustrated.
This paper focuses on
rate-independent damage in elastic bodies. Since the driving energy is nonconvex,
solutions may have jumps as a function of time, and in this situation it is known that the classical concept
of energetic solutions for rate-independent systems
may fail to accurately describe the
behavior of the system at jumps.
Therefore, we resort to the (by now well-established) vanishing viscosity approach to rate-independent modeling
and approximate the model by its viscous regularization.
In fact, the analysis of the latter PDE system presents
remarkable difficulties, due to its highly nonlinear character.
We tackle it by combining a variational approach to a class of abstract doubly nonlinear evolution equations, with
careful regularity estimates tailored to this specific system relying on a q-Laplacian type gradient regularization of the damage variable.
Hence, for the viscous problem we conclude the existence of weak solutions satisfying a
suitable energy-dissipation inequality that is the starting point for the vanishing viscosity analysis.
The latter leads to the notion of (weak) parameterized
solution to our rate-independent system,
which encompasses the influence of viscosity in the description of the jump regime.
We study the approximation of Wasserstein gradient structures
by their finite-dimensional analog. We show that simple
finite-volume discretizations of the linear Fokker-Planck
equation exhibit the recently established entropic gradient-flow
structure for reversible Markov chains. Then we reprove the
convergence of the discrete scheme in the limit of vanishing
mesh size using only the involved gradient-flow structures.
In particular, we make no use of the linearity of the equations
nor of the fact that the Fokker-Planck equation is of second order.
Spectral deferred correction methods for solving stiff ODEs are known to
converge rapidly towards the collocation limit solution on equidistant grids,
but show a much less favourable contraction on non-equidistant grids such as
Radau-IIa points. We interprete SDC methods as fixed point iterations for the
collocation system and propose new DIRK-type sweeps for stiff problems based
on purely linear algebraic considerations. Good convergence is recovered also
on non-equidistant grids. The properties of different variants are explored on
a couple of numerical examples.
The well-known structure theorem of Hadamard-Zol\'esio states that the derivative of a shape functional is a distribution on the boundary of the domain depending only on the normal perturbations of a smooth enough boundary. However a volume representation (distributed shape derivative) is more general than the boundary form and allows to work with shapes having a lower regularity.
It is customary in the shape optimization literature to assume regularity of the domains and use the boundary expression of the shape derivative for numerical algorithm. In this paper we describe the numerous advantages of the distributed shape derivative in terms of generality, easiness of computation and numerical implementation. We give several examples of numerical applications such as the inverse conductivity problem and the level set method.
For the solution of optimal control problems governed by nonlinear parabolic PDEs, methods working on the reduced objective functional are often employed to avoid a full
spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of
the adjoint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories
have to be stored. To get numerical results which are accurate enough, in many case very fine discretizations in time and space are necessary, which leads to a significant
amount of data to be stored and transmitted to mass storage. Lossy compression methods were developed to overcome the storage problem by reducing the accuracy of the stored
trajectories. The inexact data induces errors in the reduced gradient and reduced Hessian. In this paper, we analyze the influence of such a lossy trajectory compression
method on Newton-CG methods for optimal control of parabolic PDEs and design an adaptive strategy for choosing appropriate quantization tolerances.
Elementary flux mode (EM) analysis is an important tool in the constraint-based analysis of genome-scale metabolic networks. Due to the combinatorial complexity of these networks, as well as the advances in the level of detail to which they can be reconstructed, an exhaustive enumeration of all EMs is often not practical. Therefore, in recent years interest has shifted towards searching EMs with specific properties. We present a novel method that allows computing EMs containing a given set of target reactions. This generalizes previous algorithms where the set of target reactions consists of a single reaction. In the one-reaction case, our method compares favorably to the previous approaches. In addition, we present several applications of our algorithm for computing EMs containing two target reactions in genome-scale metabolic networks.
Resolving thin conducting sheets for shielding or even skin layers inside by the mesh of numerical methods like the finite element method (FEM) can be avoided by using impedance transmission conditions (ITCs). Those ITCs shall provide an accurate approximation for small sheet thicknesses $d$, where the accuracy is best possible independent of the conductivity or the frequency being small or large -- this we will call robustness. We investigate the accuracy and robustness of popular and recently developed ITCs, and propose robust ITCs which are accurate up to $O(d^2)$.
We consider multiscale Hamiltonian systems in which harmonic oscillators with several high frequencies are coupled to a slow system. It is shown that the oscillatory energy is nearly preserved over long times eps^{-N} for arbitrary N>1, where eps^{-1} is the size of the smallest high frequency. The result is uniform in the frequencies and does not require non-resonance conditions.
Sobolev stability of plane wave solutions to the cubic nonlinear Schrödinger equation on a torus
(2013)
It is shown that plane wave solutions to the cubic nonlinear Schrödinger equation on a torus behave orbitally stable under generic perturbations of the initial data that are small in a high-order Sobolev norm, over long times that extend to arbitrary negative powers of the smallness parameter. The perturbation stays small in the same Sobolev norm over such long times. The proof uses a Hamiltonian reduction and transformation and, alternatively, Birkhoff normal forms or modulated Fourier expansions in time.
Coarsening in solder alloys is a widely accepted indicator for
possible failure of joints in electronic devices.
Based on the well-established Cahn--Larch\'e model
with logarithmic chemical energy density (see, e.g., Dreyer and Mueller 2001),
we present a numerical environment for the efficient
and reliable simulation of coarsening in binary alloys.
Main features are adaptive mesh refinement based on hierarchical error estimates,
fast and reliable algebraic solution by multigrid and Schur--Newton multigrid methods,
and the quantification of the coarsening speed by the temporal growth of mean phase radii.
We provide a detailed description and a numerical assessment of the algorithm and its different components,
together with a practical application to a eutectic AgCu brazing alloy.
We investigate the convergence of an implicit Voronoi finite volume
method for reaction-diffusion problems including nonlinear diffusion
in two space dimensions. The model allows to handle heterogeneous
materials and uses the chemical potentials of the involved species as
primary variables. The numerical scheme uses boundary conforming Delaunay
meshes and preserves positivity and the dissipative property of the
continuous system. Starting from a result on the global stability of
the scheme (uniform, mesh-independent global upper and lower bounds),
we prove strong convergence of the chemical activities and their gradients
to a weak solution of the continuous problem. In order to illustrate the
preservation of qualitative properties by the numerical scheme, we present
a long-term simulation of the Michaelis-Menten-Henri system. Especially,
we investigate the decay properties of the relative free energy and the
evolution of the dissipation rate over several magnitudes of time, and
obtain experimental orders of convergence for these quantities.
In this work we present a complete algorithm for the exact computation of the guided mode band structure in photonic crystal (PhC) wave-guides. In contrast to the supercell method, the used approach does not introduce any modelling error and is hence independent of the confinement of the modes. The approach is based on Dirichlet-to-Neumann (DtN) transparent boundary conditions that yield a nonlinear eigenvalue problem. For the solution of this nonlinear eigenvalue problem we present a direct technique using Chebyshev interpolation that requires a band gap calculation of the PhC in advance. For this band gap calculation –- we introduce as a very efficient tool –- a Taylor expansion of the PhC band structure. We show that our algorithm –- like the supercell method –- converges exponentially, however, its computational costs –- in comparison to the supercell method –- only increase moderately since the size of the matrix to be inverted remains constant.
In high accuracy numerical simulations and optimal control of time-dependent processes, often both many time steps and fine spatial discretizations are needed. Adjoint gradient computation, or post-processing of simulation results, requires the storage of the solution trajectories over the whole time, if necessary together with the adaptively refined spatial grids. In this paper we discuss various techniques to reduce the memory requirements, focusing first on the storage of the solution data, which typically are double precision floating point values. We highlight advantages and disadvantages of the different approaches. Moreover, we present an algorithm for the efficient storage of adaptively refined, hierarchic grids, and the integration with the compressed storage of solution data.
A Fenchel dualization scheme for the one-step time-discretized contact problem of quasi-static
elasto-plasticity with kinematic or isotropic hardening is considered. The associated path is induced
by a combined Moreau-Yosida / Tichonov regularization of the dual problem.
The sequence of solutions to the regularized problems is shown to converge strongly to the optimal
displacement-stress-strain triple of the original elasto-plastic contact problem in the space-continuous setting.
This property relies on the density of the intersection of certain convex sets which is shown as well.
It is also argued that the mappings associated with the resulting problems are Newton- or slantly differentiable.
Consequently, each regularized subsystem can be solved mesh-independently at a local
superlinear rate of convergence. For efficiency purposes, an inexact path-following approach is proposed and
a numerical validation of the theoretical results is given.
The inverse eigenvalue problem for $T$-alternating matrix polynomials over arbitrary
algebraically closed fields of characteristic different from two is considered.
The main result shows that the necessary conditions obtained in \cite{MacMMM10} for a matrix
polynomial to be the Smith form of a $T$-alternating matrix polynomial are under mild
conditions also sufficient to be the Smith form of a $T$-alternating matrix polynomial with
invertible leading coefficient which is additionally in anti-triangular form.. In particular, this result
implies that any $T$-alternating matrix polynomial with invertible leading coefficient is
equivalent to a $T$-alternating matrix polynomial in anti-triangular form that has the
same finite and infinite elementary divisors as the original matrix polynomial.
Finally, the inverse eigenvalue problem for $T$-palindromic matrix polynomials is
considered excluding the case that both $+1$ and $-1$ are eigenvalues.
In this paper, we introduce a high order finite element (FEM) implementation using perfectly matched layer (PML) for the scattering by plasmonic structures inside layered media. The PML is proven to be very accurate and efficient by a comparative analysis with a commercial FEM software and the Multiple Multipole Program (MMP). A convergence analysis using hp-adaptive refinement inside the PML layer
shows that adaptive mesh refinement inside the PML layer is most efficient. Based on this convergence analysis an hp-strategy is proposed, which shows a remarkable error reduction for small additional computational costs.
Shielding sheets are commonly used in the protection of electronic devices. With their large aspect ratios they become a serious issue for the direct application of boundary element methods due to the occurring almost singular integrals. Impedance transmission conditions allow for boundary element formulations where only sheet mid-line has to be triangulated. We propose and analyse boundary element methods of second kind for impedance transmission conditions of two different types for the time-harmonic eddy current problem in two dimensions. In numerical experiments we study the effect on the magnetic shielding efficiency of the discretisation error due the boundary element discretisation and the modelling error due to the
impedance transmission condition.
We study the exact recovery of signals from quantized frame coefficients. Here, the basis of the quantization
is hard thresholding, and we present a simple algorithm for the recovery of reconstructable signals. The set of
non-reconstructable signals is shown to be star-shaped and symmetric with respect to the origin. Moreover, we
provide a criterion on the frame for the boundedness of this set. In this case, we also give a priori bounds.
According to the Helmholtz decomposition,
the irrotational parts of the momentum balance equations of
the incompressible Navier-Stokes equations are balanced by the
pressure gradient.
Unfortunately, nearly all mixed methods for incompressible flows
violate this fundamental property, resulting in the well-known
numerical instability of poor mass conservation.
The origin of this problem is the
lack of $L^2$-orthogonality between discretely divergence-free velocities
and irrotational vector fields.
In order to cure this,
a new variational crime using divergence-free
velocity reconstructions is proposed.
Applying lowest order Raviart-Thomas velocity reconstructions
to the nonconforming Crouzeix-Raviart element
allows to construct
a cheap flow discretization for general 2d and 3d
simplex meshes that possesses the same
advantageous robustness properties like divergence-free flow
solvers.
In the Stokes case,
optimal a-priori error estimates for the velocity gradients
and the pressure are derived. Moreover, the discrete velocity
is independent of the continuous pressure.
Several detailed linear and nonlinear
numerical examples illustrate the theoretical findings.
A basic task in signal analysis is to character-
ize data in a meaningful way for analysis and classification
purposes. Time-frequency transforms are powerful strategies
for signal decomposition, and important recent generalizations
have been achieved in the setting of frame theory. In parallel
recent developments, tools from algebraic topology, traditionally
developed in purely abstract settings, have provided new insights
in applications to data analysis. In this report, we investigate some
interactions of these tools, both theoretically and with numerical
experiments, in order to characterize signals and their frame
transforms. We explain basic concepts in persistent homology
as an important new subfield of computational topology, as well
as formulations of time-frequency analysis in frame theory. Our
objective is to use persistent homology for constructing topo-
logical signatures of signals in the context of frame theory. The
motivation is to design new classification and analysis methods by
combining the strength of frame theory as a fundamental signal
processing methodology, with persistent homology as a new tool
in data analysis.
Primal heuristics are an important component of state-of-the-art codes for mixed integer nonlinear programming (MINLP). In this article we give a compact overview of primal heuristics for MINLP that have been suggested in the literature of recent years. We sketch the fundamental concepts of different classes of heuristics and discuss specific implementations. A brief computational experiment shows that primal heuristics play a key role in achieving feasibility and finding good primal bounds within a global MINLP solver.
Enforcing solvability of a nonlinear matrix equation and estimation of multivariate ARMA time series
(2013)
The matrix equation $X+AX^{-1}A^T=B$, arising in parameter estimation of certain time series models,
is solvable only for certain values of the matrices $A,B$.
We present a numerical method to modify $A,B$ in order to make the matrix equation solvable.
Since solvability depends on the location of the eigenvalues of the palindromic matrix polynomial $\lambda^2 A+\lambda B+A^T$,
our method works by moving those eigenvalues to specified locations using first order spectral perturbation theory.
The method is heuristic but works in practice, as is supported by several compelling numerical examples.
These examples arise from parameter estimation of a common time series model, the multivariate ARMA(1,1).
We consider a general class of continuous asset price models where the drift and the volatility functions, as well as the driving Brownian motions, change at a random time. Under minimal assumptions on the random time and on the driving Brownian motions, we study the behavior of the model in all the filtrations which naturally arise in this setting, establishing martingale representation results and characterizing the validity of the NA1 and NFLVR no-arbitrage conditions.