An Alternative Fuel Refueling Station Location Model considering Detour Traffic Flows on a Highway Road System

With the development of alternative fuel (AF) vehicle technologies, studies on finding the potential location of AF refueling stations in transportation networks have received considerable attention. Due to the strong limited driving range, AF vehicles for longdistance intercity trips may require multiple refueling stops at different locations on the way to their destination, which makes the AF refueling station location problem more challenging. In this paper, we consider that AF vehicles requiring multiple refueling stops at different locations during their long-distance intercity trips are capable of making detours from their preplanned paths and selecting return paths that may be different from original paths for their round trips whenever AF refueling stations are not available along the preplanned paths. ese options mostly need to be considered when an AF refueling infrastructure is not fully developed on a highway system. To this end, we first propose an algorithm to generate alternative paths that may provide the multiple AF refueling stops between all origin/destination (OD) vertices. en, a new mixed-integer programming model is proposed to locate AF refueling stations within a preselected set of candidate sites on a directed transportation network by maximizing the coverage of traffic flows along multiple paths. We first test our mathematical model with the proposed algorithm on a classical 25-vertex network with 25 candidate sites through various scenarios that consider a different number of paths for each OD pair, deviation factors, and limited driving ranges of vehicles. en, we apply our proposed model to locate liquefied natural gas refueling stations in the state of Pennsylvania considering the construction budget. Our results show that the number of alternative paths and deviation distance available significantly affect the coverage of traffic flows at the stations as well as computational time.


Introduction
Reducing greenhouse gas (GHG) emissions in the transportation sector is one of the most vital steps in fighting against global warming in the United States (U.S.). According to the U.S. Environmental Protection Agency [1], the transportation sector generates the largest share of GHG emissions. In order to cut down tail-pipe emissions in the transportation sector, vehicles using alternative fuel (AF), such as biodiesel, hydrogen, electrical energy, and natural gas, have received significant attention in recent years because they emit less well-to-wheel GHG than that of vehicles using traditional fossil fuels, such as diesel and gasoline. Currently, 367 light-duty and 216 medium-and heavy-duty vehicles for the 2018 and 2019 model years are available in the U.S. AF vehicle market [2].
While the public interest in using AF vehicles instead of conventional vehicles has increased, the number of public refueling stations for AF vehicles is still insufficient, especially for intercity trips between urban and rural counties. Table 1 shows the 2010 census total population, number of electric charging stations, and the number of electric charging stations per 100,000 residents in the U.S. based on the Census Bureau's urban-rural classification [3,4]. In this table, counties are categorized into three groups according to their population density: mostly urban, mostly rural, and completely rural. It indicates that urban areas have more than twice as many electric charging stations than rural areas with a similar population size. One of the barriers to the invigoration and development of AF infrastructures in highway systems, which play a major role in intercity trips among urban and rural counties, is high construction cost. For example, the construction cost of one natural gas refueling station is at least $1.8 million in the Pennsylvania Turnpike [5].
Due to the sparse distribution of AF refueling stations between urban areas in highway systems, AF vehicles with a short driving range that travel long-distance intercity trips may need to use longer paths with refueling availability, including multiple refueling stops at different locations on the way to their destination, rather than their shortest paths without available refueling stations, so they can safely complete their trips. us, AF vehicle drivers may need to make detours to be able to refuel their vehicles.
Generally, road structures in highway systems are different than those in other transportation networks. First, highway roads are divided into two pathways separated by either a raised barrier or an unpaved median. In order to offer uninterrupted traffic flow, highway roads do not have any traffic intersections, and vehicles are only able to enter/leave highway systems through entrance and exit ramps. Next, highway systems have built-in service facilities where drivers are able to take a rest and refuel their vehicles. Since highway roads are physically partitioned by a barrier or an unpaved median, some built-in service facilities, called single-access stations, can only be accessed from one side of the road, while the rest, called dual-access stations, can be accessed from both sides of the road. Hwang et al. [6,7] and Ventura et al. [5] model this type of transportation systems as directed transportation networks.
In general, concerning the design of an AF refueling infrastructure along transportation networks, a number of studies allowing repeated trips between origin/destination (OD) vertices assume that vehicles make symmetric round trips traveling along preplanned (shortest) paths between the corresponding OD vertices. It also assumes that the set of candidate station locations is a subset of vertices in the network, and therefore, all the candidate locations are dual-access sites. ese assumptions imply that, if a vehicle travels from an origin to a destination on a single path by filling up at some stations, it also stops by the same stations on the return path back to the origin. ese assumptions, however, have made the refueling station location problem for AF vehicles less practical since they do not reflect the characteristics of AF vehicles well at its early stage of market.
By relaxing the assumptions listed above, this paper aims at determining more reliable locations of AF refueling stations in real-world applications based on the distinct features of AF vehicles traveling intercity trips on directed transportation networks. We first consider that some candidate sites for AF refueling stations are single-access and drivers may choose different return paths from original paths to be able to refuel their vehicles in both directions. Also, since several paths are available between ODs, AF refueling availability is considered as one of the AF vehicle drivers' top priorities when they select paths to travel in highway systems. us, we allow AF vehicles to make nonsymmetric round trips between their ODs. We first generate multiple paths between all OD pairs through a revised -th shortest path algorithm considering a maximum deviation distance. en we formulate a new mixed-integer programming model that considers a set of predetermined paths for each OD pair and a preliminary set of candidate station location sites, including single-access and dual-access sites, on a directed transportation network in which AF vehicles are able to use any of the corresponding OD paths depending on the availability of refueling service. e proposed model determines the optimal set of station locations for a given number of stations and the selected round trips for all ODs that maximize the total traffic flow covered (in round trips per time unit) by the stations. For computational experiments, we first apply the proposed model to a classical 25-vertex network with 25 candidate sites through various scenarios. We also validate our proposed model with a budget constraint to construct AF refueling stations in the state of Pennsylvania.
Our model proposed in this study is applicable to various types of AF vehicles, especially to liquefied natural gas (LNG) vehicles, with their refueling station location problems. LNG vehicles are similar to the existing long-haul vehicles powered by diesel in terms of powertrain and refueling, but LNG vehicles are known to provide economic and environmental benefits; thus, LNG vehicles are well-suited for replacing the current long-haul vehicles powered by diesel on a highway road system. For example, UPS has been working to shi their high carbon-fueled vehicles to new generation of LNG vehicles since the late 20th century and has been increasing their number [8]. In the U.S., UPS has deployed their LNG vehicles mainly in Indianapolis, Chicago, Earth City, and Nashville, and plans to use these LNG vehicles in larger areas [9]. In addition to LNG vehicles, battery-electric vehicles and hydrogen fuel cell vehicles will also be applied to our proposed model once they are fully available for long haul logistics on a highway system. e remaining of this paper is organized as follows. Section 2 reviews the related literature and shortly discusses the main distinctions of our research work over the existing studies as well. In Section 3, we first introduce the AF refueling station location problem with detour traffic flows on a highway road system and provide properties of feasible paths. Next, an algorithm that generates a set of multiple feasible paths for which the properties are satisfied on a given network is proposed. In Section 4, we provide covering conditions with candidate sites to cover round trips for each OD pair. en, we propose a mixed-integer programming model to locate a given number of AF refueling stations on a directed transportation network with the objective of maximizing the total traffic flow covered, considering multiple paths for each OD pair. In Section 5, the proposed model is tested on a well-known 25-vertex network to evaluate the effects of the number of multiple paths, maximum deviation distance, and vehicle driving range on the coverage of traffic flows. e proposed model is then validated in Section 6 with an application to the state of Pennsylvania to demonstrate its performance on a large-size problem. Lastly, we provide conclusions and discuss the future work of AF refueling station location problems in Section 7.

