Abstract

In a sparse vehicular ad hoc network, a vehicle normally employs a carry and forward approach, where it holds the message it wants to transmit until the vehicle meets other vehicles or roadside units. A number of analyses in the literature have been done to investigate the time delay when packets are being carried by vehicles on both unidirectional and bidirectional highways. However, these analyses are focusing on the delay between either two disconnected vehicles or two disconnected vehicle clusters. Furthermore, majority of the analyses only concentrate on the expected value of the end-to-end delay when the carry and forward approach is used. Using regression analysis, we establish the distribution model for the time delay between two disconnected vehicle clusters as an exponential distribution. Consequently, a distribution is newly derived to represent the number of clusters on a highway using a vehicular traffic model. From there, we are able to formulate end-to-end delay model which extends the time delay model for two disconnected vehicle clusters to multiple disconnected clusters on a unidirectional highway. The analytical results obtained from the analytical model are then validated through simulation results.

1. Introduction

In general, a vehicular ad hoc network (VANET) is formed between nodes on as-needed basis. To create a VANET, vehicles need to have wireless transceivers and computerized modules that enable the vehicles to act as network nodes. Vehicle-to-vehicle and vehicle-to-infrastructure communications using wireless short-range technologies allow vehicles to broadcast and share information that can be used in different types of applications such as active road safety applications, traffic efficiency applications, and infotainment applications [1]. In a VANET, when the distance between two vehicles is less than the communication range of the vehicles, these vehicles are able to communicate with each other using the wireless channel. Nevertheless, in a sparse vehicular network, the distance between two vehicles is usually larger than the communication range. Therefore, in this type of network, a vehicle normally employs a carry and forward approach, where it holds the message it wants to transmit until the vehicle meets other vehicles or roadside units. Consequently, the end-to-end delay in VANET is usually high. The study of the end-to-end delay in VANET can be considered as one of the most important investigations in vehicular network because extensive applications such as active safety and emergency response applications [2] require the messages to be transmitted with minimal delay.

Wu et al. [3] have determined an information propagation speed model in a VANET that is based on the carry and forward scheme. However, the derivation of the probability distribution for the time delay when using carry and forward approach or catch-up delay model is limited to only that between a cluster with informed vehicles and another cluster with uninformed vehicles, and the probability distribution of the catch-up model is not presented in the numerical results. Therefore, we are not able to ascertain the type of distribution for the catch-up model.

In this paper, we produce a closed form solution for the catch-up delay model from Wu et al. [3]. By producing a closed form solution for the catch-up time delay model between two disconnected clusters, we are able to approximate the distribution of the catch-up time delay model as an exponential distribution by using regression analysis which, to the extent of our knowledge, has not been done in the literature. Next, we investigate and derive a new distribution model to represent the number of disconnected vehicle clusters on a unidirectional highway. Using a Poisson process as our basic assumption for the arrival of vehicles and vehicle-to-vehicle (V2V) connectivity models, we formulate the distribution model for number of vehicle clusters in a certain length of a highway. Using this model, we are able to estimate the number of vehicle clusters in both sparse and dense networks. The analytical model for number of disconnected clusters is then utilized in our investigation to develop a distribution model for end-to-end delay on a highway for multiple disconnected vehicle clusters using V2V connectivity parameters that include the carry and forward mechanism. We are able to further determine the probability distribution model for the end-to-end delay between the multiple disconnected clusters as an distribution by using the exponential approximation of the catch-up delay between two disconnected clusters and the models for the multiple number of clusters. This framework enables us to assess the end-to-end delivery delay between a source and destination moving on a highway for both disconnected and well-connected vehicular networks.

Although many research studies have been carried out to incorporate the carry and forward approach in their proposed routing protocols, not many studies have been done in analyzing the end-to-end delay when the carry and forward approach is employed during the packet forwarding in VANET. Wu et al. [3] have presented an analytical study on the information propagation speed when the carry and forward approach is used in both one- and two-way highway scenarios where vehicle arrivals are based on Poisson process and the vehicle speeds are uniformly distributed in a designated range. The authors provide numerical results on information propagation speed under two network models, which are low density network and high density network.

