Owned by: China Association for Science and Technology
Sponsored by: China Coal Society
Published by: Springer Nature
About issue

The International Journal of Coal Science & Technology is a peer-reviewed open access journal. It focuses on key topics of coal scientific research and mining development, serving as a forum for scientists to present research findings and discuss challenging issues.


Coverage includes original research articles, new developments, case studies and critical reviews in all aspects of scientific and engineering research on coal, coal utilizations and coal mining. Among the broad topics receiving attention are coal geology, geochemistry, geophysics, mineralogy, and petrology; coal mining theory, technology and engineering; coal processing, utilization and conversion; coal mining environment and reclamation and related aspects.


The International Journal of Coal Science & Technology is published with China Coal Society, who also cover the publication costs so authors do not need to pay an article-processing charge.


The journal operates a single-blind peer-review system, where the reviewers are aware of the names and affiliations of the authors, but the reviewer reports provided to authors are anonymous.


  • A forum for new research findings, case studies and discussion of important challenges in coal science and mining development

  • Offers an international perspective on coal geology, coal mining, technology and engineering, coal processing, utilization and conversion, coal mining environment and reclamation and more

  • Published with the China Coal Society

Show More
Editors-in-Chief
Suping Peng, Shimin Liu
Managing Editor
Wanjie Wang
Associate Editors
Bo Hyun Kim, Dongjie Xue, Pedram Roghanchi, Wu Xiao, Zhiqiang Wu
Publishing model
Open Access. Learn about publishing OA with us
Home > Volumes and issues > Volume 10, issue 1

Reducing congestion and emissions via roadside unit deployment under mixed traffic flow

Research Article

Open Access

Published: 07 January 2023

0 Accesses

International Journal of Coal Science & Technology Volume 10, article number 1, (2023)

Abstract

It is expected that for a long time the future road traffic will be composed of both regular vehicles (RVs) and connected autonomous vehicles (CAVs). As a vehicle-to-infrastructure technology dedicated to facilitating CAV under the mixed traffic flow, roadside units (RSUs) can also improve the quality of information received by CAVs, thereby influencing the routing behavior of CAV users. This paper explores the possibility of leveraging the RSU deployment to affect the route choices of both CAVs and RVs and the adoption rate of CAVs so as to reduce the network congestion and emissions. To this end, we first establish a logit-based stochastic user equilibrium model to capture drivers’ route choice and vehicle type choice behaviors provided the RSU deployment plan is given. Particularly, CAV users’ perception error can be reduced by higher CAV penetration and denser RSUs deployed on the road due to the improved information quality. With the established equilibrium model, the RSU deployment problem is then formulated as a mathematical program with equilibrium constraints. An active-set algorithm is presented to solve the deployment problem efficiently. Numerical results suggest that an optimal RSU deployment plan can effectively drive the system towards one with lower network delay and emissions.

1.Introduction

Traffic congestion and emissions are two major negative externalities of the transportation system. Over decades, many traffic flow management strategies such as road pricing (Yin and Lawphongpanich 2006; Xu et al. 2016; Zhong et al. 2020) and tradable credit scheme (Wang et al. 2012; Gao and Sun 2014) have been designed to mitigate or internalize these externalities to improve the urban logistics efficiency and social welfare. For example, in light of the potential inconsistency between reducing congestion and reducing emissions (Nagurney 2000), Chen and Yang (2012) developed Pareto-efficient toll and rebate scheme to simultaneously minimize the congestion and emission.

In recent years, connected autonomous vehicles (CAVs) have caught great attention from the public and transportation planners due to their capability of enhancing traffic safety (Wadud 2017; Fagnant and Kockelman 2015), reducing fuel consumption (Barth and Boriboonsomsin 2009; Wadud et al. 2016), and increasing road capacity (Tientrakool et al. 2011) and driver productivity (Noruzoliaee et al. 2018), etc. Though it is tempting to replace regular vehicles (RVs) on the road with CAVs, we can envision that the road traffic will be mixed with RVs and CAVs for a long time. By leveraging the CAV flow in the mixed flow, researchers have developed a variety of new approaches for the congestion and emission management. Generally, these strategies mainly hinge on: (a) the control of CAVs’ speed profile or trajectory (Wang 2018; Yao et al. 2020; Xiong and Jiang 2021), (b) the control of CAVs’ routing behavior (Chen et al. 2020; Zhang and Nie 2018; Li et al. 2018; Sharon et al. 2018), and (c) the deployment of CAV-related infrastructure (Chen et al. 2016, 2017; Li et al. 2020; Liu et al. 2021). Specifically, controlling CAVs’ velocity, acceleration, and deceleration can achieve goals such as fuel consumption reduction (Yao et al. 2020; Xiong and Jiang 2021) and traffic flow stabilization (Wang 2018). For example, jointly optimizing traffic signal timing plans and CAVs’ trajectory can reduce the fuel consumption and CO\(_2\) emissions by 22.36% and 18.61%, respectively, in a network fully penetrated with CAVs (Yao et al. 2020). Regarding the routing behavior, CAVs’ route choices are assumed to be controllable due to their self-driving feature. Given this, Chen et al. (2020) proposed a path-control scheme to achieve the system optimum (SO), i.e., the minimum of total traffic delay, by controlling a portion of CAVs to follow the SO routing principle. It is found that the control ratio of CAVs can be very low with the implementation of the pricing scheme. The presence of CAV dedicated infrastructure, such as CAV lanes, CAV zones, and roadside units (RSUs), may affect the travel behavior of CAVs and thus the whole traffic flow, so optimizing these CAV dedicated infrastructure deployment plan can yield a reduction of traffic congestion and emissions. For example, Chen et al. (2017) optimized the deployment plan of CAV zones in a network with both RVs and CAVs to substantially reduce the social cost.

Inspired by the spirit of Chen et al. (2017), this paper aims to optimally deploy RSUs in a way to mitigate network congestion and emissions. RSUs are devices installed by the roadside to sense vehicles within their coverage areas, collect the motion information of these vehicles and communicate the collected data with CAVs on the road and the traffic management agence. In particular, RSUs can provide CAVs with information about the road traffic condition to improve their decision making in terms of choices of routes and driving profiles (Xie and Liu 2022; RGBSI 2021). Without a doubt, different RSU deployment plans may provide information with different qualities. For example, denser deployment of RSUs can collect more data, thereby offering a more comprehensive picture of the traffic on the network and enabling more reliable and robust communication; while sparse deployment implies the opposite. The information quality is expected to impact the path choice of CAV users (Chen and Du 2017; Spana et al. 2021), given that the information of higher quality is more likely to guide drivers to the fast paths and save them travel time and energy consumption (Xie and Liu 2022). As a result, this will further affect the long-term expected cost of using CAVs, which consequently impacts their adoption rate. In the meantime, the change of CAV penetration will in return influence the information quality because CAVs can also collect data in its neighborhood and share it with peer CAVs and RSUs (Xie and Liu 2022).

In this paper, RSUs will be deployed on a network where travelers make both long-term, i.e., vehicle type choice, and short-term, i.e., route choice, decisions, according to stochastic user equilibrium (SUE) principle. The information quality of CAVs is mediated by the RSU deployment plan and endogenous CAV penetration. Specifically, the impact of information quality on CAV users’ route choice is captured through the variable dispersion parameter in the path choice SUE model (Yang 1998; Lo and Szeto 2002; Yin and Yang 2003; Xie and Liu 2022), which reflects the fact that information of higher quality can increase the probability of choosing fast paths by lowering the travel time perception error. As such, different RSU deployment plans may yield different network flow distributions corresponding to different levels of network delay and emissions. Therefore, it is of critical importance to optimally design the RSU deployment plan to stimulate the change of routing behaviors of CAVs and their adoption rate in a way to mitigate the traffic congestion and emissions. To this end, we adopt the structure of Stackelberg game for the RSU optimization. Specifically, government agencies serve as the leader to design the deployment plan of RSUs, while travelers are the followers to make spontaneous reaction in terms of their choices of routes and vehicle types. Mathematically, the deployment problem can be formulated as an optimization program with equilibrium constraints (MPEC). An active-set algorithm is proposed to solve the MPEC. To the best of our knowledge, we are one of the first to explore the possibility of leveraging the RSU deployment to affect both the routing behavior of CAVs and their market penetration in a way to mitigate traffic congestion and emissions. Since the deployment of RSUs can be well combined with existing highway infrastructures, it can serve as a flexible instrument for congestion and emission management. The model and algorithm developed in this paper can be readily used to examine potential effect of a specific RSU deployment plan or to provide practical guidance to RSU deployment projects.

