HVRPTW-Dataset

Hydrogen Vehicle Routing Problem with Time Windows (HVRPTW) – Projection

1. Energy Consumption Model for HVRPTW

In the Hydrogen Vehicle Routing Problem (HVRP) literature, most studies assume that energy consumption is proportional to the traveled distance [4] , [3]). However, this simplifying assumption deviates from real-world vehicle behavior, as it ignores dynamic factors that significantly influence energy and hydrogen consumption.

To better represent real operating conditions of fuel cell electric vehicles (FCEVs), we propose an energy consumption model based on the vehicle’s longitudinal dynamics. Unlike distance-only approaches, this model incorporates several physical and operational parameters, including:

• Total vehicle mass (including empty weight and carried load)
• Road slope (grade)
• Travel speed
• Acceleration and deceleration phases
• Number of intermediate stops (traffic lights, intersections, congestion, slowdowns, etc.)
• Characteristics of the electric drivetrain, especially efficiency factors
• Energy recovery through regenerative braking
• Type of road network

---

2. Base Benchmark

Given the limited number of hydrogen refueling stations [7] and the absence of a standard benchmark for the HVRPTW, we assume that refueling occurs only at the depot. Consequently, the problem becomes close to the classical VRPTW, with the additional constraint of limited hydrogen vehicle autonomy.

The most widely used benchmark instances for the Vehicle Routing Problem with Time Windows (VRPTW) are the Solomon instances. They consist of 56 instances with a single depot and 100 customers, each defined by demand, time windows, and service times, under vehicle capacity and route duration constraints.

Key characteristics:

• Each instance is defined on a synthetic 100 × 100 grid [5]
• Speed is normalized so that 1 distance unit equals 1 time unit; thus, distance, travel time, and cost are equivalent [8])
• The depot also has a time window, and vehicles must complete routes within it [8])
• Distances between customers are computed using Euclidean distance [8])
• No explicit energy consumption model is defined in these instances, as this concept was introduced later [11])

---

2.1 Instance Classes

All instances within a class share the same customer locations.

[5] , [8]

---

2.2 Instance Types

[5]

---

2.3 Instance Structure

Example: C101_200.csv

3. Benchmark Adaptation

To use the Solomon (1987) [11]) instances in our optimization model, we adapted the benchmark by enriching it with additional parameters required to compute energy consumption on each arc.

The main objective is to transform the original benchmark, which consists of abstract Cartesian coordinates without any real geographical meaning, into a realistic road network representation. This enables the computation of hydrogen vehicle energy consumption.

---

3.1 Definition of the Projection Geographical Area

Projection City

The city of Lyon was selected as the projection area due to its role as one of the main logistics and economic hubs in France, characterized by high urban and suburban mobility demand [1]).

Its road network includes a wide variety of infrastructures, making it suitable for reproducing realistic delivery conditions.

Moreover, Lyon has a diverse topography, with hills and significant elevation changes. This altimetric heterogeneity is particularly relevant for this study, as the energy consumption model explicitly accounts for the effect of slopes on vehicle energy demand [2]).

---

Projection Area

Lyon is approximately located at:

• Latitude ≈ 45.75
• Longitude ≈ 4.83

We define a rectangular bounding box covering most of the Lyon metropolitan area, including the city center and part of the surrounding suburbs:

• Latitude range: [45.72 ; 45.80]
• Longitude range: [4.79 ; 4.90]

---

3.2 Normalization of Solomon Coordinates

Since Solomon instance coordinates are artificial, they must first be transformed into a homogeneous distribution in [0,1] to allow a controlled and consistent projection onto the Lyon geographic area.

For each instance, we compute:

<img src="images/1.png" width="200" height="100"/>

Each coordinate is then normalized as:

<img src="images/2.png" width="500" height="400"/>

---

3.3 Affine Transformation to the City of Lyon

The affine transformation resizes the dataset while preserving geometric proportions and spatial relationships between points.

For each coordinate, we apply:

<img src="images/3.png" width="400" height="300"/>

---

3.4 Correction of Projected Coordinates

After projection, some customers may end up located in unrealistic positions (e.g., parks, buildings, rivers such as the Rhône or the Saône). To address this issue, we use the OSRM (OpenStreetMap Routing Machine) API to snap each point to the nearest road segment.

For each customer (i):

<img src="images/4.png" width="100" height="50"/>

We compute:

<img src="images/5.png" width="150" height="100"/>

where:

