. This is called an initial value problem. Advection Equation (1) is being numerically solved using various one-dimensional explicit numerical schemes, for example Upwind, Downwind, FTCS and Lax-Friedrichs schemes. In these one-dimensional explicit schemes the space and time mesh sizes are restricted both by order of accuracy and numerical stability. The said schemes have been visualized with the unique schematic diagrams which are illustrated in Figures 1-4.
Having been studied carefully the space and time mesh sizes and patterns or schematic diagrams of all these schemes, another but a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. This new scheme has been visualized with the unique schematic diagram and is illustrated in Figure 5. The implication of the advection Equation (1) together with initial value problem in the context of Lighthill Whitham and Richards (LWR) traffic flow model is implemented. The LWR model is the well-known model and describes the traffic flow using a partial differential equation based on the conservation law of the vehicles in traffic     . In this study, the single lane traffic flow with constant speed and linear density-speed relationship are considered as a test case. The simulations of traffic flow model with constant speed will also be done by using the mentioned one-dimensional explicit numerical schemes including Tolesa scheme, and linear density-speed relationship will be simulated by using Tolesa scheme.
In this research work we present and discuss some of the one-dimensional explicit numerical schemes available in the literature and the newly propose scheme. Specifically, the paper is organized in three more sections that follow this Introduction. Section 2 is devoted to carry out finite difference method. The finite difference approximations of the first order hyperbolic partial differential equation using one-dimensional explicit numerical schemes are presented. Section 3 reports about macroscopic continuum traffic flow depend mainly on three quantities flux, speed and density, and present some cases for speed-density relationship. Finally section 4 deals with the numerical simulation. In this section, we present the discretization of the macroscopic continuum traffic flow model, specifically Lighthill-Whitham and Richards Traffic flow model presented in section 3, using finite difference schemes presented in Section 2.
2. Overviews of Finite Difference Method
In this study, it is mainly focused to carry out finite difference method (FDM). The FDM was introduced by Euler in 18th century and has been greatly regarded as the easiest method and widely used to solve simple geometrical problems  . The FDM is classically obtained by approximating the derivatives appearing in the partial differential equation by a Taylor expansion up to some given order which will give the order of the scheme. The application of FDM in solving a PDE is to transform a calculus problem into an algebraic problem by discretizing the continuous physical domain into a number of cells or intervals, and approximating the individual exact partial derivatives in the PDE by algebraic finite difference approximations.
Let the spatial and temporal domains be divided into N and M cells respectively. The index
represents a grid point where
temporal grids intersect. Here
. Thus, the discretization of space-time domain is assumed as
. The quantities
respectively denote the spatial and temporal step sizes. These are also known as the increments between two consecutive spatial and temporal nodes respectively.
Let the quantity
denote an approximate value of the function
at the grid point or space-time location
. That is
. Thus, the function
can now be replaced with a discrete set of point-wise approximate values
. These approximations are called the finite difference approximations. The finite difference approximations of the first order hyperbolic partial differential Equation (1) using one-dimensional explicit numerical schemes are presented in the following subsection.
One Dimensional Explicit Numerical Schemes
In this subsection, the one-dimensional explicit numerical schemes with Upwind, Downwind, FTCS and Lax-Friedrichs schemes including Tolesa scheme are presented to approximate the spatial and temporal partial derivatives of the advection Equation (1) in different views as shown in Figures 1-5. Also the orders of accuracy and numerical stability have been determined in which the space and time mesh sizes are restricted. The approximations and the schematic diagrams of the mentioned schemes are explained sequentially.
In one-dimensional explicit upwind scheme, the first order spatial and temporal partial derivatives are approximated respectively as
. In this view, the stencil is given as shown in Figure 1 and the approximations of the finite difference form of Equation (1) can be expressed as
The scheme (2) is a first order accurate in both space and time,
. This scheme is stable if the condition
In one-dimensional explicit downwind scheme, the first order spatial and temporal partial derivative are approximated respectively as
In view of this the finite difference approximations form of Equation (1) can be expressed as
Equation (3) is a first order accurate in both space and time
. This scheme is unconditionally unstable.
Figure 1. Stencil of One-dimensional explicit upwind numerical scheme.
Figure 2.Stencil of one-dimensional explicit downwind scheme.
In one-dimensional explicit FTCS scheme, as its name implies, the first order temporal and spatial partial derivatives are obtained respectively by taking a first order forward finite differencing in time and a central differencing in space. These are approximated respectively as
. Then, the approximations of the finite difference form of Equation (1) can be expressed as in Equation (4) and the schematic diagram is given in Figure 3.
Equation (4) is a first order accurate in time and second order accurate in space, i.e.,
. This scheme is unconditionally unstable.
In one-dimensional explicit Lax-Friedrichs scheme, the term
in FTCS scheme is replaced by its average value
. In Lax-Friedrich scheme the finite difference form of Equation (1) can be expressed as (5) and its stencil is viewed as in Figure 4.
Equation (5) is a first order accurate in time and second order accurate in space, i.e.,
. This scheme is conditionally stable, if the condition
Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. The one-dimensional explicit Tolesa numerical scheme is still another alternative numerical scheme to solve advection equation and can be applied to traffic flows model. The schematic diagram of Tolesa evolution scheme is shown in Figure 5 and the application of this scheme to the advection Equation (1) is straightforward. In view of this, the first order temporal and spatial partial derivatives are approximated respectively as Equations (6) and (7).
Figure 3. Stencil of one-dimensional explicit FTCS scheme.
Figure 4. Stencil of one-dimensional explicit Lax-Friedrichs scheme.
Figure 5. Stencil of one-dimensional explicit Tolesa scheme.
Let the term
on the left hand side of (6) be expressed as the average value as
In (8), the terms on the left hand side can be expanded, as in view of (5), as
Using the three expansions (8)-(10) in the two Equations (6)-(7) and on substituting them in (1), the advection equation reduces to the form as
Here in (11), the notation
is used. Thus, (11) is the newly proposed one-dimensional explicit Tolesa scheme. Equation (11) is a first order accurate in time and second order accurate in space,
Local Truncation Error Local truncation error represents the difference between an exact differential equation and its finite difference representation at a point in space and time. Local truncation error provides a basis for comparing local accuracies of various difference schemes. Accordingly the local truncation error for Tolesa scheme (11) will be
by using Taylor’s expansion. Therefore, Tolesa scheme is first-order in time and second-order in space.
The fundamental properties that every finite difference approximation of a partial differential equation should possess are consistency, convergence and stability.
Consistence: The notion of consistency addresses the problem of whether the finite difference approximation is really representing the partial differential equation. We say that a finite difference approximation is consistent with a differential equation if the finite difference equations converge to the original equation as the time and space grids are refined. Hence, if the truncation error goes to zero as time and space grids are refined we conclude that the scheme is consistent. For the explicit solution to the advection equation, the truncation error is,
, hence the Tolesa scheme is consistent with partial differential Equation (1) as long as
Stability Analysis: A finite difference scheme is stable if the scheme do not allows the growth of error in the solution with different time level. Stability analysis is a useful tool for checking validity of a given numerical scheme  . There are many approaches to analyze whether a finite difference scheme is stable or unstable. In this paper, we will consider the Von Neumann stability analysis for presented finite difference schemes. The basic idea of this analysis is given by defining the discrete Fourier transform of u as (12). Let it be assumed that the solution can be seen as eigenmodes  which at each grid point have the form
Here in (12),
is a complex number dependent on p and it works as an amplification factor; p is a real spatial wave number;
is an imaginary number. Equation (12) shows the time dependence of a single eigenmode. The differential equations are said to be stable if
Also, the Courant-Friedrichs-Lewy CFL criteria for stability say that
if and only if
. CFL is necessary condition for stability  .
Substituting (12) in (11) and solving the expression for
From Von-Neumann stability analysis, Tolesa method is stable when
After some computations, and considering different cases which satisfy equality and the CFL condition, then we decided that the
Thus, the explicit Tolesa scheme is conditionally stable (13).
Convergence: A numerical scheme is convergent if the computed solution of the discretized equation leads to the exact solution of the differential equation as the time and grid spacing lead to zero. The computed solution
must approach the exact solution
of the differential equation at any point
lead to zero while keeping
constant. In other hand, the error
satisfying the following convergence condition.
, at fixed
. Hence, the explicit Tolesa scheme is convergent; we can see this from the simulations result.
3. Lighthill-Whitham and Richards Traffic Flow Model (LWR Model)
In this section, macroscopic continuum traffic flow model is introduced and analyzed. Macroscopic continuum traffic flow depends mainly on three quantities: traffic density, traffic flow or flux and traffic velocity  . The number of vehicles on a highway per unit length is defined as traffic density and is denoted by
. The traffic flow rate or flux is defined as the number of vehicles passing through a given point x at time t and is denoted by
. In this study, highway is considered as a unidirectional roadway of finite length with no entrances and exits.
The well-known Lighthill-Whitham and Richards (LWR) model describes the traffic flow using a partial differential equation constructed based on the conservation law of the vehicles in traffic. In this model the traffic flow is represented using a first order hyperbolic partial differential equation and is put as
The flux can also be expressed in terms of the traffic density and the traffic speed as
In view of (15), the traffic flow model (14) with initial condition takes the form as
Equation (16) is called an initial value problem IVP of the macroscopic traffic flow model.
In this study, Equation (16) has been considered in two different cases depending on the speed-density relationship as constant speed, and linear speed-density relationship. Also assuming that traffic flux and speed are expressed as a function of density
Case 1: In this case Equation (16) can be expressed as (17) by considering constant speed v and the flux as a function of density
The analytical solution of the form (17) has been calculated using the method of characteristics in implicit form  as follows.
Case 2: Linear speed-density function
in 1935, Greenshields  proposed what was perhaps the first traffic flow model. According to his observations made using photographic methods, Greenshields postulated that there existed a linear relationship between speed and density. Then traffic flux
. In this case Equation (16) can be expressed as