Literature Review
In this section, we take a look at the literature relevant to this paper. First, in Subsection 2.1, we review the literature addressing -shortest path problems, and in Subsection 2.2, we examine the literature related to AF refueling station location problems. en, we organize the main distinctions of our research work by comparing with the relevant literature in Subsection 2.3.

-shortest Path Problems.
e -shortest path problem is to find shortest paths between two vertices in a given network in a nondecreasing order of length, where refers to the number of shortest paths to find. e -shortest path problem can be classified into two types according to whether paths are allowed to have cycles or required to be simple with no cycles. e first type of -shortest path problem, proposed by Hoffman and Pavley [10], allows repeated vertices along any path. Fox [11] develops an algorithm to apply this type of problem to probabilistic networks. Eppstein [12] uses the concept of binary heap data structure in a nondecreasing order of additional path length due to deviation to find the shortest paths in polynomial time. Since the experimental results of Eppstein's algorithm still take considerable time to find the shortest paths, Jiménez and Marzal [13] present a modified version of Eppstein's algorithm to improve its practical performance. Aljazzar and Leue [14] propose a directed search algorithm to search for the shortest paths between two given vertices of a network.
eir algorithm provides the same asymptotic worst-case complexity but uses less memory than Eppstein's algorithm. Liu et al. [15] propose a novel -shortest path algorithm for neural networks. Given a set of battery exchange stations, Adler et al. [16] suggest a polynomial time algorithm to solve the electric vehicle shortest-walk problem with battery exchanges considering vehicle's limited driving range. is algorithm shows the chance of extension to the -shortest path problem by transforming the original traffic network into the so called refueling shortest path network. e second type of -shortest path problem does not allow any repeated vertex along a path, and is thus called the -shortest simple (or loopless) path problem. Since this type of problem needs an additional constraint to allow only loopless paths, it turns out to be more challenging than the first type of problem [17]. Pollack [18] introduces the concept of the -shortest simple path problem and solves it by modifying Hoffman and Pavley's method [10], so as to avoid paths containing repeated vertices. Clarke et al. [19] present a branch-and-bound procedure to find the -shortest simple paths, but this method requires a significant computational effort and storage requirements in the main memory. Sakarovitch [20] first identifies several shortest paths that may contain repeated vertices by using the efficient version of Hoffman and Pavley's method [10], and then picks up -shortest simple paths among them. By applying a procedure that partitions a path into two subpaths, Yen [21] develops an efficient algorithm to find the -shortest simple path. Lawler [22] generalizes a procedure that can reduce the amount of storage required in solving the -shortest simple path problem. Katoh et al. [23] present an improved version of Yen's algorithm that solves the problem efficiently in an undirected network. e practical performance of Yen's algorithm is comparatively analyzed with other -shortest path algorithms [24,25]. An implementation of Yen's algorithm is also studied to improve its practical performance [17,26]. Hershberger and Suri [27,28] suggest the efficient replacement path algorithm for finding shortest simple paths in a directed network. Zeng et al. [29] present a heuristic -shortest path algorithm that is based on Yen's algorithm when determining an eco-friendly path that results in minimum carbon dioxide emissions from light-duty vehicles.

Alternative Fuel Refueling Station Location
Problems. e maximal covering location and the set-covering location approaches are two main streams of research for addressing AF refueling station location problems. Kuby and Lim [30] are one of the first applying a maximal covering location model to solve an AF refueling station location problem. ey introduce the flow-refueling location model (FRLM) to find the optimal location of refueling stations for AF vehicles by considering their limited driving range per refueling with the objective of maximizing the total traffic flow covered. Upchurch and Kuby [31] show that the FRLM identifies more stable locations for AF refueling stations than the -median model, especially in a statewide case study. In general, the FRLM requires a significant computational effort to pregenerate all feasible location combinations of refueling stations that allow vehicles to make round trips between ODs. Lim and Kuby [32] propose three heuristic versions for the FRLM. Kuby et al. [33] apply two of them to locate hydrogen refueling stations in Florida. Capar and Kuby [34] present a new formulation for the FRLM that skips the pregeneration of all feasible combinations on every path. Capar et al. [35] suggest an arc-cover-path-cover model that focuses on the arcs comprising each path, so as to solve the problem efficiently without the pregeneration of all In general, drivers o en deviate from their original paths of shortest travel time or distance to be able to refuel their vehicles [53]. Both the maximal covering location and the set-covering location approaches have been extended to consider driver's deviation options under a variety of situations. Kim and Kuby [54] address the deviation version of the FRLM (DFRLM) where drivers are able to deviate from their paths of shortest length between ODs, and Kim and Kuby [55] propose a heuristic algorithm for the DFRLM to solve large-scale problems. Huang et al. [56] develop a new model, called the multipath refueling location model (MPRLM), by considering multiple deviation paths between ODs. For intracity trips that require at most one refueling stop, Miralinaghi et al. [47,48] suggest deviation versions of the set-covering location problem to find potential locations of AF refueling stations. For intercity trips of large-scale problems, Yıldız et al. [57] use a branch and price algorithm, which does not need pregeneration of path generation, for the AF refueling station location problem adding the routing aspect of drivers. Most recently, Göpfert and Bock [58] and Arslan et al. [59] suggest novel projection and branch and cut methods in dealing with the deviation version of the refueling infrastructure planning, so as to extend the computational efficiency even further to solve large-size problem instances with less computational effort.

Main Distinctions of Our Research Work.
Comparing with the literature listed above, we shortly provide the main distinctions of our research work over the existing studies as follows: (i) Our problem is an extension of Hwang et al. [6] problem to consider potential deviation paths on directed transportation networks, such as highway network systems. is leads to consider (1) the mixed set of single-access and dual-access candidate sites to locate AF refueling stations, and (2) nonsymmetric round trips between ODs, where return paths are allowed to be different from original paths for refueling services in both directions. (ii) Our study is well-suited for the AF refueling station location problem, specially with LNG vehicles traveling long-distance intercity round trips. Some number of recent studies, including Miralinaghi et al. [47,48], limits this type of problem only for AF vehicles with intracity trips, which needs at most one refueling stop per trip. On the other hand, we apply covering condition procedures depending on LNG vehicle driving range, so as for LNG vehicles traveling intercity trips to allow multiple refueling stops at different locations. (iii) In our research work, paths are not fixed for every OD pair. Instead, paths are flexible to consider detour traffic flows. While Kim and Kuby [54,55] and Huang et al. [56] apply Hoffman and Pavley's algorithm [10] and Yen's algorithm [21] to take deviation paths into account, we develop a new algorithm based on Eppstein's algorithm [12] to rigorously and efficiently find shortest paths allowing repeated vertices along paths within the tolerance (i.e., maximum feasible location combinations of stations on each path for the FRLM. Jochem et al. [36] apply the Capar et al. [35] model to allocate charging stations in the German autobahn. While most AF refueling station location problems assume that the AF refueling stations can only be located at the vertices, Kuby and Lim [37] and Ventura et al. [38] consider additional candidate sites along arcs. Kweon et al. [39] extend the approach suggested in Ventura et al. [38] to locate a refueling station anywhere along a tree network to the case where a portion of drivers are willing to deviate to receive refueling service. He et al. [40] propose a bilevel model to solve the optimal locations of electric charging stations, through taking the driving range limitation of an electric vehicle, the battery charging time required, and the situation in which some electric vehicle drivers possibly charge at home into account. While Kuby and Lim [30] and the following studies solve the AF refueling station location problem using a maximal covering location problem, Wang and Lin [41] solve this problem using a set-covering model with the objective of minimizing the total cost of locating stations to cover all the traffic flow on a given transportation network. Wang and Wang [42] integrate Wang and Lin's model [41] into the classic set-covering model considering vertex-based and flow-based demands for the AF refueling station location problem. Since Wang and Lin's model [41] requires a significant computational effort to evaluate the effect of the limited vehicle driving range on the number of charging stations needed for achieving multiple origin-destination intercity travel with electric vehicles on Taiwan, MirHassani and Ebrazi [43] propose a novel approach by using the conservation of flow law, which is able to solve large-scale problems. Chung and Kwon [44] extend MirHassani and Ebrazi's model [43] to a multiperiod planning problem for allocating charging stations in the Korea Expressway. Hosseini and MirHassani [45] use the idea of MirHassani and Ebrazi's model [43] to propose a two-stage stochastic mixed integer programming model for the refueling station location problem, where the traffic flow of AF vehicles is uncertain and portable AF refueling stations are considered. Kang and Recker [46] use the idea of the set-covering problem to locate hydrogen refueling stations with the assumption that at most one refueling stop per trip is required in a city. Assuming that vehicles only require one refueling stop per trip, two refueling station location models, including the versions that consider limited capacity of refueling stations [47], and refueling demand uncertainty and driver route choice behavior [48], are developed to minimize the total cost imposed on a planner and drivers over multiple time periods. Using the Adaptive Large Neighborhood Search algorithm [49] and the Adaptive Variable Neighborhood Search algorithm [50], locations for battery swap stations and electric vehicle routes are determined to provide services with the objective of minimizing the sum of the station construction cost and routing cost. To minimize the total cost to locate electric vehicle charging stations in road networks, Gagarin and Corcoran [51] suggest a novel approach that searches for the dominating set of locations among the candidate locations whose distance is below a certain threshold from a given driver. Using a parallel computing strategy, Tran et al. [52] propose an efficient heuristic algorithm for location of AF refueling stations based on the solution of a sequence of subproblems.  [54] Kim and Kuby [55] Huang et al. [56] Miralinaghi et al. [47] Miralinaghi et al. [48] Hwang et al. [6] Hwang et al. [7] Ventura et al. [5] Deviation Allowed Allowed Allowed Allowed Allowed Allowed  [55] Huang et al. [56] Miralinaghi et al. [47] Miralinaghi et al. [48] Hwang et al. [6] Hwang et al. [7] Ventura et al. [5] Objective(s) and (2) to analyze the effects of number of alternative paths including deviation paths, maximum deviation distance, and vehicle driving range on the optimal location as well as the corresponding total traffic flow covered. e major differences between the proposed model and the existing studies that are directly relevant to ours are summarized in Table 2.

The AF Refueling Station Location Problem with Detour Traffic on a Highway Road System
In this section, we first introduce the AF refueling station location problem with detour traffic flows on a highway road system, where AF vehicles are able to make detours for refueling and to select different paths between original and return trips. Next, we provide three small instances to describe concepts of feasible paths on a directed simple network in this problem. Also, based on the examples, we derive four properties of feasible paths. en, an algorithm is presented to determine multiple paths for which the four properties are satisfied for all ODs. is paper aims to locate a given number of AF refueling stations on a directed transportation network so as to maximize the coverage of detour traffic flows by considering multiple paths between ODs. To this end, a mixed-integer programming model will be presented in Section 4.

Problem Statement.
We define the problem on a simple directed network 퐺(푉, 퐸), where 푉 = v |푖 = 1, ⋅ ⋅ ⋅ , 푛 is the set of vertices for ODs, 퐸 = v , v |for some v , v ∈ 푉 is the set of arcs having nonnegative lengths, |푉| = 푛 , and |퐸| = 푛 . e road network in this problem has neither loops nor multiple parallel road segments in the same direction between any pair of vertices.
pairs. For any 푞 v , v ∈ 푄, vehicles perform the same round trip, which can be divided into the original path 푃 v , v from origin v to destination v , and the return path are defined as the sets of original and return paths, respectively. Let 푉 푃 v , v be the set of constituent vertices in path 푃 v , v ; then, each path can be represented as a sequence of vertices such that to v , respectively, then the corresponding original and return paths are denoted as 푃 v , v and 푃 v , v . Also, we define 퐹 v , v and 퐹 v , v as the -th shortest feasible original and return paths for OD pair 푞 v , v . In Subsections 3.2 and 3.3, we will discuss in detail the properties of feasible paths.
Next, the locations of candidate sites for AF refueling stations, denoted as 푀 = 푚 |푡 = 1, ⋅ ⋅ ⋅ , 푛 , are assumed to be predetermined, where some candidate sites can only be accessed from one side of the road (i.e., single-access) and the rest can be accessed from both sides of the road (i.e., dualaccess). We define Furthermore, the distance between two candidate sites (푙) and Figure 1 shows a representation of vertex and candidate site sequences, as well as distances between vertices and candidate sites. We consider that vehicles make a complete round trip between their ODs. ey have a limited driving range under free flow conditions, denoted as , which refers to the maximum travel distance with a single refueling. Since vehicles' home locations or final destinations are generally far away from highway interchanges, we assume they have at least a half-full tank when they enter and exit the road network. From now on, we call this assumption the half-full tank assumption. e half-full tank assumption was first introduced by Kuby and Lim [30], considering that data about the actual fuel tank level of vehicles when entering and exiting a highway road system are hard to obtain or likely to be inaccurately estimated.
is assumption has been then followed by the existing literature. e half-full tank assumption originally aims at ensuring that vehicles can repeat round trips several times without running out of fuel during their round trip. at is, a vehicle accessing the last refueling station and reaching the origin/ is calculated by multiplying the length of the shortest path by a positive deviation factor, denoted as . en, a predetermined maximum number of paths, denoted as , whose length does not exceed the maximum deviation distance is nominated for each driving direction. is implies that, if a shortest distance path does not provide refueling service, then drivers can select up to alternative paths for the same OD pair depending on their refueling availability and the limited additional deviation travel distance. Table 3 summarizes the destination with at least a half-full tank is able to make the original/return trip from the origin/destination with a half-full tank and access the same station without running out of fuel. e half-full tank assumption in our study similarly makes vehicles refueled at stations positioned within a distance of 푅/2 from their origin interchanges and again at stations placed within the same distance from their destination interchanges in both the original and return trips. In this respect, we define ing three sets of candidate sites in 퐶푆 퐹 v , v , which are categorized depending on their distances from v and v : When a feasible set of AF refueling stations is located in , the corresponding paths can be covered by the refueling stations. In Section 4, we will discuss in detail the covering conditions that depend on path lengths. e half-full tank assumption can be relaxed using Hwang et al. 's [7] model to consider different fuel tank levels of vehicles at their origins and destinations when detours are available on directed transportation networks, but it is le for future research to mainly focus on addressing the refueling station location problem with deviation options on a highway road system in this study.
In general, drivers deviate from their preferred paths, e.g., the least time or shortest distance paths, in as short a distance as possible if the preferred paths do not offer refueling availability. Also, a GPS navigation system provides a certain number of routes to travelers who detour from their familiar paths. In this respect, we first consider that vehicles can deviate from their shortest path up to a maximum deviation distance, which Original and return paths between interchanges v and v 퐿 , 퐿 Sets of original and return paths Distance between two candidate sites (푙) and Set of candidate sites that are located beyond a distance of 푅/2 from both inter- Vehicle's maximum travel distance with a single refueling Positive deviation factor Predetermined maximum number of paths cannot reach v 11 with at least a half-full tank when they exit the network. For the shortest return path, vehicles that reenter the network with a half-full tank at v 11 cannot reach the AF refueling station at site 3 . us, vehicles in 푞 v 6 , v 11 would need to make detours using vertex sequence ⟩ for the original path and opposite sequence the return path. ese two deviation paths are feasible paths for 푞 v 6 , v 11 . Lastly, suppose that sites 5 and 6 are selected as AF refueling stations. Vehicles in 푞 v 12 , v 17 cannot be refueled by these two AF refueling stations if they use the shortest original and return paths. However, if vehicles make multiple cycles at vertex v 13 , then vertex sequence (6) ⟩ also provides AF refueling availability for the return path, but drivers would not choose this vertex sequence because this path is longer than the path defined by the previous vertex sequence ⟨v 17 From these observations, we can state four feasible path properties in the next subsection.

