Very large-scale neighborhood search for drone routing with energy replenishment

This thesis deals with models and methods for large scale optimization problems; in particular, we focus on decision problems arising in the context of seaport container terminals for the efficient management of terminal operations. Large-scale optimization problems are both difficult to handle and important in many concrete contexts. They usually originate from real world applications, such as telecommunication, transportation and logistics, and their combinatorial complexity often represents a major issue; therefore, optimization models are crucial to support the decision making process. In particular, column generation and branch-and-price schemes currently represent one of the most advanced and efficient exact optimization approaches to solve large scale combinatorial problems. However, the increasing size and complexity of practical problems arising in real-world applications motivates the design of new solution approaches able to tackle current optimization challenges. In this thesis, we address two complementary research streams where both methods and applications play an important role. On the one hand, we focus on the specific application of container terminals: we propose a new model for the integrated planning of operations and we provide a heuristic and an exact solution algorithm; the broader objective is to devise solution methods that can be generalized and extended to other applications and domains. On the other hand, we aim to develop new methods and algorithms for general large scale problems and, in this context, we investigate a new column generation framework that exploits the relationship between compact and extensive formulation. In particular, we focus on a class of split delivery vehicle routing problems that generalizes a large number of applications arising in the real world, such as transportation and logistics, including container terminal management. In the context of container terminals, we propose a model for the integrated planning of berth allocation and quay crane assignment: the two decision problems are usually solved hierarchically by terminal planners, whereas in the Tactical Berth Allocation Problem we optimize the two problems simultaneously. We firstly present a mixed integer programming formulation that is embedded into a two-level heuristic algorithm based on tabu search and mathematical programming techniques: our heuristic proves to be very efficient, providing good-quality solutions in a reasonable time. The problem is reformulated via Dantzig-Wolfe decomposition and solved via column generation: we propose an exact branch-and-price algorithm and our implementation, that includes state-of-the-art techniques for the master and the pricing problem, outperforms commercial solvers. Furthermore, the exact approach allows us to provide an interesting experimental comparison between hierarchical and integrated planning: computational tests confirm the added value of integration in terms of cost reduction and efficient use of resources. From a methodological point of view, this dissertation investigates a new column generation concept for difficult large scale optimization problems. In particular, we study a class of split delivery vehicle routing problems that generalizes some interesting features of Tactical Berth Allocation Problem, which are relevant also to other applications such as transportation, logistics and telecommunication. The problem, called Discrete Split Delivery Vehicle Routing Problem with Time Windows, presents two main modeling features: demand is discrete and delivered in discrete orders, opposite to the usual assumption of continuously splittable demand; the service time is dependent on the delivered quantity, opposite to the usual assumption of constant service time, regardless of the quantity. The problem is used to validate and test the new column generation approach studied in this thesis. The proposed framework, called Two-stage column generation, represents a novel contribution to recent advances in column generation: the basic idea is to simultaneously generate columns both for the compact and the extensive formulation. We propose to start solving the problem on a subset of compact formulation variables, we apply Dantzig-Wolfe decomposition and we solve the resulting master problem via column generation. At this point, profitable compact formulation variables are dynamically generated and added to the formulation according to reduced cost arguments, in the same spirit of standard column generation. The key point of our approach is that we evaluate the contribution of compact formulation variables with respect to the extensive formulation: indeed, we aim at adding compact formulation variables that are profitable for the master problem, regardless of the optimal solution of the linear relaxation of the compact formulation. We apply two-stage column generation to the Discrete Split Delivery Vehicle Routing Problem with Time Windows. Computational results show that our approach significantly reduces the number of generated columns to prove optimality of the root node. Furthermore, suboptimal compact formulation variables are detected correctly and a large number of variables is not taken into account during the solution process, thus reducing the size of the problem. However, the additional effort required by such a sophisticated approach makes the method competitive in terms of computational time only for instances of a certain difficulty. To conclude, two-stage column generation is a promising new approach and we believe that further research in this direction may contribute to solve more and more complex large scale optimization problems.

This paper discusses neighborhood search algorithms where the size of the neighborhood is "very large" with respect to the size of the input data. We concentrate on such a very large scale neighborhood (VLSN) search technique based on compounding independent moves (CIM) such as 2-opts, swaps, and insertions. We present a systematic way of creating and searching CIM neighborhoods for routing problems with side constraints. For such problems, the exact search of the CIM neighborhood becomes NP-hard. We introduce a multi-label shortest path algorithm for searching these neighborhoods heuristically. Results of a computational study on the vehicle routing problem with capacity and distance restrictions shows that CIM algorithms are very competitive approaches for solving vehicle routing problems. Overall, the solutions generated by the CIM algorithm have the best performance among the current solution methodologies in terms of percentage deviation from the best-known solutions for large-scale capacitated VRP instances.

Unmanned aerial vehicles (UAVs), or drones, are recognized for their potential to improve efficiency in last-mile delivery. Unlike the vehicle routing problem, drone route design is challenging due to several operational signatures, such as speed optimization, multi-trip operation, and energy consumption estimation. Drone energy consumption is a nonlinear function of both speed and payload. Moreover, the high speed of drones can significantly curtail the drone range, thereby limiting the efficiency of drone delivery systems. This paper addresses the trade-off between speed and flight range in a multi-trip drone routing problem with variable flight speeds (DRP–VFS). We propose a new model to specifically consider energy constraints using a nonlinear energy consumption model and treat drone speeds as decision variables. The DRP–VFS is initially formulated using mixed-integer linear programming (MILP) to minimize energy consumption. To solve large-scale instances, we propose a three-phase adaptive large neighborhood search (ALNS) algorithm and compare its performance with a commercial MIP solver. The experimental results demonstrate that the proposed method is capable of effectively identifying suboptimal solutions in practical scenarios. Furthermore, results indicate that operating drones at variable speeds leads to about 21% energy savings compared to fixed speeds, with advantages in cost savings and range extension.

