Routing and scheduling optimisation under uncertainty for engineering applications

The thesis aims to develop a viable computational approach suitable for solving large vehicle routing and scheduling optimisation problems affected by uncertainty. The modelling framework is built upon recent advances in Stochastic Optimisation, Robust Optimisation and Distributionally Robust Optimization. The utility of the methodology is presented on two classes of discrete optimisation problems: scheduling satellite communication, which is a variant of Machine Scheduling, and the Vehicle Routing Problem with Time Windows and Synchronised Visits. For each problem class, a practical engineering application is formulated using data coming from the real world. The significant size of the problem instances reinforced the need to apply a different computational approach for each problem class. Satellite communication is scheduled using a Mixed-Integer Programming solver. In contrast, the vehicle routing problem with synchronised visits is solved using a hybrid method that combines Iterated Local Search, Constraint Programming and the Guided Local Search metaheuristic. The featured application of scheduling satellite communication is the Satellite Quantum Key Distribution for a system that consists of one spacecraft placed in the Lower Earth Orbit and a network of optical ground stations located in the United Kingdom. The satellite generates cryptographic keys and transmits them to individual ground stations. Each ground station should receive the number of keys in proportion to the importance of the ground station in the network. As clouds containing water attenuate the signal, reliable scheduling needs to account for cloud cover predictions, which are naturally affected by uncertainty. A new uncertainty sets tailored for modelling uncertainty in predictions of atmospheric phenomena is the main contribution to the methodology. The uncertainty set models the evolution of uncertain parameters using a Multivariate Vector Auto-Regressive Time Series, which preserves correlations over time and space. The problem formulation employing the new uncertainty set compares favourably to a suite of alternative models adapted from the literature considering both the computational time and the cost-effectiveness of the schedule evaluated in the cloud cover conditions observed in the real world. The other contribution of the thesis in the satellite scheduling domain is the formulation of the Satellite Quantum Key Distribution problem. The proof of computational complexity and thorough performance analysis of an example Satellite Quantum Key Distribution system accompany the formulation. The Home Care Scheduling and Routing Problem, which instances are solved for the largest provider of such services in Scotland, is the application of the Vehicle Routing Problem with Time Windows and Synchronised Visits. The problem instances contain over 500 visits. Around 20% of them require two carers simultaneously. Such problem instances are well beyond the scalability limitations of the exact method and considerably larger than instances of similar problems considered in the literature. The optimisation approach proposed in the thesis found effective solutions in attractive computational time (i.e., less than 30 minutes) and the solutions reduced the total travel time threefold compared to alternative schedules computed by human planners. The Essential Riskiness Index Optimisation was incorporated into the Constraint Programming model to address uncertainty in visits' duration. Besides solving large problem instances from the real world, the solution method reproduced the majority of the best results reported in the literature and strictly improved the solutions for several instances of a well-known benchmark for the Vehicle Routing Problem with Time Windows and Synchronised Visits.The thesis aims to develop a viable computational approach suitable for solving large vehicle routing and scheduling optimisation problems affected by uncertainty. The modelling framework is built upon recent advances in Stochastic Optimisation, Robust Optimisation and Distributionally Robust Optimization. The utility of the methodology is presented on two classes of discrete optimisation problems: scheduling satellite communication, which is a variant of Machine Scheduling, and the Vehicle Routing Problem with Time Windows and Synchronised Visits. For each problem class, a practical engineering application is formulated using data coming from the real world. The significant size of the problem instances reinforced the need to apply a different computational approach for each problem class. Satellite communication is scheduled using a Mixed-Integer Programming solver. In contrast, the vehicle routing problem with synchronised visits is solved using a hybrid method that combines Iterated Local Search, Constraint Programming and the Guided Local Search metaheuristic. The featured application of scheduling satellite communication is the Satellite Quantum Key Distribution for a system that consists of one spacecraft placed in the Lower Earth Orbit and a network of optical ground stations located in the United Kingdom. The satellite generates cryptographic keys and transmits them to individual ground stations. Each ground station should receive the number of keys in proportion to the importance of the ground station in the network. As clouds containing water attenuate the signal, reliable scheduling needs to account for cloud cover predictions, which are naturally affected by uncertainty. A new uncertainty sets tailored for modelling uncertainty in predictions of atmospheric phenomena is the main contribution to the methodology. The uncertainty set models the evolution of uncertain parameters using a Multivariate Vector Auto-Regressive Time Series, which preserves correlations over time and space. The problem formulation employing the new uncertainty set compares favourably to a suite of alternative models adapted from the literature considering both the computational time and the cost-effectiveness of the schedule evaluated in the cloud cover conditions observed in the real world. The other contribution of the thesis in the satellite scheduling domain is the formulation of the Satellite Quantum Key Distribution problem. The proof of computational complexity and thorough performance analysis of an example Satellite Quantum Key Distribution system accompany the formulation. The Home Care Scheduling and Routing Problem, which instances are solved for the largest provider of such services in Scotland, is the application of the Vehicle Routing Problem with Time Windows and Synchronised Visits. The problem instances contain over 500 visits. Around 20% of them require two carers simultaneously. Such problem instances are well beyond the scalability limitations of the exact method and considerably larger than instances of similar problems considered in the literature. The optimisation approach proposed in the thesis found effective solutions in attractive computational time (i.e., less than 30 minutes) and the solutions reduced the total travel time threefold compared to alternative schedules computed by human planners. The Essential Riskiness Index Optimisation was incorporated into the Constraint Programming model to address uncertainty in visits' duration. Besides solving large problem instances from the real world, the solution method reproduced the majority of the best results reported in the literature and strictly improved the solutions for several instances of a well-known benchmark for the Vehicle Routing Problem with Time Windows and Synchronised Visits

