Multimetric Analysis of a Simulated Mixed Trafﬁc of Motorcycles and Automobiles: Flow, Energy, CO 2 , and Costs

The fleet of developing countries consists of motorcycles and cars. This heterogeneous traffic condition has its advantages and disadvantages, which results in conflicting points of view ( e.g. , motorcyclists enjoying a higher mobility while car drivers resent their decreased speed). In this paper, we corroborated the notion that traffic evaluation depends on the chosen metric ( e.g. , vehicle flow, fuel consumption, monthly costs) and the point of view (driver, rider, and policy makers). To this effect, we studied a mixed traffic condition, considering that the vehicle performance is affected by three scales: engine, vehicle, and traffic. We modeled the engine using empirical correlations of power and energy efficiency, the vehicle based on a balance of propulsive and resistive forces, and traffic with a cellular automata model. We simulated 189 traffic conditions and evaluated vehicle flow, average energy consumption, total CO 2 emission of the road, and monthly costs. We also discussed the results from the point of view of the driver, rider, and society. We concluded that the optimal condition depends both on the choice of metric and point of view, and that is not appropriate to use results from homogeneous traffic to analyze heterogeneous traffic conditions, even if both scenarios present the same total vehicle flow.


Introduction
Most of the fleets of developing countries such as Brazil, Colombia, Malaysia, Pakistan, Thailand, Philippines, and Nigeria consist of cars and motorcycles (WHO, 2015). In those countries, the use of motorcycles is motivated by the poor quality of transit, traffic congestion, and their advantageous cost of acquisition and operation, together with a higher average speed. On the other hand, motorcycle accidents are accountable for economic and social impacts (Koossalapeerom et al., 2019). The conflict between advantages and disadvantages addresses the relevance of the metric analyzed when evaluating traffic, which can be an individual one, such as travel time and direct costs; or a collective one, such as indirect costs with road infrastructure and accidents, flow of vehicles, fuel production, and CO 2 emissions (Tranter, 2012).
Traffic is a complex non-linear system affected by internal and external factors (Zegeye, De Schutter, Hellendoorn, Breunesse, and Hegyi, 2013) and it depends on instantaneous conditions such as vehicle density and the ratio of different types of vehicles. Homogeneous traffic flow models are unable to faithfully represent mixed traffic conditions, as it happens in developing countries (Agarwal, Zilske, Rao, and Nagel, 2015), in which motorcycles and passenger cars constitute the fleet and are accountable for significant emissions in major cities (Arun, Mahesh, Ramadurai, and Shiva Nagendra, 2017). Several papers claim that total vehicle flow varies for the same total density when mixed traffic occurs, increasing when there are motorcycles on the track, and decreasing when there are trucks, thus reflecting dynamic characteristics (Yang, Qiu, Yu, Sun, and Pu, 2015).
Furthermore, traffic cannot be evaluated only by its vehicle flow. Other metrics such as air pollution, CO 2 emission, direct costs, congestion, cost of changing the infrastructure, and accidents should also be taken into account (Gössling and Choi, 2015). Individual vehicle owners tend to purchase their vehicle based on flexibility and convenience. For motorcycles, low purchase and running costs are essential for the lower income population who acquire them (Koossalapeerom et al., 2019). There is also a growing concern about the global climate change, and CO 2 emissions from the transportation sector is a frequently addressed subject (Ehsani, Ahmadi, and Fadai, 2016).
This study contributes to the literature with the evaluation and discussion of different metrics such as flow, energy, CO 2, and costs, from the point of view of different stakeholders: drivers, motorcycle riders, and policy makers. Although the results from individual metrics are present in the literature, there is a lack of studies integrating the models and discussing the results. In order to evaluate those metrics, we applied a method that considers three scales: engine, vehicle, and traffic. We studied the engine using an empirical correlation of power and efficiency, the vehicle was studied modeling the propulsive and resistive forces, and the traffic with a cellular automata (CA) model, in which we simulated the traffic for 189 distinct density combinations of cars and motorcycles.

Methodology
The dynamic behavior and energy expenditure of vehicles depend on the engine, powertrain, chassis, traffic, driver/rider driving style, road/weather condition, and traffic laws. In this study, vehicle performance considers the engine, vehicle, and traffic. Other aspects are simplified, such as the driver/rider always trying to reach the maximum speed, the absence of traffic lights and intersections, and the road being a one-way straight line (Figure 1). This track model is based on the one proposed by Meng, Dai, Dong, and Zhang (2007), and it represents a single-track road for cars, with enough space for motorcycles to split and filter through them.
The engine is modeled by two empirical correlations. The cubic curve of the maximal power (P max ,eng ), Equation (1), is given by Ni and Henclewood (2008), and it is a function of engine speed (Ω). This model needs the following inputs: peak power, engine speed at peak power (Ω p ) and engine speed at peak torque (Ω t ).
The efficiency map in Equation (2) is described by Ben-Chaim, Shmerling, and Kuperman (2013), and it reports engine efficiency (η eng ) based on the engine's maximum efficiency (η eng,m ), as well as considering factors related to engine speed (µ r ) and throttle usage (µ a ). A detailed explanation regarding the integration of both models can be found in Andrade, Araújo, Santos, and Magnani (2020).
η eng (Ω,α) = n eng,m µ r µ a (2) Table 1 lists the parameters employed in this study to represent cars and motorcycles, reproducing characteristics of the usual Brazilian models. In average, Brazilian cars are lighter and less powerful than American and European cars.  The dynamic behavior of the vehicle is based on the balance of propulsive and resistive forces (Equation 3) (Cossalter, 2006;Jazar, 2014). On the left-hand side of the equation is the inertia, where m is the total mass (vehicle + driver), V is vehicle speed, and t is time. Tractive force considers αP max ,eng as the actual power of the engine, and η trans is the transmission efficiency. Throttle usage, α, is the ratio between the actual power and the maximum possible power at a specific engine speed. The braking force has β as the braking factor, µ as the friction coefficient between tire and road, and θ as the road grade. The aerodynamic drag force considers the constant k A (i.e., 1/2ρAC D ) with ρ as fluid density (air), A as the frontal area, C D as the vehicle shape, and W is the wind speed. In the equation, rolling resistance, C R , is also considered, as well as the gravitational resistance (mg), both related to road grade.
The finite difference method is used to solve the balance of forces in Equation (3). For every time step t, this Equation is employed to calculate factors α (throttling) and β (braking) to obtain the desired speed (restricted by the traffic condition, i.e., the speed and position of nearby vehicles, as well as pilot preferences, who may always try to reach the maximum allowed speed as soon as it is physically possible). Thus, those calculated α and β values (limited to the range [0-1]) are reused in Equation (3) to determine the actual instantaneous speed V (in opposition to the original desired speed). As α and β are limited (modelling the limits of the acceleration and braking powers), the actual vehicle speed will not always be equal to the desired one.
With the evolution of α (related to the throttle opening) during the vehicle movement, it is possible to calculate the vehicle energy expenditure during a time interval ∆t using Equation (4): Modelling the vehicle (Equations 3 and 4) depends on the engine model and the desired speed, which must result from the traffic model. In this study, the employed traffic model is a cellular automata (CA) model based on Meng et al. (2007). This microscopic traffic model considers the individual movement and interactions of each vehicle in heterogeneous traffic conditions. The CA model is discrete for position and time, but, even with its simplicity, it is able to reproduce complex traffic phenomena, e.g. the start-stop waves that appear in congested traffic. The CA model is also known as the NaSch model, an abbreviation from the researchers that first employed the model: Nagel and Schreckenberg, in 1992. This model was successfully improved and is still employed by several researchers (Chen and Wang, 2016;Hua, Yue, Wei, Chen, and Wang, 2020;Lan and Chang, 2005;Lan, Chiou, Lin, and Hsu, 2010;Lv, Song, Liu, and Ma, 2013; (Sean) Qian, Li, Li, Zhang, and Wang, 2017;Ruan, Zhou, Tu, Jin, and Shi, 2017;Yang, Qiu, Yu, Sun, and Pu, 2015;Zeng et al., 2021).
The modeled track has 2 000 cells with a length of 3,75 m each, totaling 7,5 km. Cars have 11 speeds from 0 to 135 km/h, and motorcycles have 5, from 0 to 54 km/h. Cars can travel on the left side of the lane and motorcycles can travel on both sides; therefore, motorcycles can split the lane with cars and other motorcycles. Additionally, this model attempts to represent the motorcycle's notable ability to weave through traffic at lower speeds ( Figure 1).
The traffic simulation has five steps, similar to the one proposed by Meng et al. (2007). Those steps are proposed both to manage the lateral and longitudinal movement and to provide safety for the vehicles, avoiding collisions. Table  2 lists the parameters and variables employed during the procedure.  To summarize, motorcycles decide if it is worth changing sides depending on traffic conditions. In the second and third steps, all vehicles accelerate or decelerate respecting the non-collision rule and their maximum speed. In the fourth step, a braking probability is introduced for all vehicles (p = 0,10), in an attempt to represent driver's distraction in real traffic. Without this step, all vehicles would move together with the same speed, as train wagons. We consider that the vehicles leaving the track return to the beginning to maintain constant densities (periodical boundary condition). In the final step, all vehicles move as previously evaluated.
Step 1 Step 2: Acceleration of the vehicles Step 3: Deceleration of the vehicles Step 4: Randomization (with the probability p) Step 5: Movement of the vehicles where n (m,c) is the number of motorcycles (n (m) ) or cars (n (c) ), and L sim is the length of the simulated track (7 500 m). A double average is applied to decrease the influence of the initial condition, one on time and other between simulations. Similarly, the average speed of the vehicles, V (m,c) [m/s], Equation (6), is given by the double average among all the vehicles on the track in the last 1 000 seconds of 30 simulations.
The flow of occupants for each type of vehicle, Q (m) [motorcycle occupant/s] and Q (c) [car occupant/s] is given by Equation (7), as the product of the occupancy rate o (m,c) [occupant/vehicle], the vehicle density ρ (m,c) , and the average speed V (m,c) . In this study, the occupancy rate for cars and motorcycles is considered equal to 1.
The total flow of occupants on the track Θ (occupant/s), Equation (8), is the sum of the flow of motorcycle and car occupants.
In this study, new metrics are introduced in relation to (Meng et al., 2007), which is possible with the use of the engine and vehicle models. The average consumed energy by vehicle by distance, E (m,c) [J/vehicle·m] is given by Equation (9).
The total consumed energy, by time and space, Ξ (J/m·s) is given by the sum of the consumption of motorcycles and cars, as provided by Equation (10). From the energy consumption, it is possible to estimate the CO 2 emissions since they are proportional. We consider the combustion of each MJ of C 8 H 18 to emit 0,07 kg of CO 2 .  Table 3 lists the Brazilian financial parameters employed in the calculation of monthly costs, C (m,c) (US$/occupant·month), Equation (11). In this Equation, the first term on the righthand side represents the monthly value relative to the vehicle acquisition, V p . To compare monthly equivalent values, the corrected number of months µ is employed, which is calculated by Equation (12) (Stoecker, 1989) and considers the monthly interest rate λ and the number of financed months κ used in the analysis. The second term of Equation (11) represents the monthly equivalent of the reselling value of the vehicle, V s . The third term calculates the operational cost with fuel, considering a fuel tariff σ (US$/J), the monthly traveled distance L m , and a factor ε to consider taxes, maintenance, and parking. The last term considers time lost in traffic, while taking the monthly wage S m into account, which the occupant earns for working H m hours.
The complete simulation can be summarized in three steps.

Results and discussion
The fundamental diagram of cars (Figure 3) presents the results of our computational algorithm reproducing Meng et al. (2007). Examining the curve for ρ (m) = 0 (no motorcycles), as car density ρ (c) increases, car flow (Equation 7) increases up to ρ (c) = 20. Car flow decreases from critical density until jam density is ρ (c) = 130. This occurs because, under low-density conditions, cars are free to move at the maximum desired speed. However, as the number of cars increases, their movement is disturbed, thus reducing flow. The maximum car flow is 2 404 cars per hour, and it happens in absence of motorcycles (ρ (m) = 0), as well as at a low car density (ρ (c) = 20). In addition, examining motorcycle influence on car flow, as the density ρ (m) increases, car flow decreases. This happens because motorcycles can travel on both lanes of the road, which disturbs car flow.  It is possible to calculate the average speed from Figures 3 and 4 using Equation (7). From the fundamental diagram, it can be seen that the speed for cars is maximum until critical density is reached; and decreases to zero at maximum density, following the pattern described by Daganzo (1994). For motorcycles, when density is low, their average speed is constant for all car densities. As expected, their average speed decreases when motorcycle and car density increases further.
In the next phase of the study, energy consumption Equation (9) is analyzed. Having performed a scale analysis on Equation (3), it can be concluded that energy consumption is predominantly a function of inertia (re-acceleration) and aerodynamic drag (speed). The force of gravity can also be crucial, but, in this study, we have a null inclination. The behavior of energy consumption for cars per km ( Figure 5) is explained by two phenomena: as the density increases, the average speed decreases, hence diminishing drag resistance and increasing the time spent for the same displacement. It also increases the frequency of accelerations. The increased energy consumption for cars until car density reaches ρ (c) = 80 is caused by the increase in the frequency of re-accelerations (drag is almost constant). When ρ (c) > 80, the decrease of car energy consumption occurs both because of the decreasing average speed and decreasing re-acceleration frequency. For very high car densities (ρ (c) > 130), another increase in consumption is noticed, which is caused by the large amount of time spent in the traffic jam. In Figure 5, it is seen that modifying the motorcycle density for the same car density affects the average car fuel consumption, that is, increasing the motorcycle density also increases the car energy consumption. Figure 6 presents the results of energy consumption for motorcycles. For lower motorcycle densities (ρ (m) = 13,3 and ρ (m) = 40,0), consumption is not affected by cars or other motorcycles in the track. For medium and higher motorcycle densities, consumption qualitatively follows the reacceleration frequency. While examining the total energy consumption (Figures 7  and 8), some trends can be noticed. First, for the same total flow, it is better to have a higher proportion of motorcycles, since motorcycles are more fuel-economic than cars. Second, in the region of higher total flow (ρ (m) = 93,3, 0 < ρ (c) < 50), increasing car density can triplicate the energy consumption (consequently, CO 2 emission) with insignificant benefit to the total vehicle flow. Another effect that is worth to discuss, it that for ρ (m) < 93,3 (Figure 7), an increase in car density tends first to increase the total flow up to its maximum value. On the other hand, for densities ρ (m) > 93,3 ( Figure 8) an increase in the number of cars always decreases the total flow for the same motorcycle density. The last metric to be analyzed is individual cost, considering the purchase and reselling of the vehicle, fuel consumption, maintenance, taxes, parking, and the cost of the time spent in traffic (Equation 11) according to the financial considerations informed in Table 3. Figure 9 presents the individual costs for cars in function of the densities of cars and motorcycles. It is possible to notice that the higher the number of cars, the higher is the monthly cost. This happens because of two reasons. First, for the same motorcycle density, as car density increases, the average speed decreases, hence increasing the time spent in traffic. The second effect is fuel consumption ( Figure 5). From zero to intermediate densities, consumption increases with car density. For higher densities, consumption decreases with car density, but it is not enough to compensate for the time lost. A higher motorcycle density increases car costs. Figure 10 presents the individual costs for motorcycles. By comparing them with the monthly car cost (Figure 9), it is possible to notice that the costs for motorcycles are always lower, since they are cheaper and their fuel consumption is lower. Cars do not strongly influence motorcycle costs for low and medium motorcycle densities.

