---
title: Flow-Refueling Location Model
url: https://www.emergentmind.com/topics/flow-refueling-location-model
type: topic
---

# Flow-Refueling Location Model

The Flow-Refueling Location Model (FRLM) is a path-based optimization model for siting refueling or charging facilities on a transportation network for vehicles with limited range. In the formulation studied in recent literature, the Fuel Refueling Location Problem is the same optimization goal and modeling framework as the FRLM introduced by Kuby and Lim (2005): candidate stations are selected so that vehicles traveling on preselected origin–destination paths, or on sets of feasible deviation paths, can complete trips under a budget on the number or cost of facilities [2409.04554]. Classical FRLM is static and typically uncapacitated, but subsequent work has extended it to deterministic path-segment MILPs, continuous candidate-site spaces, multi-path and probability-weighted coverage, station sizing, bilevel reliability requirements, robust time-expanded charging networks, and coupled transportation–power planning for electric vehicles [2109.14358].

## 1. Canonical model structure

In its classical form, FRLM is defined on a graph \(G=(N,A)\) with candidate facility sites \(N^s\subseteq N\), a set of OD flows or path flows, path-specific weights, and a vehicle range \(R\). A station-opening variable \(x_j\in\{0,1\}\) indicates whether a facility is placed at site \(j\). The central feasibility logic is that the spacing between consecutive refueling opportunities, including the path ends, must not exceed \(R\). One canonical set-covering formulation minimizes infrastructure cost while forcing every relevant subpath segment \(B_{p,k}\) to contain at least one station:
\[
\begin{array}{ll}
\text{minimize} & \displaystyle \sum_{j \in \mathcal{N}^s} o_j\, x_j \\
\text{subject to} & \displaystyle \sum_{j \in B_{p,k}} x_j \ge 1 \qquad \forall p \in \mathcal{P},\ \forall k \in \mathcal{K}_p, \\
& x_j \in \{0,1\} \qquad \forall j \in \mathcal{N}^s .
\end{array}
\]
A budget-limited variant instead maximizes covered flow:
\[
\begin{array}{ll}
\text{maximize} & \displaystyle \sum_{p \in \mathcal{P}} w_p\, y_p \\
\text{subject to} & y_p \le \displaystyle \sum_{j \in B_{p,k}} x_j \qquad \forall p \in \mathcal{P},\ \forall k \in \mathcal{K}_p, \\
& \displaystyle \sum_{j \in \mathcal{N}^s} o_j x_j \le B, \\
& x_j \in \{0,1\}\ \forall j \in \mathcal{N}^s,\quad y_p \in \{0,1\}\ \forall p \in \mathcal{P}.
\end{array}
\]
These formulations encode the classical FRLM “coverage logic” along paths by requiring every refueling-required segment to contain a facility [2209.07305].

A deterministic path-based refinement appears in the deterministic flow refuelling location problem (DFRLP), where a cyclic OD flow \(f\) travels on a predetermined path, \(x_k\) indicates facility siting, \(y_f\) indicates whether the flow is covered, and \(i_{klf}\) selects cycle segments between refueling opportunities. Its objective is
\[
\max \sum_{f\in F} v_f y_f,
\]
subject to a station-count constraint \(\sum_{k\in K} x_k=p\), a driving-range constraint
\[
\sum_{l \in L_{k f}^+} i_{k l f} t_{k l} - (1-y_f)M \le R,
\]
and linking constraints that enforce a valid origin-to-destination-and-back refueling sequence. In this formulation, explicit path segments and Big-\(M\) activation replace precomputed feasible station combinations by a single MILP [2109.14358].

## 2. Coverage logic, detours, and route realism

Classical FRLM assumes deterministic flows on preselected paths, often one preferred path per OD pair, with full initial fuel at the origin and unlimited station capacity unless capacity is modeled explicitly. A central distinction in later work is between single-path FRLM, multi-path and deviation FRLM, and models that move from flow coverage to node-based or route-based coverage. The original objective is OD- or flow-based rather than node-based: FRLM maximizes covered flow, not mere proximity of nodes to stations [2409.04554].

Deviation and multi-path generalizations enlarge the feasible route set \(R_q\) for a demand \(q\). In the probability-based multi-path refueling location model, each OD pair may use multiple deviation paths, round-trip feasibility is modeled with fuel-state variables \(B_i^{odk}\) and refueling variables \(\ell_i^{odk}\), and the objective becomes expected node coverage,
\[
\max \sum_{o\in O} p_o z_o,
\]
where \(p_o\) is the probability that origin \(o\) becomes a demand node and \(z_o\) is the fraction of destinations reachable from \(o\). Fuel conservation and capacity along a used path are enforced through Big-\(M\) equalities on consecutive arcs, and an OD is covered if at least one enumerated path is feasible. This shifts the objective from deterministic maximal coverage to maximal expected coverage [1805.06068].

