Refine
Year of publication
- 2009 (143) (remove)
Keywords
- mixed integer programming (8)
- branch-and-cut (5)
- integer programming (5)
- metastability (4)
- optimal control (4)
- BSDE (3)
- Malliavin calculus (3)
- backward stochastic differential equation (3)
- optimization software (3)
- primal heuristics (3)
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.
The goal of this paper is to show how the heat treatment of steel can be modelled in terms of a mathematical optimal
control problem. The approach is applied to laser surface hardening and the cooling of a steel slab including
mechanical effects. Finally, it is shown how the results can be utilized in industrial practice by a coupling with
machine-based control.