Abschlussarbeiten

The Steiner tree packing problem is a long studied problem in combinato-
rial optimization. In contrast to many other problems, where an enormous
progress has been made in the practical problem solving, the Steiner tree
packing problem remains very difficult. Most heuristics schemes are ineffec-
tive and even finding feasible solutions is already NP -hard. What makes this
problem special, is that in order to reach an overall optimal solution non-
optimal solutions to the underlying NP -hard Steiner tree problems must be
used. Any non-global approach to the Steiner tree packing problem is likely
to fail. Integer programming is currently the best approach for computing
optimal solutions.
The goal of this master thesis is to give a survey of models relating to the
Steiner tree packing problem from the literature. In addition, a closer look
at a model for the switchbox routing problem in VLSI-Design will be given.

When reformulating a given mixed integer program by the use of classical Dantzig-Wolfe decomposition, a subset of the constraints is partially convexified, which corresponds to implicitly adding all valid inequalities for the associated integer hull. Since these inequalities are not known, a solution of the original linear programming (LP) relaxation which is obtained by transferring an optimal basic solution of the reformulated LP relaxation is in general not basic. Hence, cutting planes which are separated using a basis like Gomory mixed integer cuts are usually not directly applicable when separating such a solution.

Nevertheless, we can use some crossover method in order to obtain a basic solution which is nearby the considered non-basic solution and generate cutting planes for the basic solution using its basis information. These cutting planes might also the solution we originally wanted to separate. So far, this problem was only considered extensively by Range, who proposed the previously described approach including a particular crossover method. We present a modified crossover method and extend this procedure by considering additional valid inequalities strengthening the original LP relaxation. Furthermore, we provide the first full implementation of a separator like this and tested it on instances of several problem classes.

In this Master Thesis we consider the problem of covering rectilinear polygons by the
minimum number of axis-parallel rectangles. This problem finds applications in the
fabrication of DNA chip arrays [Hannenhalli et al. 2002], in VLSI design, where the
rectilinear polygon is a chip that has to be covered by a huge number of rectangular
transistors.
Other applications are data compression and in particular image compression, where
large rectangular areas with the same color can be compressed into one pixel.

This work evaluates whether optimization problems resulting from modelling the optimal synthesis, design and operation of a decentralized energy supply system have an embedded structure which can be exploited by decomposition methods in solution algorithms. The objective is to determine if the the accuracy and/or the problem size in terms of number of periods of time and number of units considered may be increased, as these are limited if the branch-and-bound method combined with the simplex method is used to solve the problems. A model of the problem is formulated as a mixed-integer linear program as proposed by Yokoyama et al. (2002) and Voll (2013). The model is analyzed and two embedded structures suitable for decomposition are identified.

The first structure emphasizes the independent operation and design of every component. The second structure emphasizes the design and operation of all components and focuses on the independence of every period of time considered. The model is reformulated using the Dantzig-Wolfe decomposition principle for both proposed embedded structures. A numerical study is conducted where the synthesis, design and operation of a fictional energy supply system is optimized by both the branch-and-bound method combined with the simplex method and by the branch-and-price method. A set of instances is created for different degrees of complexity in terms of the number of units and the number of periods of time considered.

The results show that the dual bounds obtained by solving the rootnode LP relaxation can be improved in comparison to the conventional solution approach, if the reformulation emphasizing independent components is utilized. The results provide no evidence on improvements on the considered test set for the reformulation emphasizing design and operation. For the case of an optimal solution computing times required to solve the considered instances of a test set are found to be reduced by utilizing the branch-and-price method and the reformulation emphasizing components, if identical components are considered in the energy supply system in comparison to the non-commercial solver SCIP.

In der Automobilindustrie führen lange Lieferwege und das Prinzip der "Just-In-Time-Produktion" zu einem unflexibelen Produktionsplan für mehrere Wochen. Gleichzeitig ist die exakte Vorhersage von Kundenwünschen, bedingt durch das hohe Maß an Spezialisierung (Mass Customization), nahezu unmöglich. In Kombination mit laufzeittechnischen Anforderungen beschreiben diese Bedingungen ein herausforderendes Problem in der Praxis.

In dieser Abschlussarbeit wird unter Anderem ein Algorithmus vorgestellt, in welchem ganzzahlige lineare Programmierung dazu genutzt wird, die durchschnittliche Verkaufszeit eines Autos zu reduzieren. Der Algorithmus ist in der Lage, basierend auf einem existierenden Produktionplan, eine Umplanung zu berechnen, welche die strengen Lagereinschränkungen an der Fabrik und die Kundenwünsche respektiert und gleichzeitig das frühst mögliche Lieferdatum 
für den Kunden garantiert.

In dieser Abschlussarbeit soll ein Algorithmus implementiert werden, der für Schüler und Lehrer des Humboldt Gymnasiums Solingen eine zulässige Zimmerbelegung berechnet, die u.a. möglichst kostengünstig ist, also möglichst wenig Zimmer benötigt. Weitere Einschränkungen sind: - Jeder Schüler ist einem Haus und einem Zimmer zugeordnet. - Sowohl Mädchen und Jungen als auch Schüler und Lehrer schlafen getrennt. - Da in der Regel fünf Klassen auf vier Häuser verteilt werden, wird mindestens eine Klasse in mehr als einem Haus untergebracht werden müssen. Falls eine Klasse geteilt wird, soll sie auf maximal zwei Häuser verteilt werden. Eventuell wird die zu teilende Klasse vorher angegeben. - Die meisten Zimmer können mit zusätzlichen Betten ausgestattet werden. Die Anzahl der Zimmer soll minimiert werden unter der Berücksichtigung, dass die Zimmer nicht zu voll werden. Das fertige Progamm soll dem Benutzer ermöglichen Daten zu einlesen, Wünsche zu äußern (z.B. welche Klasse geteilt werden soll) und die Zimmerbelegung in einem anschaulichen Format auszuschreiben. Es werden Implementationskenntnisse in C/C++ benötigt.

Das Consecutive Dinner Problem beschreibt die Tourenplanung bzw die Gruppenfindung für teilnehmende Teams. Jede Route muss aus genau drei Stationen bestehen: Vorspeise, Hauptspeise und Nachspeise. Des Weiteren muss jedes Team genau einen Gang der eigenen Route selber zubereiten. Eine Gruppe, die sich bei einem Gang trifft, muss aus genau drei Teams bestehen, die sich nicht schon vorher bzw auch nachher nicht nochmal treffen werden. Jedes Team hat die Möglichkeit einen Wunschgang zu nennen. Die Zielfunktion des Problems minimiert einerseits die Strecke, die jedes Team im Laufe des Abends zurücklegen muss, und andererseits maximiert die Anzahl an Teams, die ihren Wunschgang zugeteilt bekommen. Mögliche Erweiterungen des Modells enthalten Differenzierungen der Teams durch Alter und/oder Geschlecht um gemischtere/homogenere Gruppenkonstellation zu erhalten. Auch könnte die Anzahl der schon teilgenommenen Veranstaltungen mit einbezogen werden. In dieser Masterarbeit soll ein IP entwickelt werden, das das oben beschriebene Problem beschreibt. Außerdem soll ein Software Tool entwickelt werden, was das Consecutive Dinner Problem löst. Daher werden Kenntnisse in der mathematischen, ganzzahligen Optimierung und Implemenationskenntnisse in C/C++ oder Java benötigt.

Die vorliegende Arbeit hat das Ziel, einen Algorithmus zu entwickeln, welcher die Knoten eines Graphen verschiedenen Clustern zuteilt. Dieser Graph wird hierbei aus einer logistischen Problemstellung des Transportes von Kraftfahrzeugen abgeleitet. Dies geschieht vor dem Hintergrund der Unterteilung eines Multi Commodity Flow Problems in mehrere kleinere Probleme, welche durch die einzelnen Cluster gebildet werden sowie durch die zwischen den Clustern verlaufenden Kanten des Graphen. Somit soll die zur Lösungsfindung des logistischen Problems benötigte Laufzeit verringert werden können. Hierfür wird zunächst ein Überblick über den derzeitigen Stand der wissenschaftlichen Diskussion zum Thema Clustering gegeben. Im Anschluss werden verschiedene Clusteringansätze hinsichtlich ihrer Anwendbarkeit auf das vorliegende Problem bewertet. Sodann sollen passende Ansätze ausgewählt und hieraus ein Algorithmus für das vorliegende Problem entwickelt, dargestellt und bewertet werden. Dieser Algorithmus wird zuletzt an verschiedenen Instanzen des der Arbeit zugrundeliegenden Problems getestet und die erzielten Ergebnisse analysiert.

The Energy Minimizing Vehicle Routing Problem (EMVRP) introduced by Kara and Yetis in 2007 is an extended version of the traditional Vehicle Routing Problem of transportation logistics. In this thesis we consider different algorithms, which heuristically solve the EMVRP in reasonable computation time. This thesis also includes a comparison between the proposed heuristics, which lays out the advantages and disadvantages of each one using different sets of test instances.

Dantzig-Wolfe (DW) reformulation is a well-known approach to produce strong dual bounds for Mixed-Integer Programs (MIPs). If a MIP has a special structure which can be exploited well by a DW reformulation the solver could be much faster on the reformulated model than on the original. Otherwise the solver may fail completely. Providing a strong DW reformulation typically requires enhanced knowledge about the underlying model structure of the MIP which makes an automatic reformulation without such knowledge difficult. Despite this problem there have been several approaches proposed which implement such an automatic DW reformulation process into a state-of the-art MIP solver, especially the Generic Column Generation (GCG) framework [1]. In such frameworks the detection of a decomposition that exploits the underlying model structure is one of the most important aspects. GCG for example uses different detectors that use heuristics and clustering algorithms to generate a set of possible decompositions. Recently supervised learning algorithms have been applied by Kruber et al. [2] to decide if a reformulation of the model should be solved (and also which if several have been detected) or if the original model should better be solved by a standard solver. These first promising results showed that machine learning could be a useful support in this decision process. Besides this single question there might be other decisions in the entire process that could also potentially benefit by the use of supervised machine learning.