A more radical generalization removes the discrete candidate-site assumption. For fleet operations on a general network, stations may be located anywhere along edges, deviations from prescribed routes are allowed, and a route \(U_t\) is covered by a point \(x\) if a vehicle can deviate at \(v_d(x)\), refuel, return to the route before the next unvisited stop, complete at least one full traversal, and return to the refueling station without violating the range \(R\) or maximum deviation distance \(D\). For off-route points, the coverage conditions are
\[
l\!\left(v_d(x), x\right) + l\!\left(x, v^-\!\left(S_{U_t}\!\left(S_{U_t}(x)\right)\right)\right) + l'\!\left(U_t\!\left(v^-\!\left(S_{U_t}\!\left(S_{U_t}(x)\right)\right), v_d(x)\right)\right) \le R,
\]
\[
l\!\left(v_d(x), x\right) \le D.
\]
This formulation treats refueling as a detour-aware coverage problem on dedicated closed routes rather than a shortest-path OD problem [2509.11441].

A frequent misconception is that FRLM is exhausted by static on-path spacing constraints. The literature shows a broader spectrum. Some models retain path coverage but allow deviations; some allow stations anywhere on edges; some replace path coverage by node-weighted “attractiveness”; and some enforce backup or reliability conditions on stations. These are related but not identical uses of the FRLM idea [2501.19000].

## 3. Formulation theory and polyhedral strength

A major theoretical development concerns the relative strength of alternative FRLM formulations. Let \(R_q\) denote the set of feasible routes for demand \(q\), \(x_j\) the siting decision, \(z_{qr}\) a path-coverage variable, and \(y_q\) an aggregated OD-coverage variable. In the disaggregated path-based formulation,
\[
\max \sum_{q\in Q}\sum_{r\in R_q} f_q z_{qr}
\]
subject to
\[
\sum_{r\in R_q} z_{qr}\le 1,\qquad
\sum_{j\in S} x_j \ge z_{qr}\quad \forall q,\forall r\in R_q,\forall S\in C_{qr},
\]
where \(C_{qr}\) is the family of refueling-required segments on route \(r\). In the aggregated formulation,
\[
\max \sum_{q\in Q} f_q y_q
\]
subject to
\[
\sum_{j\in S} x_j \ge y_q\quad \forall q,\forall S\in C_q,
\]
where \(C_q\) is built by taking unions of one segment from each route in \(R_q\). The aggregated model compresses route-level feasibility into OD-level segment families [2409.04554].

The strength comparison is expressed through OD-wise LP value functions. For any fractional \(x\),
\[
v_q^{\text{tight}}(x)\le v_q^{\text{agg}}(x)\le v_q^{\text{disagg}}(x).
\]
If \(R_q\) contains a single path, then
\[
v_q^{\text{agg}}(x)=v_q^{\text{disagg}}(x),
\]
and if that path is a single, simple path, then
\[
v_q^{\text{tight}}(x)=v_q^{\text{agg}}(x)=v_q^{\text{disagg}}(x).
\]
The interpretation given in the literature is that the set-cover polyhedron on a single simple path has a totally unimodular structure, so the LP relaxation is integral in that case. By contrast, under a cardinality budget the aggregated relaxation strictly dominates the disaggregated relaxation, and the gap can be arbitrarily large [2409.04554].

The same line of work also identifies minimal segment sets as the relevant polyhedral objects. If \(C'\subseteq C\) is the subfamily of minimal sets, then the polyhedron defined by all inequalities \(\sum_{j\in S}x_j\ge y\) for \(S\in C\) is identical to that defined by the inequalities for \(S\in C'\), and each minimal-set constraint is facet-defining in the \(0\)–\(1\) hull. Two tightening devices are emphasized under a cardinality budget: OD-specific feasible-site filtering, which replaces \(S\) by \(S\cap N_q\), and budget–noninterference facets,
\[
\sum_{j\in T} x_j \le b-w_q+\bigl(|T|-b+w_q\bigr)(1-y_q),
\]
for suitable \(T\subseteq \bar M_q\). These results clarify why aggregated branch-and-cut formulations are strong for deviation FRLM and why direct path disaggregation can yield weak LP bounds when multiple routes per OD are allowed [2409.04554].