A number of researches have developed analytical models for studying vehicular network characteristics and performance metrics [49]. However, these studies are focusing mainly on information propagation speed model, connectivity model, mobility model, and link reliability model. They do not present any probability distribution model on end-to-end delay and the information propagation speed model is normally based on the expected values of the end-to-end delivery delay. Zhang et al. in [6] have modified the information propagation speed model from Wu et al. [3] by using a traffic density for Poisson arrival model and truncated Gaussian distribution for vehicles’ speed. Nevertheless, the study done by [6] does not show any distribution model for the catch-up delay. Furthermore, the authors only include the results on the expected value for the information propagation speed for a VANET highway. The study does not show any development on the distribution model for end-to-end delay.

Wisitpongphan et al. [10] proposed a similar analytical model as [3] for VANET connectivity in a sparse network. Using empirical traffic data, the authors study and formulate VANET parameters such as interarrival time and intervehicle spacing. The authors also derive a comprehensive analytical framework that can be used to characterize a sparse vehicular network for one- and two-directional highways. Furthermore, the authors did an analysis on a parameter similar to the catch-up time in [3] which is referred to as the rehealing time. However, the study on rehealing time in [10] is focused on the two-directional highways and between adjacent vehicles. Through simulation, the authors are able to validate their analytical framework and analyze end-to-end delay for packet transmission with distance between source and destination varying from 1 to 30 km [10]. Nonetheless, the study on the end-to-end delay in [10] is based on the expected values and purely based on simulation results. There is no distribution model derived for the end-to-end parameter.

A study on end-to-end delay model is done in [11] where the authors analyze the total delay time needed by a relay to carry a packet from a source to a destination using the carry and forward system. The main goal of this study is to find the relay’s optimal location that minimizes the total delay while taking into account the effect of channel fading, path loss, and forward error correction. However, the study is based on a mobile ad hoc network scenario with only one relay between the source and destination.

In [12], the author also proposes a similar study on the end-to-end model, where the author uses the ergodic Markov chain to model the vehicle’s mobility, the exponential distribution for the initial vehicle density, and the normal distribution for the average vehicle speed. The author creates the model for vehicles that are sparsely arranged on a one-directional straight road. Using these assumptions, the author is able to obtain expressions for the exact delay time and delivery ratio. Nevertheless, the model in [12] only considered transmission between two vehicles, not between clusters of vehicles. In addition, the model does not consider the carry and forward approach during packet transmissions. Instead, the authors use -seconds-wait rule where the packets are discarded if vehicles are unable to transmit them within seconds.

3. System Model of the Time Duration for the Catch-Up Phase

Consider a scenario presented in Figure 1 where a number of vehicles independently travel along a unidirectional highway of length meters. The speed of each vehicle, , is modeled using uniform distribution over . Based on the studies done in [3, 8, 13], vehicles’ speed can be represented with uniform distribution with interval , where each vehicle freely moves at its chosen velocity. Each vehicle is then assumed to move along the highway at a constant speed,  m/s, such that the distance between the vehicle and its neighbors remains unchanged. The model assumes that vehicles arrive to the highway following a Poisson process with an average rate equal to the traffic flow rate (vehicles per unit time). Empirical studies have shown that Poisson arrival model can be used to model vehicle arrival rate in free flow phase [10] and it is commonly used model in the studies of VANETs [69]. The vehicles are then partitioned into a number of clusters. In this model, a cluster can be defined as a group of vehicles that are able to propagate messages using multihop forwarding via wireless channel. Road traffic statistics in [14] have shown that vehicles tend to travel in clusters on a highway. The clusters, which are formed in the highway, are split and merged over time due to the mobility of the vehicles. If the gap between two clusters is larger than the transmission range, , then the carry and forward strategy is used to forward messages. A vehicle is considered as an informed vehicle if the vehicle has the message that needs to be transmitted. In this scenario, we also assume that the source of the messages is located in a cluster of informed vehicles and the receiver of the messages is found in an uninformed cluster located at the end of the highway. The message is transmitted via one of two ways, through either the forward process or the catch-up process. In the forward process, the message is forwarded to other neighboring vehicles within a partition via the wireless channel, where the message rapidly propagates hop by hop until it reaches the farthest vehicle of that partition. In the catch-up process, the message travels along with the carrying vehicle until the carrying vehicle arrives within the communication range of the last uninformed vehicle in the partition ahead of it. Once the carrying vehicle is inside the partition with a group of uninformed vehicles, the message will be again propagated via forward process. Both processes alternate with each other as the message propagates along the road.