The remainder of this paper is organized as follows. In Sect. 2, we establish the network equilibrium model with mixed traffic flow and the optimization program for the RSU deployment. Section 3 presents the active-set algorithm to solve the optimization problem, followed by numerical experiments in Sect. 4. Finally, Sect. 5 concludes the paper and points out potential extensions to our work.

2.Model formulation

Consider a connected network G(VA) where V and A are the sets of nodes and links, respectively. The set of origin-destination (OD) pairs is denoted by W. For any given OD pair \(w\in W\), the travel demand \(Q^w\) is assumed fixed and deterministic, and the feasible path set is denoted by \(K^w\). Travelers may choose to purchase either a CAV or RV based on their expected travel cost or utility. We define set \(I = \{{\text{CAV}}, {\text{RV}}\}\) as the vehicle type set. Therefore, the traffic flow on the network is mixed with CAVs and RVs. Let \(f_i^k\) be the flow on path \(k\in K^w\ (w\in W)\) of vehicle type \(i\in I\) and \(x_{ia}\) the flow on link \(a\in A\) of vehicle type \(i\in I\). Then we have \(x_{ia} = \sum _{w\in W}\sum _{k\in K^w}\delta _a^kf_{i}^k\), where \(\delta _a^k\) is the path-link incidence. The expected travel time on link \(a\in A\), denoted by \(t_a\), is a strictly increasing function in the link flow \(x_a = \sum _{i\in I}x_{ia}\), i.e., \(t_a = t_a(x_a)\) and \(\text{d}t_a/\text{d}x_a>0\). The expected travel time on path k is thus given by \(T^k = \sum _{a\in A}\delta ^k_at_a\).

2.1 Stochastic user equilibrium

A logit-based stochastic user equilibrium (SUE) model will be established in this section to capture the vehicle type choices and path choices made by travelers. Travelers’ path choices are informed by both CAV penetration and the deployment plan of RSUs on the network, which in return affects the expected travel costs of CAVs and RVs. We will first establish the path choice equilibrium with a given CAV penetration rate, based on which we proceed to develop the vehicle type choice model. Table 1 presents the notation used in our methodology.

Table 1 List of notations

Notation

Definition

Sets

V

Set of nodes in the network, indexed by v

A

Set of links in the network, indexed by a

W

Set of OD pairs, indexed by w

\(K^w\)

Set of paths connecting OD pair w, indexed by k

I

Set of vehicle types, indexed by i

\(\Omega _\text{f}\)

Set of feasible path flow patterns

\(\Omega _\text{q}\)

Set of feasible vehicle type choice patterns

\(\Omega\)

Set of feasible joint path flow and vehicle type choice patterns

Parameters

\(\delta _a^k\)

Link-path incidence

\(Q^w\)

Travel demand for OD pair w

\(\theta\)

Dispersion parameter for the vehicle type choice MNL model

\(\phi _1, \phi _2\)

Parameters in the emission function

\(l_a\)

Length of link a

\(\gamma _1, \gamma _2\)

Weighting parameters adjusting the importance of two objectives

\(t_a^0\)

Free flow travel time of link a

\(B_a\)

Capacity of link a

\(\psi _1, \psi _2\)

Parameters characterizing the path choice MNL dispersion parameter for CAV users

\(\theta _{\text{RV}}^w\)

Dispersion parameter for the path choice MNL model for RV users

\(b^1_i\)

Value of time for commuting via i-type vehicle

\(b^2_i\)

Parameter that scales the vehicle price to capture other overhead cost for vehicle type i

\(b^3_i\)

Distance-based variable cost for using vehicle type i

\(\tau _i\)

Purchase price for vehicle type i

\(d^w\)

Average distance of paths for OD pair w

\(L_i\)

Average total mileage for the whole life cycle for vehicle type i

\(\bar{n}_a\)

Maximum number of RSUs allowed to deploy on link a

\(\underline{n}_a\)

Minimum number of RSUs required to deploy on link a

N

Maximum number of RSUs allowable for the network

Variables

\(x_{ia}\)

Flow on link a contributed by vehicle type i

\(x_a\)

Total flow on link a

\(f_i^{wk}\)

Flow on path k for OD pair w contributed by vehicle type i

\(t_a\)

Travel time on link a

\(T^k\)

Travel time using path k

\(q^w_i\)

Number of travelers with OD pair w who choose vehicle type i

\(P^{wk}_i\)

Probability of choosing path k for travelers with OD pair w and vehicle type choice i

\(\theta ^w_{\text{CAV}}\)

Dispersion parameter for the path choice MNL model for CAV users

\(n_a\)

Number of RSUs deployed on link a

\(\rho _a\)

RSU density on link a

\({\textbf {n}}\)

RSU deployment plan, i.e., \({\textbf {n}}=\{n_a,\forall a \in A\}\)

\(\rho ^k\)

RSU density of path k

\(\varvec{\rho }^w\)

Vector of \(\{\rho ^k: \forall k\in K^w\}\)

\(\mu ^w_i\)

Variable introduced to establish the path choice equilibrium

\(\overline{T}^w_i\)

Expected travel time for OD pair w using vehicle type i

\(C^w_i\)

Expected cost of commuting via vehicle type i

\(P^w_i\)

Probability of choosing vehicle type i for travelers with OD pair w

\(\lambda ^w\)

Variable introduced to establish the vehicle type choice equilibrium

\(\eta ^w_i, \xi ^w\)

Dual variables of the linear program associated with the VI problem

2.1.1 Path choice equilibrium

For any given OD pair \(w\in W\), let \(q_i^w\) be the number of travelers whose vehicle type choice is \(i\in I\), and \(f_i^{wk}\) the number of them who choose path \(k\in K^w\). At SUE, no traveler can reduce his/her perceived travel time by unilaterally changing his/her path choice. According to Sheffi (1984), the path choice SUE is characterized by

$$\begin{aligned} f_i^{wk} = P_i^{wk}q^w_i, \forall k \in K^w, w\in W, i\in I \end{aligned}$$
(1)

where \(P^{wk}_i\), defined as the probability of choosing path k for travelers with OD \(w\in W\) and vehicle type choice \(i\in I\), is given by the logit-based choice probability

$$\begin{aligned} P_i^{wk} = \frac{\exp {(-\theta ^w_iT^k)}}{\sum _{r\in K^w}\exp {(-\theta ^w_iT^r)}}, \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(2)

Note that the value of the dispersion parameter for CAV users \(\theta ^w_{\text{CAV}}\) decreases in the number of RSUs on the paths in \(K^w\) as well as the CAV penetration (i.e., \(q^w_{\text{CAV}}/Q^w\)), because the more RSUs and higher penetration of CAVs on the paths in \(K^w\), the more accurate information about the travel time of these paths will be collected and shared with the CAV users, and they are thus more likely to choose the path with minimum expected travel time.

Let \(n_a\) be the number of RSUs deployed on link \(a\in A\). The density of RSUs deployed on path \(k\in K^w\) is defined as the total number of RSUs deployed on path k divided by the path length, i.e., \(\rho ^k = \frac{\sum _{a\in A}\delta ^k_a n_a}{\sum _{a\in A}\delta ^k_a l_a}, k\in K^w\), where \(l_a\) is the length of link \(a\in A\). The implication here is that with a higher RSU density, there are more devices collecting and sending data per unit length of the link, thereby enabling more robust and reliable V2I communication and improving the information quality for CAV users. Without loss of generality, the dispersion parameter \(\theta ^w_{\text{CAV}}\) for CAV users with OD pair w is defined as a strictly increasing function in the RSU density on the paths belonging to \(K^w\) and the CAV penetration:

$$\begin{aligned} \theta ^w_{\text{CAV}} = g_{\theta }(\varvec{\rho }^w, q^w_{\text{CAV}}), \forall w\in W \end{aligned}$$
(3)

where \(\varvec{\rho }^w = \{\rho ^k:k\in K^w\}\) is the vector of the numbers of RSUs on the paths in \(K^w\).

We now propose an alternative condition for the path choice equilibrium. This condition, together with the one for vehicle type choice equilibrium (see Proposition 2 in Sect. 2.1.2), will facilitate the proposition of an equivalent variational inequality (VI) problem which can be easily solved.

Proposition 1

Given a RSU deployment plan \({\textbf {n}} = \{n_a: a\in A\}\), a feasible path flow pattern \({\textbf {f}}\in \Omega _{{\textbf {f}}}= \{{\textbf {f}}\ge 0:\sum _{k\in K^w}f^{wk}_i = q^w_i, \forall w\in W, i \in I\}\) is at equilibrium if and only if it satisfies