Optimised configuration of sensing elements for control and fault tolerance applied to an electro-magnetic suspension system

New technological advances and the requirements to increasingly abide by new safety laws in engineering design projects highly affects industrial products in areas such as automotive, aerospace and railway industries. The necessity arises to design reduced-cost hi-tech products with minimal complexity, optimal performance, effective parameter robustness properties, and high reliability with fault tolerance. In this context the control system design plays an important role and the impact is crucial relative to the level of cost efficiency of a product. Measurement of required information for the operation of the design control system in any product is a vital issue, and in such cases a number of sensors can be available to select from in order to achieve the desired system properties. However, for a complex engineering system a manual procedure to select the best sensor set subject to the desired system properties can be very complicated, time consuming or even impossible to achieve. This is more evident in the case of large number of sensors and the requirement to comply with optimum performance. The thesis describes a comprehensive study of sensor selection for control and fault tolerance with the particular application of an ElectroMagnetic Levitation system (being an unstable, nonlinear, safety-critical system with non-trivial control performance requirements). The particular aim of the presented work is to identify effective sensor selection frameworks subject to given system properties for controlling (with a level of fault tolerance) the MagLev suspension system. A particular objective of the work is to identify the minimum possible sensors that can be used to cover multiple sensor faults, while maintaining optimum performance with the remaining sensors. The tools employed combine modern control strategies and multiobjective constraint optimisation (for tuning purposes) methods. An important part of the work is the design and construction of a 25kg MagLev suspension to be used for experimental verification of the proposed sensor selection frameworks

A survey on low-thrust trajectory optimization approaches

In this paper, we provide a survey on available numerical approaches for solving low-thrust trajectory optimization problems. First, a general mathematical framework based on hybrid optimal control will be presented. This formulation and their elements, namely objective function, continuous and discrete state and controls, and discrete and continuous dynamics, will serve as a basis for discussion throughout the whole manuscript. Thereafter, solution approaches for classical continuous optimal control problems will be briefly introduced and their application to low-thrust trajectory optimization will be discussed. A special emphasis will be placed on the extension of the classical techniques to solve hybrid optimal control problems. Finally, an extensive review of traditional and state-of-the art methodologies and tools will be presented. They will be categorized regarding their solution approach, the objective function, the state variables, the dynamical model, and their application to planetocentric or interplanetary transfers

Analyzing SPSA approaches to solve the non-linear non-differentiable problems arising in the assisted calibration of traffic simulation models

Mathematical and simulation models of systems lay at the core of decision support systems, and their role become more critical as more complex is the system object of decision. The decision process usually encompasses the optimization of some utility function that evaluates the performance indicators that measure the impacts of the decisions. An increasing difficulty directly related to the complexity of the system arises when the associated function to be optimized is a not analytical, non-differentiable, non-linear function which can only be evaluated by simulation. Simulation-Optimization techniques are especially suited in these cases and its use is increasing in traffic models, an archetypic case of complex, dynamic systems exhibiting highly stochastic characteristics. In this approach simulation is used to evaluate the objective function, and a non-differentiable optimization technique to solve the optimization problem is used. Simultaneous Perturbation Stochastic Approximation (SPSA) is one of the most popular of these techniques. This thesis analyzes, discusses and presents computational results for the application of this technique to the calibration of a traffic simulation model of a Swedish highway section. Variants of the SPSA, replacing the usual gradient approach by a combination of normalized parameters and penalized objective function, have been proposed in this study due to an exhaustive analysis of the behavior of classical SPSA where problems arose from different magnitude variables. In this work, a varied set of Software environments have been used, combining RStudio for the analysis, Python and MATLAB for the SPSA implementation, AIMSUN as a Traffic Model Simulator, and SQLite for obtaining of simulated data and Tableau for visualizing data and results