## 4. Principal extensions of the FRLM family

The FRLM family now includes several distinct extensions that alter either the candidate-site space, the objective, the feasible-route definition, or the coupling with operational constraints.

| Variant | Main modeling addition | Source |
|---|---|---|
| DFRLP, MCFRLP, LCFRLP, CFRLP, CMCFRLP | coverage quota, location-dependent construction costs, capacity-limited stations, simultaneous sizing and partial coverage | [2109.14358] |
| Continuous-site FRLM with FDS | station can be anywhere on edges; finite dominating set \(S^*\) of endpoints \(EP\) contains an optimal solution | [2509.11441] |
| Bilevel reinforced coverage | every opened station must have another station within \(R\); lower level maximizes node attractiveness under budget | [2501.19000] |
| Probability-based multi-path refueling | node-specific probabilities \(p_o\), round-trip feasibility, and expected node coverage objective | [1805.06068] |
| Modified CFRLM with SOCP service ability | heterogeneous PEV types, time-varying demand, AC power flow, and second-order-cone service constraints | [1702.01897] |
| Diffusion-integrated FRLM | annual budgets from continuous-time adoption and supply curves, with Macro–Meso–Micro planning layers | [2402.14263] |

The capacity and sizing extensions are particularly important for EV charging. In the capacity-limited station model, \(n_k\) gives the number of charging poles at site \(k\), \(z_f\in[0,1]\) is the proportion of flow \(f\) covered, and the capacity restriction sums distance-weighted flow demand assigned to a station and bounds it by \(Cap\cdot n_k\). This allows simultaneous placement and sizing, and it permits partial coverage when capacity is insufficient. Location-dependent construction costs similarly replace a station-count budget by a budget on \(\sum_k c_k x_k\), making the number of stations endogenous [2109.14358].

The continuous-site extension replaces arbitrary discretization by an exact finite dominating set. For each route and edge, the Edge Scanning algorithm constructs refueling segments \(RS(U_t;V_a,V_b)\), extracts their endpoints, and proves that for any point \(x\in G\) there exists an endpoint \(w\in EP\) with \(T(x)\subseteq T(w)\). At least one optimal solution therefore exists in the finite endpoint set \(EP\), and the final model is a set-cover problem
\[
\min \sum_{w\in C} x_w
\quad\text{subject to}\quad
\sum_{w\in C} a_{w,t}x_w\ge 1\ \ \forall U_t\in H,
\]
with \(C\subseteq EP\) obtained by grouping identical coverage sets [2509.11441].

The bilevel reinforced-coverage model departs further from classical FRLM. Its upper level minimizes the number of opened sites under conditional covering constraints
\[
\sum_{l\in V,\ l\neq k} b_{kl}y_l \ge 1 \qquad \forall k\in V,
\]
so that every node is within range \(R\) of some station and every opened station has a distinct backup station within \(R\). The lower level then maximizes node attractiveness \(\sum_k \omega_k x_k\) under a budget, with \(x_k\) representing the number of charging devices installed at \(k\). This is a reliability-conscious covering model rather than a path-spacing FRLM in the strictest sense [2501.19000].

## 5. Exact algorithms, decomposition, and heuristics

Because FRLM variants combine siting, path feasibility, range restrictions, and often capacity or time dependence, solution methodology is a defining part of the literature. In the continuous-site setting, the Edge Scanning algorithm constructs the endpoint set \(EP\) in polynomial time. Its complexity is \(O(he(n^2+h))\), where \(h=|H|\), \(e=|E|\), and \(n=|V|\). After preprocessing, the remaining decision problem is a compact set-cover MIP on the reduced candidate set \(C\subseteq EP\) [2509.11441].

In robust charging-network planning for metropolitan taxi fleets, the strategic problem is handled by a cutting-plane master problem
\[
\min \sum_{s\in \mathcal{N}^s} o_s x_s
\]
with feasibility cuts
\[
\sum_{s\in \mathcal{I}_j} x_s \ge 1,
\]
while the operational feasibility problem is solved by branch-and-price. The column-generation master problem selects one charging-augmented route per vehicle, station capacities are enforced by route–time occupancy coefficients \(k_{vpst}\), and pricing is a resource-constrained shortest path problem on a time-expanded network with label updates
\[
C_j = C_i + \pi_i,\qquad
R_j = -r_{ij} + 
\begin{cases}
R_i, & z_i=0,\\
\min(q^{\max},R_i+\gamma(R_i)), & z_i=1.
\end{cases}
\]
The same work embeds this deterministic core in a robust framework based on adversarial sampling, with Full Scenario, \(\alpha\)-Scenario, and \(\alpha\)-Vehicle approaches [2209.07305].

For probability-based multi-path models, the mixed-integer formulation is NP-hard, and a genetic algorithm has been developed in which a binary siting vector \(x\) is the chromosome and fitness is evaluated by solving the remaining MILP with fixed \(x\). On the Sioux Falls network, the heuristic produced expected coverage values within \(0\)–\(1.9\%\) of CPLEX optima for budgets \(m=1,\dots,12\), with CPU time reductions of up to \(\sim 90\times\) at larger \(m\) [1805.06068].

When FRLM is coupled with service-level and power-flow constraints, exact solution methods move beyond MILP into conic optimization. A modified capacitated FRLM based on sub-paths is embedded in a stochastic mixed-integer second-order cone program, where station service ability is enforced by
\[
y_i^{\text{cs}} \ge \sum_{q\in\mathcal{Q}_i}\sum_{k\in\mathcal{K}} T_k\lambda_{q,k}\gamma_{q,i,k}
+\Phi^{-1}(\alpha)\sqrt{\sum_{q\in\mathcal{Q}_i}\sum_{k\in\mathcal{K}} T_k\lambda_{q,k}\gamma_{q,i,k}^2},
\]
and radial AC power flow is imposed through SOCP constraints. The resulting MISOCP is solved by branch-and-cut in a commercial conic solver, with additional path-sharing equalities used as computational strengthening [1702.01897].

## 6. Electric-vehicle planning, robustness, and current research directions

The most active contemporary use of FRLM is EV charging network planning, where static path coverage is no longer sufficient by itself. In metropolitan taxi-fleet planning, FRLM-like coverage logic is embedded in a time-expanded, state-of-charge-based operational model with station capacity, non-linear charging functions, partial recharges, and planner-selectable risk tolerance. On Munich-based instances with \(1{,}000\) vehicles and \(60\) potential charging station locations, the Full Scenario Approach yielded cost \(o=1969\), \(Z_{\text{feas}}=100\%\), and \(V_{\text{feas}}=100\%\), while allowing \(1\%\) vehicle infeasibility reduced cost by about \(20\%\) through the \(99\)-VA solutions. The same study reports that increasing battery capacities yields up to \(10\) percentage points more vehicle feasibility than equivalent increases in charging speed, and that allowing depot charging dominates both options [2209.07305].

Coupled transportation–power formulations show that once capacity and electrical feasibility are modeled jointly, the FRLM problem becomes a siting-and-sizing problem on two interdependent networks. In a 25-node highway network refined to 93 transportation nodes and coupled to a 14-bus radial distribution system, the baseline case with \(20{,}000\) PEV/day required \(28\) stations and \(1{,}169\) charging spots, whereas \(40{,}000\) PEV/day required \(45\) stations and \(2{,}340\) spots. Cases that ignored power constraints or used DC power flow produced higher post hoc unmet demand or distorted investment decisions, indicating that transport-side coverage alone is insufficient when grid congestion is material [1702.01897].

A different line of work embeds FRLM in a multi-year demand–supply framework based on innovation diffusion and fluid queues. There, annual station budgets are derived from a cumulative supply curve \(D(t)\), demand follows Bass-type adoption dynamics, and charging-location planning is performed on Macro–Meso–Micro grids. In the Chicago sketch network, capacity constraints materially altered facility counts: at the macro level with \(150\) km range, the capacity-limited version selected \(23\) sites rather than about \(10\) without capacity; at the meso level it selected \(45\) sites; and at the micro level \(91\) sites. For \(250\) km range, the corresponding capacity-limited counts were \(16\), \(6\), and \(5\). This suggests that station capacity, EV range, and spatial resolution interact strongly in medium-term planning [2402.14263].

Across these strands, several limits of classical FRLM remain visible. The original model is static, path-based, and usually uncapacitated; it abstracts from queueing, charging duration, time windows, and endogenous route reassignment. Later work therefore adds station capacities, partial coverage, service-level approximations, robust scenario analysis, multi-path routing, continuous candidate sets, or bilevel public–private structure. A plausible implication is that FRLM is best understood not as a single fixed formulation but as a modeling core—range-constrained path coverage—around which distinct operational, economic, and infrastructure couplings are built in contemporary research.

Source: https://www.emergentmind.com/topics/flow-refueling-location-model