$$\begin{aligned} \theta ^w_iT^k + \ln {f^{wk}_i} = \mu ^w_i, \forall k \in K^w, w\in W, i\in I \end{aligned}$$
(4)

where \(\mu ^w_i = \ln {\frac{q^w_i}{\sum _{r\in K^w}\exp (-\theta ^w_iT^r)}},\forall w\in W, i\in I\).

Proof

We first prove the sufficiency of the proposition. Given a feasible path flow pattern \({\textbf {f}}\in \Omega _{{\textbf {f}}}\) such that condition (4) is satisfied, and the definition of \(\mu ^w_i\), we have

$$\begin{aligned}f^{wk}_i &= \exp (\mu ^w_i)\exp (-\theta ^w_i T^k), \\&\quad \forall k \in K^w, w\in W, i \in I \end{aligned}$$
(5)
$$\begin{aligned}&q^w_i = \sum _{r\in K^w} f^{wr}_i = \exp (\mu ^w_i)\sum _{r\in K^w}\exp (-\theta ^w_i T^r), \\&\forall w\in W, i \in I \end{aligned}$$
(6)

Dividing Eq. (5) by Eq. (6), we have

$$\begin{aligned} \frac{f^{wk}_i}{q^w_i} = \frac{\exp (-\theta ^w_i T^k)}{\sum _{r\in K^w}\exp (-\theta ^w_i T^r)} \end{aligned}$$
(7)

which is identical to the equilibrium conditions (1) and (2).

We then prove the necessity. Given a feasible path flow pattern \({\textbf {f}}\in \Omega _{{\textbf {f}}}\) at equilibrium (i.e., satisfying conditions (1) and (2)), we have

$$\begin{aligned} \begin{aligned} \theta ^w_i T^k + \ln (f^{wk}_i)&= \theta ^w_iT^k + \ln (P^{wk}_iq^w_i)\\&= \theta ^w_iT^k + \ln \left( \frac{\exp (-\theta ^w_i T^k)q^w_i}{\sum _{r\in K^w}\exp (-\theta ^w_i T^r)}\right) \\&= \ln {\frac{q^w_i}{\sum _{r\in K^w}\exp (-\theta ^w_iT^r)}}, \forall k \in K^w, w\in W, i\in I \end{aligned} \end{aligned}$$
(8)

which is exactly condition (4). \(\square\)

2.1.2 Vehicle type choice equilibrium

For any given OD pair \(w\in W\), let \(C^w_i\) be the long-term expected cost of using i-type vehicle, based on which travelers make their vehicle type choice. Clearly, \(C^w_i\) depends on the daily expected travel time on paths in \(K^w\), which is given by

$$\begin{aligned} \overline{T}^w_i = \sum _{k\in K^w} P^{wk}_iT^k,\ \forall w\in W, i\in I \end{aligned}$$
(9)

Here, we define \(C^w_i\) as an increasing function of \(\overline{T}^w_i\) which captures the impact of daily expected travel time on the long-term vehicle purchasing decision. Specifically,

$$\begin{aligned} C^w_i = g_C(\overline{T}^w_i) \end{aligned}$$
(10)

Then for travelers with OD \(w\in W\), the probability of choosing i-type vehicle is

$$\begin{aligned} P^{w}_i = \frac{\exp (-\theta C^w_i)}{\sum _{j\in I}\exp (-\theta C^w_j)}, \forall w\in W, i\in I \end{aligned}$$
(11)

The logit-based SUE condition for vehicle type choice is then given by

$$\begin{aligned} q^{w}_i = Q^wP^{w}_i, \forall w \in W, i \in I \end{aligned}$$
(12)

Similar to Proposition 1 in Sect. 2.1.1, we propose an equivalent condition for vehicle type choice equilibrium.

Proposition 2

Given a RSU deployment plan \({\textbf {n}} = \{n_a: a\in A\}\), a feasible vehicle type choice pattern \({\textbf {q}}\in \Omega _{{\textbf {q}}} = \{{\textbf {q}}\ge 0:\sum _{i\in I}q^{w}_i = Q^w, \forall w\in W\}\) is at equilibrium if and only if it satisfies

$$\begin{aligned} \theta C^w_i + \ln {q^{w}_i} = \lambda ^w, \forall w \in W \end{aligned}$$
(13)

where \(\lambda ^w = \ln {\frac{Q^w}{\sum _{j\in I}\exp (-\theta C^w_j)}},\forall w\in W\).

Proof

We first prove the sufficiency of the proposition. Given a feasible vehicle type choice pattern \({\textbf {q}}\in \Omega _{{\textbf {q}}}\) such that condition (13) is satisfied, and the definition of \(\lambda ^w\), we have

$$\begin{aligned} q^{w}_i= & {} \exp (\lambda ^w)\exp (-\theta C^w_i), \forall w\in W, i \in I \end{aligned}$$
(14)
$$\begin{aligned} Q^w= & {} \sum _{j\in I} q^w_j = \exp (\lambda ^w)\sum _{j\in I}\exp (-\theta C^w_j), \forall w\in W \end{aligned}$$
(15)

Dividing Eq. (14) by Eq. (15), we have

$$\begin{aligned} \frac{q^{w}_i}{Q^w} = \frac{\exp (-\theta C^w_i)}{\sum _{j\in I}\exp (-\theta C^w_j)} \end{aligned}$$
(16)

which is identical to the equilibrium conditions (11) and (12).

We then prove the necessity. Given a feasible vehicle type choice pattern \({\textbf {q}}\in \Omega _{{\textbf {q}}}\) at equilibrium (i.e., satisfying conditions (11) and (12)), we have

$$\begin{aligned} \begin{aligned} \theta C^w_i + \ln (q^{w}_i)&= \theta C^w_i + \ln (P^{w}_iQ^w)\\&= \theta C^w_i + \ln \left( \frac{\exp (-\theta C^w_i)Q^w}{\sum _{j\in I}\exp (-\theta C^w_j)}\right) \\&= \ln {\frac{Q^w}{\sum _{j\in I}\exp (-\theta C^w_j)}}, \forall w\in W, i\in I \end{aligned} \end{aligned}$$
(17)

which is exactly condition (13). \(\square\)

To guarantee the equilibrium of the whole system, the SUE conditions in Proposition 1 and Proposition 2 must be satisfied simultaneously.

2.1.3 Equivalent VI model

In this section, we reformulate the aforementioned equilibrium conditions into a variational inequality (VI) problem which is easy to solve. The VI of interest is defined as:

To find \(({\textbf {f}}^*,{\textbf {q}}^*)\in \Omega\), such that

$$\begin{aligned}&\sum _{i\in I}\sum _{w\in W}\sum _{k\in K^w} F^{wk}_i({\textbf {f}}^*,{\textbf {q}}^*)(f^{kw}_i-f^{wk*}_i) \\&\quad +\sum _{i\in I}\sum _{w\in W} G^w_i({\textbf {f}}^*,{\textbf {q}}^*)(q^w_i-q^{w*}_i)\ge 0, \forall ({\textbf {f}},{\textbf {q}})\in \Omega \end{aligned}$$
(18)

where

$$\begin{aligned}&F^{wk}_i({\textbf {f}}^*,{\textbf {q}}^*) = \theta ^w_iT^k + \ln {f^{wk}_i}, \\&\quad \forall k \in K^w, w\in W, i \in I \end{aligned}$$
(19)

and

$$\begin{aligned}&G^w_i({\textbf {f}}^*,{\textbf {q}}^*) = \theta C^w_i + \ln {\sum _{k\in K^w}\exp (-\theta ^w_iT^k)}, \\&\quad \forall w\in W, i\in I \end{aligned}$$
(20)

The feasible set \(\Omega\) is defined by

$$\begin{aligned}&\sum _{k\in K^w}f^{wk}_i = q^w_i, \forall w\in W, i \in I \end{aligned}$$
(21)
$$\begin{aligned}&\sum _{i\in I}q^{w}_i = Q^w, \forall w\in W \end{aligned}$$
(22)
$$\begin{aligned}&f^{wk}_i \ge 0, k\in K^w, w \in W, i\in I \end{aligned}$$
(23)
$$\begin{aligned}&q^w_i \ge 0, w\in W, i\in I \end{aligned}$$
(24)

Proposition 3

A joint path flow and vehicle type choice pattern \(({\textbf {f}}^*,{\textbf {q}}^*)\in \Omega\) satisfies equilibrium conditions (4) and (13) if and only if it solves the VI problem (18).