In Zusammenarbeit mit dem Centre de recherches energetiques et municipales CREM, Martigny, Schweiz (https://www.crem.ch/).

The bachelor thesis deals with solving multi-commodity flow problems. These kinds of problems occur in different everyday situations. For example, the transport of goods in distribution networks or the movement of messages in communication networks. In addition to the considerable size of these flow problems, which often occurs, there can be further complicating constraints on nodes, edges, or goods. In the context of the work, an alternative solution approach for minimal-cost multi-commodity flow problems, which is based for the most part on the computation of k-shortest paths, is developed. Based on optimally solved single-sink single-source single-commodity flow sub-problems, whose accumulated solutions are usually an inadmissible solution of the multi-commodity flow, the developed solution approach tries to obtain a feasible solution by rerouting the flow of edges that are infeasible in the original problem. The developed approach will be evaluated in a feasibility study based on a real use case, which is a time-expanded distribution network of an automobile manufacturer with limitations on edges and goods.

This thesis aims to evaluate the boundaries of neural network training using
mixed-integer optimization. Neural networks are a powerful machine learn-
ing technique that is able to classify data non-linearly. This technique is
applied to a lot of different tasks. However, traditional approaches to neural
network training can get stuck in local minima, while algorithms exist that
solve a mixed-integer optimization problem optimally. The goal of this thesis
in the context of mixed-integer programming. We start by showing which
kind of activation functions can be expressed by a mixed-integer program
(MIP). After that variants of neural network training using mixed-integer
programming are presented and it is shown which activation functions can
be incorporated into these variants. Finally, we use a blackbox-solver to ap-
ply the resulting models to the parity function and evaluate the performance
of our approaches.

We introduce an optimization problem on graphs called Connected Vertex Clustering Problem (CVCP). The input consists of a finite graph, arbitrary linear constraints to restrict the set of feasible clusters and a linear objective function. The expected solution is a vertex clustering that optimizes the objective under the condition that each cluster induces a connected subgraph and satisfies the custom constraints. Besides partitional clustering, i.e., node partitioning, the solution may also be restricted to packings or coverings of nodes. We show that this highly configurable problem is NP-hard and propose a branch-and-price method as a solution approach. The suggested method is implemented as a framework which is capable of solving arbitrary CVCP instances. The framework can easily be extended with new features due to its plug-in architecture. This allows to exploit the characteristics of specific variants of the CVCP in order to enhance the efficiency of the solution process. We evaluate the developed framework on a districting problem for the German federal elections and on the Odd Cycle Packing Problem.

In this bachelor thesis the optimasation of the airplane boarding process is investigated. Given is a set of cabin layouts and a set of passengers with different movement speeds and amounts of carry-on luggage. To improve the overall boarding time each seat has to be assigned to a passanger and each carry-on luggage space to a luggage. The passanger satisfaction sets further requirements to the boarding process. The assignment problem is formulated as a mixed integer problem (MIP) and differentiates itself from previous works by the direct optimisation of the overall boarding time, in combination with different numbers of carry-on luggage per passanger such as limitaions of the carry-on luggage space. To investigate the assignment problem, the formulated mathematical model such as developed heuristics are implemented. The experiments with different calculations show a strong correlation between the overall boarding time and the amount of carry-on luggage to stow in the carry-on luggage space in the airplane cabin. The best solution is achieved using the heuristic \(H_{wtw}^{v2}), giving excellent boarding times in addition to a high amount of stowed carry-on luggage reaching close to maximum values.

Emissionsfreie Mobilität gewinnt immer mehr an Wichtigkeit. Für diese sind Elektrofahrzeuge unabdingbar, dadurch entstehen neue Fragen und Herausforderungen. Das iMove Projekt beschäftigt sich unter anderem mit dem „balancierten Ladeproblem (BCP)“. Hierbei wird eine Plattform entwickelt, die die Ladevorgänge von Elektrofahrzeugen optimal den Ladesäulen des Elektrische Netzes zuweisen soll. Optimalität ist bei diesem Problem zweigeteilt. Auf der einen Seite sollen die Wünsche der Elektrofahrzeugfahrer erfüllt werden. Auf der anderen Seite soll die Last auf das Elektrische Netz möglichst ausgeglichen bleiben. Hierfür wurde ein MIP entwickelt. Diese Masterarbeit beschäftigt sich mit dem Problem, dass die Elektrofahrzeugfahrer möglicherweise die Zuweisung zu einer Ladesäule der Plattform ignorieren und dieser nicht folgen. Während robuste Ansätze für diese Art der Unsicherheit nicht optimal geeignet scheinen, stellt es sich heraus das Online Ansätze eine gute Alternative bieten mit abweichenden Fahrern umzugehen. Die beste Güte erzielen Online Ansätze mit Unterstützung der Offline-Lösungen. Das Online-BCP ermöglicht es trotz abweichender Fahrer eine zulässige Lösung zu finden, allerdings zeigt sich, dass die Güte des Online-BCP noch ausbaufähig ist.

The currently ongoing energy transition increases the share of renewable energies in the electricity mix. In Germany, this is covered in particular by wind and solar energy, which are considered volatile. Since a stable power grid must guarantee a balance between the electricity produced and the electricity demanded at all times, this development presents new challenges for electricity suppliers and grid operators. One approach to counteract the fluctuation in the energy supply network is demand-side management (DSM). DSM refers to the management of demand for network-based services by customers in industry, trade and private households. Electricity suppliers benefit from DSM by increasing the reliability of supply, while electricity consumers can reduce their energy costs benefiting from financial incentives. This Master's thesis investigates, with regard to demand-side management, how optimized production planning can reduce the energy costs of a general energy-intensive production process. For this purpose, a mixed-integer linear model is developed and set up. In particular, two new possibilities for cost reduction will be modeled: on the one hand, the exploitation of volatile electricity prices by load redistribution, on the other hand, the provision of balancing power. In production planning, the various framework conditions of the production process and uncertainties in some model parameters, such as the demand for balancing power, must be taken into account. Therefore, methods from the field of robust optimization are applied as a solution for the optimization problem.

With an increasing number of man-made objects in space, the area around earth is getting increasingly polluted. Due to the chance of collisions, orbiting debris puts future missions at risk. The Kessler syndrome projects a cascading amount of debris collisions that result in an exponential increase of debris pieces. In order to avoid a catastrophic scenario, space agencies such as the European Space Agency propose missions for an active removal of space debris. Depending on the removal method, one mission aims to remove multiple debris objects. Thus, during the mission design, a combinatorial problem needs to be solved deciding the selection of debris objects and the order of space rendezvous. The problem setup shares similarities with a variant of the traveling salesman problem with city selection. The goal is to maximize the profit of the visited cities with a limited traveling budget. Since the debris orbits are exposed to perturbations, it is assumed that the costs of traveling between the objects are likely to change over time. This thesis will explore various mixed integer problem formulations for both constant costs and dynamic costs over time. For the static case, we are able to find optimal routes even for bigger debris clouds of more than 2000 objects. The proposed formulations for the dynamic case are able to solve large instances, when carefully selecting the time parameters. Moreover, we have implemented a genetic algorithm to further explore the solution space.

In this thesis, we investigate possible optimisations of the aeroplane boarding process. To accomplish this, we define and mathematically model the process, formulate it as a minimisation problem, and propose a solution procedure.We show that the problem is computationally hard in a theoretical sense and investigate the quality of heuristic approaches as well as some of its online properties.

The boarding process is of particular interest for airlines to improve profitability as well as customer satisfaction. Since boarding policy changes pose a lower implementation barrier than buying new planes or making changes to existing infrastructure, advances in aeroplane boarding research have the potential to positively impact customers and airlines relatively immediately.

The use of machine learning methods such as artificial neural networks or clustering algorithms is becoming more and more popular in the optimization community. Connecting to this trend, we propose a method that, instead of tuning mixed-integer programming solvers' behavior, facilitates testing (modified) solver behavior. Utilizing dynamic clusterings and solving mixed-integer programs, we select and visualize different diverse subsets of the test data available in the structured integer programming library strIPlib that are just as representative but exhibit a much smaller size. The generated subsets do not only comprise a strIPlib collection and a 'general-purpose' benchmark test set, but also a method to generate completely customized, diverse experiment test sets, tailored to the experiment that is to be executed with them.

Mixed-integer programming (MIP) is a common technique to model optimization problems. They can be solved with the decomposition approach of Dantzig-Wolfe reformulation (DWR) and the branch-cut-and-price (BC&P) algorithm. Many state-of-the-art software packages for solving MIPs provide scripting interfaces to interactively model and optimize these problems. We introduce PyGCGOpt, a Python interface for the generic BC&P solver GCG. The interface allows to model MIPs, to construct and visualize decompositions, and to implement solvers for pricing problems which occur from DWR. We provide example applications of our interface and show that it interacts with the solver without a significant decline in performance.

Das klassische Job-Shop-Scheduling-Problem besteht darin, eine optimale Zuordnung anfallender Aufgaben zu bestehenden Maschinen zu finden. Die Optimalität der Lösung kann unterschiedlich definiert sein, zielt in den meisten Fällen aber auf eine möglichst schnelle und/oder kostengünstige Abwicklung ab. Die vorliegende Arbeit erweitert dieses klassische Modell gemäß einer der Praxis entlehnten Problemstellung um eine dynamische Komponente. In eine schon bestehende Lösung sollen nun neu hinzukommende Aufgaben integriert werden. Die neu entstandene Situation wird separat modelliert und hinsichtlich eines eigenen Optimalitätskriteriums betrachtet.

Dantzig-Wolfe Zerlegung ist eine bekannte Methode, mit der lineare gemischt-ganzzahlige Progamme (engl. mixed integer programs), die eine gewisse Struktur haben, umformuliert und so bessere Grenzen für den optimalen Zielfunktionswert erhalten werden können. Die Dantzig-Wolfe Zerlegung wird daher in Branch-and-Bound eingesetzt, wo die Laufzeit von der Qualität der ermittelten Grenzwerte abhängt. JuMP ist eine Modellierungsprache für mathematische Optimierungsprobleme, die in der Programmiersprache Julia eingebettet ist; Coluna ist ein Framework, das JuMP erweitert und die Anwendung einer Dantzig-Wolfe Zerlegung ermöglicht. In Coluna muss jedoch der Benutzer die Problemstruktur beschreiben, welche die Anwendung einer Dantzig-Wolfe Zerlegung ermöglicht. Wir beschreiben, wie dem Benutzer diese Aufgabe abgenommen werden kann und präsentieren unsere Implementierung einer automatischen Dantzig-Wolfe Zerlegung für Coluna. Die Automatisierung der Dantzig-Wolfe Zerlegung erfolt in zwei Schritten. Zunächst erklären wir, wie eine Menge von potentiellen Strukturen für ein gegebenes lineares gemischt-ganzzahliges Programm mithilfe der Indexmengen der Bedingungen identifiziert werden kann. Im zweiten Schritt präsentieren wir verschiedene Scores aus der Literatur, welche zum Ziel haben, die für eine Dantzig-Wolfe Zerlegung potentiell beste Struktur in einer Menge von Kandidaten zu identifizieren.

The selection of the optimal decomposition for Danzig-Wolfe reformulations presents a difficult task for a generic MIP solver. The determination of the structure of the program can contain valuable information for this process. Therefore detection algorithms are used to extract the characteristics of the MIP. In this thesis a new approach is introduced to classify the MIP as a whole, in contrast to commonly used constrain and variable classification. A graph is used to summarize and analyse the extracted features of the MIP. This graph is generalised, by estimating the sets of indices, that were used for the creation of the MIP. Afterwards can the compression with other graphs, which represent different types of MIPs, result in a categorization of the examined MIP. The implementation of this approach showed promising success, when tested with multiple different types of MIPs. For well structured MIPs, the algorithm was able to reliable detect the type of the MIPs and recommend decompositions accordingly. Problems exist in the handling of small irregularities in the MIPs or structure changed by the presolving process. Therefore, further improvements of the robustness of the algorithm could enhance the classification quality.

Visual Speech Recognition, also called lipreading, is the task of interpreting speech and predicting text by only analyzing the movements of a speaker’s mouth. Lipreading is a difficult task for both humans and computers. The field of Machine Learn and its promising sub-field Deep Learn, have proven to be successful in tackling complex problems. Recent advancements have enabled lipreading systems to use deep learn models which are trained end-to-end. This thesis describes the development of such a deep learn based lipreading system which as an end goal will aid temporarily speech impaired patients in the intensive care units of hospitals in communicating more efficiently. It approaches the problem by breaking it down into two stages which are trained separately. The first stage predicts audio features from video frames and the second stage takes the predicted audio features and predicts the spoken text. The audio data that is available can be utilized for training the system. The lipreading system that was developed using this approach was able to achieve results close to the popular baseline model LipNet. Further, the collective behavior of the two stages of the model was evaluated. The presented insights serve as valuable input for future research in the broader project to which this thesis belongs.

Aufgrund der Energiewende ist es unabdingbar, dass sich der europäische Strommarkt in den nächsten Jahren deutlich verändern wird. Um diese Veränderungen vorherzusagen und darauf zu reagieren, wird das europäische Stromnetz mittels linearer Programmierung simuliert. Diese Arbeit ist in Kooperation mit der Firma Consentec GmbH entstanden, wo entsprechende Simulationen durchgeführt werden. Der Anspruch an die Genauigkeit der Simulation ist dabei in den letzten Jahren gewachsen. Konkret sollen Stromflüsse nicht mehr alleine auf der Ebene einzelner Länder oder Zonen (NTC), sondern auf Leitungsebene (flow-based) abgebildet werden. Dabei entstehen sehr viele Nebenbedingungen, von denen jedoch die allermeisten den Lösungsraum nicht verändern. Es werden verschiedene Ideen diskutiert, wie mit diesen Nebenbedingungen umgegangen werden kann. Vor allem wird ein Ansatz präsentiert, der das Herausfiltern der relevanten Nebenbedingungen deutlich beschleunigt.

Für gemischt-ganzzahlige Programme hat die Wahl der richtigen Branching-Variable einen großen Einfluss auf die Größe des resultierenden Branch-and-Bound (B&B) Baums, und damit auch auf die Laufzeit des gesamten Algorithmus. Vollständiges starkes Branching, also die Auswahl von Variablen durch vollständige Auswertung der Kindknoten die für die entsprechenden Kandidaten entstehen würden, ist gut darin, Bäume klein zu halten, wird aber wegen des hohen Berechnungsaufwand im Verhältnis zur eingesparten Baumgröße normalerweise nicht in B&B verwendet. Stattdessen wird es mit anderen Heuristiken kombiniert, was zu einigen der erfolgreichsten Kandidatenauswahlheuristiken für allgemeine Problem führt, wie zum Beispiel hybrides Pseudokosten/starkes Branching oder Zuverlässigkeitsbranching. Da die Auswertung einzelner Knoten in Branch-and-Price (B&P) wegen der durchzuführenden Spaltengenerierung im Allgemeinen länger dauert als in B&B sind sowohl die Kosten von starkem Branching als auch die Vorteile davon, einen kleinen Baum zu erhalten in B&P erhöht. Dies ändert potentiell die relative Leistung von bestehenden Auswahlheuristiken, und gibt die Möglichkeit, neue Heuristiken zu erstellen. Eine dieser neuen Heuristiken, genannt hierarchisches starkes Branching, kombiniert andere Heuristiken mit starkem Branching in einer hierarchischen Art und Weise. Wir erweitern hierarchisches starkes Branching, unter anderem indem wir es mit hybridem Pseuodokosten/starkem Branching und Zuverlässigkeitsbranching kombinieren. Danach evaluieren und vergleichen wir die Effektivität verschiedener, teilweise starkes Branching-basierten Kandidatenauswahlheuristiken, sowohl für Branching auf Originalvariablen als auch für Ryan-Foster Branching.

Die Modellierung eines Produktionsprozesses als gemischt ganzzahliges lineares Programm (MILP) erlaubt es, variable Kosten durch das Ausnutzen niedriger Strompreise zu sparen. Verkomplizierende Aspekte wie verderbliche Materialien können die Rechendauer solcher Optimierungsprobleme erheblich erhöhen. Dies führt oft zu ungenutztem Potenzial in der Lösungsqualität, weswegen effizientere Lösungsverfahren benötigt werden. Bei einer periodischen Produktionsplanung wird dasselbe MILP immer wieder gelöst, wobei sich nur die Eingangsdaten ändern. Allerdings wird die Erfahrung über Muster im Lösungsprozess meist nicht genutzt. In dieser Arbeit wird ein Produktionsplanungsproblem mit verderblichen Zwischenprodukten betrachtet. Es wird mit Hilfe von lazy constraints modelliert, d.h. verkomplizierende Ungleichungen werden separiert und während des Lösungsprozesses on-the-fly hinzugefügt. Anschließend wird gezeigt, dass ein neuronales Netz verletzte lazy constraints mit hoher Präzision vorhersagen kann, indem es Erfahrungen aus zuvor gelösten Instanzen nutzt. Das neue Modell ist bis zu 34% schneller als das schnellste konventionelle Modell und die lazy-Methode.

Eine Laser Powder Bed Fusion (LPBF) Maschine erstellt 3D Objekte durch schichtweises Aufschmelzen von Metallpuder mit Lasern. Bei diesem Prozess bieten sich mehrere Aspekte zur Optimierung an. Einerseits ist das Finden der optimalen Anzahl und Ausrichtung der Scanner für eine gegebenes Arbeitsstück von Belang. Außerdem bedarf es des bestmöglichen Pfades des Prozesskopfes über die Arbeitsfläche und letztlich muss die Arbeitslast auf die Scanner im Prozesskopf aufgeteilt werden. Diese Probleme werden als Mixed-Integer Programs modelliert und mit dem Gurobi Löser berechnet.

Die Masterarbeit befasst sich mit einem klassischen Optimierungsproblem aus der Finanzwelt, der für Banken spannenden Fragestellung der Portfolio-Optimierung. Es wird das Spannungsfeld betrachtet, dass Banken unter anderem die Maximierung ihres Zinsergebnisses anstreben, jedoch auch dazu verpflichtet sind gewisse regulatorische Vorgaben einzuhalten. Im Fokus dieser Arbeit standen einerseits das Aufsetzen einer fachlichen Konzeption des Optimierungsmodells sowie die technische Realisierung der für die Optimierung zusätzlich erarbeiteten Vorverarbeitungen. Andererseits schließt die Arbeit das Aufstellen sowie das entsprechende Umsetzen eines mathematischen Modells basierend auf ausgewählten regulatorischen Metriken ein. Es wurde ein lineares Optimierungsmodell aufgestellt, welches die Metriken NII (Net Interest Income), EVE (Economic Value of Equity) und NSFR (Net Stable Funding Ratio) sowie die Kapitalquoten und die Bilanz umfasst. Die Arbeit wurde fachlich von einem Anbieter von Softwareprodukten im Finanzwesen begleitet und die Optimierung basierend auf den Ergebnisdaten aus deren Standard-Softwarelösung für das Meldewesen von Banken aufgesetzt.

Die vorliegende Arbeit thematisiert die theoretische und praktische Analyse von gemischt ganzzahligen Programmen für orthogonale Verpackungsprobleme. Ziel der Arbeit ist es das Verpackungsproblem, welches in der Literatur als Multiple Bin Size Bin Packing Problem bekannt ist, zu untersuchen und Lösungsalgorithmen anzugeben. Insbesondere die Bedingung, dass Artikel überlappungsfrei in der Verpackung platziert werden müssen, hat in mehreren Dimensionen Probleme bereitet, wodurch Techniken aus dem eindimensionalem Fall nicht oder nur schwer erweitert werden können. Aus diesem Grund werden zunächst verschiedene Modellierungstechniken für das Einhalten dieser Bedingung an den Teilproblemen Orthogonales Packing - und Orthogonales Knapsack Problem vorgestellt und sowohl theoretisch als auch praktisch verglichen. Hierbei zeigt sich, dass kompakte Modelle, welche auf Relationen beruhen, die besten praktischen Ergebnisse liefern, obwohl andere Modellierungstechniken wie beispielsweise die Diskretisierung bessere theoretische Eigenschaften aufweisen. Die in der Praxis erfolgreichen Techniken werden daraufhin auf verschiedene Weisen für das Multiple Bin Size Bin Packing Problem erweitert und die Möglichkeiten diesbezüglich erneut in der Praxis verglichen. Der daraus resultierende Vergleich zeigt, dass die kompaktesten Modelle wiederum die besten praktischen Ergebnisse aufweisen.

Der Nord-Ostsee-Kanal oder auch Kiel Kanal genannt ist die am meisten befahrene Wasserstraße der Welt. Der Kanal ist in beide Richtungen befahrbar, an vielen Stellen aber einspurig, was die Routenplanung zu einer komplexen, bislang von Hand ausgeführten Aufgabe macht. Ich stelle in meiner Arbeit einen funktionierenden globalen Ansatz vor, der unter den bestehenden Lösungen im Schnitt mehrere Prozentpunkte an Gesamtwartezeit für die Reedereien einsparen kann. Es werden die Auswirkungen verschiedener Formulierungen diskutiert und zudem eine problemspezifische Vorverarbeitung und Stärkungen der Formulierung präsentiert. Ein neuer Ansatz versucht, gute Praxislösungen zu generieren und kann bei zusätzlicher Wartezeit von fünf Minuten pro Schiff eine Halbierung der Routenkomplexität bewirken. Schließlich stelle ich einen erweiterten Branch-&-Price Algorithmus für die Dantzig-Wolfe Reformulierung vor, der in Tests erfolgreich die MIP-Lösung einiger Instanzen übertrifft.

The ability to predict future demand accurately is of essential importance in supply chain management. This includes the prediction of sales in retail during promotion periods, where the number of sales differs greatly compared to non-promotion periods. The goal of this master thesis is to increase the accuracy of stock keeping unit (SKU) sales prediction models during promotion periods in a retail setting. Different promotional factors, together with interaction effects between different articles on promotion and consequences of the COVID-19 pandemic are analysed to achieve this goal. This introduces challenges in the form of high dimensionality and sparsity, which are met with a multiple linear regression model. Besides an accurate prediction, which is achieved by the incorporation of these features, this model allows to draw conclusions about causal relations between the influencing factors and the number of sales. The model is applied to data provide by a Dutch supermarket to evaluate its performance.

Location planning over a time horizon is a common task and therefore needs efficient solving methods. Experiments with temporal extensions to related problems indicate that branch-and-price algorithms might also be suitable for solving the temporal facility location problem (TFLP). This is supported by the reproduction of the results of two papers with related temporal problems, namely the temporal bin packing problem and the temporal knapsack problem. The knowledge obtained from these results is then used to create several TFLP instances with different parameter scenarios. The reproduction as well as the solving of the TFLP instances is performed with GCG and SCIP as branch-and-price and general-purpose solver respectively. In the experiments conducted in this thesis the instance parameter that has the biggest influence on the solving times is the opening cost, but having uncapacitated or capacitated facilities is another important factor. Comparing the results of GCG and SCIP on these scenarios leads to the conclusion that in most cases the general-purpose solver SCIP is the faster alternative. The branch-and-price solver GCG is useful for large instances with expensive opening cost and useful for some instances when the facilities are uncapacitated.

Machine learning methods and deep learning in particular have seen many use-cases for almost every domain in the last couple of years now. Among those is the domain of discrete optimization. In discrete optimization machine learning has been applied to heuristically improve some algorithms. Since similar models in discrete optimization can also be handled or processed similarly, we will utilize Graph Convolutional Network (GCN) to predict which problem class a Mixed Integer Programm (MIP) belongs to, such that suitable algorithm can be applied to a model. To be exact, the natural variable-constraint bipartite graph representation of MIPs is utilized to train a GCN in the task of classification. The trained GCN will additionally be able to represent any MIP by a fixed size feature vector, which can be used to calculate how similar different MIPs are to each other. The process of training the GCN will be documented and our final model will be evaluated by comparing its performance to its competition. In addition to this, some experiments will be conducted to get a better understanding of which features and which parts of the networks are important for good classification results.

This work explores if scheduled service network design programs can effectively support the Deutsche Post DHL Group (DP-DHL) in the tactical planning of their service networks. The DP-DHL operates service networks using trucks to transport packages between regional terminals. For cost-efficient transportation, packages are consolidated at hubs, which requires close coordination of the operated services. The underlying combinatorial optimization problem is known as the scheduled service network design (SSND) problem. SSND programs can decide which services to operate and how to transport commodities based on these services to minimize the total operational costs. In the case of DP-DHL, a variety of real-world constraints and operational options need to be considered during the program development. The proposed SSND program is solved by a price-and-branch heuristic, where services and paths for commodity flows are generated based on time-space graphs in column generation schemes. The program simplifies the DP-DHL problem to apply SSND programs known from literature to the DP-DHL case. Within the simplified framework, the currently operated service network under consideration could be optimized in terms of costs. The linear programming relaxation solution values show significant optimization potential, which could be partially translated into mixed-integer solutions. It can be challenging to fully model the real-world complexity and to design large and robust service networks with SSND programs. However, the improvements achieved indicate that SSND programs could be used as part of advanced solution procedures that can overcome modeling and runtime challenges.

Die Vorteile der Danzig-Wolfe-Reformulierung (DWR) beim Lösen von gemischt-ganzzahligen Programmen (MIPs) sind gut dokumentiert. Damit DWR gute Ergebnisse erzielt, muss jedoch eine bestimmte Struktur in der Nebenbedingungsmatrix des MIPs vorhanden sein. In den meisten heutigen Solvern muss der Benutzer solche Strukturen für den Solver bereitstellen. Es gibt in der Literatur einige Vorschläge zur Automatisierung der Strukturerkennung, von denen jedoch einer eine sehr wünschenswerte Prämisse hat, da er keine Benutzereingaben erfordert, d.h. keine Kenntnisse über die Problemstruktur voraussetzt. Der in der Arbeit "Structure detection in mixed-integer programs" von Khaniyev et al. vorgeschlagene Algorithmus erkennt eine bordered block diagonal (BBD) Struktur. Er schlägt eine Metrik für die Güte solcher BBD-Strukturen vor und implementiert einen greedy, auf community detection basierenden Ansatz. Diesen Algorithmus haben wir im GCG Solver implementiert, so dass wir die wichtigsten Designentscheidungen nachvollziehen sowie die Performance testen konnten.

Extracting special structures from coefficient matrices is a known practice to solve linear programs more efficiently. One of these extraction problems is the detection of maximum embedded reflected network matrices (DMERN). That is finding a maximal sub matrix in our coefficient matrix such that by multiplying a number of its rows with −1 this sub matrix contains in each column at most one -1 entry and at most one +1 entry, while all other entries are zero. The goal of this bachelors thesis is to give an introduction into the DMERN problem and provide a collection and explanation of solution approaches that can be found in the literature.

This thesis implements, documents and analyzes two different image-based approaches to detect structural similarity among mixed integer programs. The first approach detects structural similarity among mixed integer programs using the constraint coefficient matrix. Therefore, the constraint coefficient matrix is permuted using decomposition techniques of the Generic Column Generation solver (short GCG) and is visualized as an image. Several instances of strIPlib are used to create these images and then passed on to train a convolutional autoencoder. The trained convolutional autoencoder can compute feature vectors that represent latent structural features of the mixed integer program. Thereafter, the feature vectors are used to measure similarity among mixed integer programs. The aforementioned approach is extended by adding the left-hand and right-hand side vectors and the objective coefficient vector to the constraint coefficient matrix and thus creating more insightful images. These images are easily created using the developed visualization tool for decompositions in PyGCGOpt. This procedure generates our second approach to detect structural similarity among mixed integer programs. Finally, the analysis of both approaches shows that both approaches are equally well suited to detect structural similarity among mixed integer programs.

State estimation is responsible for the stable and efficient operation of electric power system and plays a vital role in power system monitoring and control. The increasing shift towards use of renewable energy sources (RES) in electric power system has introduced bi-directional flow of power in the electric power system, one from the transmission system and other from the renewable energy resources which can violate the network operating constraints which can further lead to the network failure. Therefore, state estimation in distribution networks have become very critical. Distribution System State Estimation (DSSE) algorithms compute the state variables of an electrical network complex voltages (magnitude and phase angle) using input measurements of electric bus active and reactive power injections, branch power flows and bus voltages magnitudes taken from the system. The main challenge for implementation of state estimation in distribution network is the limited number of measurements. System is generally unobservable (partially observable), therefore an attempt to solve this problem of lacking real time measurements is the employment of so-called pseudo-measurements to ensure the observability of the system and a reliable estimation. Therefore, the validity of state estimation depends on the accuracy of generated pseudo-measurements. The focal point of this thesis is to device an alternative approach to generate and model pseudo-measurement in reference to distribution system state estimation. In the proposed approach, Artificial neural networks (ANNs) are used to generate pseudo-measurements where few real measurements are used in conjunction with typical load and generation profiles. The error between the load and generation profiles (target output of ANN) and ANN output is modeled through Gaussian Mixture Model (GMM). The state estimation of three different networks are computed using the pseudo-measurements generated using the proposed methodology and compared with the actual state estimation values to validate the quality of state estimation results using the predicted pseudo-measurements. We have also demonstrated our proposed approach of pseudo-measurement modeling for state estimation on network with topological changes.

Energy grids are the backbone of every developed country. Currently, the grid topology is considered to be fixed. Topological changes are mostly executed in case of overloads or violations of safety margins. Thus, primary measures to control the power grid are market-based ones, like actively steering power plants. Recent research also identified the potential of optimizing the power flow including topological measures, like transmission element switching. The mathematical complexity of this formulation prohibits its practical usage. Currently, results are either obtained too slowly or violate critical security requirements using heavily simplified models. This work presents implementations of the traditional problem formulation and their restrictions. Alternative solution methods like nearest neighbor and reinforcement learning are implemented and compared against traditional optimization methods.

In this work, we explore combinatorial attacks on (binarized) neural networks. The "attacks" are small changes to the input data, such that the network misclassifies the input. To construct such combinatorial attacks, mixed-integer programs are used, instead of the gradient of the network, which is usually used. In the first part of the thesis, we reproduce the Iprop attack proposed by Gupta et al. in "Combinatorial Attacks on Binarized Neural Networks". Binarized neural networks are hereby aritficial neural networks with weights restricted to 0 and 1. We compare our results with those reported in the paper and analyze challenges in the implementation of the proposed method. In the second part, we extend the attacks to continuous neural networks and compare LP-based attacks with gradient-based attacks. We answer the question in which scenarios the LP-based approaches offer better performance and evaluate the strength of our formulation. Furthermore, we analyze whether successful "attacks" can be transferred to other network architectures and if LP-based attacks generalize better than gradient-based attacks.

In the course of planning railway transportation systems multiple problems have to be solved. Typically, those problems are solved sequentially in different stages, but since the planning stages are highly dependent on each other, it is beneficial to follow a more integrated approach. Inspired by the problem statement of the informatiCup 2022 "Abfahrt!", the integrated problem of creating a feasible train timetable with an associated scheduled passenger routing that minimizes the overall delay of all passengers will be investigated. This means, we direct trains and passengers individually through a public transportation network (PTN), while respecting all capacity constraints at all times. This problem turns out to be strongly NP-hard.  We propose different solution approaches based on integer linear programming and implement them to analyze their performance on examples. As a result, the approach based on a space-time graph representing actions on the original PTN outperforms the other presented concepts.

Diese Arbeit untersucht das Verhalten von Dantzig-Wolfe Dekompositionen mit Hinblick auf die Klassifizierung von Restriktionen. Anstatt für die Untersuchung alle möglichen Teilmengen der Restriktionen als Dekomposition zu betrachten, erlaubt das Klassifizieren eine Möglichkeit deutlich weniger Dekompositionen zu betrachten. Somit sind Analysen auf größeren Instanzen möglich. Es wird geprüft, ob man anhand einer Zuordnung der Restriktionen in Klassen eine Regel für das Verhalten der dualen Schranken ableiten kann. Zudem wird analysiert, welchen Einfluss die Klassen auf das Lösen der einzelnen Dekompositionen haben, insbesondere auf die Lösungsdauer.

Collection and Dispute Management, a part of the broader Order to Cash process, is responsible for collecting the credit back from the customer and improving the cash flow. Bayer's collections process suffers a severe overdue problem, owing to delayed payment by the customers. The current reactive nature of the collection process poses a significant challenge for reducing the average delinquent days. Hence a proactive collection strategy with a specific collection action for each customer is the need of the hour. This thesis focuses on understanding the business problem and providing an analytical solution for the collection team of one of the country. The solution employs Machine Learning via two modules: Invoice prediction and Days to Pay prediction. Finally, the Predictive Worklist module is the decision model that ensembles the results from the other two modules and the behavioral characteristics of the customer to create a recommended action for each customer. A detailed analysis of the problem and performance of those modules are evaluated. The proposed model serves as base for further extension of the project for other countries.

The accurate determination of the measurement uncertainty is important in the production environment to avoid economic losses due to erroneous decisions caused by measurement errors. Prerequisite for determining the measurement uncertainty is the description of the measurement process with the help of a model. Artificial neural networks (ANN) have the potential to create a model that validly describes the measurement process. However, there are no guidelines for the selection of the hyperparameters and topology of the ANN that are used in this context. In this master’s thesis, a method for optimized hyperparameter and topology selection for model building in measurement uncertainty determination using ANN was developed. The optimization of the hyperparameters and topology of ANN was first formally described as a mathematical optimization problem. To solve the problem, different algorithms (including classical algorithms, genetic algorithms) were identified and adapted for the problem of modeling in the determination of measurement uncertainty. Finally, the methods were evaluated using several measurement data sets.

In this Master Thesis, we investigate Mixed Integer Programs (MIPs) for optimally solving the problem of covering rectilinear polygons with a minimum number of axis parallel rectangles. Since larger instances involve millions of possible rectangles inducing long solving times, it might be necessary to choose a MIP formulation that can be solved with column generation. This allows us to start with a basic feasible solution and dynamically add new rectangles improving the current solution. In general, solutions to this NP-hard problem find applications in the fabrication of DNA chip arrays, in VLSI design, and image compression. In particular, aiming at image compression, other shapes like ellipses will be tested to improve the compression rate. Additionally, we investigate an approach that allows the model to select the shapes on its own to achieve even higher compression rates. Furthermore, the MIPs can be extended to allow for a lossy image compression. This thesis will elaborate and implement the introduced approaches to test them for their computational feasibility.

The nesting problem belongs to the class of cutting and packing problems. A two-dimensional cutting problem is called a nesting problem when the two-dimensional objects that need to be cut are non-rectangular. The objective is usually to minimize waste by utilizing the available material in the most efficient manner. Nesting problems are often encountered in manufacturing industries such as wood, steel, and textile industry. The geometric properties of irregularly shaped objects make solving nesting problems in practice very challenging. In this thesis, we cover the existing approaches to solving nesting problems as well as the geometric tools utilized by these approaches. Additionally, we explore a methodology combining two different strategies, namely, a heuristic algorithm and a mixed-integer programming model.

Reduced cost variable fixing is a classic technique for reducing the problem size when solving mixed-integer program. It allows fixing variables whose reduced costs are higher than the absolute gap. Reformulating problems with Dantzig-Wolfe decomposition and solving them with branch-price-and-cut algorithms often leads to tight dual bounds and smaller absolute gaps, but application of RCVF is not straightforward in these algorithms. Fixing the variables in the reformulation is not very promising but fixing variables of the original problem is not immediately possible, as no reduced costs are known for them. This thesis explores various ways to compute reduced costs of original variables and therefore apply RCVF in BP&C both from a theoretical and computational perspective.

In dieser Arbeit untersuchen wir verschiedene Relaxationen für ganzzahlige lineare Programme und betrachten dabei insbesondere deren Anwendung auf das Stable Set Problem. Wir analysieren Dantzig-Wolfe- Reformulierungen, die Corner-Relaxation und deren Variationen, die Sherali-Adams-Relaxation und die Lovasz-Schrijver-Relaxation anhand folgender Gesichtspunkte: Wie genau ist die Relaxation? Kann man die konvexe Hülle bestimmen? Wie hoch ist der Zeitaufwand? Wie hängen diese Faktoren möglicherweise von der Struktur des Graphen ab? Kann man das Problem auf ein Teilproblem reduzieren? Kann man verschiedene Verfahren kombinieren? Kann man die Verfahren auch auf andere Probleme anwenden?

Texts written in natural language are an unstructured data source that is hard for machines to understand. The amount of text in the world wide web is growing every minute. To deal with this huge number of unstructured data automated text analysis is crucial. Natural Language Processing (NLP) is part of artificial intelligence that makes natural language texts comprehensible for machines. Natural language processing (NLP) techniques can contribute to this problem by offering automated means to do preprocessing, text classification, feature extraction, and topic modelling. Social networks, like Facebook or Twitter, Application like Amazon or eBay are a phenomenon that has recently transformed numerous aspects of our lives. The impact of these platforms is no longer limited to entertainment or personal presentation of individuals. Indeed, they have formed a new sheer of business, with some emerging companies focusing primarily on this platform and others, more traditional ones, expanding their marketing efforts there. My intense research is being performed in order to efficiently categorize, filter, detect text and reviews on the customer-generated content on social networks and application reviews. One of the key segments of this effort is sentiment analysis and goal is to create different models on sentiment of the posts, comments, reviews or other forms of reactions users generate in relation to a company or a product. On the basis of that checking computational performance and results.

Um die Frage zu beantworteten, ob sich die Investition in ein bestehendes Produktionssystem lohnt, muss man beide Konfigurationen anhand von Kenngrößen miteinander vergleichen. Dazu wird die Produktionsstruktur als Maschinen, welche über Puffer miteinander verkettetet sind, abstrahiert. Dabei fangen die Puffer Maschinenausfälle ab, damit nicht weitere Arbeitsschritte auf Grund von Materialmangel oder mangelnder Kapazität im Output stillstehen. Ein Puffer kann zum Beispiel ein Förderband sein, auf dem sich Zwischenprodukte aufstauen können. Desweitern soll das Modell auch konvergente und divergente Produktionsströme abbilden, das heißt Maschinen können mit mehr als zwei Maschinen verbunden sein und somit zum Beispiel eine Endmontage beschreiben. Zur Bewertung verschiedener Konfigurationen kann man die Ausbringung als Kenngröße verwenden, welche die Produktionsmenge bezüglich einer Zeiteinheit beschreibt. Zudem ist der Wirkungsgrad, also die Ausbringung im Vergleich zur kleinsten Produktionsrate einer Maschine, interessant. Zur Messung dieser Werte in der Realität, müsste allerdings die Produktionsstruktur umgebaut werden, was mit erheblichen Kosten verbunden. Somit muss auf eine Simulation oder ein numerisches Verfahren zurückgriffen werden. In dieser Arbeit werden Algorithmen, zur Bewertung von konvergenten und divergenten Produktionsstrukturen aus der Literatur erläutert. Der Input ist ein Netzwerk aus Maschinen, wobei diese durch Kennzahlen zur Produktionsrate, Ausfallzeiten und Reparaturdauer beschrieben werden und die Puffer über ihre Kapazität definiert sind. Des Weiteren wird der Fokus auf der Anwendung dieser Algorithmen in der Praxis liegen, insbesondere im Vergleich zur Simulation. Dazu wird ein Algorithmus prototypisch umgesetzt. Zu Beginn der Arbeit werden die Ansätze der verschiedenen Algorithmen aus der Literatur erläutert. Im Anschluss wird die Vorgehensweise der Algorithmen erklärt, dabei wird der mathematische Hintergrund beleuchtet. Der letzte Abschnitt beschäftigt sich mit der Güte der Bewertung des implementierten Algorithmus und mit seiner Laufzeit.

Recent times have seen the rise of a research current which aims to solve classical Combinatorial Optimization (CO) problems through the help of Machine Learning (ML), and specifically in this work Deep Reinforcement Learning (RL) methods. One of the problems where these two exciting fields meet one another tackles the question of how to guide the branching process within a Branch and Bound algorithm, where sequential decisions must be done in order to reach an optimal solution. To this extent, the goal  of the current work is to gain insight of the state of the art methodologies dealing with such a problem, and through an implementation and an experimental study, to provide an comprehensive view of the impact of employing RL to aid the solving of CO problems.

Aktuell sind Klimaschutzmaßnahmen in vielen Sektoren von großem gesellschaftlichem und politischem Interesse, um den anthropogenen Treibhauseffekt einzudämmen. Der Verkehrssektor nimmt dabei eine entscheidende Rolle ein, da dieser zu den weltweit größten Verursachern von Treibhausgasen zählt. Im Rahmen der Mobilitätswende wird angestrebt, in den kommenden Jahren die bisher überwiegend verbreiteten Verbrennungsmotoren durch Fahrzeuge mit hybridem oder vollelektrischem Antrieb weitestgehend zu ersetzen. Damit einher geht der Ausbau der nötigen Ladeinfrastruktur und Entscheidungen bezüglich der Standortwahl und Größe von Ladestationen. In dieser Masterarbeit wird das Problem basierend auf dem Temporal Facility Location Problem mit Methoden der diskreten Optimierung mathematisch modelliert. Dieses stellt eine zeitliche Erweiterung des klassischen Standortproblems dar, in der Kunden jeweils nur für ein gegebenes Zeitintervall an einen Standort angebunden werden. Die neue Restriktion lautet, dass zu keinem Zeitpunkt mögliche Kapazitäten der Standorte überschritten werden dürfen. In einer detaillierten theoretischen Analyse werden zunächst verschiedene Varianten des grundlegenden deterministischen Temporal Facility Location Problems untersucht. Neben einer jeweiligen Komplexitätsanalyse werden beispielsweise Untersuchungen bezüglich der Approximierbarkeit durchgeführt sowie verschiedene mathematische Formulierungen hergeleitet und verglichen. Daraufhin werden die Probleme in ihren Online-Varianten untersucht und Resultate bezüglich der Kompetitivität vorgestellt. Anschließend werden realitätsnahe Charakteristika des Anwendungsbereichs in die Problemdefinition integriert, indem verschiedene Modifizierungen vorgenommen werden, und erneut eine detaillierte theoretische Analyse durchgeführt.

In this thesis we take a look at the variety of graph spanner problems in the current literature and some exact mixed integer programming approaches to solve these. We present five spanner types used in problem definitions, the path formulation and multicommodity flow formulation we use to solve them. We run an experiment using the Erdős-Rényi random graph model as well as the Watts-Strogatz small world model to analyze the behavior of different combinations from spanner types, formulations, and graph models. We find that pairwise and sourcewise spanners prove slower to solve, while subsetwise spanners reveal themselves to be fastest to solve of the spanner types. The random graph model turns out to provide faster running times than the small-word model and the multicommodity flow formulation provides faster model building and subsequently faster overall running times than the path formulation. As part of the experiment we generate a library of input instances for the spanner problems and graph models. Since randomness is involved in the generation of the instances, this library helps in the proper analysis and comparison of solution approaches and easily allows for further expansion on it.

Mixed integer programming is a powerful tool to model and solve many modern optimization problems like network design, frequency assignments, or chip verification. As solving mixed integer programs, or MIPs in short, is NP-hard, many modern solvers rely on primal heuristics to assist the solving process of the branch-and-cut algorithm. These heuristics have the ability to shorten the solving process by using the given information to find high-quality feasible solutions and are therefore very important when solving MIPs of higher complexity. While modern solvers, such as the open-source solver SCIP, have already a wide range of implemented heuristics, there is a great interest in heuristics with new approaches due to increasing MIP complexity. This thesis is going to describe the new heuristic "Feasibility Jump", introduced by Luteberget et al., and integrate it into SCIP. This heuristic tries to find a high-quality feasible solution by iteratively setting one variable to a value where the most constraints are satisfied. Evaluation of the impact feasibility jump has on the solving process of SCIP yields positive results regarding solving time or primal integral while having no significant effect on the solving process of MIPs where it does not find a solution.

Knowledge about structure in an Integer Linear Program can help to improve the solving process, such as a block-diagonal structure by enabling a Dantzig-Wolfe decomposition. While GCG, a solver based on SCIP, can automatically detect structures in Integer Linear Programs, other commercial and non-commercial solvers require the user to enter the structure information manually to utilize the benefits of a Dantzig-Wolfe-Decomposition. Even though one of the main factors for a performant decomposition is to find a good structure within a problem, we are interested in the quality of the decomposition algorithms themselves. We compare and evaluate which solver performs best on certain problem instances providing all with the same information about the problem itself and its structure.

Dantzig-Wolfe decomposition is a popular technique for solving (mixed integer) linear programs where a large subset of constraints can be divided into independent subprograms. In such a program's coefficient matrix, the rows and columns corresponding to the subprograms can be reordered to form a block-diagonal. A similar structure is the staircase, where neighboring blocks are linked by columns or rows. Staircases can appear in coefficient matrices of linear programs modeled over a time horizon where each period is linked to its predecessor and successor. If the periods are linked by columns, Dantzig-Wolfe decomposition can still be applied by utilizing Lagrangian decomposition. To apply it, the periods' constraints and the existence of a staircase structure either must be known beforehand or a staircase must be identified in the coefficient matrix. Since the constraints are likely in the wrong order to identify the structure by just looking at the matrix, algorithms must be used for the detection. In this thesis, an algorithm is presented that was designed to detect staircases by joining rows to blocks via agglomerative hierarchical clustering. It was tested on five temporal (mixed) integer programs and its results were compared to those of two detection algorithms that have already been part of the solver Generic Column Generation. The proposed algorithm was the best staircase detector overall but still showed some flaws, the major one being its high computation time. Nevertheless, it demonstrates that block-diagonal detectors can be utilized for the detection of staircases. At the same time, one of the two existing staircase detectors demonstrated potential. For two of the temporal programs, staircase structures were used for solving instances from the literature. In one of them, they had a very negative impact on the solving times. In the other, the differences in visual quality of various staircases had no obvious impact.

This thesis presents a new web-based platform for the structured Integer Program library, strIPlib, a collection of 21,000 mixed-integer program (MIP) instances with an exploitable problem structure. Each instance is annotated with more than 150 metadata attributes that describe its structure, characteristics, and the algorithmic behavior of the Generic Column Generation (GCG) software while solving it. The metadata is beneficial for compiling test sets to evaluate the algorithmic components of solvers in isolation or together by selecting instances that closely match the objective. However, the current website that provides access to the instances and metadata has limitations in presentation and functionality, such as a lack of in-depth categorization of instances or missing download and filter functions. These limitations hinder researchers from finding instances that match their research question and accessing the raw linear program files. The new user interface displays instances hierarchically, reflecting problem types and contributors. Each of the 150 metadata attributes has been made filterable with graphs showing their distribution and correlation. After filtering, a test set containing instance and structure files can be compiled and downloaded. The integration of GCG reduces the entry barrier and allows researchers to detect and visualize problem structures on the website with one click. Due to the maintainable design, the platform is prepared to integrate other existing projects, such as graphics for visual evaluation of GCG runtime data or a tool for reducing the number of instances after filtering while maintaining the overall diversity.

Abstract: Cluster analysis of statistical data sets attempts to identify similarity structures in data sets in order to be able to combine individual statistical units into groups in a meaningful way. Data sets can be very large and contain features that are either redundant or irrelevant for grouping. The selection of the relevant features is a fundamental problem of cluster analysis. In the paper [Benati et al., 2018] 2 mixed integer models are presented, which calculate the problem of the selection of relevant features. Test results that have been done in the paper solving these models with CPLEX show that both formulations have promising results with well separated clusters, but have their limits for problems where clusters overlap or features are not unique. In this bachelor thesis the two mixed integer programs were solved with the original data sets using the framework GCG [Gamrath and Lübbecke, 2010], which follows a branch price and cut approach and were compared with the results from the original paper. Results show that although the underlying structure of the problem is promising, solving the instances with GCG is slower by orders of magnitude.

Power transmission grids face increasing stress due to changes in power production and consumption patterns. The shift towards renewable, usually less steady energy sources or the deployment of "smart" devices have their contribution to this process. Currently, not all degrees of freedom in the power grid operation are exhausted; thus, the operation is sub-optimal concerning congestion and costs. While the power dispatch is optimized, this is not the case for topology measures, mainly due to computational burdens. This thesis analyzes the issue of using topology measure optimization integrated into an optimal power flow problem, the so-called Security Constrained DC Optimal Transmission Switching (SC-DC-OTS) problem. An integrated problem formulation of a Power Transfer Distribution Factor (PTDF) based SC-DC-OTS is provided together with its implementation. Further, an offline algorithm is proposed to find topology measure candidates efficiently, thereby reducing the solution space of the SC-DC-OTS model, the so called Power Set Poisoning Depth-First Search (PSPDFS) algorithm. In parallel, the implementation explores improvements to the ergonomics of modeling tasks. Four simulation studies are conducted to evaluate the functioning and performance of the model implementation, as well as the potential of the PSPDFS algorithm. Considering topology measures and security constraints leads to a vast amount of side conditions in the model, and thus, considering all potential topology measures is not manageable. The simulations show that in some cases, PSPDFS can considerably reduce the number of topology measures.

Today, machine learning methods based on statistical learning theory are enjoying growing interest in both research and applications. Simultaneously, operations research, in particular, linear programming, experiences a renewed surge in popularity, providing solutions that are guaranteed to be optimal for a wide range of problems. As part of the efforts made to integrate those fields for mutual benefit, recent studies have attempted to apply machine learning to solutions and non-solutions of linear programs. The goal here is to obtain a linear programming instance from such data that describes the problem as accurately as possible. In practice, data, i.e., solutions and non-solutions, can often be collected independently of and much more easily than modeling the problem. However, since modeling a problem as a linear program needs domain as well as methodological knowledge, skilled experts are hard to find and thus also expensive. As a step in the direction of automatic modeling, learning linear programs could therefore be of great interest for everyone applying linear programming. Therefore, in this thesis, an exploratory study of a learning algorithm is conducted with respect to the following questions: Can reasonably good results be obtained by the algorithm on synthetic linear programs obtained from different sampling methods? Can the results obtained be transferred to real-world known instances? Of what quality are the obtained solution instances? And what properties must therefore be met by the examples used? The core of this work is an implementation of the algorithm to verify and investigate these and other aspects through various experiments, and compare it to benchmark algorithms.

In this thesis, we take a look at tackling the challenge of solving complex combinatorial optimisation problems using a combination of machine learning and combinatorial optimisation techniques. As modern-day problems become increasingly intricate, it is essential to develop efficient and effective solutions to address them. We focus on the task of predicting shortest paths in the context of the Warcraft Shortest Paths dataset, a scenario where traditional machine learning approaches may not be sufficient for accurate predictions. To begin with, we look at an overview of the fundamental concepts of machine learning and combinatorial optimisation and investigate various architectural approaches, including machine learning-based architectures, two-stage approaches utilising both machine learning and combinatorial optimisation separately, and combined pipelines. To evaluate the effectiveness of these approaches, I apply them to the Warcraft Shortest Paths dataset and analyse their performance using relevant metrics such as loss and gap. The main contribution to this topic includes a comprehensive comparison of these three approaches applied to the same dataset, offering insights into their advantages, disadvantages, and potential improvements. This analysis also includes a brief discussion regarding the importance of the metrics used for evaluating different architectures, seen through the lens of the optimality gap and losses.

Facility layout design, the efficient arrangement of functional units inside a facility, is a crucial part of factory planning. Despite extensive research on the topic, mathematical optimization approaches are still underutilized in practice, partly due to differences between the facility layout problem and facility layout design. This study addresses both aspects, beginning with the development of an efficient mixed integer programming (MIP) model. The effect of different modeling choices is examined, and a new linear area approximation via non-equidistant support points is introduced. The improved formulation allows operators to set a maximum area deviation and distributes the deviation equally. Much of the problem’s computational complexity comes from the non-overlapping constraints. Two different sets of constraints are presented and compared. Additionally, the importance of a two-stage approach is established, and two solution approaches are proposed. A collection of common test instances is introduced based on which a set of computational tests are performed. The tests show that the MIP approach can efficiently solve small to medium-sized instances with up to 14 functional units. Test instances with 20 to 60 functional units were solved more efficiently by the heuristic. In addition to theoretical test instances, this research includes a real-world planning scenario from the industry. Both solution approaches were successfully adapted to solve the industry test instance. However, the heuristic tends to generate mostly infeasible layouts. This raises the question of whether the heuristic can provide any benefit for solving the facility layout design in practice. Despite this, the overall approach made it possible to improve both the mathematical optimality of the solution and also the applicability to the real-world scenario. Nevertheless, the newly generated layouts exhibit less structural transparency than the original layout, and further fine-tuning and acceptance efforts may be required to ensure the adoption of the solution.

Binarized Neural Networks (BNNs) are feed-forward neural networks with fully connected layers that purely use binary weights. The evaluation is highly efficient, which makes it attractive to use them in low-power settings since their performance is close to the performance of linear neural networks. However, robust networks are important, especially in settings that are safety-critical. So-called attacks on neural networks are small perturbation vectors that are added to the input that lead to different (pre-defined) results. Existing methods to generate attacks mostly use gradient-based methods, but the binary, non-differentiable nature of the BNNs renders gradient-based methods impossible. One more recent approach already makes use of mixed integer linear programming (MILP) models. In this body of work, we reproduce the results of existing approaches that make use of MILP models. We come up with new formulations and we show that they are stronger performance- and results-wise. Furthermore, we show that a reformulation of these models into a column generation approach does not seem to be practical and demonstrate an MILP model to alter the network to be robust against the previously generated attacks.

Many challenges in the industry and economy, e.g. production planning, vehicle routing, and territorial partitioning, can be modeled and solved with Mixed Integer Programming (MIP). MIP is a powerful modeling technique for mathematical optimization and feasibility problems. Solving an MIP problem is generally NP-hard. MIP solvers like Solving Constraint Integer Programs (SCIP) and Gurobi are based on the branch and bound algorithm. MIP solvers use a technique called primal heuristics to find feasible solutions, which can lead to the early pruning of the branch and bound tree. This approach can accelerate the solution process and, in some cases, even enables the discovery of an optimal solution. Examples of primal heuristics are rounding heuristics, diving heuristics, and large neighborhood search. A variant of the branch and bound algorithm is the branch and price algorithm, which exploits a reformulation of the original MIP via Dantzig-Wolfe Decomposition (DWD) to solve the optimization problem. Until recently, implementations of the branch and price algorithm used to be specific to a particular problem and its structure. Current state-of-the-art research and development enables a generic way to reformulate a given MIP and the corresponding structure into the DWD and solve the problem with branch and price automatically. An example of such a solver is Generic Column Generation (GCG), which is based on SCIP. It automatically detects, reformulates and solves a given MIP with branch and price. Primal heuristics, when applied in branch and price algorithms, may not fully exploit the reformulation’s structure in their generic form, or they might be customized for specific problems. Recent research explores the adaptation of generic heuristics to the branch and price algorithm in column generation-based MIP solvers. The thesis will investigate the effectiveness of using primal heuristics in mixed integer programming solvers based on the branch and price algorithm. This involves testing different primal heuristics to determine which ones result in the most efficient and accurate solving process, as well as analyzing the computational performance of the solvers when employing these heuristics. The goal of this research is to provide insight into the potential benefits and drawbacks of using generic primal heuristics in MIP solvers and to ultimately improve their effectiveness in solving complex optimization problems.

Dantzig-Wolfe decomposition is an algorithm that considers an integer program and convexifies a subset of constraints. Solving the resulting master problem leads to potentially stronger dual bounds than the respective linear programming relaxation. As the subset can be chosen arbitrarily, this includes the trivial cases of convexifying no and all constraints, resulting in a weakest and strongest reformulation, respectively. The strength of such reformulations is expressed through the associated dual bound and can impact the performance of corresponding algorithms. We repeat available measures on the strength of DW reformulations from the literature. Next we discuss computational results of experiments on the strength presented among others by Bastubbe, Lübbecke and Witt (2018). They observed that some classes of constraints, as defined by MIPLIB, have less to no effect on the dual bounds of DW decomposition when they are reformulated without other classes of constraints. This situation is called ‘level 0’ and means that only one particular constraint class is convexified. These observations are proven for the set-partitioning, set-packing, set-covering, cardinality and invariant knapsack classes. We can show that they never improve the dual bound when they are chosen to be convexified without other types of constraints. Next the classical edge formulation for the stable set problem is analyzed regarding its strength. We characterize weakest and strongest possible Dantzig–Wolfe reformulations and present a detailed example that illustrates the results.

The hub location problem is a mathematical optimization problem of determining the optimal locations for hubs with the objective of minimizing overall costs by routing traffic between customers through determined hubs. Hub location problems are in general NP-hard. A possible heuristic to receive solutions for large-scale instances is to contract the network of such an instance, drastically reducing the computation time needed for solving the resulting small-scale instance. The solution of the contracted network can then be rewritten to the original network, leading to a feasible solution of the large-scale instance. In this thesis, we explore several contraction and rewrite methods and compare the quality of their solutions.

Die Lösung von NP-vollständigen Problemen – wie die Erfüllbarkeitsüberprüfung von SAT-Formeln (SAT-Probleme)– kann mithilfe von Heuristiken beschleunigt werden. Außerdem ist es beim Lösen von SAT-Problemen ein verbreiteter Ansatz, das Hauptproblem in Teilprobleme aufzuteilen (divide-and- conquer). Dabei ist Konfliktanalyse, in den Teilproblemen, die einen Konflikt enthalten, ein essenzieller Bestandteil aktueller SAT-Solver. In diesen werden Informationen in Form neuer Konfliktklauseln vom Solver abgeleitet. Diese Form der Klauseln wird verwendet, um bekannte Konflikte festzuhalten und in folgende Berechnungen miteinzubeziehen. Darüber hinaus sind SAT-Problemen ein spezielles Fall vom gemischten ganzzahligen Optimierungsproblemen (MIPs). Zudem ist der Branch-and-Bound (B&B)-Algorithmus, der auch ein divide-and-conquer Prinzip folgt, ein essenzieller Bestandteil der aktuellen MIP-Solver. Allerdings werden dabei aufgetretene Konflikte ignoriert. In der Arbeit „Conflict analysis in mixed integer programming“ wurde ein Ansatz entwickelt, bei dem auch das Teilproblem, das einen Konflikt enthält, untersucht wird. Dazu wurde eine Methode vorgestellt, die die Methoden der Konfliktanalyse von SAT-Solvern auf die Lösungsverfahren von MIP überträgt. Dieses Generalisieren der Konfliktanalyse vom SAT-Solving auf MIP-Solving hat in den zum Testen verwendeten MIP-Instanzen zum einen eine Beschleunigung des MIP-Solvings von unzulässigen Problemen gezeigt, zum anderen wurde der Suchbaum bei lösbaren Problemen zwar verkleinert, aber die Gesamtlaufzeit hat sich vergrößert. Des Weiteren existieren mehrere Strategien zum Teilen eines gegebenen Optimierungsproblems mit B&B, um einen Lösungsbaum zu erstellen. Sowohl die Rechenintensität bzw. die gesamte Laufzeit als auch die Größe des resultierenden Lösungsbaums sind entscheidend, um eine bestimmte Strategie auszuwählen. Demzufolge ist es ein Aspekt ML in B&B zu integrieren, diejenigen Teil der Strategien, die zwar rechenintensiver sind, aber einen Lösungsbaum erstellen können, der weniger Knoten enthält und somit zu schnelleren Laufzeit als anderen Strategien führen kann, durch schnelle Approximation zu ersetzen. Der Schwerpunkt dieser Arbeit ist es, an die Idee der Konfliktanalyse beim Lösen von MIP aufzusetzen und zu untersuchen, wie eine Methode zur Anwendung von Machine-Learning (ML) in B&B Perspektiven zur Anwendung von ML für die Konfliktanalyse in MIPs liefert und einen Ansatz zur Adoption dieser Method präsentiert, damit die Konfliktanalyse von MIPs mithilfe von ML durchgeführt werden kann.

Laser Powder Bed Fusion is an additive manufacturing technique for 3D laser printers. Hereby the 3D object gets divided into multiple 2D layers. The 2D object on the layer then gets completed through Laser Powder Bed Fusion. We look at a special variant with a moveable print head and multiple lasers in the print head. Which creates a new challenge in how to find a good pathing to finish the 2D layer fast. We create an integrated model of the path planning and then try to make the whole model more efficient, through changing the objective function parameters and different methods like Dantzig-Wolfe decomposition. When it comes to optimizing the idle time of the lasers we created a model for object with up to 13.720 scan vectors

This master thesis integrates the component bound branching rule, originally proposed by Vanderbeck et al. and later reformulated by Desrosiers et al., into the branch-price-and-cut solver GCG. This rule, similarly to Vanderbeck's generic branching scheme, exclusively operates within the Dantzig-Wolfe reformulated problem, where branching decisions generally have no corresponding actions in the original formulation. The current GCG framework requires modifications for such branching rules, especially within the pricing loop, as seen in Vanderbeck's method implementation. These rules also fail to utilize enhancements like dual value stabilization. A significant contribution of this thesis is the enhancement of the GCG architecture to facilitate the seamless integration of new branching rules that operate solely on the reformulated problem. This allows these rules to benefit from current and future advancements in the branch-price-and-cut framework, including dual value stabilization, without necessitating alterations to the implementation of the branching rule itself. The thesis proposes an interface to manage constraints in the master problem that lack counterparts in the original formulation. These constraints require specific modifications to the pricing problems to ensure their validity in the master. The generic mastercut interface, tightly integrated into the GCG solver, spans the pricing loop, column generation, and dual value stabilization. Enhancements to the existing branching rule interface complement this integration, enabling effective utilization of the generic mastercuts. Finally, the component bound branching rule will be implemented using this new interface and evaluated on a set of benchmark instances. Its performance will be benchmarked against the existing Vanderbeck branching rule, offering a practical comparison of both approaches.

In the process of solving mixed-integer linear programs state-of-the-art solvers construct a so- called conflict graph. The conflict graph stores information about pairs of assignments to binary variables that lead to infeasibility. Commonly such structures are used to obtain clique cuts during separation. However, all inequalities valid for the stable set polytope of the conflict graph can be used to strengthen the linear relaxation of the original problem. In this thesis, we will experimentally investigate the potential strength of the conflict graph. We do so by applying a Dantzig-Wolfe reformulation with convexification. Further, we also investigate the effect of obtaining the whole conflict graph by a probing procedure, as the conflict graph is usually only partially discovered. Here we notice that for some types of problems we can substantially improve the dual bound by considering the whole conflict graph. Additionally, we evaluate the strength of a binary relaxation that is obtained when each constraint is relaxed such that all non-binary variables are eliminated. This shows us that in some cases we do need to consider non-binary variables to obtain all edges. We also apply convexification on separated cuts and see that this often does not lead to a substantial improvement. At last, we discuss further ideas for experiments and to improve techniques based on the conflict graph.

In linear optimization, staircase structures can be found in the coefficient matrix of a mixed integer program by the permutation of its rows and columns. The structure consists of independent blocks along the diagonal, where two consecutive blocks may be connected by variables. Structures in the coefficient matrix of large scale linear programs can be exploited to solve the program more efficiently by procedures like Dantzig-Wolfe decomposition. For this, the structures first have to be detected. However, there are different properties in a staircase structure to look for, as the number of stairs, the number of linking variables and the stair heights. Therefore it is difficult to compare existing detection methods and results. This work gives and overview of different staircase structure variants that highlight different properties of the structure. Different methods for detecting these structures are considered, namely the Longest Shortest Path algorithm, the Clustering Heuristic, and a MIP formulation. The efficacy of these methods is demonstrated through computational experiments and a score for evaluating the quality of staircase decompositions is proposed.

Due to continuously growing metropolitan regions, the requirements for ground- based transport systems are increasing. With increasing electrification, there is a potential of sustainable air mobility systems to relieve existing transport systems and to complement the existing transport networks. This potential of Urban Air Mobility (UAM) is based, among other aircraft, on electrically operated vertical take-off and landing (eVTOL) air taxis. In addition to the concepts, operational strategies have already been investigated in previous research work. For this pur- pose, agent-based approaches provide a relevant basis for the simulation of realistic demand data. In this context, the pooling of passengers into the aircraft with dif- ferent passenger capacities, ranges and loading times was insufficiently analyzed, so that a comparison of the market potentials could be distorted. For this purpose, two potential solutions were investigated: the transition to a new agent-based modeling (ABM) platform and the development of the post-processing tool to the existing de- mand data generated by a simulation in MATSim. The current state of the scientific literature on ABM approaches shows that the platforms OPUS and RePast4Py can be seen as alternatives to MATSim, but require a lot of effort for an implementation. Consequently, a linear optimization model is relevant as a post-processing tool and should be modeled in the context of current operations research (OR) approaches. For this purpose, a second literature research was made in order to compare similar past models and to learn the approaches of the operation research from them. Using the OR methods, the mathematical model was then created and fully described. The development of the mathematical model includes the following aspects: changing of an aircraft by passengers, maintaining passenger capacities and ranges of an aircraft, tracking the battery charge state, dynamic demand. Subsequently, the model was validated and some parameter studies were made. The results of the parameter stu- dies show that passenger pooling is best possible with large passenger flows, but is limited in the context of UAM. However, possible pooling can be optimal and inclu- de the Aurora PAV and Joby S4 vehicles as a total fleet. Evaluations of the runtime show limitations of the application on overall scenarios of UAM and motivate further research with regard to algorithmic optimization.

The thesis investigates many heuristics used to speed up to the pricing in Column Generation by applying these in GCG with different implementations and values on a representative test set. Partial pricing, multiple pricing and heuristic solving of subproblems have been extensively tested, yielding some interesting takeaways on how and when (not) to use them. A newly implemented random subproblem selection technique for partial pricing gave a considerable speedup compared to the other implemented selections in GCG. Also, acceleration techniques like the column pool have been checked for their effectiveness, revealing that it is regularly not beneficial. Furthermore, a finding regarding cutting stock instances formulated as set-covering is reported, speeding up multiple instances extraordinarily. Last but not least, it is reported which configurations worked particularly well for each problem type, and indications are given which heuristics are worth modifying.

Transporting goods plays a crucial role in global trade and e-commerce. This introduces the need for complex transportation networks, which must be efficient and cost-effective. Hub Location Problems model hub-and-spoke networks with the most common objective being the minimization of the total costs of the network. The Service Network Design and Hub Location Problem is a variant of Hub Location Problems that models trucks and their routes on the hub-to-hub connections. However, this problem is challenging to solve with an exact solver for larger instance sizes in a reasonable timeframe. Therefore, we propose a General Variable Neighborhood Search heuristic with the goal of obtaining high-quality solutions for larger instance sizes of the problem in a reasonable timeframe. Our heuristic is able to improve most initial solutions of the problem in a short timeframe, leading to near-optimal, sometimes optimal, solutions for smaller instances. For larger size instances, we were able to obtain solutions beating those of the exact solver for a set maximum duration in a fraction of the time.

Solder Paste Printing (SPP) is an important process in surface mount technology (SMT) for electronics manufacturing and there are two critical steps involved the process:- (i) exchanging the process tools, namely stencils at regular intervals of time during the process and (ii) blocking/replacing the stencils when they cannot be used anymore, both to ensure high quality printed circuit boards (PCB) after the process. However, the major challenges involved here are to perform these steps in such a way that the quality of prints is not affected and at the same time the stencils are utilized to the right extent, and finally the costs or risks of product scrap or tool exchanges are minimized. Hence, it is a crucial step to identify the feasible idle time after which the stencils can be exchanged during the SPP process, and block and replace them to avoid risks during future production. In an attempt to reduce product scrap rates and internal defect costs (IDC) and increase global competitiveness in future PCB manufacturing, this thesis aims to explore and identify the best decision-making approaches for the solder paste printing process based on historical data and determine a feasible idle time for tool exchange and an end-of-life period of stencils through data analytics techniques.

In computational experiments, Bastubbe, Lübbecke, and Witt observed that for a variety of MIP instances the spectra of dual bounds resulting from all their respective Dantzig-Wolfe reformulations contain only a few distinct values. In contrast, we prove that the polyhedra from reformulations are asymptotically almost surely different in the case of the stable set problem on complete graphs. Nevertheless, our experiments illustrate that even in this case there are only a few distinct dual bounds. As a byproduct, we provide a complete classification of when two reformulations of this problem differ, extending previous theorems of Lübbecke and Witt that classify the weakest and strongest reformulations of stable set problems. We also present two general results on the strength of Dantzig-Wolfe reformulation, which generalize Geoffrion's necessary condition.

In this thesis, we look into possible effects of presolving algorithms on mixed-integer program decompositions, with a specific focus on column generation using the GCG solver. We systematically analyze the effects of a range of presolvers on the generation of problem decompositions and their selection through GCG’s implemented scoring method, as well as on the runtime of the selected decomposition and its semantic structure. Based on this, we argue that certain presolvers are decomposition-preserving, meaning they will make it possible to re-detect a decomposition after presolving, that was detectable before presolving as well. Among the presolvers examined, most show only minimal or negligible effects on decomposition detection, with some such as trivial and tworowbnd presolvers demonstrating significant impacts on runtimes for specific instances. Notably, presolvers like convertinttobin and boundshift were found to preserve decompositions, where decompositions with similar structure and semantics were detectable before and after presolving. Further, our analysis of presolvers domcol and sparsify showed a larger impact on the amount of overall generated decompositions as well as for the possibility of aggregating blocks into similar pricing problems within the detected decompositions. This is especially notable for instances of the map labeling problem and the container loading problem. We conclude, that the results of our experiments, even though interesting in certain problem categories and niches, are based on a limited dataset. As demonstrated in our analysis, the majority of instances remained unaffected in terms of decomposition detection, fueling an opportunity for further research to look into, e.g. larger and more diverse testsets including problems from the MIPlib.

In order to minimize generation costs in power grids, linear optimization can be used. Because of the occurrence of Braess’ paradox in power grids, the optimal transmission switching problem (OTS) optimizes over decision variables indicating whether a branch is switched on or off. Modeling the internal structure of substations is called optimal transmission switching with substation reconfiguration (OTS-SR). This thesis aims at evaluating several formulations and different algorithmic approaches for solving OTS-SR more efficiently regarding the solving time. Several alterations to the formulation of OTS-SR, two Dantzig-Wolfe reformulations and two Benders decompositions are proposed and evaluated using GCG, Coluna, SCIP and Gurobi. Results indicate that not only the introduction of transmission switching makes the OTS problem difficult to solve but also modeling the power flow according to Ohm’s law. Furthermore, it seems to be challenging to make use of the internal structure of the substations in order to solve the OTS-SR problem more efficiently.

When solving a mathematical optimization problem with the GCG solver, in general a model with unknown structure and semantics is given. GCG attempts to detect a structure in order to find a good decomposition for a Dantzig-Wolfe reformulation. If a problem is modeled using the modeling language Pyomo, more information about the model and its structure is available. In this thesis several approaches how this model information can be used to improve the solving process with GCG are analyzed. For example, in Pyomo the different types of constraints and variables and their index sets are directly known, so it is not necessary that GCG tries to detect the constraint and variable classes. Another idea is to represent a generalized version of the model as a graph. This graph can be used to determine the type of the given problem. Depending on the problem type, a predefined decomposition can be applied. All of the tested methods improved the running time when solving some of the test instances. However, there are cases for all methods where they perform poorly and significantly increase the running time.

Branch-and-bound is a widely used algorithm to solve discrete optimization problems, such as integer linear programs. In this thesis, we will explore existing methods to visualize the solving trees that the algorithm employs for program optimization. We will define the balance of a branch-and-bound tree as an index indicating which children tend to make more progress in completing the solving process by making a greater impact on the optimality gap of the bounds. Subsequently, we will compare the balance of trees generated by different solvers and settings, examining their potential impact on the algorithm’s performance.

The thesis aims to address the critical issue of fault detection and diagnosis in large-scale process plants, focusing on minimizing production downtimes and maintaining operational efficiency. My approach involves researching and evaluating machine learning (ML) algorithms such as Mahalanobis Distance, Isolation Forest, and Auto Encoder, etc. Specifically, the motive is to explore innovative techniques such as topological methods and transformer-based anomaly detection. The thesis will focus on two primary areas: A. Research and Identification of Additional ML Models: Conduct a thorough literature review to identify cutting-edge ML methods that could be applied to anomaly detection in large-scale process plants. Focus on innovative approaches such as topological methods and transformer-based anomaly detection. B. Comparative Analysis and Evaluation: Implement the selected ML models and integrate them into an existing framework. Evaluate the performance of the new models against the existing ones using real plant data, assessing their accuracy, robustness, and practical applicability.

The increasing integration of decentralized and renewable energy sources has transformed the low voltage power distribution system into a bidirectional grid. Grid monitoring is required to prevent thermal overload. However, due to high investment costs for measurement devices and the time-consuming installation, distribution system operators (DSOs) often rely on estimation techniques. One possible method of generating pseudo-measurements based on weather data, calendar data and metadata (e.g. the number of electric vehicles in a low-voltage feeder) does not consider historical power usage, making it less effective for accurate power estimation. This thesis explores the application of machine learning techniques in load forecasting to enhance the accuracy of active power estimation in the low voltage grid. By comparing the performance of new models against an existing pseudo-measurement model, the study evaluates the precision of load forecasts at the level of individual feeders and across different forecast horizons. Data for this thesis were collected in Baden-Württemberg, Germany, aiming to provide the DSO with more reliable tools for preventing grid failures.

With the increasing prevalence and application of large language models (LLMs) in problem-solving across a wide range of domains, from routine tasks to complex business applications, this thesis aims to evaluate the capability of LLMs in formulating combinatorial optimization (CO) models, such as those relating to machine scheduling and vehicle routing, using (mixed)-linear integer programming. Our experiments focus on some of the most widely used LLMs, including ChatGPT, Google Gemini, and Microsoft Copilot. For each CO model under consideration, we assess the accuracy and correctness of the model formulations generated by these LLMs. To provide a comprehensive analysis, we employ a range of prompting strategies, from well-defined problem statements to more ambiguous descriptions, to assess the ability of models in processing and interpreting inputs of varying complexity. A larger objective of this thesis is to gauge the reliability of LLMs in supporting operations research practitioners—particularly those with limited experience—in developing optimization models that are used to address complex, real-world problems.

When solving an integer program with a branch-and-price algorithm, cutting planes can be used to strengthen a node’s relaxation. These cutting planes can either be formulated on the original variables or on the master variables. Whenever a cutting plane is added to a node’s reduced master problem, the objectives of the associated pricing problems must be modified to ensure that the column generation still works as intended. If the cutting plane was formulated on the master variables, this requires the addition of new variables and constraints to the pricing problems. Subset-Row inequalities are a class of cutting planes which were originally designed for vehicle routing problems with time windows (VRPTW) and can be viewed as a sub-class of Chvatal-Gomory inequalities. As part of this thesis, a separation plugin which generates subset-row inequalities for a given master problem was added to the branch-price-and-cut solver GCG and its impact on the solving process was evaluated on a set of integer program instances.

Machines can be either in an on or off state. Switching this state frequently is often not possible due to technical limitations, e.g. due to a cool-down or start-up process. In addition, many state changes may cause damage to machines, increase pollution emissions, or cause other side effects. When we use linear programming we can add min-up or min-down constraints to avoid this behavior. These constraints set a minimum uptime and a minimum downtime before the machine is allowed to switch its state again. This can be used, for example, to schedule power plants that have a very energy-intensive start-up process. Therefore, the operator of these power plants may prefer to keep them activated for some consecutive time steps before shutting them down again. This type of model modification adds complexity and often increases the time and the memory required for solving the model. Therefore, we create a constraint handler that aims to allow users to easily add min-up and min-down constraints to their models and to solve them more efficiently. This constraint handler uses several algorithmic techniques that aim to enhance the performance when solving mixed integer linear programming problems with min-up and min-down constraints. We evaluate the efficiency of this constraint handler by comparing several metrics, such as the solving time, when solving with and without the constraint handler. We also implement an algorithm that exploits the structure of the formulation of the min-up and min-down constraints to apply the algorithmic techniques automatically, without additional effort by the user.

Symmetry in mixed-integer programming can lead to redundant computation and increased solving times, especially when problems are decomposed into independent blocks. While symmetry detection is well-studied for general MIPs, less attention has been given to detecting symmetric blocks in automatically generated decompositions. This thesis compares the current approach of GCG—where blocks are compared in pairs and the symmetry detection graph contains only two blocks at a time—to a single detection graph approach, where all blocks are checked together. Results show that pairwise pre-checks are highly effective for quickly filtering out instances where no symmetry detection is needed. For the actual symmetry detection, using a single graph is more efficient, especially when leveraging the graph pre-processor SASSY. The methods are implemented and tested in GCG, and the results indicate that a combination of both approaches achieves the best performance

PyGCGOpt is a Python scripting interface for GCG, a decomposition-based solver that applies a branch-price-and-cut algorithm to solve mixed-integer programs. This thesis extends the PyGCGOpt interface by integrating cutting planes and branching rules into the framework. These enhancements improve the prototyping capabilities of PyGCGOpt and make it more accessible, thereby facilitating research on branch-price-and-cut algorithms. In addition, a computational study is conducted to evaluate the effectiveness of cutting planes within the context of Dantzig-Wolfe reformulations. The results offer insights into the impact of cutting planes on original variables in a branch-price-and-cut setting.

Recently, Machine Learning (ML) techniques have been increasingly leveraged in the context of Column Generation (CG). That is, researchers try to exploit recent advances in ML to accelerate the CG algorithm, e.g., by implementing dual stabilization techniques, pricing heuristics, and column selection procedures using ML algorithms. With respect to partial pricing, the purpose of this thesis is to investigate whether ML can be employed to speed up the pricing phase of the CG algorithm through a pricing problem selection. That is, this thesis investigates for MIPs with many distinct pricing problems (blocks), whether an ML model can be trained that predicts which pricing problems will return valuable columns in the present pricing round and which ones will only generate unworthy columns (if any). In this context, different metrics (scores) are defined that aim at quantifying the benefit that solving a pricing problem in a pricing round yielded. Based on data collected from a broad test set of different MIP models, the suitability of the defined scores to serve as targets is examined by analyzing for different MIP instances whether notable variations between their blocks can be seen according to the scores. Ultimately, (ML-) models are trained with the help of the defined scores and tested for various MIP instances in a partial pricing setup. Thus, it is evaluated whether such an (ML-) model supported partial pricing approach based on the defined scores yields notable performance improvements or not.

For an economically and ecologically sustainable operation of energy systems, both the plant design and its optimal operation are crucial. In recent years, the increasing integration of renewable energies and the combination of various energy sources and storage technologies have led to significantly more complex multi-energy systems (MES). When using mixed-integer linear programming to optimize such systems, standard solvers often reach their computational limits. In this work, we investigate decomposition methods, in particular, the Dantzig-Wolfe reformulation, as a way to tackle these challenges in the operational optimization of industrial MES. To this end, we identify common modeling structures that arise in MES optimization and construct various example systems exhibiting these structures. We develop several decomposition approaches, and evaluate in numerical studies whether the generic branch-price-and-cut solver GCG offers advantages over conventional branch-and-cut solvers, especially for complex instances with a large number of components or extended time horizons. To ensure broad applicability beyond specific case studies, we also propose a generic mathematical framework capable of expressing and applying the developed decomposition approaches.

Bei vielen Logistiknetzwerken kommt es im Laufe der Zeit immer wieder zu Veränderungen der Ausgangsbedingungen, wie zum Beispiel veränderten Kosten oder Transportbedarfen, sodass eine Anpassung an die neuen Gegebenheiten notwendig wird. Bereits kleine Veränderungen eines vorhandenen Netzwerkes können großen Aufwand und Kosten in der Praxis bedeuten. Darum ist eine Einschränkung für die Anpassung eines solchen Systems von Interesse, um diese auch anwendbar und praktikabel zu gestalten. In dieser Bachelorarbeit sollen verschiedene Methoden untersucht werden, ein schon vorhandenes System möglichst gut an die veränderten Bedingungen anzupassen. Hierbei sollen nur beschränkt viele Änderungen zugelassen werden und trotzdem ein neues Netzwerk entstehen, welches im Vergleich zu einer optimalen Anpassung möglichst gut ist. Dazu werden im Folgenden verschiedene Metriken für die Distanz zwischen zwei Lösungen entwickelt und verglichen, indem gemischt-ganzzahlige Programme ausgearbeitet, implementiert und für verschiedene Test-Daten gelöst werden. Ferner sollen die erhaltenen Lösungen hinsichtlich des Abstands zueinander, der Lösungsgüte, der Rechenzeit und ihrer Veränderung am Netzwerk analysiert werden.