Next Article in Journal
A Lean Scheduling Framework for Underground Mines Based on Short Interval Control
Previous Article in Journal
Optimal Comfortable Load Schedule for Home Energy Management Including Photovoltaic and Battery Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Level Full Factorial Design Approach for the Analysis of Multi-Lane Highway Section under Saturated and Unsaturated Traffic Flow Conditions

1
Department of Civil Engineering, College of Engineering, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
2
College of Transportation Engineering, Tongji University, 4800 Cao’an Highway, Shanghai 201804, China
3
Department of Civil and Environmental Engineering, Faculty of Engineering—Rabigh Branch, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Department of Civil Engineering, College of Engineering, Al-Baha University, Al-Baha P.O.Box 1988, Saudi Arabia
5
Department of Architecture Design, College of Architecture and Planning, Sudan University of Science and Technology (SUST), P.O. Box 407, Khartoum 11111, Sudan
6
Department of Civil and Resource Engineering, Technical University of Kenya, Haile Sellasie Avenue, Nairobi P.O. Box 52428-00200, Kenya
7
Department of Civil and Construction Engineering, University of Nairobi, Harry Thuku Road, Nairobi P.O. Box 30197-00100, Kenya
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(12), 9194; https://doi.org/10.3390/su15129194
Submission received: 15 January 2023 / Revised: 28 March 2023 / Accepted: 31 March 2023 / Published: 7 June 2023

Abstract

:
Oversaturation of highways occurs due to their inadequate assessment and design. In this paper, we propose both a mathematical queuing model and a Discrete-Event Simulation (DES) framework based on Newell’s triangular flow-density relationship for the performance analysis of a multi-lane highway section. The proposed framework is a finite capacity queuing system, which captures an increase in the flow with the vehicle density up to the capacity of the section in an unsaturated condition and a decrease in the flow in the case of a saturated condition, depicting the actual traffic conditions on the highway section. First, the Birth–Death Process is used to build the mathematical queuing model (BDP), and the average number of vehicles (average queue length) and blocking probability on the highway section are estimated. Then, the accuracy of the mathematical queuing model is verified by the proposed DES framework. The “significance and effects” of different design factors are evaluated using the two-level full factorial design technique. The analysis of the experimental results reveals that the length of the highway section and the number of lanes are the most significant factors affecting the average queue length and blocking probability, while the jam density only has a significant effect on the average queue length and does not affect the blocking probability. In case of a two-way interaction, the combined effect of the “length-lanes” significantly affects the average queue length. In the end, a multiple-factor linear regression model is also developed for the prediction of the average number of vehicles on the highway section based on the design factors.

1. Introduction