Proof

Let \(({\textbf {f}}^*,{\textbf {q}}^*)\) be a solution to the VI problem (18). Note that there are a few variables dependent on \(({\textbf {f}},{\textbf {q}})\), such as \(\theta ^w_i\) and \(T^k\), etc. For simplicity, we use asterisk to indicate their values at \(({\textbf {f}},{\textbf {q}}) = ({\textbf {f}}^*,{\textbf {q}}^*)\). For example, we use \(\theta ^{w*}_i\) to denote \(\theta ^w_i({\textbf {f}}^*,{\textbf {q}}^*)\).

Since \(({\textbf {f}}^*,{\textbf {q}}^*)\) solves VI problem (18), it must satisfies the optimality condition of the VI problem. That is, there exists a vector \((\varvec{\mu }^*,\varvec{\lambda }^*,\varvec{\sigma }^*,\varvec{\pi }^*)\) such that

$$\begin{aligned}&F^{wk}_i({\textbf {f}}^*,{\textbf {q}}^*) - \mu ^{w*}_i -\sigma ^{wk*}_i = 0, \\&\quad \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(25)
$$\begin{aligned}&G^w_i({\textbf {f}}^*,{\textbf {q}}^*) + \mu ^{w*}_i - \lambda ^{w*}-\pi ^w_i = 0, \forall w\in W,i\in I \end{aligned}$$
(26)
$$\begin{aligned}&f^{wk*}_i\sigma ^{wk*}_i = 0, \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(27)
$$\begin{aligned}&q^{w*}_i\pi ^{w*}_i = 0, \forall w\in W, i\in I \end{aligned}$$
(28)
$$\begin{aligned}&\sigma ^{wk*}_i\ge 0, \pi ^{w*}_i\ge 0, \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(29)

Eqs. (21)–(24).

Note that if \(f^{wk}_i\) approaches to 0, then \(F^{wk}_i\) will approach to negative infinity, so (18) will not hold for all \(({\textbf {f,q}})\) and\(({\textbf {f}}^*,{\textbf {q}}^*)\in \Omega\). Therefore, we must have \(f^{wk}_i>0\). From (27), we have \(\sigma ^{wk*}_i=0\ (\forall k\in K^w, w\in W, i\in I)\). Using this and Eq. (25), we have

$$\begin{aligned}&f^{wk*} = \exp (\mu ^{w*}_i)\exp (-\theta ^{w*}_iT^{k*}), \\&\quad \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(30)

By summing Eq. (30) over \(k\in K^w\) and using Eq. (21), we have

$$\begin{aligned}&q^{w*}_i = \sum _{k\in K^w} f^{wk*}_i = \exp (\mu ^{w*}_i)\sum _{k\in K^w}\exp (-\theta ^{w*}_iT^{k*}), \\&\quad \forall w\in W, i\in I \end{aligned}$$
(31)

From Eqs. (30) and (31), we immediately have the equilibrium flow pattern condition

$$\begin{aligned} \frac{f^{wk*}_i}{q^{w*}_i} = \frac{\exp (-\theta ^{w*}_i T^{k*})}{\sum _{r\in K^w}\exp (-\theta ^{w*}_i T^{r*})}, \forall w\in W, i\in I \end{aligned}$$
(32)

and the alternative expression for multipliers \(\mu ^{w*}_i\):

$$\begin{aligned} \mu ^{w*}_i = \ln {\frac{q^{w*}_i}{\sum _{r\in K^w}\exp (-\theta ^{w*}_i T^{r*})}}, \forall w\in W, i\in I \end{aligned}$$
(33)

which is identical with the condition in Proposition 1.

The condition for vehicle choice equilibrium can be derived in similar manner from the optimality condition of the VI problem. Since \(f^{wk}_i>0\), we must have \(q^w_i>0\). Using Eqs. (26), (28) and (33), we have

$$\begin{aligned} q^{w*}_i = \exp (\lambda ^{w*})\exp (-\theta C^{w*}_i), \forall w\in W, i \in I \end{aligned}$$
(34)

Summing Eq. (34) over \(i\in I\) and using Eq. (22), we have

$$\begin{aligned} Q^{w} = \sum _{i\in I}q^{w*}_i = \exp (\lambda ^{w*})\sum _{j\in I}\exp (-\theta C^{w*}_j), \forall w\in W \end{aligned}$$
(35)

From Eqs. (34) and (35), we have the vehicle type choice equilibrium condition and the alternative expression for multipliers \(\lambda ^{w*}\):

$$\begin{aligned} \frac{q^{w*}_i}{Q^w}= & {} \frac{\exp (-\theta C^{w*}_i)}{\sum _{j\in I}\exp (-\theta C^{w*}_j)}, \forall w\in W, i\in I \end{aligned}$$
(36)
$$\begin{aligned} \lambda ^{w*}= & {} \ln {\frac{Q^w}{\sum _{j\in I}\exp (-\theta C^{w*}_j)}},\forall w\in W \end{aligned}$$
(37)

which are identical to the condition in Proposition 2. Therefore, any \(({\textbf {f}}^*,{\textbf {q}}^*)\) that solves the VI problem (18) must satisfy both SUE conditions in Propositions 1 and 2.

Conversely, given an equilibrium path flow and vehicle type choice pattern \(({\textbf {f}}^*,{\textbf {q}}^*)\), we can find the optimal multipliers \(\mu _i^{w*}\) and \(\lambda ^{w*}\) using Eqs. (33) and (37). Since the optimality condition is both necessary and sufficient for VI problems with convex polyhedral feasible set (Facchinei and Pang 2003), \(({\textbf {f}}^*,{\textbf {q}}^*)\) is also a solution to VI problem (18). \(\square\)

2.2 RSU deployment model

The goal of the optimization model is to decide a RSU deployment plan to minimize the total network delay and emissions. To measure the emissions of the traffic flow, we adopt the following emission function based on average link speed (Wallace et al. 1984; Chen and Yang 2012; Yin and Lawphongpanich 2006):

$$\begin{aligned} e_{a}(x_a) = 0.2038t_a(x_a)\exp \left( \frac{0.7962l_a}{t_a(x_a)}\right) ,\forall a\in A \end{aligned}$$
(38)

where \(e_a(x_a)\) is the volume of CO emissions (g/veh) on link \(a\in A\) when the link flow is \(x_a\). \(l_a\) is the length (km) of link a and the link travel time \(t_a(x_a)\) is in the unit of km/min.

Now we consider the problem of minimizing the total network delay and the emissions by optimizing the deployment of RSUs on the network. Let \({\textbf {n}}=\{n_a:\forall a\in A\}\) be the RSU deployment plan. It parameterizes the mixed flow traffic equilibrium and its equivalent VI problem (18). For any given RSU deployment plan \({\textbf {n}}\), let \(SOL(\Omega ,{\textbf {F(n)}},{\textbf {G(n)}})\) denote the solution set of VI problem (18). The optimization problem for mitigating delay and emissions via RSU deployment (named by RSUD) is then defined as a mathematical program with equilibrium constraints (MPEC):

$$\begin{aligned}&\text {(RSUD)}\ \min _{{\textbf {n}}}\ \gamma _1\sum _{a\in A} t_a(x_a)x_a + \gamma _2\sum _{a\in A} e_a(x_a)x_a \end{aligned}$$
(39)
$$\begin{aligned}&s.t.\ \sum _{a\in A} n_a \le N \end{aligned}$$
(40)
$$\begin{aligned}&\underline{n}_a\le n_a \le \bar{n}_a, \forall a\in A \end{aligned}$$
(41)
$$\begin{aligned}&x_a = \sum _{i\in I}\sum _{w\in W}\sum _{k\in K^w} \delta ^k_af^{wk}_i, \forall a\in A \end{aligned}$$
(42)
$$\begin{aligned}&n_a\in Z_+, \forall a\in A \end{aligned}$$
(43)
$$\begin{aligned}&({\textbf {f}},{\textbf {q}})\in SOL(\Omega ,{\textbf {F(n)}},{\textbf {G(n)}}) \end{aligned}$$
(44)

The objective (39) is to minimize the weighted sum of the network delay and emissions. Constraint (40) means that the total number of RSUs deployed on the network is no more than N due to limited budget. Constraint (41) is the RSU requirement for each link. It requires that the number of RSUs deployed on link \(a\in A\) is no less than \(\underline{n}_a\) and no more than \(\bar{n}_a\), where \(\underline{n}_a\) and \(\bar{n}_a\) are predetermined parameters whose values are related to the link length, the coverage of each RSU and the desired system performance in terms of communication reliability and robustness, etc. Constraint (42) is the definitional expression of link flow \(x_a\). Constraint (43) requires \(n_a\) to be integers. Constraint (44) is the equilibrium constraint which ensures the objective is minimized under equilibrium.