The term in Figure 1 is the time duration for the catch-up process, where packets are being carried by the carrying vehicle until the vehicle is able to forward the packets via wireless transmission to the last uninformed vehicle in the partition ahead of it. Although the open form cumulative distribution function (CDF) of has been extensively studied and derived in [3], as shown in (1), the authors do not present the CDF or probability density function (PDF) of in their numerical results. Based on the assumptions, expressions, notations, and model parameters provided by the authors in [3], we are able to produce the CDF and PDF of via numerical integration. Notations and model parameters needed for the derivation of distribution are shown in Notations in this analysis.

The cumulative distribution function (CDF) of as presented in [3] is shown in

Based on (1), we conclude that we need to derive closed form solutions for , , and , which are not presented in [3].

3.1. The Derivation of Closed Form Solution for CDF and PDF of Message Propagation Distance,

Let denote the distance traveled by a first vehicle in the front most informed cluster after passing a random location, , during the time interval (refer to Figure 1). The CDF of is expressed aswhere is given as

Let in (3) denote the speed of the source vehicle located at the location at time, . Let denote the speed of the vehicle and let denote the arrival time of a vehicle at location after .

Based on the assumption of a Poisson arrival, is defined as the number of vehicles arriving at the highway during and expressed as

Using the Poisson process theorem [15, 16], given that vehicles have passed the location between time , the arrival times at which the events occur, considered as unordered random variables, are distributed independently and uniformly in the interval . Therefore, with this theorem, can be presented as . With the assumptions that and are independent and identically distributed (i.i.d.) and uniformly distributed at the intervals and , respectively, (3) can be expressed as

Let denote the vehicles’ speed and it is uniformly distributed at interval . Therefore, can be expressed as

The probability in (5) denotes the probability of the distance traveled by vehicles after passing the location at the speed between the interval and at the time in the interval . With the assumption that the random variable is independent and identically distributed (i.i.d.) and uniformly distributed in the interval , the probability can be formulated as

Using (4), (6), and (7), we are able to solve (2), which is the CDF of and is shown in where Next, we take the derivative of to derive the PDF of , which is displayed in where

3.2. The Derivation of Closed Form Solution for Cumulative Distribution Function of : Distance that the Partition Tail Moves

Let denote the distance traveled by the last vehicle in a uninformed cluster that is in front of an informed cluster during the time interval . The CDF of is expressed in where

In (13), is defined as the speed of the uninformed vehicle at location at time and denote the number of vehicles that pass location . As shown in Figure 1, we conclude there are no other uninformed vehicles passing location after the last uninformed vehicle during a catch-up phrase, and hence the earliest time for the last uninformed vehicle to pass location is , with the condition that the vehicle is in location at the time 0 and the vehicle speed is in the interval . Using the same properties of independent and identically distributed (i.i.d) for and as in Section 3.1, (13) can be rewritten as

Next, given that , under the properties of the Poisson process [15], the jump epochs have the same distribution as the order statistics corresponding to independent variables uniformly distributed on the interval . Therefore, (4) has to be rewritten according to the the new distribution of , which is shown in

In (14), is expressed as

Let denote the probability of the distance traveled by vehicles after passing location in the interval of with a speed between . Therefore, (7) can be rewritten as

Finally, we substitute (15), (16), and (17) in (12) to solve , which is shown in where

3.3. Derivation of PDF for : The Distribution of Disconnected Vehicles Gap

Let denote the gap between two neighboring vehicles and let denote the gap between two neighboring disconnected vehicles. According to Wisitpongphan et al. in [10], the spacing between two vehicles can be expressed by an exponential distribution and the validity of assumption for VANETs has been confirmed by the empirical measurement reported in [10]. Therefore, the gap between two neighboring vehicles is denoted aswhere the parameter , which is the vehicle density, can be estimated as .

Two neighboring vehicles are considered disconnected if the gap between the vehicles is larger than communication range, . If the spacing between two neighboring vehicles follows an exponential distribution, then the derivation of the PDF of is as follows:

3.4. Derivation of CDF and PDF of via Numerical Integration and Approximation of Distribution

It may be noted that the CDF of given in (1) does not have a closed form solution but has to be evaluated via numerical integration. Plots in Figure 2 present the probability of distribution against time in seconds for different values of flow rate.