Traffic congestion on a highway network is a condition that appears when the volume of traffic on a highway reaches or exceeds its capacity at a particular segment. The issue of traffic congestion is a real problem with the ever-increasing population and use of vehicles, and it draws the interest of academic researchers and traffic engineers. Due to long traffic jams and queues, the lack of productive time has an adverse impact on socioeconomic costs. Studying the interactions between cars, drivers, and the environment (roads, traffic control systems, etc.), with the goal of designing the ideal highway network and ensuring efficient traffic flow as well as minimal traffic congestion issues, is a highly challenging problem for researchers. The smooth movements of vehicles, smooth operation of traffic, and safety on the highways are the main concerns affecting the recognition of traffic and the highway system.
Scientific studies started in the 1930s to study traffic issues in order to understand and help in the prevention and resolution of traffic congestion problems. Several speed–density–flow models were developed over the years and categorized as single-regime or multi-regime models. Single-regime models assume a continuous relationship between speed, density, and traffic flow [1,2,3], while the relationship is discontinuous in multi-regime models depending on the level of density. Similarly, Lighthill and Whitham presented one of the most popular macroscopic traffic models, which was based on the fluid dynamic theory [4]. Treating the traffic flow as a 1-D compressible fluid, they studied the traffic jam as a shockwave. Prigogine presented the gas-kinetic model based on the Boltzmann equation [5]. Using the premise of a delayed adaptation of velocity, Newell proposed the microscopic, optimum velocity model [6].
Several queuing models were developed by different researchers for the analysis and design of transportation systems due to the inherent characteristics of urban transit systems, namely, the connection between the service facilities (walkways, stairs, ticket machines at transit stations, etc.) and the flow of entities (pedestrians) [7,8,9,10,11,12,13,14,15]. Similarly, Jain and MacGregor Smith [16] established state-dependent queuing models for simulating vehicle traffic flows and revealed that they are more realistic in extremely congested conditions. Afterwards, several other researchers carried out performance assessments of multi-lane highways based on the queuing theory and simulation models [17,18,19,20,21]. However, these queuing and simulation modeling models neglected the bottleneck and spillback phenomenon on highways and in transit terminal facilities. In the spillback effect, traffic congestion propagates upstream, which leads to the delay of and speed reduction by passengers as well as vehicular traffic, which is the true saturation condition. Therefore, in this research, we describe the realistic propagation of vehicular traffic congestion by the insertion of different phenomena such as queuing, the spillback effect from downstream traffic, and dissipation along the highway. Here, we develop both a mathematical queuing model and a simulation model, which incorporates Newell’s kinetic wave model [22]. Newell’s flow model is basically based on a triangular flow–density ( q K ) relationship, which is shown by the fundamental diagram in Figure 1.
This basic diagram highlights a few key features, which are the highway capacity ( C h ) based on the jam density (density occurring at zero speed) ( K j ) , number of lanes ( N ) and length of the highway ( L ) , critical density ( K c r ) (the density and speed experienced during peak operations), and maximum flow ( q max ) . Within this triangular relationship, there are two traffic conditions, i.e., a congested traffic flow condition and an uncongested traffic flow condition. In the case of the congested traffic condition ( K < K c r ) , the vehicles’ speed (wave pace) ( V ) is the slope of ( q K ) . When the traffic is not congested, low flow densities result in a monotonic increase in the traffic flow. The vehicles continue to move freely until the maximum flow is reached. At a critical density, the traffic state changes from free flow traffic to congested traffic. Beyond that point, in the case of the congested condition ( K > K c r ) , the space between vehicles becomes less as the density rises, which ultimately slows down the traffic. Since queues take up space, there is spillback in this case, and a bottleneck queue propagates upstream on the highway.
Based on the above analysis, our research is divided into three stages. First, we propose a mathematical queuing model for the vehicular traffic flow on a multi-lane highway section based on Newell’s triangular flow–density relationship to consider both the free-flow and congested conditions. The proposed queuing model is used to estimate the average number of vehicles on the multi-lane highway section. Secondly, for the verification of the accuracy of the proposed mathematical queuing model, a Discrete-Event Simulation (DES) framework is also developed in SimEvents®. Third, in order to assess the “significance and effect” of various design factors, a statistical two-level full factorial design approach is implemented in MINTAB software [23,24]. A regression model is also developed in the end to predict the average number of vehicles on the highway section based on various design factors.

2. Description of Multi-Lane Highway Section as a Queuing System

A section of highway is shown in Figure 2, in which vehicles occupy spaces as they enter the section. The available spaces on the highway section act as “servers”, and they become busy as vehicles enter the multi-lane highway segment.
The occupancy of these spaces causes an increase in the lane density ( K ) . When all the vacant spaces are occupied, the vehicles’ movement ceases, and the traffic flow stops. Jam density refers to the lane density associated with the jam condition on a highway segment ( K j ) . If a highway segment has a length ( L ) and number of lanes ( N ) , and the jam density ( K j ) of the section is known, the capacity ( C h ) of the multi-lane highway section can be calculated by using Equation (1).
C h = K j L N
Due to the limited capacity ( C h ) of the highway section, the number of empty spaces reduces as the number of vehicles ( n ) increases. The density ( K ) increases with the ultimate reduction in vehicular speed ( V ) . Thus, a multi-lane highway section is described as a state-dependent multi-server finite-capacity queuing system. The vehicular arrival rate can be described by X , and the state-dependent service rate of the multi-lane highway section based on the flow–density model is represented by Y ( n ) q - k . The available spaces on the multi-lane highway section act as servers. The capacity of the multi-lane highway section is the sum of all the servers, i.e., the number of available spaces. Therefore, the general form to represent a finite-capacity multi-lane highway section is a X / Y ( n ) q - k / C h / C h queuing system.

3. Mathematical Modeling Based on Birth–Death Process (BDP)