Four Properties of Feasible Paths.
We derive four properties of feasible paths on a directed transportation network. Based on these four properties, an existing -th shortest path algorithm will be modified in the next subsection to generate multiple feasible paths for all OD pairs. Since we consider round trips, relevant notation and parameters. Note that data about the actual values of deviation factor and maximum number of paths for each OD pair are difficult to obtain or likely to be inaccurately predicted. Besides, different drivers would have different standards for the proper values of and , as well as different vehicle driving ranges . us, this study does not fix their values. Instead, in Section 5 we change their values for a given number of refueling stations to demonstrate the coupled effects of these parameters on the coverage of OD traffic flows.

ree Small
Instances. e concept of feasible path on a directed transportation network is illustrated with three small disconnected networks with seventeen vertices v 1 , ⋅ ⋅ ⋅ , v 17 , two dual-access candidate sites 푚 1 , 푚 3 , and five singleaccess candidate sites 푚 2 , 푚 4 , 푚 5 , 푚 6 , 푚 7 in Figure 2. Suppose we are trying to determine whether vehicles with 푅 = 100 can perform round trips without running out of fuel.
First, since vehicles are considered to have at least a half-full tank at theirs ODs, for 푞 v 2 , v 4 , shortest paths 푃 푆 v 2 , v 4 and 푃 푆 v 4 , v 2 are feasible because the corresponding trips can be covered by placing a single refueling station at site 1 . However, for 푞 v 1 , v 3 , since candidate sites are available in neither 푃 푆 v 1 , v 3 nor 푃 푆 v 3 , v 1 , drivers need to deviate from these paths to receive refueling service in their round trips. If a refueling station is located at site 1 , then vertex sequence for the original path and opposite sequence  of 푞 v 1 , v 3 in Figure 2, a feasible vertex sequence , v 3 = v (4) ⟩ passes through a refueling station at dual-access site 1 twice, but in opposite directions. In this respect, prior to checking whether paths go over candidate sites multiple times, we treat a dual-access candidate site as two distinct single-access sites. From this observation, we derive the second property for candidate sites of feasible paths.
is an original path that satisfies Property 1 and passes through a single-access candidate site 푚 1 ∈ 푀 twice: en, there must exist an original path, denoted as 퐹 2 v , v , that satisfies Property 1 and goes over site 1 only Proof. Let (푎) and (푏) be the -th and -th candidate sites in By eliminating all consecutive candidate sites from (푎+1) to (푏) and the corresponding arcs in (푏) are placed at the same site When drivers decide to make detours from their preplanned paths for refueling, they would be reluctant to travel along unnecessary paths to reach available refueling stations. For instance, let us take an original path of 푞 v 1 , v 3 from the previous network in Figure 2. Recall that a vertex sequence ⟨v 1 e candidate site sequence for this second path is ⟨푚 1 = 푚 (0) , 푚 2 = 푚 (1) , 푚 1 = 푚 (2) ⟩, which enables successful trips from v 1 to v 3 without running out of fuel. However, drivers would not travel further along the subpath of the second path, ⟨v 4 = v (2) , v 5 = v (3) , v 4 = v (4) ⟩ because they can reach v 3 without refueling at site 2 . ey would return at v 4 a er visiting the refueling station at site 1 , i.e., they would use the first feasible path, ⟨v 1 the four properties are applicable to both feasible original and return paths. e first property represents that feasible paths have a proper sequence of candidate sites to be selected for AF refueling stations which cover trips in feasible paths. Property 1. Let 퐹 v , v be the -th shortest feasible original path from v to v that contains feasible candidate sites for AF refueling stations to be able to cover trips in 푃 v , v , and one of the following two cases must be satisfied: there is a single candidate site available for a refueling station on 퐹 v , v , and by the half-full tank assumption, this site is located within a distance of 푅/2 from v and v to ensure that vehicles on 퐹 v , v make a successful trip from v to v . en, by definition of 푆푆 퐹 v , v and 퐸푆 퐹 v , v , this site should be located in the common In case (b), there are at least two candidate sites available In practice, vehicles do not need to visit the same refueling station located at a single-access site several times. For example, let us consider two return path alternatives of 푞 v 14 , v 16 in Figure 2. While a vertex sequence ⟨v 16 (5) , v 14 = v (6) ⟩ goes over site 6 twice, vehicles can also use a vertex sequence ⟨v 16 16 to v 14 with a single visit to site 6 . at is, a shorter feasible path is always preferred to a longer feasible path that goes through a single-access candidate multiple times. On the other hand, it is unavoidable for vehicles to revisit the same refueling station when the station is placed at a dual-access candidate site. From one of the original paths then we say that 퐹 2 v , v is dominated by 퐹 1 v , v . is implies that 퐹 2 v , v is unnecessary to be considered in our proposed model. In practice, drivers would use only For example, let us take the previous example of 푞 v 6 , v 11 in Figure 2. Recall that en, we can find another feasible original path, which is 퐹 2 v 6 , v 11 = ⟨v 6 ey have the same sequence of candidate sites, but vehicles may prefer Furthermore, in order to clarify the concept of dominated feasible paths, let us see another example of 푞 v 14 , v 16 in Figure 2. For the return path, we have two feasible return paths, we have five possible sets of candidate sites to cover the trip, i.e., 6 , 푚 6 , 푚 5 , 푚 7 , 푚 6 , 푚 7 , 푚 5 , and 푚 7 , 푚 6 , 푚 5 . In case of 퐹 2 v 16 , v 14 , 푚 7 , 푚 5 is the only feasible set. at is, In this respect, Property 4 shows that the shortest feasible paths can dominate other feasible paths when the candidate sites of the shortest feasible paths include the candidate sites of the other paths.
Proof. By definition of dominated paths, we need to show that any subset of candidate sites to cover the trips in 퐹 2 v , v is also able to cover the trips in is a set of candidate sites that covers trips in then we can construct a shorter feasible path whose instead of the second one for the original path of 푞 v 1 , v 3 . Similarly, drivers can consider two possible paths for the return path of 푞 v 12 , v 17 . While vehicles using a vertex sequence ⟨v 17 (5) , v 12 = v (6) ⟩ require each refueling stop at both sites 5 and 6 , vehicles in a vertex sequence ⟨v 17 (2) , v 12 = v (3) ⟩ need only one refueling stop at site 6 . Since the length of the latter path is shorter, vehicles in the return path of 푞 v 12 , v 17 would use the second vertex sequence instead of the first one. Based on these observations, we conclude that drivers would not make detours along a vertex sequence which requires unnecessary travel to reach available refueling stations. From now on, we call a site in unnecessary subpaths an irrelevant candidate site, so that feasible paths can be represented by vertex sequences without irrelevant candidate sites. Property 3 provides conditions under which there exist irrelevant candidate sites in paths.
when one of the following conditions is satisfied: Proof. In case (a), suppose that is irrelevant because vehicles covered by a single refueling station at site (0) can complete trips from v to v . Similarly, it can be shown when implies that any additional refueling stop is unnecessary between (푎) and (푏) . en, by Conditions (1), (2), and (3), if any subset of candidate sites that covers the trips in 퐹 2 v , v is also able to cover the trips in 퐹 1 v , v , path of 푃 v , v , for 푘 = 1, ⋅ ⋅ ⋅ , 퐾, generated through the minheap data structure, is called a candidate path from v to v and is denoted as 퐹 ὔ v , v .
When generating 퐹 ὔ v , v , we determine whether or not 퐹 ὔ v , v is included in the set of SNDF paths for 푞 v , v by Properties 1, 2, 3, and 4. First, 퐹 ὔ v , v needs to have AF refueling availability by Property 1. If either then 퐹 ὔ v , v has AF refueling availability. Second, if 퐹 ὔ v , v repeats some single-access candidate sites, then the path is excluded from the set of SNDF paths because we also construct a similar feasible path without repeating these candidate sites whose length is shorter than 푑 퐹 ὔ v , v by Property 2. Next, regarding Properties 3 and 4, we compare is removed from the set of SNDF paths for 푞 v , v because the associated OD pair already has shortest feasible paths which dominate 퐹 ὔ v , v by Property 4. Until reaching the maximum number, , of SNDF paths, Eppstein's algorithm updates a binary heap that includes sets of sidetracked arcs to find the next paths. However, the -th SNDF path algorithm can even stop before reaching if the sum of sidetracked arcs' δ v v , v on the root of the binary heap becomes greater than the value of 훽 ⋅ 푑 푃 v , v . is represents that the length of a newly constructed path cannot exceed a maximum deviation distance of drivers. A er determining at most SNDF paths from interchanges v 's to a certain interchange v , we similarly apply the above procedures to other interchanges v 's. Regarding 2퐾 SNDF paths in 푞 v , v , Algorithm 1 describes the pseudocode of the -th SNDF path algorithm for all OD pairs.
We prove below that at most 2퐾 SNDF paths (i.e., up to original paths and another up to SNDF return paths) for each OD pair can be determined by the -th SNDF path algorithm.
then it also contradicts that 퐹 1 v , v is the shortest feasible path. us, 퐹푆 퐹 2 v , v can cover trips in 퐹 1 v , v by Property 1. ☐

-th Shortest Nondominated Feasible Path Algorithm.
In this subsection, we describe the -th shortest nondominated feasible (SNDF) path algorithm to generate multiple paths that (1) offer AF refueling availability, (2) contain unique singleaccess candidate sites, (3) do not have irrelevant candidate sites, and (4) are not dominated by other paths for all OD pairs. We extend the -th shortest path algorithm developed by Eppstein [12] to first construct candidate paths between ODs for SNDF paths, and then eliminate some of the paths that violate the four properties discussed on the previous subsection. Finally, given a deviation factor and a maximum number of SNDF paths, the -th SNDF path algorithm determines up to SNDF paths from v to v and another up to SNDF paths from v to v . e first step of the -th SNDF path algorithm is to transform a simple directed network 퐺(푉, 퐸) with the set of vertices and the set of arcs into an expanded directed network, denoted as 퐺 ὔ 푉 ὔ , 퐸 ὔ , by transforming each dual-access candidate site in into two distinct single-access candidate sites, one in each side of the road, and setting all candidate sites in as vertices in ὔ . en, the maximum cardinalities of ὔ and ὔ can be written as 儨 儨 儨 儨 푉 ὔ儨 儨 儨 儨 = 푛 + 푛 + 2푛 and 儨 儨 儨 儨 퐸 ὔ儨 儨 儨 儨 = 푛 + 2푛 + 4푛 , respectively, where and refer to the number of single-and dual-access candidate sites. Next, for interchange v ∈ 푉 ὔ , v denotes a shortest path arborescence rooted at vertex v containing the shortest path to each v in 푉 ὔ \ v . e shortest path arborescence v for interchange v in 퐺 ὔ 푉 ὔ , 퐸 ὔ can be constructed using Dijkstra's algorithm in time 푂 儨 儨 儨 儨 퐸 ὔ儨 儨 儨 儨 + 儨 儨 儨 儨 푉 ὔ儨 儨 儨 儨 log 儨 儨 儨 儨 푉 ὔ儨 儨 儨 儨 . Let v denote the set of arcs composing v . If we select arc v , v ∈ 퐸 ὔ \ 퐸 푇 v and follow the shortest path from v to v instead of taking the shortest path in v , then the additional distance from v , denoted as δ v v , v , is computed as: ese arcs aside from 푃 v , v do not belong to v , and Eppstein's algorithm calls these arcs the sidetracks of 푃 v , v . e sidetracks are used to generate the k-shortest path of 푃 v , v . For 푃 v , v , a special tree-based data structure, called the min-heap data structure, in which it is a completely binary tree with data structure in a nondecreasing order of values of sidetracks, is applied to build shortest paths of 푃 v , v using the depth first search algorithm in time log . e k-th shortest for 푘 = 1, ⋅ ⋅ ⋅ , 퐾 − 1, then 퐹 ὔ v , v is not included in the set of SNDF paths for 푞 v , v . Otherwise, the algorithm assigns 퐹 ὔ v , v to the -th SNDF path 퐹 v , v and increases the value of by one. If 푘 ≤ 퐾, a next temporary path can be constructed by adding new sequences of sidetracked arcs into v . e additional distance of the next path can be calculated as is the set of newly added arcs and If the additional travel distance does not exceed 훽 ⋅ 푑 푃 v , v , then the above processes are repeated until 푘 > 퐾. Otherwise, the algorithm moves to the next interchange v in the shortest path arborescence v . A er determining at most SNDF paths, we repeat the algorithm by building a new shortest path arborescence for other interchanges. erefore, for each OD pair, this algorithm determines up to SNDF paths from origin to destination and up to other SNDF paths from destination to origin, which leads to at most 2퐾 SNDF paths. ☐ Since the original network has single-and dual-access candidate sites for AF refueling stations, vehicles can select different return paths from their SNDF original paths when they come back to their origin interchanges. is means that nonsymmetric round trips can be considered for each OD pair in our proposed model. In practice, drivers would make different detours from their SNDF original paths if the original paths do not provide AF refueling service in the opposite direction.

Repeat
Set 푘 = 0. Construct a min binary heap data structure, denoted as 퐻 (v ,v ) , with sequences of sidetracked arcs using the corresponding values of Note that the binary heap initially has the empty sequence of sidetracked arcs because the empty sequence of sidetrack arcs represents 푃 v , v .

Repeat
Construct a candidate path 퐹 ὔ v , v by inserting a sequence of sidetracked arcs on the root of If 퐹 ὔ v , v is satisfied with Properties 1, 2, 3, and 4, then define respectively, and set 푘 = 푘 + 1. Remove the current sequence of sidetracked arcs on the root of 퐻 (v ,v ) , and insert new sequences of sidetracked arcs from the min-heap data structure to the binary heap 퐻 (v ,v ) by maintaining the min-heap property. (iii) Type 2 -th SNDF original path:

Coverage for Type 1 -th SNDF Path.
In case of Type 1 -th SNDF original paths, since the distance between two candidate sites in 푆푆 퐹 (1) is at most , vehicles can travel from v to v by visiting two refueling stations at the sites. Sometimes, if there exists a single candidate site in the common area of 푆푆 퐹 (1) then locating a refueling station at this common area is sufficient to cover trips on Type 1 -th SNDF original paths.
us, for proper AF refueling availability on Type 1 -th SNDF original paths, AF refueling stations should be located at either one common or two distinct candidate sites in 푆푆 퐹 (1)

Figure 3 depicts trips in
Similar covering conditions are also applied to Type 1 -th SNDF return paths.

Coverage for Type 2 -th SNDF Path. When two candidate sites in each
to locate AF refueling stations on Type 2 SNDF original paths, the distance between the two sites can be longer than because the length of Type 2 -th SNDF original paths is greater than . at is, if vehicles at a station in 푆푆 퐹 (2) 푘 v 푖 , v 푗 cannot reach the other stations in 퐸푆 퐹 (2) 푘 v 푖 , v 푗 due to the limited driving range, then additional refueling stations at candidate sites in 퐼푆 퐹 (2) 푘 v 푖 , v 푗 should be located to ensure that the distance between consecutive stations is less than or equal to . In order to consider this covering condition for Type 2 -th SNDF original paths, we define identification coefficients that indicate whether a vehicle that is refueled at a station located in a candidate site can reach another station located at a subsequent candidate site to receive refueling services. For candidate sites are determined as follows: = 1 indicates that site 푚 ( ὔ ) follows (푙) along the -th path and their distance is less than or equal to . In an example in Figure 4, we have for Properties 3 and 4. ese processes are repeated until either feasible paths are generated or the length of a newly constructed path does not Next, we move to other interchanges v 's, ∀v ∈ 푉(퐺), 푖 ̸ = 푗, which requires time 푂 |푉| 퐾|푀| 2 + 儨 儨 儨 儨 퐸 耠儨 儨 儨 儨 퐾 log 퐾 . In summary, generating SNDF paths from other interchanges to a certain interchange v requires time Finally, since the number of interchanges for ODs is |푉|, the total computational complexity of the algorithm takes time original and return paths in all OD pairs. ☐

Model Formulation
Due to the limited driving range and the amount of fuel remaining at ODs, AF refueling stations should be properly placed at candidate sites along paths to ensure vehicles safely travel from one point to another without running out of fuel.
In other words, in order to cover trips on 퐹 v , v , when 푑 퐹 v , v > 푅, at least one AF refueling station should be located at each candidate site in 푆푆 퐹 v , v and 퐸푆 퐹 v , v , and then locating additional stations at candidate sites in 퐼푆 퐹 v , v are required if the distance between the two candidate sites in 푆푆 퐹 v , v and 퐸푆 퐹 v , v is greater than . On the other hand, when 푑 퐹 v , v ≤ 푅, since the distance between the two candidate sites in 푆푆 퐹 v , v and 퐸푆 퐹 v , v is less than or equal to , additional AF refueling stations are unnecessary to be considered for AF refueling stations between the sites in 푆푆 퐹 v , v and 퐸푆 퐹 v , v . Similar covering conditions are applied to the return path. From this observation, the following types of -th SNDF path are defined according to (i) Type 1 -th SNDF original path: (ii) Type 1 -th SNDF return path: subject to In this formulation, the objective function (10) maximizes the traffic flow that is covered by the set of AF refueling stations and 푑 푖 푗푘 푚 (2) , 푚 (3) ≤ 푅, then the corresponding identification coefficients for site (3) are set as 훼 푘, ese values mean that a refueling station must be located at either site (1) or (2) , so as to reach site (3) in 퐹 (2) 푘 v 푖 , v 푗 . Note that Hwang et al. [6] apply a similar procedure that uses identification coefficients to setup coverage restrictions to a unique path for each OD pair. For Type 2 -th SNDF return paths, 푈푆 퐹 (2) 푘 v 푗 , v 푖 is similarly defined to determine identification coefficients.

Mixed-Integer Programming Model.
Given a predetermined maximum number of feasible paths and predetermined locations of candidate sites for AF refueling stations on a directed transportation network, a mixed-integer programming model is formulated to select the optimal station locations that maximize the traffic flow covered (in round trips per time unit) along multiple feasible original and return paths between ODs. We first introduce two types of parameters and three types of decision variables as follows:

Parameters.
푓 v , v : average traffic flow in round trips per time unit between interchanges v and v , 푘,푚 ,푚 ὔ : identification coefficient for sites and ὔ for Type 2 -th SNDF path.

Decision Variables.
en, we propose the following a mixed-integer programming model: by the current set of refueling stations, 0, otherwise, covered by the current set of refueling stations, 0, otherwise.  (17) and (18), respectively. ese constraint sets ensure vehicles on a directed transportation network to make nonsymmetric round trips if necessary for refueling service. Constraint (19) allows exactly AF refueling stations to be located at the predetermined candidate site locations. Lastly, the three kinds of binary decision variables are defined in Constraint sets (20) to (23). e main stakeholders of the proposed model that maximizes the covered traffic flows can be shareholders, investors, employers and employees, and any other private entities whose stake is directly or indirectly tied to this objective function, assuming that the revenue is proportional to the traffic flow located at the selected candidate sites on the transportation network. Constraint sets (11) and (12) guarantee that trips along SNDF original paths of Types 1 and 2 receive refueling service from one station in set of candidate sites 푆푆 퐹 v , v and another one in set of candidate sites 퐸푆 퐹 v , v . e -th SNDF return paths have similar refueling conditions in Constraint sets (13) and (14). Next, coefficients in Constraint set (15) are determined by identification coefficients 푘,푚 ,푚 ὔ for 푚 푡 ∈ 퐶푆 퐹 (2) otherwise, 훼 푘,푚 ,푚 ὔ = 0. is constraint set guarantees that, if a refueling station is located at ὔ in 푈푆 퐹 (2) 푘 v 푖 , v 푗 , then one of sites in 퐶푆 퐹 (2) 푘 v 푖 , v 푗 should be selected to locate another station which would provide refueling for vehicles to reach the station at ὔ. Similarly, Constraint set (16) is applied to   time because it is directly related to the problem size. In the -th SNDF path algorithm, 2퐾 SNDF original and return paths are generated for all OD pairs, so that the number of decision variables and constraints in the mathematical model for this optimization application largely increase as increases. However, the results show that, the longer the limited driving range, the shorter the computational time.
is can be explained by the relationship between the limited driving range and Type 2 SNDF paths. at is, if the limited driving range is sufficiently large, the paths that satisfy the four properties of feasible paths are quickly determined and the number of constraints for Type 2 SNDF paths decreases, so that the computational time also decreases.