In general, Figure 2 indicates, with the increase in vehicle flow rates, the distance between vehicles decreases and therefore, decreases the catch-up delay. The -axis in Figure 2 denotes the catch-up delay between two disconnected clusters of vehicles. A large value of indicates that the data packets are carried by an informed vehicle in a catch-up phase most of the time and a small value of indicates that the catch-up phase occurs in a short time and the data packets are transmitted using forwarding phase most of the time. As shown in Figure 2(a), in a sparse network, the catch-up phase occurs at a high delay. At  veh/hr, the frontmost vehicle of an informed cluster of vehicles has a probability of 70% to catch up with the last vehicle in an uninformed cluster of vehicles at a time delay larger than or equal to 200 seconds. However, in a dense network with  veh/hr, the catch-up process is highly likely to happen at approximately 100 seconds or less.

Figure 2 shows that the distribution has the shape of an exponential distribution. Therefore, we use exponential regression analysis to approximate the distribution with an exponential distribution using the following equation [17, 18]:where and are constant called the model regression coefficients.

An exponential regression analysis is performed by applying the logarithm to the base of of both sides of (22). Subsequently, (22) can be written by

By substituting , a linear regression analysis equation, which is expressed as , can be rewritten aswhere the regression coefficients and are expressed aswhere

There are two main parameters in regression analysis that can indicate that an exponential distribution is a good fit for distribution [17, 18].

(1) Sum of Square Errors (SSE). In general, this parameter measures the difference between data points and an estimation model with a value closer to zero to indicate a good fit.

Coefficient of Determination (). This parameter indicates how well data points fit an approximation curve with a value approaching one to demonstrate a good fit.

Table 1 shows the output of the exponential regression analysis for distributions for vehicle flow rates 360, 1080, 1800, and 3600 veh/hr. For each of the traffic flow rates, the exponential regression yields a high value for parameter, which is approaching one. Table 1 also indicates that the parameter SSE yields values that are very close to zero. With the parameter yielding values close to one and the SSE values producing values near to zero in the exponential regression analysis, we can ascertain that our distribution can be approximated with an exponential distribution.

Figure 3 displays the exact distribution function of with its approximation counterpart, given in (22) for the respective vehicle flow rate. From these figures, we can establish the high accuracy between the exact distribution and its approximation, and we determine that the exact distribution function of can be approximated using an exponential distribution expression.

4. Derivation of the End-to-End Time Delay for Unidirectional Highway

Using Figure 1 as an example, total end-to-end delay from a source vehicle to a destination vehicle, , can be expressed aswhere is forwarding time in cluster and is catching up time from cluster to cluster .

It is assumed that low vehicle density in a vehicular network causes the communication range to become smaller than the average intervehicle gap. Therefore, we ascertain that the message transmission time can be approximated using entirely the vehicle movement while ignoring the message transmission time within a cluster, that is, , as it is very small, since the packets are transferred via wireless channel within the cluster, resulting in . We consider that this assumption is plausible based on the current researches on performance analysis on IEEE 802.11p technology [1921]. These researches show that the average transmission delay between two neighboring vehicles via wireless channel takes on the values of milliseconds as opposed to the transmission delay via carry and forward strategy, which has been discussed in Section 3.4.

Therefore, from the above assumption, (27) can be rewritten asGeneralizing (28) for clusters, we are able to derive as

Let denote the number of vehicle clusters on a highway with the assumption that is a random variable that is independent of ’s where derivation will be explained later in Section 4.1. Hence, we can find the conditional PDF of given that using the conditional characteristic function of . From (28), the conditional characteristic function of given that can be expressed as

The conditional PDF of can be found by taking the inverse transform of . However, since the distribution of in (1) is found using numerical integration, a closed form solution is not feasible for (30). In Section 3.4, using exponential regression analysis, we have determined that the catch-up time () distribution in Figure 2 can be approximated with an exponential distribution. Henceforth, can be denoted with the following expression:where the value of should be found using minimum mean square error (MMSE) between the exact distribution function of in (1) and its approximation. Using the approximation of as an exponential distribution, the conditional PDF of given that can be found using characteristic function expressions of and sums of ’s. The characteristic function of can be expressed asFrom (32), we can find the expression for the characteristic function of given that :

Thus, the conditional PDF of can be found using inverse Laplace transform of the conditional characteristic function of in (33) which is expressed in

From (34), we ascertain that the conditional PDF of given that follows an distribution.

4.1. Analysis on the Distribution of Number of Clusters,

To find the PDF of , we have to derive the distribution model for the number of clusters, . Based on the message propagation scenario in Figure 1 of Section 3, we consider a unidirectional highway of length meters. The source vehicle is located in the first informed cluster and the destination vehicle is located at the end of the highway. Therefore, a message from the source vehicle has to be propagated over multiple clusters of vehicles in order to be transmitted to the destination vehicle. In addition, we assume that the vehicles enter the highway according to a Poisson process with traffic flow rate of . Therefore, by employing the Poisson process assumption, we consider that the number of clusters can be modeled by using the Poisson distribution. Subsequently, by applying the vehicular network analytical framework provided in [3], we are able to formulate a distribution model for number of clusters for a unidirectional highway of length .

Let denote number of clusters in a unidirectional highway and its derivation is based on vehicular network parameters provided in [3].

4.1.1. Average Distance during Forwarding Phase,

Let represent the distance traveled by messages during a forwarding phase where the expected value is given aswhere is the transmission range and is the average cluster size and it is expressed aswhere is the average number of vehicle gaps in a cluster and is the average gap between two connected neighboring vehicles.

Next, we derive the following closed form expressions for and , respectively: Therefore, using (37) and (38), (35) can be solved as follows:

4.1.2. Average Distance during Catch-Up Phase,

Let denote distance traveled by messages during a catch-up phase and the expected value is expressed asIt may be noted that the closed form expression is not feasible for (40) since is solved using numerical integration. Therefore, (40) is solved through numerical integration using the exact distribution of in (1).

Finally, the probability mass function for can be expressed aswhere is the length of the highway, is the expected value of the forwarding distance, and is the expected value of the catch-up distance.

4.2. Derivation of Distribution Based on and Distributions

It is stated in [15, 16] that the Erlang distribution is obtained by the -fold convolution of independent exponential distribution. In Section 4, we have shown that the distribution of is the summation of the catch-up delay, , for vehicle clusters. In addition, we establish that is a random variable that exhibits the properties of an exponential distribution.

Therefore, let denote the sum of the for and let denote the number of vehicle clusters on a unidirectional highway. Based on the Law of Total Probability, the PDF of can be formulated aswhere

The PDF of , as shown in (42), exhibits properties of an Erlang distribution, which will be proved in Section 5.1.

5. Results and Analysis for Distribution of Total Catch-Up Time in One-Way Street

5.1. Numerical Results

In this section, we present some pertinent numerical results regarding the analysis done in this chapter. Figure 4 shows the probability mass functions for from (41) and for traffic flow rate and 2520 veh/hr. It can be seen in Figure 4 that, at  veh/hr, the peak of the distributions occurs at and at  veh/hr; the peak of the distribution is at . Consequently, the plots in Figure 4 demonstrate that as the traffic density in a highway increases, clusters in the highway merge to form a larger cluster, thereby reducing number of clusters in the highway.

In addition, there is a correlation between number of clusters and traffic density where the number of clusters decreases as traffic density increases. The relationship is shown in Figure 5, which presents the average distance traveled by messages during forward phase, , and catch-up phase, . Figure 5 shows that the value of also increases with the increment of traffic flow rate, confirming that as the traffic density increases, messages are mostly transmitted via wireless channel rather than being carried by vehicles. From Figure 5, we ascertain that the average cluster length increases as vehicle traffic flow rate increases; the intervehicle spacing reduces until the gap is less than the transmission range.

Furthermore, the trend in the plots of and in Figure 5 derived from our analysis in Section 4.1 shows an exact match with the numerical results of the original and from [3]. The same trend displayed in Figure 5 and the original result in [3] validated the accuracy of our work on the distribution of from [3]. Figure 6 displayed the information propagation speed based on the distribution, which further validates the accuracy of our analysis as Figure 6 shows similar trend with the information propagation speed from [3]. Therefore, we are able to arrive with the same conclusion as [3] in which higher vehicle density leads to a larger partition size and shorter intercluster distance and henceforth reduces the catch-up time. Figure 6 shows a sharp increase as the vehicle flow rate increases and the propagation speed is shown much faster than the vehicle movement.