In this section, we discuss the mathematical queuing model’s formulation for a highway section in which the arrival process of vehicles is assumed to be based on Poisson’s distribution, and the state-dependent service process of the highway section is based on exponential distribution. The uni-directional flow of vehicular traffic is considered in this research. To formulate the M / M ( n ) q - k / C h / C h model for the performance analysis of the highway section based on the state-dependent Newell’s flow–density model, we employ the Birth–Death (BD) Process [25].
In the case of Newell’s flow–density model, there are two conditions, i.e., a saturated condition and an unsaturated condition. Therefore, it is necessary to determine a transition point that demarcates the unsaturated condition from the saturated condition. Let n c r be the traffic volume below which the condition is unsaturated ( n < n c r ) and above which the condition is saturated ( n n c r ) . For a finite capacity M / M ( n ) q - k / C h / C h queuing system, the state transition diagram is shown in Figure 3.
The vertical dotted line demarcates the saturated condition and the unsaturated condition. An infinitesimal generator matrix with birth rate λ and death rate μ ( n ) can be created from the state transition diagram, based on the BDP, which is given by Equation (2).
( λ 0 λ 0 μ 1 ( λ 1 + μ 1 ) λ 1 2 μ 1 ( λ 2 + 2 μ 2 ) λ 3 3 μ 1 ( λ 3 + 3 μ 3 ) λ 4 4 μ 1 ( λ 4 + 4 μ 4 ) λ n n μ 1 ( λ n + n μ n ) λ c C μ 1 ( λ c + C μ c ) )
Let π be the state probability, i.e., the number of vehicles on the highway section.
π Q = 0
Using Equations (2) and (3), the steady-state (equilibrium) equations can be obtained as
λ 0 π 0 + μ 1 π 1 = 0
π 1 = ( λ 0 / μ 1 ) π 0
λ i 1 π i 1 ( λ i + μ i ) π i + μ i + 1 π i + 1 = 0
π i + 1 = ( λ i + μ i μ i + 1 ) π i ( λ i 1 μ i + 1 ) π i 1
Using the principle equation, making appropriate substitutions, and using the mean of recursion, all the state probabilities can be expressed in terms of zero-state probability. The zero-state probability can be obtained by the normalization condition, i.e., n = 0 C π n = 1 .
π 0 = [ 1 + i = 1 C λ n ( i ! ) i = 1 C u i ] 1
The vehicles’ flow rate on any highway section is related to the speed of the vehicles on the section, the number of vehicles on the section, and the length of the section, which is given by Equation (9):
μ = n V L
However, there are two conditions, i.e., the unsaturated and saturated traffic conditions. The traffic flow rate can be expressed in terms of both the unsaturated and saturated flow rates, as shown by Equation (10).
μ n = { n V L N n < n c r w ( C n ) L N n n c r
With the substitution shown by Equation (10), we obtain all the state probabilities for both the saturated and unsaturated conditions, as shown by Equation (11).
π n = { λ n n = 1 n c r 1 u n π o n < n c r λ n n = 1 n c r 1 u n n = n c r C u n π o n > n c r
In order to obtain the performance measures, such as the blocking probability and the average number of vehicles on the highway section, we can use Equations (12) and (13), which are based on the state probabilities from Equation (11).
π C = λ C n = 1 n c r 1 u n n = n c r C u n π o
E [ N ] = i = 1 C i π i 1

4. Discrete-Event Simulation Architecture of Multi-Lane Highway Section

To ensure that our suggested queuing model is accurate for the multi-lane highway section, the DES model is developed in the SimEvents® toolbox of the MATLAB programming environment, as shown in Figure 4.
The DES model consists of two phases, i.e., (1) average vehicles’ arrival phase; (2) flow–density based on the service phase of the highway. The major blocksets of the SimEvents® toolbox for the DES model are the FIFO_Queue blocks, which assign queueing spaces to the vehicles, the Server blocks, which temporarily store the passengers, and the Start and Read Timers, which display the average dwelling time of the vehicles on the highway section. The Level-2 S-function block computes and updates the state-dependent speed of the vehicles as the number of vehicles increases on the section. The Constant blocks are used to input values for various parameters in the DES model, such as arrival rate, length of the section, number of lanes, etc. The Display block displays the output, i.e., average number of vehicles on the highway section.
As discussed earlier, the highway section has a finite capacity (shown in Figure 2) in which the arriving demand cannot overcome the overall capacity of the highway section. Therefore, during the vehicles’ arriving phases, it must be assured that the number of arrivals stays less than or equal to the capacity ( C h = K j L N ) of the highway section. To execute this scenario in the DES environment, the vehicles’ arrival from the upstream section stays in the FIFO_Queue block before entering the Server block. The vehicles are subsequently sent to the downstream section. The MATLAB® Function block is used to perform the following functions:
  • It assesses the number of vehicles arriving at the highway section with its capacity ( C h = K j L N ) . This block prevents the entry of entities into the FIFO_Queue block when the highway section is at capacity n = C h .
  • This block computes the blocking probability as the highway jams, i.e., operates at capacity. When the entrance to the FIFO_Queue block is blocked for vehicles, the blocked vehicles are simultaneously activated and registered through the output switch block’s second entity port (OUT2). Based on the blocking probability, the average number of vehicles is estimated from the FIFO_Queue block of the DES model.

5. Computational Experiments

In this section, first the accuracy of the mathematical queuing model is verified by the proposed simulation model. The two-level full factorial design approach is used to statistically analyze the effect of various factors on the average number of vehicles on the highway section.

5.1. Verification of Proposed M / M ( n ) q - k / C h / C h Model

Before conducting the two-level full factorial design of the experiment for analyzing the significance of various factors, it is necessary to first verify the accuracy proposed by M / M ( n ) q - k / C h / C h . For verification purposes, we perform the experiments on a simple highway section, where the number of lanes is equal to one and two and the length of the section is equal to 0.10 mi and 0.2 mi with different arrival rates, which are according to the actual conditions on the collector road sections in Islamabad, Pakistan. The average number of vehicles (average queue length) ( E N ) on the highway section, blocking probability (PC), and average vehicles’ dwelling time (EW) are obtained from both the mathematical queuing and DES models. The results are presented in Table 1.
During simulation experiments, the mean values of 30 replications are recorded and tabulated. Each simulation test is performed with a simulation time of 50,000 time units to achieve outputs at steady-state conditions. The 95% confidence interval approximations for the distribution of the mean, which follow the normal distribution, are also recorded in Table 1. All the computational experiments are carried out on a PC with Intel® Core™ i5-4570 CPU@ 3.20 GHz and 8 GB of RAM under a Windows® operating system. The mean computational times (CPU times in minutes) is also recorded for all the experiments. As shown in Table 1, both the mathematical and simulation models show a clear consistency. Therefore, the models can be used for the performance assessment of a multi-lane highway section under saturated and unsaturated conditions.

5.2. Factorial Design Approach

The two-level full factorial design approach is conducted to evaluate the significance and effect of four different factors such as the length of the highway section, the number of lanes of the highway section, the average vehicles’ arrival rate, and the jam density on the average number of vehicles (mean queue length) ( E N ) on the highway section. The average number of vehicles on the highway section is obtained from our proposed mathematical queuing and simulation models. Table 2 displays the four factors along with their low and high levels and units. The numbers in parenthesis (−1 or +1) refer to the coded values for the factors under study at their low (−1) and high (+1) levels, respectively.
Table 3 displays the main effect estimates as well as the interaction effect estimates. The difference between the average response of a factor at a high level and that of a factor at a low level is referred to as the effect. The interaction effect between two factors (say, A*B) is defined as the mean difference between the effects of the “length of the highway section” at a high level for the “number of lanes” and the effect of the “length of the highway section” at a low level for the “number of lanes”.
MINITAB-19 software was used to calculate the full model regression coefficients in terms of their coded values [14]. The effect estimates of the factors can be obtained by multiplying the regression coefficient by a factor of 2. Table 4 shows the analysis of variance (ANOVA) of up to two factors’ interactions. An ANOVA with higher order interactions of three or four variables results in an ill model, so it is deserted due to the sparsity of effect principle.
It asserts that the majority of processes are regulated by a few major factors and a few low-order interactions.

5.2.1. Effects of the Factors

From Figure 5, it can be observed that the effect estimates of the length of the highway section, number of lanes, and jam density are positive (positive slope), which means that they have a positive effect on the average number of vehicles (average queue length) on the highway section.
The length plot shows a relatively high positive slope, but the jam density plot has a somewhat positive slope. This suggests that increasing the amounts of these parameters reduces the average queue length. This shows that with a slight increase in the length of the highway section, i.e., from 0.1 mi to 0.2 mi, the average queue length also increases, when keeping all other factors constant.
Similarly, the average vehicles’ arrival rate is negative, which means that it has a negative impact on the queue length and the increase in this factor causes an increase in the average queue length, which is quite obvious in all cases. A cumulative normal probability plot of the main and interaction effects of the factors at α = 0.05 is shown in Figure 6.
It is obvious that none of the components are in a straight line. Outlier variables, shown by dark red squares, mostly concern estimating the average queue length and are thought to be the most dominating elements. The factors that are insignificant tend to fall along a straight line on this plot, whereas the significant main effects and interaction have nonzero means and do not lie along the straight line. We can note that the average vehicle arrival rate is not a potential outlier; consequently, the influence of the average vehicle arrival rate explored in this study (for this specific highway section) is not significant compared to the length of the highway section, which is a potential outlier. This could be explained by the fact that the average vehicle arrival rate range (1500–2000 veh/h) covered in this study is rather small. The significance of factors can also be assessed by using a Pareto chart, as shown in Figure 7.
The factors with a horizontal bar extended beyond the dashed vertical line illustrate the significance of the factors.

5.2.2. Interaction Effects of the Factors

An interaction is the failure of one factor to produce the same effect on the response at the different levels of another factor. The interaction between two factors is said to occur when a change in the values of one variable alters the effect on another factor. This implies that insignificant factor interactions produce similar trends in response to the different levels of another factor. The two-way interaction (two-factor interaction) effect is shown in Figure 8, and an insignificant two-factor interaction (B*D, p-value = 0.657) and one significant two-factor interaction (A*B, p-value = 0.012) at α = 0.05 are presented.
It is clear that when the length of the highway section is increased, while keeping the number of lanes at a low level, the average queue length (average number of vehicles) increases. Similarly, with an increase in the length of the highway section, while keeping the number of lanes at a higher level, the queue length again increases.

5.2.3. Regression Model

In order to predict future observations, the regression analysis shows the statistical relationship between one or more independent factors and the response factor. In this research, by considering only the main (linear) and interaction terms (two-factor interaction), a multiple-factor linear regression model is obtained for the prediction of the average vehicles’ queue length on the highway section, which is given by the following equation:
EN = 22.668 + 8.742 Length + 7.221 Number of Lanes − 0.682 Arrival Rate + 2.578 Jam Density + 3.309 Length × Number of Lanes + 0.704 Length × Arrival Rate + 1.924 Length × Jam Density − 0.817 Number of Lanes × Arrival Rate+ 0.403 Number of Lanes × Jam Density − 0.727 Arrival Rate × Jam Density
The coefficient of the determination R 2 of the regression model is 0.976, and the adjusted R 2 is 0.930. R 2 = 93% in the model shows that more than 93% of the variability of the responses can be explained by the factors of regression model.

6. Conclusions and Recommendations

In this research, we proposed a mathematical queuing and simulation model for the performance evaluation of a multi-lane highway section under congested and uncontested conditions using Newell’s triangular flow–density relationship. This model is more realistic as it considered different traffic flow conditions. The two-level-full-factorial design approach is used to determine the significant and insignificant design parameters such as the length of the highway section, number of lanes, average vehicles’ arrival rate, and jam density. It was observed that the length of the highway section, number of lanes, and jam density have a positive effect, which shows that increasing the level of these factors results in an increased value for the average queue length. The jam density has a mild positive slope, which shows that it has a slight significance compared to the other highly significant factors. Additionally, it was found that the majority of the highest-order interactions are insignificant, which supports the sparsity of effect principle. Only the two-way interaction of the “length-lanes” is a significant interaction affecting the average queue length on the highway section. A regression model to predict the average queue length from the two-level full factorial design is developed, which takes into account various design factors and their interactions.
The effects of the four factors studied are limited to a specific range of values. More research with a broader range at both low and high levels is recommended. Other different sites can be considered for the performance evaluation based on the full factorial design approach. Additionally, the proposed model could be modified slightly and used for the analysis of freeways. With a few minor adjustments, the proposed model may also be used to analyze passenger flow in airport and metro station corridors.

Author Contributions

Conceptualization, A.K. and C.M.M.; formal analysis, A.K. and H.A.; investigation, S.A.; methodology, A.K. and H.A.; project administration, R.A.; resources, H.A. and A.E.; software, H.A. and A.A.; supervision, S.A.; visualization, R.A. and A.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicaple.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

The researchers would like to acknowledge the Deanship of Scientific Research, Taif University, for funding this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Greenshields, B.D. A study in highway capacity. Highw. Res. Board Proc. 1935, 14, 448–477. [Google Scholar]
  2. Greenberg, H. An analysis of traffic flow. Oper. Res. 1959, 7, 79–85. [Google Scholar] [CrossRef]
  3. Underwood, R.T. Speed, volume, and density relationship: Quality and theory of traffic flow, yale bureau of highway traffic. New Haven Conn. 2004, 141–188. [Google Scholar]
  4. Lighthill, M.J.; Whitham, G.B. On kinematic waves II. A theory of traffic flow on long crowded roads. Proc. R. Soc. London. Ser. A Math. Phys. Sci. 1955, 229, 317–345. [Google Scholar] [CrossRef]
  5. Prigogine, I. A Boltzmann-like approach to the statistical theory of traffic flow. In Theory of Traffic Flow; Elsevier Publishing: Amsterdam, The Netherlands, 1961. [Google Scholar]
  6. Newell, G.F. Nonlinear effects in the dynamics of car following. Oper. Res. 1961, 9, 209–229. [Google Scholar] [CrossRef]
  7. Khattak, A.; Yangsheng, J. Modeling of subway stations circulation facilities as state-dependent queuing network based on phase-type distribution. In Proceedings of the 2016 IEEE International Conference on Intelligent Transportation Engineering (ICITE), Singapore, 20–22 August 2016; pp. 133–138. [Google Scholar]
  8. Jiang, Y.; Zhu, J.; Hu, L.; Lin, X.; Khattak, A. AG/G(n)/C/C state-dependent simulation model for metro station corridor width design. J. Adv. Transp. 2016, 50, 273–295. [Google Scholar] [CrossRef]
  9. Zhu, J.; Hu, L.; Jiang, Y.; Khattak, A. Circulation network design for urban rail transit station using a PH(n)/PH(n)/C/C queuing network model. Eur. J. Oper. Res. 2017, 260, 1043–1068. [Google Scholar] [CrossRef]
  10. Ni, W.; Yu, J.; Cai, H.; Bai, M.; Wu, B. Study on Optimization of Passenger Flow at a Metro Station Based on AnyLogic—Case Study of Youfangqiao Station of Nanjing Metro Line 2. Complex Syst. Model. Simul. 2021, 1, 242–252. [Google Scholar] [CrossRef]
  11. Hu, L.; Jiang, Y.; Zhu, J.; Chen, Y. A PH/PH(n)/C/C state-dependent queuing model for metro station corridor width design. Eur. J. Oper. Res. 2015, 240, 109–126. [Google Scholar] [CrossRef]
  12. Jiang, Y.; Lin, X. Simulation and optimization of the ticket vending machine configuration in metro stations based on anylogic software. In ICTE 2013: Safety, Speediness, Intelligence, Low-Carbon, Innovation, Proceedings of the Fourth International Conference on Transportation Engineering, Chengdu, China, 19–20 October 2013; ASCE: Reston, VA, USA, 2013; pp. 754–760. [Google Scholar]
  13. Khattak, A.; Yangsheng, J.; Abid, M.M. Optimal configuration of the metro rail transit station service facilities by integrated simulation-optimization method using passengers’ flow fluctuation. Arab. J. Sci. Eng. 2018, 43, 5499–5516. [Google Scholar] [CrossRef]
  14. Xu, X.Y.; Liu, J.; Li, H.Y.; Hu, J.Q. Analysis of subway station capacity with the use of queueing theory. Transp. Res. Part C Emerg. Technol. 2014, 38, 28–43. [Google Scholar] [CrossRef]
  15. Jiang, Y.; Luo, N. Urban rail transit station automatic fare collection configuration research and simulation. In ICTE 2013: Safety, Speediness, Intelligence, Low-Carbon, Innovation, Proceedings of the Fourth International Conference on Transportation Engineering, Chengdu, China, 19–20 October 2013; ASCE: Reston, VA, USA, 2013; pp. 1142–1149. [Google Scholar]
  16. Jain, R.; Smith, J.M. Modeling vehicular traffic flow using M/G/C/C state dependent queueing models. Transp. Sci. 1997, 31, 324–336. [Google Scholar] [CrossRef]
  17. Heidemann, D. A queueing theory approach to speed-flow-density relationships. In Transportation and Traffic Theory, Proceedings of the 13th International Symposium on Transportation and Traffic Theory, Lyon, France, 24–26 July 1996; National Cademies: Washington, DC, USA, 1996. [Google Scholar]
  18. Vandaele, N.; Van Woensel, T.; Verbruggen, A. A queueing based traffic flow model. Transp. Res. Part D Transp. Environ. 2000, 5, 121–135. [Google Scholar] [CrossRef] [Green Version]
  19. Khattak, A.; Hussain, A.; Ibrahim, F. A PHD-DES Framework for the Performance Assessment of Multi-lane Highways Under Random Traffic Flow. Arab. J. Sci. Eng. 2019, 44, 8615–8628. [Google Scholar] [CrossRef]
  20. Khattak, A. Level of Service of Multi-Lane highway Section under fluctuated traffic flow condition: A state-dependent Discrete-Event Simulation framework. In Proceedings of the 2019 5th International Conference on Transportation Information and Safety (ICTIS), Liverpool, UK, 14–17 July 2019; pp. 1245–1250. [Google Scholar] [CrossRef]
  21. Van Woensel, T.; Vandaele, N. Empirical validation of a queueing approach to uninterrupted traffic flows. 4OR 2006, 4, 59–72. [Google Scholar] [CrossRef]
  22. Newell, G.F. A simplified theory of kinematic waves in highway traffic, part I: General theory. Transp. Res. Part B Methodol. 1993, 27, 281–287. [Google Scholar] [CrossRef]
  23. Ryan, T.P.; Morgan, J.P. Modern experimental design. J. Stat. Theory Pract. 2007, 1, 501–506. [Google Scholar] [CrossRef]
  24. Minitab, Inc. Minitab 16 Statistical Software; Computer Software; Minitab, Inc.: State College, PA, USA, 2010; Available online: www.minitab.com (accessed on 3 December 2021).
  25. Stewart, W.J. Probability, Markov Chains, Queues, and Simulation: The Mathematical Basis of Performance Modeling; Princeton University Press: Princeton, NJ, USA, 2009. [Google Scholar]
Figure 1. Newell’s flow–density relationship.
Figure 1. Newell’s flow–density relationship.
Sustainability 15 09194 g001
Figure 2. Queuing system representation of multi-lane highway section.
Figure 2. Queuing system representation of multi-lane highway section.
Sustainability 15 09194 g002
Figure 3. State transition diagram for M / M ( n ) q - k / C h / C h queuing model.
Figure 3. State transition diagram for M / M ( n ) q - k / C h / C h queuing model.
Sustainability 15 09194 g003
Figure 4. Discrete-Event Simulation architecture of multi-lane highway section.
Figure 4. Discrete-Event Simulation architecture of multi-lane highway section.
Sustainability 15 09194 g004
Figure 5. Main effects plot for mean number of vehicles.
Figure 5. Main effects plot for mean number of vehicles.
Sustainability 15 09194 g005
Figure 6. Normal plot of the standardized effects.
Figure 6. Normal plot of the standardized effects.
Sustainability 15 09194 g006
Figure 7. Pareto chart of the standardized effects.
Figure 7. Pareto chart of the standardized effects.
Sustainability 15 09194 g007
Figure 8. Interaction plot for EN.
Figure 8. Interaction plot for EN.
Sustainability 15 09194 g008
Table 1. Verification of proposed mathematical queuing model.
Table 1. Verification of proposed mathematical queuing model.
Length: 0.2 mi
Lanes: 1
Jam Density: 80 veh/mi/ln
Speed: 60 mi/h
Length: 0.1 mi
Lanes: 1
Jam Density: 90 veh/mi/ln
Speed: 60 mi/h
Length: 0.2 mi
Lanes: 2
Jam Density: 100 veh/mi/ln
Speed: 90 mi/h
Proposed
M / M ( n ) q - k / C h / C h
Model
Proposed
DES Model
(95% Confidence Interval)
Proposed
M / M ( n ) q - k / C h / C h
Model
Proposed
DES Model
(95% Confidence Interval)
Proposed
M / M ( n ) q - k / C h / C h
Model
Proposed
DES Model
(95% Confidence Interval)
λ = 2100 veh/h λ = 1525 veh/h λ = 1720 veh/h
EN
(veh)
15.8516.33
(15.98, 16.85)
7.256.89
(6.11, 7.22)
39.9237.97
(36.61,38.22)
PC0.880.862
(0.853, 0.871)
0.540.52
(0.511, 6.04)
0.930.93
(0.923, 0.947)
EW
(s)
28.9130.13
(29.11, 31.87)
17.1216.30
(15.98, 17.25)
83.5483.54
(81.97, 84.88)
Table 2. Levels of factors studied.
Table 2. Levels of factors studied.
FactorsAbb.UnitsLow LevelHigh Level
LengthAMiles0.1 (−1)0.2 (+1)
Number of LanesB-1 (−1)2 (+1)
Vehicles’ Arrival RateCVeh/hr1500 (−1)2000 (+1)
Jam DensityDVeh/mile/lane90 (−1)120 (+1)
Table 3. Coefficients and main and interaction effects’ estimates.
Table 3. Coefficients and main and interaction effects’ estimates.
TermEffectCoefSE CoefT-Valuep-ValueVIF
1-Way Interaction
Constant 22.6680.85526.500.000
  • Length
17.4848.7420.85510.220.0001.00
  • Number of Lanes
14.4417.2210.8558.440.0001.00
  • Arrival Rate
−1.364−0.6820.855−0.800.4611.00
  • Jam Density
5.1562.5780.8553.010.0301.00
2-Way Interaction
  • Length*Number of Lanes
6.6193.3090.8553.870.0121.00
  • Length*Arrival Rate
1.4090.7040.8550.820.4481.00
  • Length*Jam Density
3.8491.9240.8552.250.0741.00
  • Number of Lanes*Arrival Rate
−1.634−0.8170.855−0.960.3831.00
  • Number of Lanes*Jam Density
0.8060.4030.8550.470.6571.00
  • Arrival Rate*Jam Density
−1.454−0.7270.855−0.850.4341.00
Table 4. Analysis of variance.
Table 4. Analysis of variance.
SourceDFAdj SSAdj MSF-Valuep-Value
MODEL102434.86243.4920.800.002
One-Way Interaction42170.71542.6846.370.000
  • Length
11222.731222.73104.470.000
  • Number of Lanes
1834.20834.2071.280.000
  • Arrival Rate
17.447.440.640.461
  • Jam Density
1106.35106.359.090.030
Two-Way Interactions6264.1544.033.760.084
  • Length*Number of Lanes
1175.23175.2314.970.012
  • Length*Arrival Rate
17.947.940.680.448
  • Length*Jam Density
159.2559.255.060.074
  • Number of Lanes*Arrival Rate
110.6810.680.910.383
  • Number of Lanes*Jam Density
12.602.600.220.657
  • Arrival Rate*Jam Density
18.458.450.720.434
  • Error
558.5211.70
Total152493.38
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Almujibah, H.; Khattak, A.; Alotaibi, S.; Alahmadi, R.; Elhassan, A.; Alshahri, A.; Matara, C.M. Two-Level Full Factorial Design Approach for the Analysis of Multi-Lane Highway Section under Saturated and Unsaturated Traffic Flow Conditions. Sustainability 2023, 15, 9194. https://doi.org/10.3390/su15129194

AMA Style

Almujibah H, Khattak A, Alotaibi S, Alahmadi R, Elhassan A, Alshahri A, Matara CM. Two-Level Full Factorial Design Approach for the Analysis of Multi-Lane Highway Section under Saturated and Unsaturated Traffic Flow Conditions. Sustainability. 2023; 15(12):9194. https://doi.org/10.3390/su15129194

Chicago/Turabian Style

Almujibah, Hamad, Afaq Khattak, Saleh Alotaibi, Raed Alahmadi, Adil Elhassan, Abdullah Alshahri, and Caroline Mongina Matara. 2023. "Two-Level Full Factorial Design Approach for the Analysis of Multi-Lane Highway Section under Saturated and Unsaturated Traffic Flow Conditions" Sustainability 15, no. 12: 9194. https://doi.org/10.3390/su15129194

APA Style

Almujibah, H., Khattak, A., Alotaibi, S., Alahmadi, R., Elhassan, A., Alshahri, A., & Matara, C. M. (2023). Two-Level Full Factorial Design Approach for the Analysis of Multi-Lane Highway Section under Saturated and Unsaturated Traffic Flow Conditions. Sustainability, 15(12), 9194. https://doi.org/10.3390/su15129194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop