I. INTRODUCTION
In cities around the world, road transport is a major source of greenhouse gas (GHG) emissions and air pollution, affecting the environment and human health, respectively. From 2022 to 2050, global energy consumption for road transport is projected to increase by 28 %(1). In Ecuador, the energy balance calculated for 2021 by the Ministry of Energy and Mines estimates that transport with approximately 49 % is the sector that showed the highest energy consumption. In fact, diesel and gasoline were the energy sources with the highest demand, with annual growth rates of up to 25 %(2). As a result, the sector with the highest GHG emissions in Ecuador in 2021 was transport with 50.7 %. Hence, the energy and environmental problems produced by road transport cannot be ignored.
Fuel consumption and emissions of vehicles under real traffic conditions are the result of the interaction of several factors and conditions(3). In general, these conditions can be categorized into: (I) vehicle design (for example, vehicle type, powertrain configuration and technology, fuel type and aftertreatment system technology). (II) driver characteristics. (III) travel conditions (for example, passenger load and use of auxiliary systems). (IV) traffic flow conditions (for example, congestion). (V) road conditions (for example, slope level and road type). (VI) environmental conditions (for example, temperature). And (VII) geographical conditions (for example, altitude). Hence, unlike laboratory testing with controlled conditions, evaluation under real traffic conditions is a strategy that the scientific community is currently applying to obtain real-world vehicle performance.
In recent times, several studies have reported that vehicles show significant differences between the emission levels measured experimentally in real traffic conditions and the levels obtained in vehicle type-approval processes with controlled operating conditions. This has led to additional and rapid changes in enforcement in the United States(4), China and the European Union(5). In fact, Europe already requires the use of portable emission measurement systems (PEMS) and on-board diagnostic (OBD) systems to verify vehicle emissions, both for type-approval and real-traffic testing(6). In Latin America, Brazil has adopted OBD requirements like the Euro V standard. In the case of OBD systems, they have been designed to monitor the performance of the engine and aftertreatment system components of the vehicles. Consequently, OBD systems have become a key component of vehicle inspection and maintenance programs, as well as for verifying that a vehicle meets regulatory compliance limits(7). In general, PEMS measurements and OBD data collection are the two strategies commonly applied worldwide to assess the performance of vehicles and engines in real-traffic conditions. However, in Latin America PEMS equipment is not widely available. Hence, the OBD approach becomes a valid strategy to be applied in studies developed in countries such as Ecuador.
In parallel, to minimize discrepancies between type-approval emission levels and those in real traffic, the European Union has introduced the Regulation (EU) 2016/427(8), which includes a testing procedure called ‘Real Driving Emissions’ (RDE) as a complement to the globally harmonized light-duty vehicle test cycles (WLTC)(9). RDE test method is based on the use of PEMS to monitor on-road emission levels of light-duty vehicles. In addition, the routes defined in the RDE tests need to meet certain characteristics: environmental, dynamic, driving behaviour and equipment accuracy, among others described in the Regulation (EU) 2017/115(10). Therefore, the design of a test route according to the RDE regulation would be a strategy to improve the reliability of the data obtained in a study.
Previous work has focused on assessing the fuel efficiency and emission levels of vehicles using four approaches: (a) emission modelling and vehicle simulators, (b) vehicle specific power (VSP) based assessments, (c) real traffic testing with PEMS, and (d) engine map development. Some relevant work for each approach can be mentioned below. In the case of emission modelling, Wen et al.(11)recently used a random forest model to map real-time vehicle emissions on a motorway in Chengdu (China). Regarding the VSP approach, this was used by He et al.(12)to develop and statistically compare VSP distributions as a function of speed and vehicle typology in Beijing. In the case of experimental PEMS data, Bishop et al.(13)used it to identify discrepancies in fuel efficiency and nitrogen oxide (NOx) emissions between homologation processes and real traffic for various diesel and gasoline vehicle models. Regarding engine maps, Mera et al.(14)combined engine operation data and the VSP approach to reduce emission estimation error. It can also be mentioned that Castresana et al.(15)recently investigated the robustness of artificial neural networks to predict multiple performance and emission parameters using information from a diesel engine map for marine applications. In general, the engine mapping approach characterizes the fuel consumption and emission performance of a vehicle’s engine for different engine speed and torque ranges(16). In addition, this resulting information is very useful as input data for emission models and vehicle simulators.
Based on the literature review, three gaps were identified concerning engine maps: (I) maps have not been developed using only data from the vehicle’s OBD system. (II) transient maps have not been constructed from a vehicle operating under a test route designed according to RDE regulations. (III) in Latin America, the engine map approach has not been applied to assess vehicle performance. Hence, this study focuses on the development of engine maps of fuel consumption and carbon dioxide (CO2) emissions for a vehicle operating under real traffic conditions on an RDE test route in the city of Ibarra, Ecuador.
There are three research questions posed for this study:
Is it possible to develop contour engine maps using information from the OBD system of a light-duty vehicle operating in real traffic?
What are the most used torque and RPM ranges in engine maps of a vehicle operating on an RDE test route?
How do fuel efficiency and CO2emission levels vary for different engine map zones?
What are the average fuel consumption and CO2emission factors of a light-duty vehicle for the motorway, suburban and urban sections of the RDE route?
II. METHODOLOGY
A. Test vehicle
A representative passenger vehicle powered by gasoline of the ecuadorian vehicle fleet was selected to analyze the performance of an engine in real traffic conditions. In fact, in the last decade, this vehicle has been one of the most used models for the taxi service in several cities of Ecuador(17). The tested vehicle was a Chevrolet Aveo Emotion, equipped with an engine of 1.6 liters. In general, the characteristics of the test vehicle are representative of passenger capacity, weight, engine torque and power and gearbox type for the automotive segment of the country.Table Iprovides additional information of the tested vehicle. Moreover, it is important to note that the tested vehicle underwent an exhaustive preventive and corrective maintenance program to guarantee the results of this study. In fact, this maintenance ensured that the vehicle successfully passed the mandatory annual technical inspection required by ecuadorian regulations. Consequently, the evaluated vehicle was in optimal technical condition before the execution of the experimental campaign.
B. Test route
In this study, an experimental campaign was conducted along a test route designed in the city of Ibarra, Ecuador. This city is the capital of the province of Imbabura and has a population of approximately 200 000 inhabitants. In terms of vehicle fleet, Ibarra has experienced a significant growth of approximately 6 % per year in recent years(18), which has generated challenges in terms of traffic and mobility. In this study, the test route was designed under the regulation (EU) 2017/1151(10), which establishes different kinematic parameters, such as distance by route section, vehicle speeds, cumulative altitudes, and total trip duration. Hence, this study was carried out on an RDE test route (with approximately 53 km) including urban (0-60 km/h), suburban (60-90 km/h), and motorway (>90 km/h) sections, as shown inFig. 1. However, it is important to mention that this study is not intended to represent an RDE homologation process to evaluate the fuel consumption and CO2emissions of the vehicle. In addition, details of the test route are provided inTable II.
C. Instrumentation
In general, the study required instruments for measuring the following three types of information: (I) vehicle kinematics (for example, position and speed), (II) vehicle and engine operating parameters (for example, torque and speed engine values), and (III) instantaneous fuel consumption data. A global positioning system (GPS) logger device, called ‘GL-770’, was installed in the test vehicle to record latitude, longitude, altitude, and vehicle speed with a frequency of 1Hz. It is important to note that this device is non-intrusive, resulting in the collection of vehicle kinematic information without affecting the regular operation of the tested vehicle operation. To record engine operating parameters and the vehicle fuel consumption data, an OBD interface device (called ELM Electronics 327) was connected to a mobile phone application (called Torque Pro). The ELM 327 was connected to the OBD2 diagnostic port to read the engine operating parameters in real time from the engine control unit (ECU). At the same time, these parameters were sent via Bluetooth from the ELM 327 to the Torque Pro mobile application to be recorded. In this study, the paid version of Torque Pro was used to set the 1 Hz sampling time of the data. The three instruments used for field data collection are shown infigure 2.
D. Experimental campaign
As mentioned above, a Chevrolet Aveo vehicle equipped with GPS and OBD interface devices was evaluated on a test route in the city of Ibarra. The experimental campaign of this study was conducted in February 2023. The campaign for the tested vehicle included both peak (7:00 and 13:00) and non-peak (10:00 and 15:00) periods. The vehicle was driven normally on the test route, respecting the traffic signals and speed limits established in the city. During the experimental campaign, the vehicle was powered by a commercially available type of gasoline called ‘Extra’ in Ecuador. This gasoline is characterised by an octane number of 85. The air conditioning of the vehicle was turned off to avoid additional loads on the engine operation during the experimental campaign. Before the start of each test, the vehicle was warmed up for 15 minutes until the engine reached normal operating temperature. To ensure the repeatability of the tests, the same driver drove the vehicle during all tests. Regarding the passenger load, the campaign was conducted with two people: a driver and a staff technician. The experimental campaign included six entire trips on the RDE route. These tests were conducted during peak periods (7:00 and 13:00) and off-peak periods (10:00 and 15:00). This strategy was designed to ensure that the results obtained accurately reflected vehicle performance at different levels of congestion. The average duration of each entire trip (urban, suburban and motorway road sections) on the test route was approximately 1.7 hours, and the average trip speed of the vehicle was approximately 35 km/h.Table IIIshows an overview of the operating conditions for the tested vehicle by road section in the experimental campaign.
E. Pre-processing real-world vehicle data
In this study, the applied preprocessing method was developed using R Studio software(19). As the measurements of the experimental campaign were recorded using two devices with different initial recording times, a preprocessing step was developed to synchronise the different signals from the two devices through the vehicle speed. This speed was recorded second-by-second using both the GPS position data and through the OBD system of the vehicle. It is important to mention that the speed profiles were smoothed using a moving window filter(20). Overall, the data pre-processing stage ensures to guarantee of the results presented in this study.
F. Method for developing engine maps based on OBD data
Three variables are necessary in order to develop an engine map graph. Typically, torque and speed are used as the primary factors for the horizontal and vertical axes of the maps. Meanwhile, the third variable (Brake Thermal Efficiency (BTE), fuel consumption, brake-specific fuel consumption (BSFC) and emission rates) served as the parameter defining the theme of the map. In this study, the method applied to develop the engine maps included the following two stages: (1) generate grid engine maps by engine speed and torque ranges with averaged values of the third variable, and (2) develop two-dimensional contour engine maps based on grid data. These two stages were performed using R Studio(19).
Initially, the grid engine maps were developed following the methodology proposed by Rosero et al.(21). Initially, grids were generated based on torque and speed engine intervals, and the remaining engine operating data (referred to as the third variable) were assigned to their corresponding grid. Subsequently, the grouped data were averaged to obtain a single value by grid (speed-torque). Outliers were filtered by examining the relative frequency and standard deviation of each data group. Finally, the averaged values were discretized and assigned specific colours for each data interval when plotted on the engine maps.
Subsequently, to plot the two-dimensional contour engine maps, this study utilized the RStudio software(19)with an R package called ‘ggplot2’ R(22). This package includes a function named ‘geom_contour_filled()’, which generates contour plots using z variable observations (for example, fuel consumption, and emission rates) that should be previously grouped into interval combinations of two variables, x (engine speed) and y (torque). Based on this requirement, the data matrices resulting from the grid engine were used as input data to generate the contour engine maps. Finally, to develop the Contour-plot engine maps, a divergent colour scale was applied for each defined interval of the third variable. As a result of these two processing stages, this study obtained contour maps of BTE, BSFC, brake-specific CO2emissions, fuel consumption rates and CO2rates.
III. CALCULATION METHODS
A. Effective power and torque of the engine
The OBD system provided engine load values ( L e obd ) as a percentage relative to the maximum effective power ( 𝑃 𝑒 𝑚𝑎𝑥 ) that the engine could produce. Based on the technical specifications of the tested vehicle manufacturer, the 𝑃 𝑒 𝑚𝑎𝑥 assumed for this study was 76 kW. To determine the actual effective engine power ( 𝑃 𝑒 ) in kW at any given moment, equation 1 was utilized.
The instantaneous engine torque values (𝜏) expressed in N.m, were calculated by multiplying the effective power and the engine speed.
B. Fuel energy efficiency and CO 2 emission factors
The grid engine map generation process involved the calculation of the second-by-second BTE and BSFC values expressed in (%) and (𝑔/𝑘𝑊ℎ), respectively, for the tested engine. These calculations were performed using equation 2 and 3, as follows:
Here, 𝐹𝐶 is the instantaneous fuel consumption expressed in (g/s), ?? 𝑡ℎ represents the instantaneous thermal power (kW) resulting from the fuel being injected, while 𝐶 𝑉 𝑓𝑢𝑒𝑙 is the calorific value of the gasoline used (46.3 𝑀𝐽/𝑘𝑔)(23).
In this study, based on the stoichiometry of gasoline combustion, the CO2emission rate (𝐸 𝑅 𝐶 𝑂 2 ) was calculated using equation 4 as follows:
Here, since the fuel utilized is gasoline (C8H18), the assumed values were 18 for hydrogen (H) and 8 for carbon (C)(24).
To assess the performances of the tested engine for different road sections (𝑖), the average values of BTE (%), BSFC (g/kWh), energy-specific CO2emission factor (g/kWh), as well as fuel consumption and CO2emission rates (g/s), were estimated using equations 5, 6, 7, 8, 9, as follows:
Here, T is the number of second-by-second test data.
Moreover, to evaluate the vehicle performance for different road sections, the average distance-specific fuel consumption 𝐹𝐶 𝑑 and CO2 emission 𝐸 𝐹 𝐶 𝑂 2 𝑑 factors, expressed in (g/km), were estimated using equations 10 and 11, as follows:
IV. RESULTS AND DISCUSSION
The results provided in this section address the research questions and are organized into three topics: (a) typical engine operation patterns by torque and speed for different road sections, (b) contour engine maps of BTE, fuel efficiency, and CO2 emissions, and (c) the overall fuel consumption and emission factors.
A. Typical engine operation patterns by torque and speed
Figure 3presents the patterns of engine operation by torque and speed for motorway, suburban, and urban roads. Grey circles represent the cumulative relative frequency (CRF) of engine operation points calculated for each torque and speed interval, while yellow shaded areas identify the highest engine operation patterns. In the case of the motorway section, the engine operated approximately 85 % of the time in a speed range of 2 500-3 500 rev/min, which are half-load intervals considering that the maximum engine power is reached at 5 800 rev/min. On the suburban section, engine performance was heterogeneous, operating approximately 95 % between low and high (75-100 N·m) torque ranges and low and medium (3 000-3 500 rev/min) speed ranges. In the case of the urban section, approximately 50 % of engine operation was between the lowest torque and lowest speed ranges. Based on Table II, these low-load operation results for urban roads are related to the fact that the vehicle operated with approximately 30 % at idle conditions and an average travel speed of less than 20 km/h. Overall, these results demonstrate that the type of road section significantly influences the torque and speed ranges that a vehicle engine can operate in real traffic conditions. This is consistent with that reported by Bishop et al.(25), who also found that the engine load patterns can vary with the type of road section, especially due to large changes in vehicle kinematics. In fact, Bishop et al. found that a light passenger vehicle engine ran approximately 35 % in low (90-110 N·m) torque ranges and low (1 000-1 500 rev/min) speed ranges. This value of 35 % operation at low load conditions is slightly higher than the 30 % found in the present study. Moreover, Tsiakmakis et al.(26)reported that the average travel speed and 𝑣∗𝑎 factors for light passenger vehicles (LDVs) vary significantly with the type of road section. In fact, Tsiakmakis et al. found that between motorway and urban sections, the average trip speed was reduced approximately 3.5 times (from 67 to 19 km/h). Whereas, in the present study, this speed was strongly reduced by 5 times (from 92 to 18 km/h).
B. Contour engine maps based on OBD data
1. Fuel energy efficiency and brake-specific CO2 emission maps
Figure 4illustrates the contour engine maps for BTE (a), BSFC (b), and brake-specific CO2emissions (c). Infigure 4. a) the areas with the highest BTE (light green) values correspond to operating ranges with medium torque and engine speed, while those with the lowest efficiency (blue) are in areas of low torque demand. Additionally, at high torque and speed ranges, the engine shows values above 25 %. Infigure 4.b) and 4.c) it can be seen that the contours with the lowest fuel efficiency and brake-specific CO2emissions (respectively) are in low torque ranges, which are common for an engine running in idle mode.Figure 4. b) shows that the largest contour corresponds to a range of BSFC with values between 175 and 250 g/kWh. This range is typical for engines mounted on light vehicles(27).Figure 4.c) shows that the lowest CO2emissions (light green) areas with values below 600 g CO2/kWh are in the mid-range of torque and speeds. Overall, these results indicate that the lowest fuel consumption and CO2emissions correspond to medium torque and engine speed ranges. Similar patterns were found by previous studies that evaluated engine performance of engines based on the engine mapping approach(13),(25). In fact, Bishop et al.(13)found that between medium and high load conditions, the brake-specific CO2 emissions of an engine were reduced by approximately 2.5 times (from 500 to 200 g CO2/kWh). Moreover, these three figures show that in idle mode, which is typical of urban operation conditions, the engine has a lower fuel efficiency in terms of energy.
2. Fuel consumption and emission rates maps
Figure 5.a) and 5.b) present contour engine maps for fuel consumption and CO2emissions rates, respectively. Overall, these figures have similar sharply increasing trends (from light green to blue areas) of strong increase. In fact,Figure 5.a) and 5.b) show that the engine under maximum power operating conditions with an increase in torque and speed, fuel consumption and CO2emissions rates reached the highest values (blue) of 3.5 and 10.0 g/s, respectively. In contrast, as expected, when there is a lower demand for engine power, the fuel consumption and emission rates have the lowest values (light green). Although these green areas show lower fuel consumption expressed in terms of rates, based on Figure 5, it can be stated that these same areas are the least efficient in energy terms. The resultant increasing trends are consistent with previous works for LDV’s(13),(14),(28)and heavy-duty vehicles (HDV’s)(21),(29), and emissions models(30). For instance, Mera et al.(13)found that between low and high load conditions, the CO2emissions rates of an LDV engine increased approximately 20 times (from 0.32 to 7.2 g CO2/s). Similarly, in the present study, the CO2emissions rates were increased by about 14 times (from 0.69 to 10 g CO2/s). It is important to note that the information contained infigure 5, expressed in terms of (g/s), can be used as input data to develop an empirical emission model. Previously, the information from the maps must be arranged in the form of a matrix by torque and engine speed intervals(21),(31). The model can predict second-by-second fuel consumption and CO2emission values of the vehicle for any given operating scenario. Hence, future work will focus on using the engine map matrix resulting from this study for predicting fuel consumption and emissions of a vehicle by using torque and speed signals from engine operation as input data.
C. Overview of vehicle engine performance
Table IVshows the overall engine performance of the tested vehicle for motorway, suburban, and urban road sections. The performance of the vehicle is expressed in terms of energy-emission factors, average emission rates and distance-emission factors. In the case of BTE, the average value was 27 % for all sections, although in urban a BTE of approximately 25 % was found. This result is consistent withfigure 3and4, which showed that in urban conditions, the engine frequently ran in low torque and speed ranges, where there was lower fuel efficiency. In the case of energy-emission factors, this study found an average value of BSFC and CO2of approximately 275 and 716 g/kWh, respectively. The minimum BSFC was for the motorway section with 250.18 g/kWh. In the case of the CO2emission factor in terms of distance, the average value was approximately 190 g/km. Comparing urban and motorway road sections, CO2emission factors of the vehicle increase by up to 50 %. The average value of 190 g CO2/km found in the present study is higher than the range of 140-150 g CO2/km reported by Zurita et al.(32), who used standardized driving cycles (for example, New European Driving Cycle (NEDC) and Federal Test Procedure (FTP) 75) to evaluate the performance of a Chevrolet Aveo under controlled laboratory conditions. This demonstrates the importance of a vehicle being evaluated with driving cycles that reflect the actual kinematic conditions of the spatial domain in which the vehicle operates(33). Otherwise, there may be a discrepancy between the values obtained in type-approval processes and those obtained in real traffic conditions. Overall, these results demonstrate that the vehicle and its engine performed better in highway conditions. Hence, to improve engine performance, especially in urban conditions, the use of start-stop technologies in engines is a valid strategy(28). This leads to reduced engine running times in low load ranges (idle conditions), which are characterized by low fuel efficiency as this study has shown.
V. CONCLUSION
Under real traffic conditions a vehicle interacts with the road, other vehicles, and traffic control devices. The level of traffic influences driving patterns and, consequently, this could affect travel times, fuel consumption and emission levels of vehicles. In this context, this study aimed to create contour engine maps of fuel consumption and CO2emissions from a vehicle operating under real-world traffic conditions. A typical gasoline-powered light vehicle from the ecuadorian car fleet was chosen for the experimental campaign developed on an RDE route. This study presents a novel methodology for constructing two-dimensional engine contour maps using real-time recording of engine and vehicle parameters from the OBD system.
The results showed that the engine of the tested vehicle operated with a BTE of 27 %, a BSFC of 275 g/kWh, and a CO2energy-emission factor of 716 g/kWh under real traffic conditions. Furthermore, in terms of distance, the CO2emission factor for the tested vehicle was approximately 190 g/km. Overall, these results demonstrate two aspects: (I) the tested vehicle and its engine had lower fuel efficiency for urban conditions, which were characterized by frequent engine operation at low-load conditions; and (II) it was demonstrated that the OBD approach is a valid method to evaluate the performance of a vehicle in real traffic conditions. Hence, the OBD approach is a valid alternative for Latin American countries, where PEMS equipment is not readily available. Finally, it is important to note that the resulting engine maps of this study include valuable information that can be used as input data for vehicle simulators, as well as for developing emission inventories in cities.