Drone routing with energy function: Formulation and exact algorithm

Partial-route inequalities for the multi-vehicle routing problem with stochastic demands

The electric vehicle routing problem (EVRP) is constrained by the limited driving range and time-consuming energy replenishment. These characteristics shift the focus of the EVRP from simple path optimization to an integrated optimization of routing and energy replenishment. Consequently, the energy replenishment strategy becomes a critical determinant of the feasibility and economic viability of EVRP solutions. This paper presents a systematic literature review structured around a core classification of replenishment strategies. The strategies are categorized into two primary modes: charging and battery swapping. This framework addresses common gaps in existing research, such as imprecise strategy definitions and fragmented analyses. For the charging strategy, we establish a three-dimensional classification framework, which comprises the charging function, charging policy, and charging station type. Within this context, wireless charging is considered as a special method of energy replenishment. The battery swapping strategy relies on battery swapping stations (BSSs): the EVRP with BSSs (EVRP–BSSs) and the BSS location-routing problem with electric vehicle (BSS–EV–LRP). Our review identifies several limitations in the current body of research. These include an imbalance between modeling accuracy and computational efficiency, insufficient coverage of diverse operational scenarios, and a superficial integration of emerging technologies. Furthermore, many studies lack a multi-stakeholder perspective that considers collaborative solutions. Future research should prioritize addressing these gaps. Key directions include developing effective methods for solving nonlinear charging functions and expanding research into more specialized scenarios. Additionally, there is a need to improve collaborative algorithms for battery swapping and to develop shared BSS models that serve multiple enterprises. The strategy-driven framework proposed here offers a clear reference for modeling and scenario adaptation in future EVRP studies.

In this article, we propose a branch‐and‐price‐and‐cut (BPC) algorithm to exactly solve the manpower routing problem with synchronization constraints (MRPSC). Compared with the classical vehicle routing problems (VRPs), the defining characteristic of the MRPSC is that multiple workers are required to work together and start at the same time to carry out a job, that is, the routes of the scheduling subjects are dependent. The incorporation of the synchronization constraints increases the difficulty of the MRPSC significantly and makes the existing VRP exact algorithm inapplicable. Although there are many types of valid inequalities for the VRP or its variants, so far we can only adapt the infeasible path elimination inequality and the weak clique inequality to handle the synchronization constraints in our BPC algorithm. The experimental results at the root node of the branch‐and‐bound tree show that the employed inequalities can effectively improve the lower bound of the problem. Compared with ILOG CPLEX, our BPC algorithm managed to find optimal solutions for more test instances within 1 hour. © 2016 Wiley Periodicals, Inc. Naval Research Logistics 63: 138–171, 2016

The aquaculture service vessel routing problem with time dependent travel times and synchronization constraints

On a Pickup to Delivery Drone Routing Problem: Models and algorithms

This paper considers the vehicle routing problem with soft time windows, a challenging routing problem where customers' time windows may be violated at a certain cost. The vehicle routing problem with soft time windows has a lexicographic objective function, aimed at minimizing first the number of routes, then the number of violated time windows, and finally the total routing distance. We present a multistage very large-scale neighborhood search for this problem. Each stage corresponds to a variable neighborhood descent over a parameterizable very large-scale neighborhood. These neighborhoods contain an exponential number of neighbors, as opposed to classical local search neighborhoods. Often, searching very large-scale neighborhoods can produce local optima of a higher quality than polynomial-sized neighborhoods can. Furthermore, we use a sophisticated heuristic to determine service start times allowing us to minimize the number of violated time windows. We test our approach on a number of different problem types, and compare the results to the relevant state of the art. The experimental results show that our algorithm improves best known solutions on 53% of the most studied instances. Many of these improvements stem from a reduction of the number of vehicles, a critical objective in vehicle routing problems.

An optimization algorithm for the workforce management in a retail chain

Optimization based models and algorithms for routing and scheduling with time window constraints: Koskosidis, Ioannis Alkiviadis, Ph.D. Princeton University, 1988, 282pp.

This paper proposes a novel admission and routing scheme which takes into account arbitrarily assigned priorities for network flows. The presented approach leverages the centralized Software Defined Networking (SDN) capabilities in order to do so. Exact and heuristic approaches to the stated Priority Flow Admission and Routing (PFAR) problem are provided. The exact approach which provides an optimal solution is based on Integer Linear Programming (ILP). Given the potentially long running time required to find an exact and optimal solution, a heuristic approach is proposed; this approach is based on Genetic Algorithms (GAs). In order to effectively estimate the performance of the proposed approaches, a simulator that is capable of generating semi-random network topologies and flows has been developed. Experimental results for large problem instances (up 50 network nodes and thousands of network flows), show that: i) an optimal solution can be often found in few seconds (even milliseconds), and ii) the heuristic approach yields close-to-optimal solutions (approximately 95% of the optimal) in a fixed amount of time; these experimental results demonstrate the pertinence of the proposed approaches.

A Contribution to the Reentrant Flow-Shop Scheduling Problem

Application of Genetic Algorithm methodologies in fuel bundle burnup optimization of Pressurized Heavy Water Reactor