Effect of the Number of SNDF Paths.
We first analyze the effect of the different number of SNDF paths on the coverage of OD traffic flows for each deviation factor. e three graphs in Figure 6 show the coverage of traffic flows (in %) with 훽 = 0%, 50%, and 100% of 푑 푃 v , v as we increase the number of refueling stations from 푝 = 1 to 25. Note that we fix 푅 = 20 in this subsection.
First, when 훽 = 0%, the coverage is not significantly affected by the values of , but a slight improvement is obtained between 퐾 = 1 and 퐾 = 3, 5 when the number of refueling stations increases from 푝 = 3 to 15. For example, when 푝 = 8 and 퐾 = 1, the first SNDF return path of is not covered. In the case of 퐾 = 3 and 5 with 훽 = 0%, the mixed-integer programming model locates one refueling station at site 14 and then 퐹 3 v 13 , v 3 = ⟨v 13 (2) , v 3 = v (3) ⟩, which is the second SNDF return path of 푞 v 3 , v 13 , is covered by the one refueling station. is indicates that some OD pairs have multiple shortest paths which satisfy the properties of feasible paths, so that all these paths can be selected in our model. Note that 퐾 = 3 is enough for multiple SNDF paths to be considered in the network when 훽 = 0% because there is no difference in coverage between 퐾 = 3 and 퐾 = 5 from 푝 = 1 to 25. Second, as we allow drivers to take paths with long detour distances, the effect of the value of on the coverage of traffic flows by a given number of refueling stations is significantly noticeable. In particular, the coverage improvement between 퐾 = 1 and 퐾 = 5 ranges from 0% to 8.31% when 훽 = 50%, while the improvement increases up to 13.07% when 훽 = 100%. e -th SNDF path algorithm can generate more SNDF paths under 훽 = 100%, so that more paths can be considered to locate the optimal sets of refueling stations in the mixed-integer programming model. en, more traffic flows using 2퐾 SNDF original and return paths can be covered by a given number of stations. Note that the coverage improvement between different number of SNDF paths fades away slowly as the refueling network becomes mature. Lastly, in all cases, OD pairs are not fully covered when we locate refueling stations at all candidate sites. For example, considering 퐾 = 5 with 훽 = 100% for both original and return paths, which is the maximum deviation factor in the experiments, 70.58% of all traffic flows can be refueled with 푝 = 25. is is because some refueling stations at single-access candidate sites can provide refueling service only to traffic covered by a given number of stations; that is, the more traffic flow covered by the stations, the higher revenue the stakeholders expect. If we assume that the construction cost does not significantly differ by station type or region, then Constraint (19) can play a role as the budget constraint. is assumption, however, would not be always true in practice [60]. So, we can relax this assumption by replacing Constraint (19) by the following constraint: where refers to the capital cost for building a refueling station in candidate site 푚 ∈ 푀, and the refueling infrastructure budget. Constraint (24) allows different capital costs by station type (e.g., single-access vs. dual-access) or by region and forces the model to determine candidate sites for refueling stations within the budget. is constraint is used to build AF refueling stations in a statewide network to showcase that the proposed model is well-suited for solving practical problems in Section 6.

Computational Experiments
e proposed mixed-integer programming model with the -th SNDF path algorithm is applied to a well-known directed transportation network with 25 vertices, 43 arcs, and 300 OD pairs provided by Simchi-Levi and Berman [61]. Figure 5 shows the arc lengths next to each arc. It also displays 25 randomly selected candidate station locations, including 10 single-and 15 dual-access sites, to allocate the AF refueling stations in the test network. A singleaccess candidate site is symbolized by a triangular shape on one side of the arc, and a dual-access candidate site is depicted by a diamond shape in the middle of the arc. To analyze the effects of number of SNDF paths, maximum deviation distance, and vehicle driving range on the OD flow coverage, we consider three route options, 퐾 = 1, 3, and 5, three deviation factors, 훽 = 0%, 50%, and 100%, and three driving ranges, 푅 = 15, 20, and 30.
All computational experiments in Section 5 were conducted on an Intel i5 2.2 GHZ Dual-Core laptop with 12 GB RAM. Our computational experiments have three procedures. In the first procedure that is denoted as P1, to determine 2퐾 SNDF original and return paths for all OD pairs, we ran the -th SNDF path algorithm with the three limited driving ranges (i.e., 푅 = 15, 20, and 30) and the three deviation factors (i.e., 훽 = 0%, 50%, and 100%) using MATLAB R2013a. Note that the same deviation factors were applied to both original and return paths. en, in the second procedure, denoted as P2, the mixed-integer programming model with the SNDF paths and set of candidate sites was programmed in MATLAB to generate CPLEX format files that follow the syntax rules of CPLEX. In the last procedure, denoted as P3, all the problems were solved by version 12.4 of CPLEX. Table 4 provides the computational times for the three solution procedures above, P1, P2, and P3. In general, the overall computational time increases as the number of SNDF paths increases and the deviation factor is larger. Specifically, the large number of SNDF paths affects the computational deviation factor values of the vehicles that deviate from their shortest paths, denoted as , when the number of refueling stations ranges from 푝 = 1 to 25. Note that we fix 푅 = 20 and consider 퐾 = 1, 3, and 5 in the table. For example, given 푝 = 10, when vehicles are allowed to deviate to alternative paths whose lengths are at most 50% longer than those of the corresponding shortest paths and 퐾 = 1, the covered detour flows account for only 12.47% of the total covered flows and their average deviation factor is 19.99% of their shortest path lengths. First, regardless of the value of , we observe that the value of generally increases from 푝 = 1 to 25 when we consider only one SNDF path between ODs. is represents that the deviation factor largely affects the coverage if vehicles have few options to select their route regardless of the number of refueling stations. However, with 퐾 = 3 and 5, the value of flows in the same driving direction. Table 5 summarizes the coverage of traffic flows with 훽 = 0%, 50%, and 100% for 푝 = 1, ⋅ ⋅ ⋅ , 25 when we allow a different number of SNDF paths for each OD pair.

Effect of the Maximum Deviation Distance.
is subsection analyzes how much traffic flow makes detours from their preplanned paths for refueling and how long is the average deviation distance of covered detour traffic flows. For each OD pair, vehicles are able to deviate from shortest original path 푃 v , v and shortest return path 푃 v , v up to 푑 푃 v , v 1 + 훽 and 푑 푃 v , v 1 + 훽 , respectively, in order to reach AF refueling stations. Table 6 shows the proportions of the covered detour traffic flows with respect to the total flow covered, denoted as , and the average for 훽 = 50%, 푝 = 22 for 훽 = 100%), and then the coverage difference is steady. is implies that, if an additional refueling station is added to a network, the coverge of vehicles with short driving ranges can increase significantly when a long deviation distance is allowable. Next, when 푅 = 30, the coverage difference between 훽 = 0% and 훽 = 50% is insignificant. However, if vehicles are able to deviate further from their preplanned paths, i.e., 훽 = 100%, the trends of improvement fluctuate according to the number of refueling stations, especially for 푝 = 1-8. is shows that a significant coverage improvement can be expected by adding one refueling station at a time on a network where the infrastructure of AF refueling stations is insufficient if vehicles with a long driving range are allowed to make adequate detours. Table 7 summarizes comparisons of optimal coverages of no-detour versus detour flows and their coverage differences for 푅 = 15, 30 when 퐾 = 3.
tends to increase when refueling stations are scarce, while the value of starts decreasing a er the traffic flows are mostly covered. ese patterns imply that the effect of on the traffic flow covered gradually diminishes as the number of refueling stations increases, especially when vehicles have multiple path options between OD pairs. Next, regarding , we have an extreme case where the value of drops rapidly when we consider 훽 = 100% and the number of refueling stations is small. However, in most cases, the value of does not change significantly as the coverage of traffic flows increases. is result can be explained by the predetermined locations of candidate sites. In other words, even though we allow trips along many feasible paths with long detour distances between OD pairs, the coverage proportion of detour flows is usually small with respect to all traffic flows because the candidate sites are predetermined before establishing AF refueling stations.    [5]. e OD flows between counties are estimated by a simple gravity spatial interaction model [62]. We use the Freight Analysis Framework Version 4 long distance truck volume for the year of 2045 in the state of Pennsylvania, which was developed by the Federal Highway Administration