• <img src="images/6.png" width="20" style="display:inline-block; vertical-align:middle;"> : a point on the road network
• <img src="images/7.png" width="50" style="display:inline-block; vertical-align:middle;"> : distance between the two points
• <img src="images/8.png" width="50" style="display:inline-block; vertical-align:middle;"> : argument minimizing the distance

The final corrected position is:

<img src="images/9.png" width="150" height="70"/>

These snapped coordinates become the official customer locations.

Finally, the original Solomon instance coordinates are replaced with these projected coordinates to generate the new adapted datasets.

---

Example

C101_200_coords.csv

3.5 Benchmark Enrichment

The projection of customers onto a real road network provides geographically consistent coordinates within the urban environment of Lyon. However, this information alone is not sufficient to compute hydrogen vehicle energy consumption.

For this reason, the Solomon benchmark was enriched by generating a second file (JSON format) containing detailed characteristics of each arc connecting two customers.

For each node pair (i, j), the following information is computed and stored:

---

---

Acceleration Model

Due to the lack of experimental data on acceleration and deceleration profiles for light-duty fuel cell electric vehicles (FCEVs), we use values reported for battery electric vehicles (BEVs) with similar characteristics. This substitution is justified by the close architectural similarity between both technologies:

• A BEV uses an electric motor, battery pack, and controller, while a FCEV (or fuel cell electric vehicle) uses the same electric drivetrain architecture, except that energy is supplied by a hydrogen fuel cell instead of a rechargeable battery [6] , [9].
• Longitudinal dynamics follow Newton’s second law regardless of powertrain type, where the traction force required to overcome resistances and produce acceleration is determined by the vehicle mass (m) and motor torque, which depends on wheel radius [10].

---

Acceleration and Deceleration Profiles

[10]

---

8. References

<a id="ref-1"></a> [1] Métropole de Lyon. (2024). Document d’Orientations sur la Logistique des biens et des services (DOLB&S). Direction des Mobilités.
https://www.grandlyon.com/fileadmin/user_upload/media/pdf/deplacements/orientation-logistique-urbaine.pdf

<a id="ref-2"></a> [2] Métropole de Lyon. (2020). Plan de Protection de l'Atmosphère et orientations topographiques urbaines / Aire de mise en valeur de l’architecture et du patrimoine.
https://www.grandlyon.com/fileadmin/user_upload/media/pdf/urbanisme/sites-patrimoniaux/spr_avap_croix_rousse_diagnostic.pdf

<a id="ref-3"></a> [3] Abibou, M., et al. (2024). The Hydrogen Vehicle Routing Problem with Time Windows: Formulations and distance-proportional approximations. European Journal of Operational Research.

<a id="ref-4"></a> [4] Abibou, M., et al. (2025). Optimizing green logistics: Models and metaheuristics for hydrogen fleet deployment. Computers & Operations Research, 164, 106-121.

<a id="ref-5"></a> [5] Gendreau, M., & Tarantilis, C. D. (2010). Managing Bio-inspired and Metaheuristic Algorithms for Vehicle Routing Problems. In Managing Green Logistics, Springer, pp. 115-142.

<a id="ref-6"></a> [6] Gondal, I. A., et al. (2018). Comparative Review of Energy Drivetrains: Fuel Cell Electric Vehicles versus Battery Electric Vehicles. International Journal of Hydrogen Energy, 43(11), 5912-5929.

<a id="ref-7"></a> [7] Guo, X., et al. (2025). Strategic Planning and Location Optimization for Hydrogen Refueling Infrastructure under Autonomy Constraints. Transportation Research Part E: Logistics and Transportation Review, 182, 103-124.

<a id="ref-8"></a> [8] Kumar, R., et al. (2012). On the Equivalence of Distance, Time, and Energy in Normalized Network Space for Routing Problems. Journal of Operational Research Society, 63(8), 1089-1102.

<a id="ref-9"></a> [9] Parikh, S., et al. (2023). Performance Evaluation and Architectural Similarities of Electric Powertrains: BEV vs. FCEV. IEEE Transactions on Transportation Electrification, 9(3), 3120-3135.

<a id="ref-10"></a> [10] Settey, T., et al. (2021). Application of Longitudinal Vehicle Dynamics and Newton’s Second Law to Powertrain Efficiency Assessment. Transport, 36(4), 289-301.

<a id="ref-11"></a> [11] Solomon, M. M. (1987). Algorithms for the Vehicle Routing and Scheduling Problems with Time Window Constraints. Operations Research, 35(2), 254-265.Constraints*. Operations Research, 35(2), 254-265.