Figure 7 displays the PDF of distribution. The probability of is plotted against the time delay in seconds for different vehicles flow rates. As shown in Figure 7, the PDF of exhibits similarities to an distribution, which confirms our analysis that the PDF of follows distribution given that .

5.2. Simulation Results

In this section, we present simulation results in Figures 9(a)9(c) for our proposed analytical model using the network simulator NS-2 [22, 23]. Our simulation scenario, which is displayed in Figure 8, is based on a one-directional highway with the length of 15 km. The highway is assumed to have multiple one-directional lanes, where vehicles can overtake each other without changing their lane or maneuvering. Vehicles are generated using a Poisson process with flow rates of 360, 1080, 1800, and 3600 veh/hour. Each vehicle is assigned a random speed based on a uniform distribution between the intervals  m/s and  m/s and the assigned speed does not change over the simulation time. We perform the simulation for 1200 seconds and repeat the simulation for 1000 iterations. Since NS-2 is built to simulate a network environment, we configure media access control (MAC) and physical (PHY) layers in NS-2 to retain [3] assumptions of ideal MAC and PHY layers for the model so that the simulation is executed under ideal communication channel. The packets are generated using Poisson traffic with mean of 0.1 seconds and the transmission range is set to 250 meters. In addition, we configure the source vehicle to be the only informed vehicle at time and to be located at position ; and the destination vehicle is the first vehicle to pass location at time .

Figure 9 displays the comparison between the numerical results of the analysis with the simulation results. Figure 9(a) presents the numerical and simulation results for the PDF of from (42). The results are plotted as a function of the time delay with different values of traffic flow rates. From Figure 9(a), at  veh/hr, the simulation result shows that the peak of the distribution occurs at approximately seconds, whereas, at  veh/hr, the distribution’s peak is at seconds. It can be concluded that as the traffic flow rate increases, the end-to-end delay decreases since the traffic density on the highway will be large enough to allow packet transmission via multihop forwarding on the wireless channel.

In addition, the plots for distribution in Figure 9(a) show that the numerical and simulation plots of are very close to each other and exhibit the same trend as the distribution model. From Figure 9, we ascertain that the small discrepancy between the analytical plots of distribution and the plots of simulation is caused by the configuration of NS-2. We conclude that even though NS-2 has been configured with the ideal communication channel, MAC and PHY conditions in NS-2 still affect the simulation results. Furthermore, we ascertain that the use of Poisson traffic for data packets generation in the simulation causes a small deviation between analytical and simulation results.

6. Conclusion

In this section, we propose an analytical framework for the end-to-end delay model for a vehicular network on a unidirectional highway by extending the catch-up delay model between two adjacent vehicle clusters to multiple vehicle clusters as well as using traffic characteristics models to determine the distribution model for number of clusters. We approximate the distribution of the catch-up delay model between two disconnected clusters using the exponential regression analysis to derive the catch-up delay model for multiple clusters. Using the approximation, we establish that the catch-up delay model for multiple clusters follows an distribution. We also validated our analytical results through simulation. In this study, we are able to confirm that the carry and forward process often occurs in a low density network and the process has caused the packets to be delivered at high end-to-end delay. Nonetheless, the study is done without considering a number of real-world communication aspects such as channel fading and contention issues at media access layer which can be considered as the research future works.

Notations

The length of a road
: Gap between two neighboring vehicles
: Gap between two neighboring connected vehicles
: Gap between two neighboring disconnected vehicles
: Traffic flow rate (vehicles/unit time)
: Vehicles density (vehicles/unit distance)
: Minimum mean square error between the exact distribution of and an exponential distribution
: Number of vehicles arriving at the highway during interval
: Number of vehicle clusters
: Characteristic function of a random variable
: Vehicle radio range
: Time duration of a catch-up phase
: Time duration of a forward phase
:The sum of multiple catch-up times,
: Average speed of vehicle , ; a random variable in the interval
: The distance traveled by messages during a forwarding phase
: The distance traveled by messages during a catch-up phase
: Message propagation distance during
: Distance that the partition tail moves during
: The size of a cluster.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank Universiti Teknikal Malaysia Melaka (UTeM) and Memorial University of Newfoundland for their continuous financial support in making this research possible.