Stochastic perturbation of the Lighthill–Whitham–Richards model via the method of stochastic characteristics

In this paper we apply the method of stochastic characteristics to a Lighthill–Whitham–Richards model. The stochastic perturbation can be seen as errors in measurement of the traffic density. For concrete examples we solve the equation perturbed by a standard Brownian motion and the geometric Brownian motion without drift.


Introduction
Many traffic flow models go back to scalar conservation laws of generally non-linear type, i.e.
see e.g. [4]. There are many different results, making the area of scalar conservation laws to an area of very active research in mathematics, see e.g. [3,5,15,16,21,23], to name just a few of those. The function u in (1) describes the density of vehicles on a road and thus has values on the compact set [0, 1]. A conservation law is derived under the assumption, that the time propagation of a mass on a certain interval is only affected by the flux at the boundary of the interval. Hence one often chooses f (u) = u · v(u), where v is the Eulerian velocity of the traffic. There are many different models for traffic flow discussed in the literature where the most famous ones are Aw-Rascle-Zhang model [1,2,24] and Lighthill-Whitham-Richards model, which uses f (u) = u · (1u), i.e. a velocity depending linearly on the density. The scalar conservation law (1) now reads u t + (1 -2u) · u x = 0. ( The flux function in the Lighthill-Whitham-Richards model is in a relatively good agreement with traffic measurements, see [19]. The main problem is that measurements show that data points are quite accurate for low and high densities but are noisy around the maximum point. The flux function is hence rather given by where H is a function vanishing at u = 0 and u = 1 and M t is a suitable nice enough stochastic process. Plugging this in the above conservation law (1) yields The random perturbation is chosen in such a way that H(u) x controls the dependence of the noise with respect to the density in space. The stochastic process itself give a temporal noise, in order to represent fluctuation in time. Of course with the use of a random field the model may be more general. However, the method of characteristics can not be applied in the same way as in this article.
The choice of the Stratonovich integral is due to the existence of the chain rule in the same form as in the deterministic case. For Itô integral a correction term is needed, which leads in some cases to non-existing solutions.
General results for such stochastic conservation laws can be found in e.g. [3, 6-8, 10-14, 20]. The results base on Itô-and Stratonovich type noise in finite and infinite dimensional spaced as well as rough noise and going much further than the model studied in this short article. We hence can of course not give a comprehensive list and refer to the above works and their references.
Industrial application and especially the study of traffic flows are nowadays interested in shocks and nonlinearities. Such shocks occure if the corresponding characteristic intersect, see [4]. In the method of stochastic characteristics, which we treat in this article, this is exactly the same. Our focus in this paper, however, is to show how stochastically perturbed conservation laws behave until the time of shocks. The investigation of shocks would be an interesting study for the future.
In this manuscript, we use the method of stochastic characteristics to solve such equations for the stochastically perturbed Lighthill-Whitham-Richards model explicitly for different cases, where the driving process is given by a Brownian motion or a geometrical Brownian motion. This work can be seen as the starting point to the investigation of different stochastically perturbed hyperbolic equations in direct applications.
The main advantage of the applied method is that the solutions are given in an explicit expression. This makes tedious and highly involved numerical simulation unnecessary.
The motivation of the method is given by traffic flow equations, however, it can be applied to other stochastically perturbed conservation laws.

Prelimiaries
We look at some examples for the Lighthill-Whitham-Richards model for different initial conditions. Due to the underlying model the initial condition describes the density of our traffic problem at time t = 0 and at position x ∈ [0, 1].

Lemma 2.1 Consider the following partial differential equation on
where g(x) is a smooth function. Let (ξ t , η t ) be the solutions to the so called characteristic equations given by ⎧ ⎨ Hence we obtain Then the solution to (4) is given by u( Proof The proof is an direct consequence of Chap. 3 in [9], in particular [9, § 3.2, Theorem 2] with an application of the inverse mapping theorem [ Due to Lemma 2.1 we obtain with the initial condition g(x) = 1x: Thus the solution of the above PDE (6) is given by If we change the initial condition to be g(x) = 1x 2 , we obtain Thus the corresponding solution is equal to One can easily verify that (10) solves indeed the PDE (6).
The main advantage of this method is the precise expression of a solution to a PDEprovided that the corresponding initial condition g(x) and coefficient functions are explicitly given. Due to this fact and for a better comparison between the deterministic and stochastic case we present a collection of solutions in Appendix A.
Along the characteristics the solution remains constant. In the case of the traffic problem and under the considered initial conditions the characteristics never cross each other which means that no shocks appear and hence the solutions are global. As written in the introduction we will study the perturbed case (3) for H(u) = 0 and for M t to be the standard Brownian motion (Bm) as well as the so called geometric Brownian motion defined in the following way.
where W t is a Brownian motion. Hence the geometric Brownian motion without drift is given by

Definition 2.4
Let W t be a standard one-dimensional Brownian motion on a complete separable probability space ( , F, P, F t ), with right-continuous filtration (F t ) t>0 . Then we define for any smooth function The integrals are given in the sense of Stratonovich.
Based on these definitions we are able to apply the so called method of stochastic characteristics to the PDEs as (4) but perturbed by Brownian motion respectively geometric Brownian motion. Since we consider partial differential equations with perturbations by (geometric) Brownian motions we get an ω -dependence in the solutions. The idea of the method is nearly the same as before: now we fix ω ∈ and transform a stochastic partial differential equation (SPDE) into a system of stochastic differential equations (SDEs), solve it and determine the solution to the original SPDE by using stopping times. Hence the precisely determined solutions are given for almost all ω and all space and time variables (x, t) up to a certain stopping time denoted by σ (x). In contrast to the deterministic case we will introduce in the following the method of stochastic characteristics in a more detailed way. Based on Definition 2.4 a perturbed Lighthill-Whitham-Richards model (3) is equivalent to the Cauchy problem Therefore the solution to equation (12) is denoted by u(x, t, ω), but for short notation we only write u(x, t). Suppose u is a solution to (12) and at least one-times continuously differentiable with respect to space and time for fixed ω ∈ . Furthermore, we assume that there exists a curve ξ s (r) which maps the point r ∈ to a point of a neighborhood in at time s. Additionally, we assume ξ 0 (x) = x for all x ∈ [0, 1] as the initial condition. Due to these assumptions we consider and define the following functions, now for fixed ω, r ∈ [0, 1] and s ∈ [0, T]: In the next step we combine (12) with equations (13) and obtain d dt By similar calculations as in [9, § 3.2.1, equation (11)] we get The above stochastic differential equations (SCE) are called stochastic characteristic equations, for a more detailed description and proofs, see also [17]. Given a point x ∈ [0, 1] and assuming that there exist unique solutions to (SCE) starting from x at time t = 0, these solutions solve the corresponding integral equation with initial function g: Let us assume that the solutions (ξ t (x), η t (x), χ t (x)) exist up to a stopping time T(x). As mentioned above we have to work on different stopping times based on the following definition of an explosion time. In the case of our solutions (ξ t , η t , χ t ) this yields to the following definitions of stopping times.

Definition 2.6 Let T(x) be the infimum of all explosion times of the solutions (ξ t , η t , χ t ).
Then we define for all x, y ∈ [0, 1] the stopping times where Dξ t denotes the Jacobian matrix.
Now let the inverse process ξ -1 t of ξ t exist up to some stopping time σ (x). Then we define for almost all ω and for all (x, t) with t < σ (x, ω) the solution Detailed derivations and introductions can be found in [22,Chap. 3]. Now we are able to solve different SPDEs concerning the Lighthill-Whitham-Richards model by using the method of stochastic characteristics.

Application & representation
In particular industrial applications the equations are much more involved than in the examples this played in here. Of course to use our method the characteristics have to computed explicitly. After that standard ODE or SDE methods can be used, if needed.
Based on the flow rate function H(u) and the continuity equation the most natural choice of the drift term is H(u) = uu 2 . In a first step we perturb the Lighthill-Whitham-Richards model by a standard Brownian motion. Hence we consider u(x, 0) = g(x). (15) By using direct computation one can show that the corresponding stochastic characteristic equations are given for almost all ω and all (x, t) up to a stopping time σ (x) by Due to the linearity in the space derivative u x the solution η t (x) = g(x) is always valid. Therefore we receive the solution At this point we compare the characteristics in the deterministic case with the corresponding perturbed one, see Fig. 1. As the initial condition we use here g(x) = 1x.  In this case of g(x) = 1x there exists obviously a process ξ -1 t , such that the inverse property is fulfilled for almost all ω and all (x, t) up to stopping time σ (x), i.e.
The solution to the considered SPDE (15) is given for almost all ω and all (x, t) up to stopping time σ (x) by which looks similar to the deterministic solution (8). Due to the explicit expression of the solution we are able to visualize a sample path easily, see Fig. 2. As introduced in Definition 2.6 the stopping time can be determined explicitly in this example by The perturbation by a geometric Brownian motion as given in Definition 2.3 is in this case straightforward. According to Definition 2.4 we practically can replace the Brownian motion W t by exp(-t/2 + W t ) -1. Let us consider By an application of the method of stochastic characteristics we finally get the precise solution for almost all ω and (x, t) up to a stopping time σ (x) by where we can use the classical l'Hospital argument. The corresponding stopping time is equal to In Fig. 3 we display one sample path with initial condition 1x 2 perturbed by the term Due to this approach we have to verify that the equations (17) as well as (19) really solve the underlying problems. For the sake of simplicity these necessary but lengthy calculation can be found in the Appendix C for equation (17).

Conclusions and discussion
The method of stochastic characteristics can be used effectively to solve a stochastically perturbed Lighthill-Whitham-Richards model. The solutions are explicitly given up to a stopping time in closed form. Numerical simulations based on these models can hence been implemented straightforward. However one has to be careful, that the intersection of characteristics due to stochastic perturbation can lead to solutions which are only defined on a smaller time interval than the non-perturbed ones. On the other hand, it may be also possible, that the stochastic perturbations increase the time interval where solutions are defined. An example for a solution which is ill-defined due to intersecting characteristics can be seen in Fig. 4.
A collection for different examples of stochastic perturbations can be found in Appendix B. Note that with the considered perturbations measurement errors can be modelled effectively. This could be of high interest for more complicated traffic flow models.
well as for different diffusion terms H(u). Taking into account that these might not model the original traffic flow problem perfectly, the approach of the method of stochastic characteristics will give explicit solutions. Firstly we perturb the equation by standard Brownian motion.
• The solution to the equation is given for almost all ω and all (x, t) up to a certain stopping time by • The solution to the equation is given for almost all ω and all (x, t) up to a certain stopping time by • The solution to the equation is given for almost all ω and all (x, t) up to a certain stopping time by .
Replacing the standard Brownian motion by the geometric Brownian motion without drift we are able to determine also explicit solutions to different SPDEs.
• The solution to the equation is given for almost all ω and all (x, t) up to a certain stopping time by • The solution to the equation is given for almost all ω and all (x, t) up to a certain stopping time by e -2 t 0 exp(exp(-s/2 + W s )) ds .
Proof In a first step we determine the partial derivatives u t and u x by using •dW t dt =Ẇ t . We obtain du dt = d dt Finally we have to verify that u t + (1 -2u)u x + (1 -2u)u xẆt = 0.