3.Solution algorithm

The (RSUD) model mentioned above is an integer program with equilibrium constraints which is difficult to solve due to the non-convex feasible region and the violation of Magasarian-Fromovitz constraint qualification (MFCQ). In order to obtain a promising solution to the (RSUD), we apply the active-set algorithm proposed by Zhang et al. (2009). It solves the problem by iteratively solving a series of simple subproblems and eventually reaching a strongly stationary point.

First, we introduce the gap function-based method (Aghassi et al. 2006) to solve the VI traffic assignment problem, which will constitute a subproblem incorporated in the active-set algorithm. The gap function-based method transforms a VI problem with polyhedral constraints into an equivalent nonlinear mathematical program and has been successfully used to solve many similar problems in literature (He et al. 2013; Chen et al. 2016, 2016). Specifically, the solution to VI problem (18) can be obtained by solving the following mathematical program (MP-1).

$$\begin{aligned}&\text {(MP-1)}\ \min _{{\textbf {f}},{\textbf {q}},\varvec{\eta },\varvec{\zeta }}\ \sum _{i\in I}\sum _{w\in W}\sum _{k\in K^w} F^{w,k}_i({\textbf {f}},{\textbf {q}})f^{wk}_i \\&\quad + \sum _{i\in i}\sum _{w\in W}G^w_i({\textbf {f}},{\textbf {q}})q^w_i - \sum _{w\in W} Q^w\xi ^w \end{aligned}$$
(45)
$$\begin{aligned} s.t.,\ (3),(9),(10),(21)-(24),(42), \text{and} \end{aligned}$$
$$\begin{aligned}&\theta ^w_iT^k + \ln {f^{wk}_i}\ge \eta ^w_i, \forall k\in K^w w\in W, i\in I \end{aligned}$$
(46)
$$\begin{aligned}&\theta C^w_i + \ln {\sum _{k\in K^w}\exp (-\theta ^w_iT^k)}\ge \xi ^w-\eta ^w_i, \\&\quad \forall w\in W, i\in I \end{aligned}$$
(47)

The objective function (45) is the gap between the primal and dual value of linear program (LP) associated with VI problem (18). Constraints (46) and (47) ensure the dual feasibility of the LP. Note that constraint (46) can be recast as \(f^{wk}_i\ge \exp (\eta ^w_i-\theta ^{wk}_iT^k)\), which renders any point with \(f^{wk}_i=0\) infeasible. (MP-1) is a continuous nonlinear program which in our study is solved using the nonlinear optimization solver CONOPT.Footnote 1

The second subproblem is defined as a nonlinear mathematical program (MP-2).

$$\begin{aligned} \text {(MP-2)}\ \min _{{\textbf {y,f,q}},\varvec{\mu ,\lambda }}\ \gamma _1\sum _{a\in A} t_a(x_a)x_a + \gamma _2\sum _{a\in A} e_a(x_a)x_a \end{aligned}$$
(48)
$$\begin{aligned} s.t.,\ (3),(9),(10),(21)-(24),(42), \text{and} \end{aligned}$$
$$\begin{aligned}&\theta ^w_iT^k + \ln {f^{wk}_i} - \mu ^w_i = 0, \forall k\in K^w, w\in W, i\in I \end{aligned}$$
(49)
$$\begin{aligned}&\theta C^w_i + \ln {q^w_i} - \lambda ^w_i = 0, \forall w\in W, i\in I \end{aligned}$$
(50)
$$\begin{aligned}&n_a = \sum _{m=1}^{M_a} 2^{m-1}y_{am}, \forall a\in A \end{aligned}$$
(51)
$$\begin{aligned}&y_{am} = 0, \forall (a,m)\in \Pi _0 \end{aligned}$$
(52)
$$\begin{aligned}&y_{am} = 1, \forall (a,m)\in \Pi _1 \end{aligned}$$
(53)

Constraints (49) and (50) are the equilibrium conditions (or the KKT conditions of VI problem (18)). (MP-2) shares the same objective with (RSUD) to minimize the delay and emissions under equilibrium constraints. However, every decision variable \(n_a\) in (MP-2) is fixed by \(M_a\) binary variable \(y_{am}\) via constraints (51)–(53), where \(M_a\) is the minimum number of binary variables required to represent \(n_a\) up to the level of \(\bar{n}_a\). The two index sets \(\Pi _0\) and \(\Pi _1\) are introduced to indicate whether \(y_{am}\) is set as 0 (if \((a,m)\in \Pi _0\)) or set as 1 (if \((a,m)\in \Pi _1\)). Note that since all binary variables \(y_{am}\) involved in (MP-2) are fixed using equality constraints (52) and (53). (MP-2) is treated as a continuous optimization problem. The goal of solving (MP-2) is to obtain the multipliers associated with constraints (52) and (53), which will provide guidance for us to adjust the two index sets. Solution obtained for (MP-1) will serve as an initial feasible solution to (MP-2) to facilitate the solution process. In (MP-2), the 0-1 variables \(y_{am}\)’s are fixed as either 1 or 0 by constraints (52) and (53). Hence, it can be regarded as a continuous nonlinear program. When (MP-2) is solved via CONOPT, the multipliers associated with constraints (52) and (53) will also be derived.

After we obtain the multipliers \(\alpha _{am}\) and \(\beta _{am}\) associated with constraints (52) and (53), respectively, the third subproblem (KP) plays its role to generate a potential adjustment to the deployment plan, i.e., the two index sets.

$$\begin{aligned}&\text {(KP)}\ \min _{{\textbf {u,v}}}\ \sum _{(a,m)\in \Pi _0} \alpha _{am}u_{am} - \sum _{(a,m)\in \Pi _1}\beta _{am}v_{am} \end{aligned}$$
(54)
$$\begin{aligned}&s.t.\ \sum _{(a,m)\in \Pi _0}2^{m-1}u_{am} - \sum _{(a,m)\in \Pi _1}2^{m-1}v_{am}\le N \\&\quad - \sum _{(a,m)\in \Pi _1}2^{m-1} \end{aligned}$$
(55)
$$\begin{aligned}&\underline{n}_a \le \sum _{m:(a,m)\in \Pi _1}2^{m-1} + \sum _{m:(a,m)\in \Pi _0}2^{m-1}u_{am} \\&\quad - \sum _{m:(a,m)\in \Pi _1}2^{m-1}v_{am}\le \bar{n}_a, \forall a\in A \end{aligned}$$
(56)
$$\begin{aligned}&\sum _{(a,m)\in \Pi _0} \alpha _{am}u_{am} - \sum _{(a,m)\in \Pi _1}\beta _{am}v_{am}\ge \kappa \end{aligned}$$
(57)
$$\begin{aligned}&u_{am}\in \{0,1\}, \forall (a,m)\in \Pi _0 \end{aligned}$$
(58)
$$\begin{aligned}&v_{am}\in \{0,1\}, \forall (a,m)\in \Pi _1 \end{aligned}$$
(59)

where \(u_{am}\) and \(v_{am}\) are binary decision variables to specify whether (am) should be shifted from set \(\Pi _0\) to \(\Pi _1\), and (am) shifted from set \(\Pi _1\) to \(\Pi _0\), respectively.

The objective function of (KP) measures the potential reduction resulted from the adjustment \(({\textbf {u,v}})\). Constraints (55) and (56) ensure that the adjustment will not violate the budget and link requirement constraints. (57) is used to rule out adjustments that fail to actually reduce the objective of (MP-2). Constraints (58) and (59) indicate the binary nature of the decision variables.