Case Study: The State of Pennsylvania
In this section, we apply the proposed model with different capital costs to build LNG refueling stations for long haul logistics on a highway system by their station type (singleaccess vs. dual-access) into the state of Pennsylvania with 2,211 OD pairs, so as to demonstrate the performance in a large-size problem. From the Topologically Integrated Geographic Encoding and Referencing (TIGER) database, we first make a map of Pennsylvania with center points of 67 counties and centroid-to-centroid routes, as shown in Figure 8. en, we stations, construction of additional single-access stations can be beneficial to trucks making detours between ODs.
Next, we observe that we can spend a large amount of budget effectively when trucks have more route options between ODs. If we have more than $120 M and trucks with 푅 = 300 have three SNDF paths (i.e., 퐾 = 3), more dual-access LNG refueling stations are selected instead of single-access LNG refueling stations comparing to the case of trucks with a single SNDF path (i.e., 퐾 = 1). It means that since trucks can have many alternative paths to make round trips between ODs, the proposed model with the budget constraint effectively selects more dual-access LNG refueling stations on their routes and improves the optimal coverage for a given budget. e case of trucks with 푅 = 600 also shows a similar situation.
Next, as shown in Figure 9, the marginal optimal coverage resulting from additional refueling stations generally decreases with the amount of construction budget available. In other words, initially, a large percentage of the optimal coverage can be achieved with a relatively small construction budget. A er that, the rate of increase of optimal coverage declines as the construction budget increases. One of the interesting things here is that the optimal coverage gap between no-detour and detour flows, however, increases as the construction budget increases. For example, for 푅 = 300, when the budget is initially $10 M, the optimal coverage gap between no-detour (훽 = 0%) and detour flows (훽 = 50%) with 퐾 = 3 is 7.91%; in 2016 [63]. Also, to consider a conservative range for current and older models of LNG trucks with a single fuel tank [64] and improving technologies of dual fuel tank systems [65], we set 푅 = 300 and 600 for LNG trucks. Lastly, we conduct this experiment by using Intel(R) i7-6700 K CPU 4.00 GHz and 16 GB RAM desktop PC. Table 8 provides comparisons of number of selected stations and optimal coverage for no-detour versus 50-percent detour flows with 퐾 = 1 and 3 when 푅 = 300 and 600. Also, Figure 9 shows trade-off between budget and optimal coverage for the two LNG truck driving ranges (푅 = 300, 600) and three route options (훽 = 0%, 훽 = 50% with 퐾 = 1, and 훽 = 50% with 퐾 = 3). We apply our coverage maximization model with the budget constraint to the state of Pennsylvania using the budget range from $10 M to $200 M in $10 M increments.
First, we easily identify the effect of detour flows (both and ) on the optimal coverage for different budgets. For example, the smallest budget of $10 M allows only construction of two dual-access LNG refueling stations with the optimal coverages of 19.99% and 21.14% for 푅 = 300 and 600, respectively, when trucks merely use shortest paths between each OD. When we allow trucks to make detours between ODs, the corresponding coverages increase to 27.91% and 29.40% for the both limited driving ranges. is implies, when we do not have enough budget to construct dual-access LNG refueling    with 10 single-and 15 dual-access candidate sites, and (ii) a statewide network of Pennsylvania state containing 2,211 OD pairs with 50 single-and 50 dual-access candidate sites. For our experiment on the classical 25-vertex network, we have analyzed the effects of number of SNDF paths, deviation distance, and vehicle driving range on the coverage of no-detour and detour flows for a given number of AF refueling stations. Computational experiments have shown that the number of SNDF paths and maximum deviation factor largely affect the coverage of traffic flows, especially when the network does not have enough AF refueling stations to cover most of the traffic flows. Next, for our experiment on the case study of Pennsylvania state, we have also analyzed the effects of detour flows, station type (single-access vs. dual-access), and longer driving range for different construction budget levels for the LNG refueling infrastructure to present more insights into our society. For future research, our proposed model can be applied to a variety of real-world problems to locate AF refueling stations in highways. To reduce the complexity of the -th SNDF path algorithm, we can also attempt to develop a more efficient version of the -th SNDF path algorithm for large-scale networks. Upchurch et al. [66] consider available fueling capacity of AF refueling stations to cover traffic flows. Similarly, our proposed model can be extended to contemplate the limited capacity of AF refueling stations. Lastly, we generalize our solution approach to address the continuous version of the AF refueling station location problem where stations can be located anywhere along the network. As we have discussed the effect of predetermined candidate sites on the coverage in Subsection 5.2, if we were able to locate stations anywhere along the network, the suboptimality of the optimal solutions of the discrete version of the problem could be significantly improved.
Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.
when the budget later becomes $200 M, this gap increases to 11.63%. is result implies that, even if the LNG refueling infrastructure becomes mature (e.g., the budget is around $200 M or more), the effect of detour flows on the optimal coverage would be even more noticeable. e next interesting thing is that the effect of the different number of SNDF paths on the optimal coverage declines as the construction budget increases. at is, when LNG trucks are allowed to deviate from their shortest path up to 50% of 푑 푃 v , v , the optimal coverage gap between detour flows with 퐾 = 1 and detour flows with 퐾 = 3 declines as the construction budget increases. For example, for 푅 = 300, at the budget of initial $10 M, this optimal coverage gap is 5.96%; later at the budget of $200 M, the gap declines to 1.85%. It supports our result analysis discussed in Subsection 5.1 on the effect of the number of SNDF paths. Another interesting one here is that the effect of driving range of 푅 = 600 on the optimal coverage also declines as the construction budget increases, especially when considering detour options. In other words, comparing Figures 9(a) and 9(b) show that the optimal coverage gap between 푅 = 300 and 600 for detour flows with 퐾 = 3 declines overall as the budget increases. For example, at the budget of initial $10 M, this optimal coverage gap is 1.49%; later at the budget greater than or equal to $120 M, this gap declines to less than or equal to 0.15%. us, if the LNG refueling infrastructure becomes mature (e.g., if the budget is around $120 M or more) and if drivers allow reasonable detour options, using LNG trucks with driving range 푅 = 600 would not lead to a significant improvement on the LNG refueling coverage.
Lastly, the last row of Table 8 presents the computational times for P1, P2, and P3 when we run the case study of the state of Pennsylvania. As the computational times of the small instances in Section 5 are shown, the large number of SNDF paths increases the overall computational time. Also, we verify the relationship between the limited driving range and the computational time; that is, the longer the limited driving range, the shorter the computational time because, when the limited driving range is long, the paths that meet the four properties of feasible paths are determined quickly and the number of constraints for Type 2 SNDF paths decreases.

Conclusions and Future Research
Most path-based demand models available in the literature assume that AF vehicles travel along a unique path between ODs to complete their round trips. In this research work, we have considered that, if preplanned paths for each OD pair cannot offer any AF refueling service, then alternative paths can be used by AF vehicles to receive refueling service. With the four feasible path properties derived in Subsection 3.3, we have proposed an algorithm that generates SNDF paths in each path direction for all OD pairs. A new mixed-integer programming model has been also proposed to optimally locate AF refueling stations at candidate sites with the objective of maximizing the coverage of traffic flows using multiple SNDF paths on a directed transportation network. In order to test our proposed model with the -th SNDF path algorithm, we have run two experiments: (i) a classical 25-vertex network containing 300 OD pairs