Conclusions
According to this analysis, the best scenario can be determined from each point of view. For car drivers, considering the objective of travelling at a higher average speed, the best scenario occurs in the absence of motorcycles (ρ (m) = 0) and with a low car density (ρ (c) < 20). Considering individual costs, car drivers would rather have an empty road, without other cars (ρ (c) = 0) or motorcycles (ρ (m) = 0). For motorcycle riders, from the point of view of average speed, a low motorcycle density (ρ (m) ≤ 40) would be compelling because they would travel through traffic at maximum speed. The same condition also is interesting financially, since for motorcycles the cost due to lost time is the more relevant.
For society, the metrics should be others. From the point of view of vehicle flow, it would be interesting to have few cars (ρ (c) = 15) and an intermediate presence of motorcycles (ρ (m) = 93,3), with a total density of ρ = 108,3, as seen in the fundamental diagram for vehicles ( Figure 4). Considering CO 2 emissions as proportional to energy consumption (Figures 7 and 8) for a determined flow, a higher motorcycle proportion is responsible for lower the consumption on the track, since motorcycles are more fueleconomic. This study concludes that, for the same total flow, the relative density of cars and motorcycles can have a huge impact on the total CO 2 emissions of the road.
As hypothesized in the introduction, we conclude that the best traffic conditions depend on the metric (e.g., average speed, flow, CO 2 , and costs) and on the point of view (individual or social). It is shown that that the percentage of motorcycles and cars strongly influences the metrics. These results inform that it is not feasible to use results from one study (e.g., richer countries with homogeneous traffic conditions and a small motorcycle fleet) to analyze another one (e.g., developing countries with a mixed traffic of cars of cars and motorcycles), even though both may have the same total flow of vehicles.