ScienceDocs Optimization Engineering Team and Partner Case Studies

The ScienceDocs Optimization Engineering team brings together specialists in operations research, artificial intelligence, and advanced data analytics. These case studies demonstrate how optimization modeling, machine learning, and AI-driven decision tools can significantly improve operational performance across industries.

Optimization Engineer/Internal Case Studies

Case Study #1 - A data-driven AI model to predict potential late deliveries

warehouse optimization

Figure 1: Pipeline of AI model development.

One of the key operational challenges that many companies face is minimizing late deliveries to customers. Reliable and timely delivery is essential in many industries, as delays can disrupt customer operations and negatively affect long-term relationships. The goal of this project is to develop a data-driven approach to understand the root causes of delivery delays and identify opportunities to improve delivery performance.


The project begins with a comprehensive analysis of the company’s historical delivery data (Figure 1). This includes exploratory data analysis to identify patterns and relationships among variables, such as shipment characteristics, order attributes, customer locations, transportation conditions, and delivery performance. Visualization techniques and statistical summaries are used to uncover key drivers of late deliveries and highlight operational bottlenecks within the supply chain (Figure 2).


Building on these insights, predictive models are developed to estimate the likelihood of delivery delays for future shipments. Several classification models, including classification tree, Random Forest, Logistic Regression, and GridSearch are trained and evaluated to identify the most effective approach for predicting late deliveries. The results enable the company to proactively identify high-risk shipments and implement preventive actions. In addition, the analysis estimates the financial impact of delivery delays and the economic value of improved prediction accuracy. Overall, the project provides companies with actionable insights to improve delivery reliability and prioritize operational improvements. Click here to observe the correlation in exploratory data analysis (EDA).

This Project was performed by Dr. Mohammad



Case Study #2 - Warehouse Optimization

storage assignment algorithms Warehouse order picking is one of the most time-consuming processes in logistics operations, often accounting for more than half of total warehouse labor time. This case study focused on improving storage assignment strategies in order to significantly reduce picking time and improve operational efficiency.

The project began with extensive data analysis, including pattern recognition of orders, inventory levels, and part relationships. Pareto analysis and correlation analysis were used to identify meaningful relationships among parts, shipments, and customer orders. These insights allowed the development of clustering models to identify families of parts that are frequently ordered together.

Multiple analytical methods were applied, including Principal Component Analysis (PCA), Singular Value Decomposition (SVD), and integer programming models to identify optimal part clusters. Storage assignment algorithms were then developed to position these clusters strategically throughout the warehouse in order to minimize total travel distance while also reducing aisle congestion between workers.

The optimized storage design significantly improved operational efficiency, reducing order picking time by more than 40% compared to the existing layout.

This Project was performed by Dr. Mohammad


Case Study #3 - Optimizing Warehouse Sorter Reconfiguration

warehouse optimizationOverview

A large-scale warehouse fulfillment operation needed to improve how packages were sorted and consolidated for outbound shipping. Packages were grouped by destination and carrier (e.g., UPS, USPS, FedEx) and routed into containers—either large boxes or pallets—before being loaded onto scheduled trucks throughout the day.

The core objective was to maximize container utilization to reduce the need for additional trucks, while avoiding sorter congestion caused by insufficient chute allocation for high-demand destinations.

The Challenge

The warehouse sorter system consisted of a conveyor loop with multiple chutes, each feeding into a designated container. The system operated under several constraints:

  • Packages remained on the conveyor until they reached their assigned chute
  • If a container was full, packages could miss their chute while a manual swap (taking several minutes) was performed
  • The conveyor loop was divided into two “sides,” limiting routing flexibility
  • Some facilities included a linear extension sorter for low-demand destinations, accessible only via a connecting chute
  • Demand fluctuated continuously based on outbound truck schedules

These dynamics created a complex, real-time balancing problem:

  • Too few chutes assigned to a destination → congestion and missed drops
  • Too many chutes → underutilized containers and unnecessary trucks

Traditional rule-of-thumb approaches varied by site and were not sufficient for consistent, system-wide optimization.

The Solution

An integer programming optimization model was developed to dynamically reconfigure chute assignments in response to changing demand.

Key elements of the solution included:

  • Incremental optimization runs every 30 minutes, using updated package forecasts
  • Input from the warehouse control system, including:
    • Current sorter configuration
    • Forecasted package volumes by destination
  • Deployment as a Java library integrated into the main control system
  • Use of the Xpress optimization engine to solve the model

To meet strict operational constraints, the model was designed for high-speed execution:

  • Required runtime: ≤ 60 seconds
  • Achieved runtime: < 10 seconds

This performance improvement was enabled by recognizing that:

  • Chutes on each side of the loop could be treated as generic resource pools
  • Once location-level allocation was determined, assigning specific chutes became a much simpler downstream problem

Implementation & Outputs

The optimization engine returned multiple solution views to support both operators and developers, including:

  • Recommended chute-to-destination assignments
  • Configuration changes from the current state
  • Diagnostic views for evaluating system performance

This allowed for both operational decision-making and continuous refinement of the model.

Results & Impact

The deployment of this optimization framework produced measurable operational improvements:

  • Reduced reliance on additional outbound trucks through better container utilization
  • Improved sorter throughput by minimizing congestion and missed chute events
  • Standardized decision-making across hundreds of warehouse facilities
  • Replaced inconsistent, local heuristics with a scalable, data-driven optimization approach

Takeaway

This project demonstrates how real-time optimization—when tightly integrated into operational systems—can transform complex logistics environments. By combining fast-solving integer programming with system-aware modeling, the solution delivered both efficiency gains and enterprise-wide consistency at scale.

Key Ideas (Quick Takeaways)

  • Real-time optimization (≤10 sec solve time) is critical for live warehouse systems
  • Treating physical constraints as abstract resources can dramatically speed up models
  • Frequent re-optimization (every 30 min) enables responsiveness to demand shifts
  • Standardization across sites is a major hidden value driver
  • Optimization reduces both cost (fewer trucks) and operational friction (less congestion)

This Project was performed by Dr. Michael

Case Study #4 - Supply Chain Network Optimization Under Uncertainty

Overview

Supply chain network design is the strategic design of the physical and decision structure through which products move from supply sources to end customers. In practice, it addresses questions such as which suppliers, plants, warehouses, distribution centers, or hubs should be used, how the network footprint should be configured, which customer markets should be served from which facilities, and how material and product flows should move across the network to meet service requirements at the lowest practical total landed cost.  See the details of this Supply Chain Network Optimization Case Study.

Case Study #5: Online Retail Inventory Placement Optimization

Key Ideas

  • Online retailers need to find cost-effective inventory placement plans that balance variable costs and fixed costs, a challenging problem due to its large scale and computational complexity.
  • We design a large-scale optimization framework that solves the problem with a three-step process: aggregating items into clusters, solving the cluster-level problem with column generation, and disaggregating cluster-level solutions into item-level placement plans with a performance guarantee.
  • The framework delivers value for businesses through scale (handling millions of items), speed (solving problems 18.77 times faster than the previous industry standard), and savings (reducing the placement cost by 1.7%).

Overview

In online retail, customer orders can be shipped from any fulfillment center, near or far, at different costs. In contrast, traditional brick-and-mortar retail customers are served primarily by the closest physical store. The flexibility of online retail shipping makes it possible to achieve lower costs and higher service levels. However, it also requires all the fulfillment centers to be considered together as a whole rather than in isolation, which increases computational complexity and calls for new inventory optimization techniques.

This case study focuses on the strategic decision of inventory placement, the problem of selecting fulfillment centers in which to place items. We define a placement plan as a chosen subset of fulfillment centers for placing each item, along with an estimate of the inventory levels.

The goal is to determine a cost-minimizing placement plan for each item, while satisfying customer demand requirements and fulfillment center capacity constraints. The objective is to minimize the sum of two types of costs: variable costs and fixed costs.

  • Variable costs are proportional to the volume of inventory handled. We focus on the outbound shipping cost, but there may also be other variable costs, such as inbound shipping and receiving costs, or the cost of labor for picking and packing items before outbound shipment.
  • Fixed costs are charged when an item is placed at a fulfillment center, regardless of its inventory level. This can represent actual monetary costs, such as a flat-rate inbound delivery cost charged by upstream suppliers for every separate delivery location. Fixed costs also incentivize items to be consolidated in fewer fulfillment centers to reduce logistical overhead, a concept known as “sparsity”. This concept is particularly important for fulfillment centers that have limited containers for carrying unique items, as well as for low-demand items that have few units and should not be split across too many locations.

A good placement plan balances the trade-off between variable costs and fixed costs. On the one hand, placing items in more fulfillment centers can achieve a lower variable cost, as the items can be placed closer to the customers. On the other hand, consolidating items in fewer fulfillment centers incurs a lower fixed cost. Finding a placement plan that achieves the “happy medium” can significantly reduce operational costs, which is crucial for online retail businesses because they often have a low profit margin.

The Challenge

The difficulty of the inventory placement problem arises from its scale and computational complexity in various aspects:

  • Large number of items: An online retailer may offer thousands or millions of items in their catalog.
  • Large fulfillment network: There may be dozens or hundreds of fulfillment centers, serving hundreds of customer demand regions.
  • Fixed costs: Fulfillment fixed costs introduce discrete (binary) variables into the optimization formulation, increasing the computational complexity of the problem.
  • Capacity constraints: Fulfillment center capacity constraints create coupling between items because they share the same space. This requires all items to be optimized jointly rather than independently.

To demonstrate the impracticality of a naive solution approach with off-the-shelf software, we solve a small sample problem with a network of 88 fulfillment centers, 99 customer demand regions, and 1,000 items. Using Gurobi, an optimality gap of 6-7% is reached after 5 hours on a personal laptop, and the gap does not improve much after an additional 5 hours. For problems with more items, even constructing the formulation becomes prohibitive in terms of computer memory. Thus, a more efficient approach is needed.

The Solution

We develop a large-scale optimization framework that consists of three main steps (Figure 1):

  1. Item aggregation. First, the items are aggregated into clusters based on similarity in demand patterns. In our numerical experiments, we use k-means clustering based on the Euclidean distance between weighted vectors of item attributes. In practice, similarity measures and cluster procedures should be chosen according to the business context. For instance, existing item categorization may naturally lend itself to some partition. Also, if historical demand data identifies certain items that are likely to be purchased together, it may be a good idea to aggregate them in the same cluster, so that they are predisposed to being shipped from the same fulfillment center in one package.
  2. Column generation. Next, the cluster-level optimization problem is solved efficiently using the technique of column generation, a specialized, highly efficient technique for quickly solving large problems. We develop a column-based reformulation wherein a column for a cluster is a placement plan in which each customer region is fully served by a single fulfillment center. Intuitively, the columns can be thought of as a set of “bases” that span the solution space of the original cluster-level problem. By only generating and retaining low-cost columns as we go, the reformulated problem can be solved quickly.
  3. Cluster disaggregation. Finally, cluster-level solutions are disaggregated into item-level placement plans. Since each cluster contains many items, there is flexibility to perform disaggregation in a way that does not deviate too much from the cluster-level plan.

Mathematically, these three steps can be viewed as a series of approximations that break down the scale and computational complexity of the original item-level placement problem. To ensure that the final approximate solution has a near-optimal cost, we derive mathematical bounds on the optimality gap that can be easily calculated and provide guarantees on the solution quality.

Users have the ability to choose the number of clusters depending on their computational resources. Having more clusters means a more fine-grained approximation and therefore a smaller optimality gap, but it would also require more time and computer memory to solve.

Figure 1

large scale optimization framework

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Implementation and Outputs

We demonstrate the effectiveness of this solution framework with a sample problem that consists of 1,000,000 items, to be placed in a network of 88 warehouses serving 98 customer regions in the U.S. (Figure 2). Our program is implemented in Python and Gurobi.

Figure 2

sample placement problem

For this sample problem, our large-scale computation framework produces a placement plan in 3.54 hours, with a 0.74% optimality gap (that is, the cost is no more than 0.74% higher than the optimal cost).

We compare our results with a sequential placement heuristic similar to the status quo algorithm used by our industry partner. The heuristic is to place one item at a time, according to some user-specified order. After an item is placed, the remaining fulfillment center capacities are given as input to the placement problem of the next item. The main advantage of the sequential placement heuristic is its simplicity. By placing one item at a time, it avoids the complexity of jointly optimizing a large number of items. The main drawback, however, is that sequentially solving 1,000,000 instances of a simple problem can still take a long time, and that the myopic approach can result in an arbitrarily suboptimal cost.

For the same sample problem, the sequential placement heuristic produces a placement plan with a 2.45% optimality gap in 66.45 hours.

Compared to the status quo heuristic, our solution is 1.7% lower in cost while being 18.77 times faster. (Figure 3)

Figure 3

performance comparison

Results and Business Impact

Our large-scale optimization framework helps online retail companies stay competitive by delivering inventory placement plans with scale, speed, and savings.

  • Scale: The solution framework is able to handle large fulfillment networks and millions of items. Even at this scale, it still accurately models a combination of variable and fixed costs, as well as capacity constraints, allowing it to capture different types of realistic business considerations.
  • Speed: The solution framework solves the problem in a few hours rather than a few days. This means that decision-makers can re-run the optimization frequently as updated inputs are received, allowing them to respond to changing conditions quickly.
  • Savings: The solution framework outperforms status quo heuristics in producing a low-cost placement plan, resulting in cost savings. Moreover, it also produces a performance guarantee (optimality gap) that gives businesses the peace of mind that the placement plan has a near-optimal cost.

 

Takeaway

Cost-minimizing inventory placement is crucial for reducing the operational cost of online retail businesses. The problem is challenging to solve with off-the-shelf software due to its scale and computational complexity. Our large-scale optimization framework produces low-cost placement plans quickly by utilizing the three-step technique of item aggregation, column generation, and cluster disaggregation. For a large, realistic sample problem, our framework produces a placement plan in 3.54 hours that is within 0.74% of the optimal cost, outperforming the status quo heuristic by being 1.7% lower in cost and 18.77 times faster. Our framework delivers scale, speed and savings for online retail businesses.

This case study was performed by Dr. Chen

Case Study #6 - Revealing Near-Optimal Strategic Alternatives with Modeling to Generate Alternatives (MGA)

Overview
Optimization models are often used to identify a single best solution under a fixed set of assumptions. In many planning problems, however, decision-makers need more than one mathematically optimal answer. What is often needed is a structured way to identify other materially different strategies that remain close to optimal, especially when important real-world considerations are not fully represented in the model.

This case study centers on the use of Modeling to Generate Alternatives (MGA) to expose that hidden decision space. MGA was embedded within a large-scale optimization framework to identify solutions with similar objective values but meaningfully different structures. The objective was not merely to compute a least-cost plan, but to determine how much flexibility existed near the optimum and which strategic tradeoffs remained available once a hard system target was imposed.

The Challenge
Traditional scenario analysis can vary inputs and assumptions, but it does not fully address structural uncertainty. A model may produce one optimal configuration even when several other near-optimal configurations are also plausible and potentially more attractive once external factors are considered. Those factors may include implementation risk, stakeholder preference, technology acceptance, policy constraints, or operational practicality.

MGA was used specifically because conventional scenario analysis and Monte Carlo methods do not fully account for this type of structural uncertainty. It was framed as a way to search the near-optimal solution space rather than stopping at the single least-cost answer. That framing made it possible to answer questions such as:

  • Are there multiple fundamentally different ways to satisfy the same system target?
  • How narrow or broad is the near-optimal decision space?
  • Which choices are essential across many alternatives, and which are substitutable?
  • When does the model indicate flexibility, and when does it indicate a constrained path?

The Method
The baseline optimization first identified a least-cost solution under a stringent system-wide target. MGA was then applied by adding a cost slack and re-solving the model in a way that intentionally searched for solutions that were as different as possible in decision space while still remaining near-optimal in cost. The original objective was converted into a constraint with an allowed deviation above the optimum, and the search was redirected through updated objective coefficients so that a large set of alternative solutions could be explored systematically.

In this implementation, random weights were assigned to activity variables across repeated runs, producing a broad exploration of the suboptimal decision space. Ex-post screening was then used to identify the solutions that were maximally different from one another. This approach was adopted after more conventional Hop-Skip-Jump MGA did not produce sufficiently different solutions in such a high-dimensional model.

Why use MGA?
A single optimum can create a false sense of certainty. In this work, near-optimal but structurally different solutions were treated as decision-relevant because externally imposed preferences or constraints may cause one of those solutions to become preferable in practice. MGA therefore served as a way to reveal alternatives that a conventional least-cost optimization would leave hidden.

The degree of diversity across MGA solutions was also informative in itself. A diverse set of near-optimal alternatives indicated that multiple pathways could satisfy the same target at similar cost. Highly similar alternatives, by contrast, indicated a narrower and more constrained decision space.

Results
The MGA analysis was run over 200 iterations across multiple slack values to characterize how much diversity existed in the near-optimal solution space. The central result was methodological as much as numerical: the feasible decision space changed sharply as even small amounts of slack were introduced. At very small slack, the model reproduced solutions that were nearly identical to the original optimum. As slack increased, the range of technology combinations broadened substantially, revealing multiple structurally different but still near-optimal pathways.

Figure 1 below shows 2050 technology-specific capacities across 200 MGA iterations for six slack values. At 0.01% slack, no major alternative emerged and the generated solutions remained nearly identical to the baseline optimum. At 0.05% and 0.1% slack, visibly different solutions began to appear, especially in utility-scale solar PV, nuclear, and combustion turbine capacities. As slack increased further to 1%, 2%, and 5%, the explored solution space widened substantially and a broader set of technology configurations became feasible within a near-optimal cost range.

Figure 1

Modeling to Generate Alternatives (MGA)

This is the central value of MGA. The method did not simply produce another scenario. It showed how rapidly the structure of the solution space changed once a small departure from strict optimality was allowed, and it made those tradeoffs visible in a way that a single least-cost run could not.

Transferability Beyond Energy
Although the original application was in a large-scale system planning model, the method is broadly transferable. The same logic applies wherever stakeholders need more than one answer from an optimization model.

Examples include:

  • supply chain network design
  • production planning
  • scheduling and workforce design
  • facility location
  • inventory strategy
  • routing and service design
  • capital planning under uncertainty

In these settings, MGA can help convert a model from a single-answer engine into a more useful decision-support tool. It can reveal whether several near-optimal strategies exist with different operational implications, implementation burdens, or risk profiles.

Takeaway
This case illustrates how optimization can be used not only to identify a best solution, but also to map the credible alternatives around it. MGA is particularly valuable when decision-makers care about robustness, strategic flexibility, and practical choice, not just objective-function optimality. By exposing near-optimal but structurally different options, the method makes hidden decision space visible and supports stronger decisions under uncertainty.

 

This Project was performed by Dr. Hadi.

Case Study # 7 - Breaking the Complexity Barrier in Multi-Component Maintenance Optimization

Overview

In capital-intensive industries, maintenance optimization is a high-stakes combinatorial puzzle. While grouping repairs saves on massive setup costs— for example ranging from $5,000 to $100,000 per hour in production losses in chemical plants—the resulting mathematical models are typically too complex to solve. I am among the very first to bridge this gap by formulating these challenges as multi-stage stochastic integer programs. By leveraging advanced decomposition methods and proving specific structural properties of the problem, I developed algorithms that increased the solvable instance size by 100x while keeping the optimality gap within 10%.

Challenge

Multi-component maintenance optimization is a notoriously complex field because it requires the seamless integration of stochastic failure processes with high-stakes combinatorial grouping problems. Also, the problem is challenging to solve due to decision-dependent (endogenous) uncertainties in the parameters.

More specifically:

  • Modelling Gap: Currently, there is a lack of general, efficient models or methods available to manage this type of endogenous uncertainty in multi-component systems.
  • Computational Gap: For a system with just 8 components, the number of potential scenarios becomes so vast that industry-standard solvers like CPLEX typically run out of memory or time before they can even approach an optimum. Most existing research is forced to rely on highly simplified assumptions—such as limiting analysis to two-component systems—or small-scale simulations that fail to capture the high-dimensional complexity of real-world industrial environments.

Figure 1. Exponential and decision-dependent scenario tree

Methods and Results

I addressed these challenges by developing a multi-stage stochastic programming model for multi-component system maintenance to minimize the total cost. I approximate this multi-stage model with a novel two-stage stochastic linear integer model in a rolling horizon.

To solve this two-stage stochastic model, we discussed solution methods in the following maintenance scenarios:

Time-Based Maintenance (TBM): Progressive Hedging & Shifting Windows

For systems relying on historical failure distributions, I implemented a heuristic-based Progressive-Hedging algorithm to find the high quality solution

  • Decomposition Strategy: I evaluated standard Benders Decomposition and the Integer L-shaped method, but found them inefficient for this problem; Benders cuts often fail due to the lack of strong duality in integer programs, and the L-shaped method requires excessive cuts for every feasible solution.
  • Proposal novel Heuristic: I developed a specialized heuristic within the Progressive Hedging framework that utilizes a shifting window to group maintenance activities in polynomial time. This allows the model to scale to a large number of components while maintaining an optimality gap of 10%.
  • Benchmark comparison: I further compared the proposed algorithms with two widely adopted benchmark algorithms when the solver computational time becomes impractical. Through an extensive numerical study, the proposed algorithm can achieve more than 40% of costing saving than the two benchmarks

Table 1. Cost savings of the proposed algorithm (Algorithm 4) compared to two benchmarks

heuristic-based Progressive-Hedging algorithm

Condition-Based Maintenance (CBM): Structural Partitioning

For systems with real-time sensor data, I reduced problem complexity based on the problem structural properties. I proved that the problem can be converted to a special set partitioning problem which finds the optimal maintenance partition by moving components between “do-nothing” and “maintenance” sets.  This bypasses the need for exhaustive scenario tree searches, delivering optimal results near-instantaneously. Table 1 shows the computational time  (in seconds) of the proposal algorithm compared to CPLEX w.r.t. different number of components n.

Table 2. the computational time  (in seconds) of the proposal algorithm compared to CPLEX w.r.t. different number of components n.

Condition-Based Maintenance (CBM): Structural Partitioning

Real-World Impact

Case 1: Wind Turbine Blade Degradation (3 blades, 10-year horizon)

  • Cost Parameters: PM = $200k/blade, CM = $600k-$1M/blade, Setup = $130k
  • Result: Considering economic dependence altered maintenance timing decisions, enabling high-value grouping opportunities. Potential savings: $260k-$390k per maintenance event through grouping.

Case 2: Crude Oil Pipeline Corrosion (17 pipelines, 5-year horizon)

  • Cost Parameters: PR = $5/pipe, CR = $20/pipe, Setup = $200
  • Result: 47 maintenance decisions altered when economic dependence was considered. With high setup cost relative to component costs, grouping effects dominate. Bootstrap analysis confirms robustness: mean cost $770.57 ± $260.5.

Broader Applications

The proposed solution methodology—using two-stage stochastic programming to handle endogenous parameter uncertainties—extends naturally to classical OR problems such as vehicle routing with stochastic demands, multi-product inventory control, and production scheduling with machine changeovers. The structural properties and decomposition techniques developed for maintenance optimization apply directly to these problem classes, where shared setup costs dominate decision-making and create similar grouping trade-offs.

Summary & Key Takeaways

Our work achieves a computational breakthrough by increasing solvable system sizes from 10–15 components to 100+ components, with 2,000×+ speedups over commercial solvers. The approach delivers 30–50% cost savings in real applications by explicitly modeling economic dependence—a factor ignored by standard methods. By combining stochastic programming with proven structural properties, the framework provides both exact solutions (for condition-based systems) and near-optimal solutions within a 10% optimality gap (for time-based systems), making advanced optimization immediately accessible to practitioners using standard reliability models.

This case study was performed by Dr. Zhu.

Case Study #8 - Collaborative Truck-Multi-Drone Pollution Routing for Cost and Emission Reduction in Last-Mile Logistics

Overview

Last-mile delivery is the most expensive (53% of total costs), time consuming and polluting (30% of total emission) segment of the shipping process. With the rapid expansion of e-commerce and urbanization, demand is expected to intensify further — potentially leading to a 21% increase in congestion and a 32% increase in emissions in the top 100 global cities by 2030. To address this, hybrid truck-drone delivery systems have emerged as a particularly promising solution. In these systems, the truck serves as a mobile depot that both delivers parcels and acts as a launch/retrieval station for drones, enabling collaborative operations that significantly outperform traditional truck-only routing in time, cost, and sustainability.

The Challenge

  • Traditional routing approaches fail to address the following real-world complexities in hybrid truck-drone delivery systems:
  • Variable truck speeds caused by traffic, terrain, and time-of-day conditions, which dramatically affect fuel consumption and pollution levels
  • Simultaneous coordination of one truck with multiple drones performing pickup-and-delivery tasks
  • Dynamic synchronization of truck movements, drone launch/recovery points, and customer time windows
  • Large-scale optimization requiring joint decisions on routing, customer assignment, and vehicle synchronization
  • High interdependence between all system components, where small inefficiencies can significantly impact overall cost, completion time, service reliability, and environmental performance

The Solution

ScienceDocs Optimization Engineering Consultant Dr. Arash developed the Truck Multi-Drone Pollution Routing Problem with Pickup and Delivery under Variable Truck Speeds (TMD-PRP-PDVS). Using mixed-integer linear programming and advanced metaheuristics, the model simultaneously optimizes:

  • Truck routes and variable speeds
  • Drone launch and recovery locations
  • Assignment of deliveries between truck and drones
  • Total cost (fuel, drone energy, driver labor)

To ensure both accuracy and scalability, smaller problem instances were solved to optimality using Gurobi, providing strong benchmarks. For large-scale, real-world scenarios, an Adaptive Large Neighborhood Search (ALNS) framework was developed with problem-specific destroy-and-repair operators. Additionally, a Speed-Level Adjustment for Truck (SLAT) mechanism was developed and embedded within the ALNS to dynamically adjust truck speeds during the search process. Sensitivity analyses were conducted across key parameters including drone fleet size, drone energy capacity, and customer density, to evaluate robustness and quantify the impact of operational decisions under different environments.

Key Visuals from the Study

truck-multi-drone system

Figure 1: Schematic diagram of the collaborative truck-multi-drone system, showing how drones are launched and recovered from the moving truck while serving multiple customers.

 

Real-world route visualizations

Figure 2: Real-world route visualizations on city maps ((a) Phoenix, USA (LDCS) and (b) Bengaluru, India (EDCS)), illustrating the optimized synchronized paths of the truck and two drones under variable speed conditions and different cost structures. Values next to the arcs are truck speed levels (mph).

 

Cost and emission comparison

Figure 3: Cost and emission comparison for an instance, demonstrating the significant reductions achieved by the proposed model versus traditional truck-only routing. Values next to the arcs are truck speed levels (mph).

Results & Impact

The optimized solution produced clear, measurable improvements across different cost environments. Under both labor-dominant cost settings (LDCS)—typical in Western economies with high driver wages—and energy-dominant cost settings (EDCS) —common in regions with higher fuel costs—the hybrid truck–drone system consistently outperformed traditional truck-only operations.

  • Substantial reduction in total operational cost (up to ~38–43% cost saving), with similar benefits observed in both LDCS and EDCS environments.
  • Additional cost improvements from intelligent truck speed optimization (≈7–10% on average), adapting differently across cost structures (faster speeds in LDCS vs. fuel-efficient speeds in EDCS)
  • Better on-time performance even under realistic road condition variability through synchronized truck–drone coordination.
  • Scalable performance for real fleet sizes (validated up to 100 customers), where collaborative savings remain significant across scenarios.
  • Sensitivity studies show that increasing drone fleet size and battery capacity amplifies collaborative savings, while speed‑optimization consistently provides 5–10% extra savings across different fleet sizes, energy capacities, and service area sizes.

The model was tested on real-world city maps of Phoenix (USA) and Bengaluru (India). The results from both cities fully verified the overall findings, demonstrating that the collaborative truck-multi-drone approach with variable speed optimization delivers robust cost and emission reductions regardless of geographic and cost-structure differences.

Takeaway

This work demonstrates how advanced mathematical modeling and decision support systems can transform last-mile logistics — delivering networks that are simultaneously cheaper, faster, and greener. By integrating pollution-aware routing with collaborative truck-drone operations, the framework provides logistics operators with a practical, scalable tool for sustainable and cost-effective operations.

A key insight is that truck speed adjustment and drone utilization are the primary drivers of cost savings. In LDCS, the optimal strategy is to allow the truck to travel at higher speeds, which significantly reduces labor costs. In contrast, in EDCS, lower truck speeds combined with longer routes and greater drone utilization yield the greatest overall savings. The model shows that jointly optimizing routing, variable truck speeds, and drone deployment rather than treating these decisions independently, is essential for achieving maximum efficiency and resilience in complex, dynamic delivery networks.

Figures and results adapted from: Alizadeh, A., Srinivas, S., Noble, J., & Otto, A. (2026). Collaborative truck multi-drone pollution routing problem with pickup and delivery under variable truck speeds: a cost minimization approach. International Journal of Production Research, 1-33.

This case study was performed by Dr. Arash.

Search ScienceDocs

Find grants, services, and resources across our site.