The overall procedure of the active-set algorithm applied to solve (RSUD) is described as follows. The flow chart of the algorithm is provided in Fig. 1.

  1. (1)

    Step 0: Set \(\zeta\) = 1, \(\Pi _0^1=\{(a,m): a\in A, m=1,2,\ldots ,M_a\}\), and \(\Pi _1^1=\emptyset\)

  2. (2)

    Step 1: Solve (MP-1) with RSU deployment plan \({\textbf {n}}^\zeta\) determined by \({\textbf {y}}^\zeta\) compatible with \((\Pi _0^\zeta ,\Pi _1^\zeta )\). Then use the solution to (MP-1) as the initial feasible solution for solving (MP-2) with \((\Pi _0^\zeta ,\Pi _0^\zeta )\), and obtain objective value \(Z^\zeta\), the solution \(({\textbf {y}}^\zeta ,{\textbf {f}}^\zeta ,{\textbf {q}}^\zeta , \varvec{\mu }^\zeta , \varvec{\lambda }^\zeta )\) and the multipliers \(\alpha ^\zeta _{am}\) and \(\beta ^\zeta _{am}\) associated with constraints (52) and (53), respectively.

  3. (3)

    Step 2: Set \(\kappa =-\infty\), and adjust the active sets \((\Pi _0^\zeta ,\Pi _1^\zeta )\) as follows:

    1. (a)

      Solve (KP) with multiplies \(\alpha ^\zeta _{am}\) and \(\beta ^\zeta _{am}\), and obtain the solution \((\hat{{\textbf {u}}},\hat{{\textbf {v}}})\). If the objective value is 0, stop and \(({\textbf {y}}^\zeta ,{\textbf {f}}^\zeta ,{\textbf {q}}^\zeta , \varvec{\mu }^\zeta , \varvec{\lambda }^\zeta )\) is the best solution. Otherwise, go to Step 2b.

    2. (b)

      Set

      1. i.

        \(\chi\) = \(\sum _{(a,m)\in \Pi _0} \alpha ^\zeta _{am}\hat{u}_{am} - \sum _{(a,m)\in \Pi _1}\beta ^\zeta _{am}\hat{v}_{am}\)

      2. ii.

        \(\hat{\Pi }_0 = \left( \Pi _0^\zeta -\left\{ (a,m)\in \Pi _0^\zeta :\hat{u}_{am}=1\right\} \right) \cup \left\{ (a,m)\in \Pi _1^\zeta :\hat{v}_{am}=1\right\}\)

      3. iii.

        \(\hat{\Pi }_1 = \left( \Pi _1^\zeta -\left\{ (a,m)\in \Pi _1^\zeta :\hat{v}_{am}=1\right\} \right) \cup \left\{ (a,m)\in \Pi _0^\zeta :\hat{u}_{am}=1\right\}\)

    3. (c)

      Solve (MP-1) with a deployment plan \({\textbf {y}}^\zeta\) compatible with \((\hat{\Pi }_0^\zeta ,\hat{\Pi }_1^\zeta )\). If the objective value \(Z'\) is less than \(Z^\zeta\), go to Step 2d because the pair \((\hat{\Pi }_0^\zeta ,\hat{\Pi }_1^\zeta )\) leads to a reduction in the objective. Otherwise, set \(\kappa = \chi + \varepsilon\), where \(\varepsilon >0\) is sufficiently small, and return to Step 2a.

    4. (d)

      Set \(\Pi _0^{\zeta +1}=\hat{\Pi }_0^\zeta ,\Pi _1^{\zeta +1}=\hat{\Pi }_1^\zeta\) and \(\zeta = \zeta +1\). Go to Step 1.

Fig. 1
figure 1

Flow chart of the solution algorithm

It has been proved that the active-set algorithm terminates after a finite number of iterations and that the solution obtained by the algorithm must be a strongly stationary point. Interested readers are recommended to refer to Zhang et al. (2009) for more details.

4.Numerical examples

In this section, we will conduct numerical experiments to validate the performance of our proposed model and algorithm. The numerical experiments are based on the standard Nguyen-Dupuis network (see Fig. 2). There are four OD pairs whose demand is given in Table 2. The paths connecting the four OD pairs are listed in Table 7 in the Appendix.

Fig. 2
figure 2

Nguyen-Dupuis network (Nguyen and Dupuis 1984)

Table 2 Travel demand of each OD pair

 No.

OD pair

Demand (veh/h)

1

1->2

400

2

1->3

800

3

4->2

600

4

4->3

200

We use the Bureau of Public Roads (BPR) function to measure the flow-dependent travel time on each link, i.e., the link travel time is

$$\begin{aligned} t_a(x_a) = t^0_a\left( 1+0.15\left( \frac{x_a}{B_a}\right) ^4\right) , \forall a\in A \end{aligned}$$
(60)

where \(t^0_a\) represents the free-flow travel time on link a and \(B_a\) is the capacity of the link. Their values are presented in Table 8 in the Appendix.

In Xie and Liu (2022), the dispersion parameter \(\theta _{\text{CAV}}^w\) for CAV users is assumed as a linear function in the penetration rate of CAV. Here, we adopt the same assumption, defining \(\theta ^w_{\text{CAV}} = g_{\theta }(\varvec{\rho }^w, q^w_{\text{CAV}})\) as a linear function in both the CAV penetration (\(\frac{q^w_{\text{CAV}}}{Q^w}\)) and the average RSU density of paths in \(K^w\) (\(\frac{1}{|K^w|}\sum _{k\in K^w}\rho ^k\)). Specifically,

$$\begin{aligned} \theta ^w_{\text{CAV}} = \theta ^w_{\text{RV}}+\psi _1 \frac{q^w_{\text{CAV}}}{Q^w} + \frac{\psi _2}{|K^w|}\sum _{k\in K^w}\rho ^k , \forall w\in W \end{aligned}$$
(61)

where \(\psi _1\) and \(\psi _2\) are coefficients that gauge the contribution of CAV penetration and RSU density to the information quality. The first term is the dispersion parameter for RV drivers to capture the information quality received by RV drivers. Equation (61) implies that the information quality of CAVs is higher than that of RVs due to the CAV technology and the assistance from RSUs.

For the function of long-term expected cost \(C^w_i\), we use the same formulation in Xie and Liu (2022), i.e.,

$$\begin{aligned} C^w_i = b^1_i \overline{T}^w_i + \frac{b^2_i\tau _i}{L_i}d^w + b^3_id^w, \forall w\in W, i\in I \end{aligned}$$
(62)

where \(b_i^1\) is the value of time when driving i-type vehicle, so the first term is the average travel time cost per trip. \(\tau _i\) is the vehicle purchase price and \(b^2_i\) is a scaling factor applied to the purchase price to capture other overhead costs (e.g., the purchase tax and the fee for granting a license plate, etc.). \(L_i\) is the maximum mileage that the vehicle can run in its life cycle. \(d^w\) is the average length of the paths in \(K^w\). The second term is the fixed cost of owning a vehicle averaged to a single trip. With \(b^3_i\) denoting the operational cost (e.g., the refuel cost and maintenance cost, etc.) per unit distance, the third term represents the operational cost per trip. The setting of these parameters are summarized in Table 9 in the Appendix. They are basically set in line with the settings in Xie and Liu (2022), with some of them scaled due to unit conversion.

The link requirement parameters \(\underline{n}_a\) and \(\bar{n}_a\) are set according to the settings in Li et al. (2020), in which the coverage of a RSU varies from 0.25 km to 1 km. Accordingly, the maximum RSU density is 4 RSUs per kilometer and we do not impose a lower bound on the number of RSUs to be deployment on each link, i.e., \(\underline{n}_a=0\). The specific value for \(\underline{n}_a\) and \(\bar{n}_a\) for each link is shown in Table 8 in the Appendix. Using this setting, the maximum number of RSU required for the network is 383. We thus set the RSU deployment budget N as 200 to make sure that the numerical example can reflect the impact of limited budget. Finally, we attribute equal weights to the two objectives, i.e., \(\gamma _1 = \gamma _2=1\).

Table 3 Validation of the path choice equilibrium

OD

Path

RV

CAV

\(\theta ^w_iT^k+\ln {f^wk_i}\)

\(\mu ^w_i\)

\(\theta ^w_iT^k+\ln {f^wk_i}\)

\(\mu ^w_i\)

1

1

3.63

3.63

5.93

5.93

2

3.63

5.93

3

3.63

5.93

4

3.63

5.93

5

3.63

5.93

6

3.63

5.93

7

3.63

5.93

8

3.63

5.93

2

1

4.74

4.74

7.65

7.65

2

4.74

7.65

3

4.74

7.65

4

4.74

7.65

5

4.74

7.65

6

4.74

7.65

3

1

4.55

4.55

7.24

7.24

2

4.55

7.24

3

4.55

7.24

4

4.55

7.24

5

4.55

7.24

4

1

3.33

3.33

5.89

5.89

2

3.33

5.89

3

3.33

5.89

4

3.33

5.89

5

3.33

5.89

\(\mu ^w_i=\ln {q^w_i}-\ln {\sum _{r\in K^w}\exp (-\theta ^w_iT^r)}\)

Table 4 Validation of the vehicle type choice equilibrium

OD

Vehicle type

\(\theta C^w_i + \ln {q^w_i}\)

\(\lambda ^w\)

1

RV

7.37

7.37

CAV

7.37

2

RV

8.44

8.44

CAV

8.44

3

RV

7.92

7.92

CAV

7.92

4

RV

6.65

6.65

CAV

6.65

\(\lambda ^w = \ln {Q^w}-\ln {\sum _{j\in I}\exp (-\theta C^w_j)}\)

To verify the correctness of our solution algorithm, we validate whether the flow and vehicle type choice pattern obtained from the algorithm is at equilibrium. The validation is accomplished through checking whether the solution satisfies the equilibrium conditions proposed in Proposition 1 and 2. The results are presented in Tables 3 and 4. As we can see from Table 3, for each OD pair w and vehicle type i, the expression \(\theta ^w_iT^k+\ln {f^wk_i}\) has the same value for all path \(k\in K^w\), and it equals \(\mu ^w_i=\ln {q^w_i}-\ln {\sum _{r\in K^w}\exp (-\theta ^w_iT^r)}\), which means that flow pattern is at equilibrium according to Proposition 1. Likewise, from Table 4 we can see that for each OD pair w, \(\theta C^w_i + \ln {q^w_i}\) has the same value for both vehicle types, which equals \(\lambda ^w = \ln {Q^w}-\ln {\sum _{j\in I}\exp (-\theta C^w_j)}\). Therefore, we can claim that the vehicle type choice is at equilibrium by Proposition 2.

Fig. 3
figure 3

Convergence of the solution algorithm

Figure 3 shows the convergence behavior of the active-set algorithm. In our example, the algorithm terminates after a handful of iterations, which agrees with the theoretical property of the algorithm (Zhang et al. 2009). Moreover, it returns a considerably better solution than the initial solution, which must be a strongly stationary point of the origin problem (Zhang et al. 2009).

Table 5 Link performance before and after RSU deployment

Link

Before deployment

After deployment

Flow (veh/h)

Travel time (h)

Emissions (g/veh)

Flow (veh/h)

Travel time (h)

Emissions (g/veh)

Deployment (\(n_a\))

1

748.76

0.80

10.42

747.86

0.79

10.38

16

2

451.24

0.73

9.86

452.14

0.74

9.91

16

3

507.36

1.08

14.11

484.60

0.93

12.20

21

4

292.64

0.34

5.43

315.40

0.39

6.00

0

5

769.65

0.23

3.05

749.68

0.21

2.84

7

6

486.47

0.20

3.46

482.78

0.20

3.44

0

7

759.77

0.15

2.37

758.08

0.15

2.37

12

8

382.70

0.40

6.23

350.56

0.34

5.62

30

9

280.52

0.10

1.84

298.82

0.11

1.90

0

10

479.25

0.30

4.59

459.26

0.27

4.32

21

11

358.93

0.16

2.98

392.00

0.16

3.01

0

12

492.21

0.18

3.41

483.58

0.18

3.40

23

13

286.91

0.25

3.98

314.60

0.29

4.48

0

14

874.91

0.44

6.02

834.14

0.38

5.30

14

15

641.07

0.62

8.47

608.00

0.53

7.39

21

16

713.09

0.77

10.23

685.40

0.68

9.09

19

17

372.82

0.33

4.73

358.96

0.30

4.37

0

18

78.41

0.23

4.54

93.18

0.23

4.54

0

19

286.91

0.30

4.87

314.60

0.35

5.48

0

Delay

4058.85 h

3725.90 h (− 8.20%)

Emissions

57.33 kg/h

53.24 kg/h (− 7.13%)

% CAV

55.35%

55.89%

In Table 5, we present the vehicle flow, travel time and hourly emissions on each link before and after the optimal RSU deployment plan. It is assumed that there is no RSU installed on the network before the deployment. We can see that the optimal deployment plan has yield a significantly different pattern on the link flow, link travel time and link emissions. Specifically, such a different pattern has resulted in a considerable reduction on the total network delay (8.20%) and the emissions (7.13%), as shown in Table 5. It indicates that it is viable to effectively leverage the flow through RSU deployment to reduce both total system delay and emissions. Compared to the total delay and emissions, the increase in CAV penetration due to RSU deployment is not significant (0.54%), implying that the reduction of the former mainly roots in the change of path choices of CAV users but not the CAV adoption rate.

Table 6 Path flow and travel time before and after the deployment

OD

Path

Before deployment

After deployment

RV flow (veh/h)

CAV flow (veh/h)

Total flow (veh/h)

Travel time (h)

RV flow (veh/h)

CAV flow (veh/h)

Total flow (veh/h)

Travel time (h)

1

1

28.51

49.90

78.41

1.12

26.86

66.32

93.18

1.13

2

25.69

37.16

62.85

1.47

24.38

40.40

64.78

1.45

3

21.10

21.29

42.39

2.13

20.76

17.75

38.51

1.99

4

18.76

15.27

34.03

2.52

18.96

11.14

30.10

2.29

5

26.00

38.44

64.44

1.43

24.64

42.67

67.31

1.42

6

21.35

22.03

43.38

2.09

20.98

18.74

39.73

1.95

7

18.98

15.80

34.78

2.48

19.16

11.77

30.93

2.26

8

20.39

19.34

39.73

2.24

20.17

15.30

35.47

2.09

2

1

60.11

69.88

129.99

2.28

60.28

68.25

128.53

2.14

2

53.45

50.12

103.56

2.67

55.04

42.00

97.05

2.44

3

60.84

72.28

133.12

2.24

60.93

72.26

133.19

2.10

4

54.09

51.84

105.93

2.63

55.64

44.47

100.10

2.40

5

58.10

63.47

121.57

2.39

58.56

58.47

117.03

2.23

6

75.04

130.79

205.83

1.54

70.20

153.90

224.10

1.63

3

1

61.69

91.54

153.23

1.72

59.72

107.01

166.73

1.55

2

50.66

52.39

103.05

2.37

50.85

44.02

94.87

2.09

3

45.04

37.54

82.59

2.76

46.43

26.63

73.07

2.39

4

48.38

45.99

94.37

2.53

48.88

35.35

84.23

2.22

5

64.21

102.54

166.76

1.58

60.98

120.12

181.10

1.48

4

1

14.01

13.31

27.32

2.53

14.31

10.12

24.44

2.23

2

12.45

9.37

21.82

2.92

13.07

6.24

19.31

2.54

3

13.38

11.60

24.98

2.68

13.76

8.20

21.95

2.37

4

17.75

27.05

44.81

1.74

17.17

26.63

43.80

1.63

5

22.93

58.15

81.08

0.88

20.58

69.92

90.50

1.02

Note: Boldface: shortest path

Table 6 presents the path flow and travel time before and after the optimal deployment of RSUs. The shortest path between each OD pair is in boldface. As we can see, the four shortest paths carry significantly more flow after the deployment, and three of them have a higher travel time than before. Moreover, all non-shortest paths become less congested after the deployment. This suggests that the RSU deployment can reduce travel time differences among paths between the same OD pair. It is interesting to note that though the travel time among different paths become more even, the flow distribution of RVs and CAVs exhibits different changes. For each OD pair, the RV flow is distributed more evenly on the paths after the deployment while the CAV flow tends to more concentrated. This is because RSUs provide more accurate information for CAV users, so they are more likely to select the shortest path despite of the reduced travel time difference among paths. Therefore, the shortest path for each OD pair attracts more CAVs after the deployment of RSUs. On the other hand, since RV drivers can hardly benefit from the RSU deployment in terms of information quality, their path choice will become more random as the travel time for different paths gets closer to each other.

Fig. 4
figure 4

Network delay and emissions under different deployment budgets

Finally, we conduct a sensitivity analysis to figure out the resulted network delay and emissions under different budgets. The result is shown in Fig. 4. As expected, when we allow more RSUs to be deployed on the network, the network delay and emissions can be further reduced. However, the benefit from increasing deployment budget tends to diminish with the increase of budget. Therefore, before implementing any RSU deployment project aiming at reducing congestion and emissions, it is necessary for the practitioners to conduct careful cost-benefit analysis in order to well balance the investment and return, thereby attaining a desirable return-to-investment ratio.

5.Conclusions

Motivated by the fact that the deployment of RSUs may improve the quality of information received by CAVs, thereby influencing CAVs’ routing behaviors and adoption rate, we investigate the optimal deployment plan of RSUs under the mixed traffic flow to reduce network congestion and emissions. First, we develop an SUE model to describe the equilibrium path flow and vehicle type choice pattern given a fixed RSU deployment plan. With the established equilibrium model, the deployment of RSUs is modeled as a mathematical program with equilibrium constraints. The active-set algorithm is adopted to solve the deployment model efficiently, in which the equilibrium conditions are reformulated as a VI problem and solved with gap function-based method. Numerical experiments are conducted to demonstrate the viability of deploying RSUs to achieve a desirable flow pattern with less congestion and emissions. It is also shown that RSU deployment have different impacts on RV and CAV flows. The sensitivity analysis suggests that the marginal benefit from RSU deployment may decrease with the budget for RSU deployment.

There are a few possible extensions to this work. For example, models can be extended to integrate other management strategies such as congestion/emission pricing and tradable credit to further regulate CAV and RV users’ routing behaviors. Moreover, since RSUs can provide information in a real-time manner, CAV users may make their route choice dynamically. Therefore, future research may extend the deployment model of RSUs by considering the dynamic traffic assignment in the lower level problem.

References

[1] Aghassi M, Bertsimas D, Perakis G (2006) Solving asymmetric variational inequalities via convex optimization. Oper Res Lett 34(5):481–490
[2] Barth M, Boriboonsomsin K (2009) Energy and emissions impacts of a freeway-based dynamic eco-driving system. Transp Res D Transp Environ 14(6):400–410
[3] Chen S, Du L (2017) Simulation study of the impact of local real-time traffic information provision strategy in connected vehicle systems. Int J Transp Sci Technol 6(4):229–239
[4] Chen L, Yang H (2012) Managing congestion and emissions in road networks with tolls and rebates. Transp Res B Methodol 46(8):933–948
[5] Chen Z, He F, Yin Y (2016) Optimal deployment of charging lanes for electric vehicles in transportation networks. Transp Res B Methodol 91:344–365
[6] Chen Z, He F, Zhang L, Yin Y (2016) Optimal deployment of autonomous vehicle lanes with endogenous market penetration. Transp Res C Emerg Technol 72:143–156
[7] Chen Z, He F, Yin Y, Du Y (2017) Optimal design of autonomous vehicle zones in transportation networks. Transp Res B Methodol 99:44–61
[8] Chen Z, Lin X, Yin Y, Li M (2020) Path controlling of automated vehicles for system optimum on transportation networks with heterogeneous traffic stream. Transp Res C Emerg Technol 110:312–329
[9] Facchinei F, Pang J-S (2003) Finite-dimensional variational inequalities and complementarity problems. Springer, New York
[10] Fagnant DJ, Kockelman K (2015) Preparing a nation for autonomous vehicles: opportunities, barriers and policy recommendations. Transp Res A Policy Pract 77:167–181
[11] Gao G, Sun H (2014) Internalizing congestion and emissions externality on road networks with tradable credits. Procedia Soc Behav Sci 138:214–222
[12] He F, Wu D, Yin Y, Guan Y (2013) Optimal deployment of public charging stations for plug-in hybrid electric vehicles. Transp Res B Methodol 47:87–101
[13] Li R, Liu X, Nie YM (2018) Managing partially automated network traffic flow: efficiency vs. stability. Transp Res B Methodol 114:300–324
[14] Li Y, Chen Z, Yin Y, Peeta S (2020) Deployment of roadside units to overcome connectivity gap in transportation networks with mixed traffic. Transp Res C Emerg Technol 111:496–512
[15] Liu Z, Chen Z, He Y, Song Z (2021) Network user equilibrium problems with infrastructure-enabled autonomy. Transp Res B Methodol 154:207–241
[16] Lo HK, Szeto WY (2002) A methodology for sustainable traveler information services. Transp Res B Methodol 36(2):113–130
[17] Nagurney A (2000) Congested urban transportation networks and emission paradoxes. Transp Res D Transp Environ 5(2):145–151
[18] Nguyen S, Dupuis C (1984) An efficient method for computing traffic equilibria in networks with asymmetric transportation costs. Transp Sci 18(2):185–202
[19] Noruzoliaee M, Zou B, Liu Y (2018) Roads in transition: Integrated modeling of a manufacturer-traveler-infrastructure system in a mixed autonomous/human driving environment. Transp Res C Emerg Technol 90:307–333
[20] RGBSI (2021) 7 types of vehicle connectivity. https://blog.rgbsi.com/7-types-of-vehicle-connectivity#:~:text=Vehicle%20to%20cloud%20(V2C)%20communication,Redundancy%20to%20DSRC%20communication
[21] Sharon G, Albert M, Rambha T, Boyles S, Stone P (2018) Traffic optimization for a mixture of self-interested and compliant agents. In: Proceedings of the AAAI conference on artificial intelligence, vol 32
[22] Sheffi Y (1984) Urban Transportation Networks: Equilibrium Analysis with Mathematical Programming Methods. https://books.google.com.hk/books?id=zx1PAAAAMAAJ
[23] Spana S, Du L, Yin Y (2021) Strategic information perturbation for an online in-vehicle coordinated routing mechanism for connected vehicles under mixed-strategy congestion game. IEEE Trans Intell Transp Syst 23(5):4541–4555
[24] Tientrakool P, Ho Y-C, Maxemchuk NF (2011) Highway capacity benefits from using vehicle-to-vehicle communication and sensors for collision avoidance. In: 2011 IEEE vehicular technology conference (VTC Fall). IEEE, pp 1–5
[25] Wadud Z (2017) Fully automated vehicles: a cost of ownership analysis to inform early adoption. Transp Res A Policy Pract 101:163–176
[26] Wadud Z, MacKenzie D, Leiby P (2016) Help or hindrance? The travel, energy and carbon impacts of highly automated vehicles. Transp Res A Policy Pract 86:1–18
[27] Wallace CE, Courage K, Reaves D, Schoene G, Euler G (1984) Transyt-7f user’s manual. Technical report
[28] Wang M (2018) Infrastructure assisted adaptive driving to stabilise heterogeneous vehicle strings. Transp Res C Emerg Technol 91:276–295
[29] Wang X, Yang H, Zhu D, Li C (2012) Tradable travel credits for congestion management with heterogeneous users. Transp Res E Logist Transp Rev 48(2):426–437
[30] Xie T, Liu Y (2022) Impact of connected and autonomous vehicle technology on market penetration and route choices. Transp Res C Emerg Technol 139:103646
[31] Xiong B-K, Jiang R (2021) Speed advice for connected vehicles at an isolated signalized intersection in a mixed traffic flow considering stochasticity of human driven vehicles. IEEE Trans Intell Transp Syst
[32] Xu W, Fan M, Xu H (2016) Road pricing for a multi-modal transportation network with emission considerations. Syst Eng Theory Pract
[33] Yang H (1998) Multiple equilibrium behaviors and advanced traveler information systems with endogenous market penetration. Transp Res B Methodol 32(3):205–218
[34] Yao Z, Zhao B, Yuan T, Jiang H, Jiang Y (2020) Reducing gasoline consumption in mixed connected automated vehicles environment: a joint optimization framework for traffic signals and vehicle trajectory. J Clean Prod 265:121836
[35] Yin Y, Lawphongpanich S (2006) Internalizing emission externality on road networks. Transp Res D Transp Environ 11(4):292–301
[36] Yin Y, Yang H (2003) Simultaneous determination of the equilibrium market penetration and compliance rate of advanced traveler information systems. Transp Res A Policy Pract 37(2):165–181
[37] Zhang K, Nie YM (2018) Mitigating the impact of selfish routing: an optimal-ratio control scheme (ORCS) inspired by autonomous driving. Transp Res C Emerg Technol 87:75–90
[38] Zhang L, Lawphongpanich S, Yin Y (2009) An active-set algorithm for discrete network design problems. In: Transportation and traffic theory 2009: golden jubilee. Springer, pp 283–300
[39] Zhong R, Xu R, Sumalee A, Ou S, Chen Z (2020) Pricing environmental externality in traffic networks mixed with fuel vehicles and electric vehicles. IEEE Trans Intell Transp Syst 22(9):5535–5554

About this article

Cite this article

Liu, Y., Chen, Z., Gong, S. et al. Reducing congestion and emissions via roadside unit deployment under mixed traffic flow.Int J Coal Sci Technol 10, 1 (2023).
  • Received

    30 April 2022

  • Accepted

    21 October 2022

  • DOI

    https://doi.org/10.1007/s40789-022-00557-2

  • Share this article

    Copy to clipboard

For Authors

Explore