Viscous scalar conservation law with stochastic forcing: strong solution and invariant measure
Sofiane Martel (SIMSMART), Julien Reygner (CERMICS)

TL;DR
This paper studies one-dimensional viscous scalar conservation laws with stochastic forcing, proving existence and uniqueness of strong solutions and invariant measures under broad conditions, including degenerate noise.
Contribution
It establishes the first rigorous results on strong solutions and invariant measures for viscous scalar conservation laws with degenerate stochastic forcing.
Findings
Existence and uniqueness of strong solutions
Existence and uniqueness of invariant measures
Results hold under degenerate noise conditions
Abstract
We are interested in viscous scalar conservation laws with a white-in-time but spatially correlated stochastic forcing. The equation is assumed to be one-dimensional and periodic in the space variable, and its flux function to be locally Lipschitz continuous and have at most polynomial growth. Neither the flux nor the noise need to be non-degenerate. In a first part, we show the existence and uniqueness of a global solution in a strong sense. In a second part, we establish the existence and uniqueness of an invariant measure for this strong solution.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Viscous scalar conservation law with stochastic forcing: strong solution and invariant measure
Sofiane Martel
INRIA Rennes - Bretagne Atlantique, 35042 Rennes, France.
and
Julien Reygner
Université Paris-Est, CERMICS (ENPC), 77455 Marne-la-Vallée, France.
Abstract.
We are interested in viscous scalar conservation laws with a white-in-time but spatially correlated stochastic forcing. The equation is assumed to be one-dimensional and periodic in the space variable, and its flux function to be locally Lipschitz continuous and have at most polynomial growth. Neither the flux nor the noise need to be non-degenerate. In a first part, we show the existence and uniqueness of a global solution in a strong sense. In a second part, we establish the existence and uniqueness of an invariant measure for this strong solution.
Key words and phrases:
Stochastic conservation laws, Invariant measure
2010 Mathematics Subject Classification:
35A01,35R60,60H15
This work is partially supported by the French National Research Agency (ANR) under the programs ANR-17-CE40-0030 - EFI - Entropy, flows, inequalities and QuAMProcs.
1. Introduction
1.1. Stochastic viscous scalar conservation law
We are interested in the existence, uniqueness, regularity and large time behaviour of solutions of the following viscous scalar conservation law with additive and time-independent stochastic forcing
[TABLE]
where , , is a family of independent Brownian motions. Here, denotes the one-dimensional torus , meaning that the sought solution is periodic in space. The flux function is assumed to satisfy the following set of conditions.
Assumption 1** (on the flux function).**
The function is on , its first derivative has at most polynomial growth:
[TABLE]
and its second derivative is locally Lipschitz continuous on .
The parameter is the viscosity coefficient. In order to present our assumptions on the family of functions , , which describe the spatial correlation of the stochastic forcing of (1), we first introduce some notation. For any , we denote by the subset of functions such that
[TABLE]
The norm induced on is denoted by . For any integer , we denote by the intersection of the Sobolev space with . Equipped with the norm
[TABLE]
and the associated scalar product , it is a separable Hilbert space. On the one-dimensional torus, the Poincaré inequality implies that and . Actually, the following stronger inequality holds: if , then and for all ,
[TABLE]
The spaces generalise to the class of fractional Sobolev spaces , where , which will be defined in Section 2.1. We may now state:
Assumption 2** (on the noise functions).**
For all , and
[TABLE]
Let be a probability space, equipped with a normal filtration in the sense of [10, Section 3.3], on which is a family of independent Brownian motions. Under Assumption 2, the series converges in , for any , towards an -valued Wiener process with respect to the filtration , defined in the sense of [10, Section 4.2], with the trace class covariance operator given by
[TABLE]
Thus, almost surely, is continuous in and for all , the process is a real-valued Wiener process with variance
[TABLE]
1.2. Main results and previous works
First, we are interested in the well-posedness in the strong sense of Equation (1). In particular, we look for solutions that admit at least a second spatial derivative in order to give a classical meaning to the viscous term, in the sense of the following definition:
Definition 1** (Strong solution to (1)).**
Let . Under Assumptions 1 and 2, a strong solution to Equation (1) with initial condition is an -adapted process with values in such that, almost surely:
- (1)
the mapping is continuous from to ; 2. (2)
for all , the following equality holds:
[TABLE]
In the above definition, the first condition ensures that the time integral in Equation (7) is a well-defined Bochner integral in . For a careful introduction of the general concepts of random variables and stochastic processes in Hilbert spaces, the reader is referred to the third and fourth chapters of the reference book [10].
Our first result is the following:
Theorem 1** (Well-posedness).**
Let . Under Assumptions 1 and 2, there exists a unique strong solution to Equation (1) with initial condition . Moreover, the solution depends continuously on initial data in the following sense: if is a sequence of satisfying
[TABLE]
then, denoting by the family of associated solutions, for any , we have almost surely
[TABLE]
Similar results have already been established: the case where the flux is strictly convex is treated in [4, Appendix A], and the case where is globally Lipschitz continuous is treated in [22]. Furthermore, the case of mild solutions (in spaces) has been looked at in [21]. Here, no global Lipschitz continuity assumption nor restrictions on the convexity of the flux function are made. We can also point out that the well-posedness of stochastically forced conservations laws in the inviscid case (i.e. when ) has been under a great deal of investigation in the recent years. In this "hyperbolic" framework, the appearance of shocks prevents the solutions to be smooth enough to be considered in a strong sense as in our present work. Therefore, the study of entropic solutions [19] or kinetic solutions [13, 20] to the SPDE have been the two main approaches, both of which rely on a vanishing viscosity argument: the entropic or kinetic solution is sought as the limit of its viscous approximation as the viscosity coefficient tends to [math].
More recent works concern the Burgers equation with stochastic transport noise in the viscous and inviscid cases [1], or the spatial regularity for solutions of the viscous Burgers equation with additive noise [23]. A natural extension of our works would be to consider a viscous conservation law with multiplicative noise or even, as in [1], a transport noise.
Let denote the set of continuous and bounded functions from to . As a consequence of Theorem 1, we can define a family of functionals on by writing
[TABLE]
where the notation indicates that the random variable is the solution to (1) at time starting from the initial condition .
Corollary 1**.**
Under Assumptions 1 and 2, the family is a Feller semigroup and the process is a strong Markov process in with semigroup .
Proof.
The uniqueness of a strong solution and the fact that, for all , the processes and have the same distribution, ensure that is a semigroup, and therefore that is a Markov process. The Feller property is a straightforward consequence of the result of continuous dependence on initial conditions given in Theorem 1, whereas it is a classical result that the strong Markov property of follows from the Feller property of (see for instance the proof of [7, Theorem 16.21]). ∎
Let denote the Borel -algebra of the metric space , and refer to the set of Borel probability measures on . The Markov property allows us to extend the notion of strong solution to (1) by considering not only a deterministic initial condition but any -measurable random variable on . In this perspective, we define the dual semigroup of by
[TABLE]
In particular, is the law of when is distributed according to .
Definition 2** (Invariant measure).**
We say that a probability measure is an invariant measure for the semigroup (or equivalently for the process ) if and only if
[TABLE]
Theorem 2** (Existence, uniqueness and estimates on the invariant measure).**
Under Assumptions 1 and 2, the process solution to the SPDE (1) admits a unique invariant measure . Besides, if is distributed according to , then and, for all , .
A few similar results exist in the literature. Da Prato, Debussche and Temam [9] have studied the viscous Burgers equation (which corresponds to the flux function ) perturbed by an additive space-time white noise whereas Da Prato and Gatarek [28] studied the same equation but with a multiplicative white noise. Both showed the well-posedness of the equation as well as the existence of an invariant measure. These results are moreover put in a much detailed context in the two reference books [10, 11]. Boritchev [3, 4, 5] showed the existence and uniqueness of an invariant measure for the viscous generalised Burgers equation (which corresponds to the case of strictly convex flux function) perturbed by a white-in-time and spatially correlated noise. E, Khanin, Mazel and Sinai [18] showed the existence and uniqueness of an invariant measure for the inviscid Burgers equation with a white-in-time and spatially correlated noise. Debussche and Vovelle [14] generalised this last result by extending it to non-degenerate flux functions (roughly speaking, there is no non-negligible subset of on which is linear). Besides, the fact that these results from [18, 14] also hold when makes them quite powerful: it shows indeed that the presence of a viscous term is not a necessary condition for the solution to be stationary. On this topic, we refer the reader to a recent nicely detailed survey by Chen and Pang [8].
The stochastic Burgers equation is mainly studied as a one-dimensional model for turbulence. By showing a stable behaviour at large times, this model manages, to some extent, to fit the predicitions of Kolmogorov’s "K41" theory about the universal properties of a turbulent flow [25, 24]. Whether it is modelled by the Burgers equation or a by more general process such as Equation (1), turbulence is then described through the statistics of some particular small-scale quantities in the stationary state [16, 17]. Sharp estimates were given by Boritchev for these small-scale quantities [4], which were furthermore shown to be independent of the viscosity coefficient. One of the purposes of this paper is to lay the groundwork for the numerical analysis of Equation (1). In a companion paper [6], we introduce a finite-volume approximation of (1) which allows to approximate the invariant measure . Generating random variables with distribution shall eventually lead us to compute said small-scale quantities and analyse the development of turbulence in the model established by Equation (1).
1.3. Outline of the article
The proofs of Theorems 1 and 2 are respectively detailed in Sections 2 and 3.
2. Well-posedness and regularity
This section is dedicated to the proof of Theorem 1. This proof is decomposed as follows. In Subsection 2.1, we introduce a weaker formulation of Equation (1), the so-called mild formulation. In Subsection 2.2, we show that Equation (1) is well-posed locally in time both in the mild and in the strong sense. In Subsection 2.3, we give higher bounds for the Lebesgue and Sobolev norms of this local solution. Eventually, these estimates allow us to extend the local solution to a global-in-time solution, and thus to prove Theorem 1 in Subsection 2.4. In the sequel, some results (Propositions 1, 2, 3 and 4) are either standard or mild adaptations of results which are proved elsewhere. We omit their proof here and refer to Subsection 2.2.5 in [26] for details.
2.1. Mild formulation of (1)
In this subsection, we collect preliminary results which shall enable us to provide a mild formulation of Equation (1), for which we prove the existence and uniqueness of a solution on a small interval.
2.1.1. Fractional Sobolev spaces
For all , let us define , and , . The family is a complete orthogonal basis of such that, for all , is on and . With respect to this basis, we define the fractional Sobolev space , for any , as the space of functions such that
[TABLE]
We take from [4, Appendice A] the following proposition and adapt it to our case of a flux function satisfying Assumption 1:
Proposition 1**.**
Under Assumption 1, for any , the mapping
[TABLE]
is bounded on bounded subsets of . Moreover, when or , it is Lipschitz continuous on bounded subsets of .
By virtue of Proposition 1, for all , we denote by and two finite constants such that:
- •
for all such that , ;
- •
for all such that , .
2.1.2. Heat kernel
Let us denote by the semigroup generated by the operator :
[TABLE]
Some of its properties are gathered in the following proposition.
Proposition 2** (Properties of the heat kernel).**
The semigroup satisfies the following properties.
- (1)
For any , for any , for any , and ; besides, the mapping is continuous on . 2. (2)
For all , there exists a constant such that
[TABLE] 3. (3)
For any , and , the process belongs to .
2.1.3. Stochastic convolution and mild formulation of (1)
Let be a normal filtration on the probability space and be a -Wiener process in with respect to this filtration. Given that the orthonormal basis of the space satisfies , the family is an orthonormal basis of . We set
[TABLE]
so that by (6), is a real-valued Brownian motion with variance . Next, we write
[TABLE]
Proposition 3**.**
Under Assumption 2, for all , the series
[TABLE]
converges in , and its sum defines an -adapted, -valued process almost surely continuous.
The process is called the stochastic convolution associated to the -Wiener process .
In the sequel, we let be a -stopping time, almost surely finite. We shall say that a process is -adapted if for all , the random variable is -measurable.
Definition 3** (Local mild solution).**
Let be an -measurable, -valued random variable. Under Assumptions 1 and 2, a (local) mild solution to the SPDE
[TABLE]
on is an -valued, -adapted process such that, almost surely:
- (1)
the mapping is continuous on ; 2. (2)
for all ,
[TABLE]
The combination of Propositions 1 and 2 ensures that all terms of the identity (11) are well-defined.
We now clarify the relationship between the notions of mild and strong solutions.
Proposition 4** (Mild and strong solutions).**
Under the assumptions of Definition 3, let be a mild solution to (10) on . If , then:
- (1)
for all , and the mapping is continuous on ; 2. (2)
for all ,
[TABLE]
Conversely, any -valued, -adapted process satisfying these two conditions almost surely is a mild solution to (10) on .
2.1.4. Existence and uniqueness of a mild solution on a small interval
For any integer , let us define
[TABLE]
where we recall that the constant is defined in Proposition 2, the constants and are defined after Proposition 1, and the constant is defined in (2).
Notice that , almost surely.
In the spirit of [9, 4], we obtain the existence and uniqueness of a mild solution to (10) on the "small" interval by a fixed-point argument.
Lemma 1** (Local existence and uniqueness).**
Let and be two -measurable random variables taking values respectively in and such that . Furthermore, let us set . Then, under Assumptions 1 and 2, there is a unique mild solution to (10) on .
Proof.
Let us introduce the random set
[TABLE]
Thanks to Propositions 2 and 3, we may define the random operator by
[TABLE]
and notice that any satisfies Equation (11) if and only if .
We first write, for some and for any ,
[TABLE]
On the one hand, by the first assertion of Proposition 2, ; on the other hand, we know thanks to the second assertion of Proposition 2 that
[TABLE]
furthermore, thanks to Proposition 1, if , then is bounded in uniformly in time, i.e. for all , . Thus,
[TABLE]
By definition of , it follows that whenever .
We now take . Then, for any ,
[TABLE]
where we have used the same arguments as above. Using now the Lipschitz continuity result in Proposition 1 and the definition of , we get for all ,
[TABLE]
meaning that is a contraction mapping on , which is complete. Then, by the Banach fixed-point theorem, admits a unique fixed point in . To show that this solution to Equation (11) is unique among all the -valued continuous processes, let us first notice that our choice of implies
[TABLE]
Assume that there is another solution of (11) not belonging almost surely to . Then we have with positive probability
[TABLE]
This means that the double inequality holds on some non-negligible event. On this event, the fixed-point argument also holds in the set
[TABLE]
which is formally a subset of . Thus, by uniqueness of the fixed point, we have and in particular , which is absurd. As a consequence, is the only -valued process with continuous trajectories satisfying Equation (11) on .
Finally, let and define the sequence of processes , by . It is clear from the definition of the operator and from Proposition 3 that each process is -adapted. On the other hand, the Banach fixed-point theorem asserts that almost surely, the sequence converges to in . As a consequence, for any , the sequence of -measurable random variables converges almost surely to , which makes this limit also -measurable. Thus, the process is -adapted. ∎
2.2. Construction of a maximal solution to (1)
In this subsection, we use the notions introduced in Subsection 2.1 to prove the following existence and uniqueness result for (1).
Lemma 2** (Existence and uniqueness result of a maximal solution to (1)).**
Under Assumptions 1 and 2, for any , there exists a pair such that:
- (1)
for any -stopping time such that almost surely, and , is the unique mild solution to (1) on ; 2. (2)
almost surely, or .
The random time is called the explosion time and the process is called the maximal solution to (1).
Proof.
Let . Let . By Lemma 1, Equation (1) possesses a unique mild solution on , where . We now define the filtration by
[TABLE]
and recall that the process defined by is a -Wiener process with respect to . Therefore, applying Lemma 1 again with this -Wiener process, and initial condition and , we obtain a mild solution of on , where . It is then easily checked that defining and for any , we obtain a unique mild solution to Equation (1) on .
We now proceed by induction and set for all ,
[TABLE]
where at each iteration we use Lemma 1 to extend the process to the unique mild solution of Equation (1) on . It is then clear that satisfies the first assertion of Lemma 2.
Since the sequence of integers is nondecreasing, if and only if there exists and such that, for all , . Hence, we can write
[TABLE]
However, by the strong Markov property, for any , the random variables , , are independent and identically distributed, and by the definition of , they are almost surely positive. As a consequence, by Borel’s 0-1 law,
[TABLE]
As the countable union of negligible events is still negligible, we get
[TABLE]
This implies that almost surely, if then , so that , which is the wanted result. ∎
2.3. Estimates on the maximal solution
Let . Let be the maximal solution to Equation (1) given by Lemma 2. By Proposition 4, is a continuous -valued process. Besides, Lemma 2 allows us to define, for any , the stopping time
[TABLE]
which always satisfies . In the sequel, we shall prove that , which shall imply that , almost surely.
Lemma 3**.**
Under Assumptions 1 and 2, for any and for all , we have:
[TABLE]
Moreover, there exist two constants depending only on , and such that
[TABLE]
Proof.
Let . We want to apply Itô’s formula on to the -valued process with the function . Since this process writes
[TABLE]
with , the standard formulation of Itô’s formula in Hilbert spaces [10, Theorem 4.32] requires at least to be continuous on , which is not the case for here. Hence, we shall proceed to approximate with a sequence of smooth functions , , apply Itô’s formula to the functions and then take the limit .
Step 1. Approximation of the -norm. Let be a function from to such that and whose support is contained in the interval . For any , we set the regularised Heaviside function and its antiderivative
[TABLE]
We now define a truncated -norm by setting
[TABLE]
The first differential and the second differential have the following expressions: ,
[TABLE]
[TABLE]
Step 2. Itô’s formula. First, let us notice that the process can be seen as an -valued -Wiener process where the operator has covariance
[TABLE]
Indeed, Assumption 2 ensures that and . We now have
[TABLE]
so that we can apply Itô’s formula [10, Theorem 4.32] for the real-valued process , which leads to
[TABLE]
Since the -norm of is bounded uniformly in time, the third term of the right-hand side is a square integrable martingale [10, Theorem 4.27]. Thus, for , integrating in time up to and taking the expectation, we get
[TABLE]
Step 3. Passing . We want now to pass to the limit . Regarding the left-hand side in the above equation, the family of functions is non-decreasing with respect to , so that the monotone convergence theorem yields
[TABLE]
For the flux term, we have almost surely, for all and for all , . Furthermore,
[TABLE]
Thus, the dominated convergence theorem applies and yields
[TABLE]
We now integrate by parts the viscous term:
[TABLE]
and this last integrand is dominated uniformly in by , where . Furthermore, thanks to (16), we have
[TABLE]
Thus, we get from the dominated convergence theorem,
[TABLE]
With similar computations, for the noise term, we have
[TABLE]
and
[TABLE]
Letting go to in (19), (20), (21) and (22), we get
[TABLE]
It turns out that the flux term disappears:
[TABLE]
where is an antiderivative of . As regards the noise coefficients, we have
[TABLE]
thanks to (3) and (4). As a consequence, we get from (23) the inequality
[TABLE]
Rewriting the integrand in the left-hand side, we get
[TABLE]
Since has a zero space average and is continuous in space (because it belongs to ), almost surely the function vanishes somewhere on the torus. Thus, we can apply the Poincaré inequality on the left-hand side which leads, after multiplying by on both sides, to the inequality
[TABLE]
For , we get
[TABLE]
and the claimed result for arbitrary follows by induction and from the inequalities and . ∎
Remark 1**.**
By Jensen’s inequality, the bound (18) also holds for any real number .
Lemma 4**.**
Under Assumptions 1 and 2, there exist two constants depending only on , , and , such that for all and all ,
[TABLE]
Proof.
We want to apply Itô’s formula to the squared -norm of the process . As for the proof of Lemma 3, we proceed by truncation of this function.
Step 1. Approximation of the -norm. We set
[TABLE]
The first differential and the second differential have the following expressions: ,
[TABLE]
[TABLE]
Step 2. Itô’s formula. Itô’s formula applied to yields almost surely and for all ,
[TABLE]
We first check that the third term of the right-hand side is a square-integrable martingale:
[TABLE]
Thus, taking the expectation, the stochastic integral disappears and we get
[TABLE]
On one hand, we can rewrite the viscous term as follows:
[TABLE]
On the other hand, applying Young’s inequality on the flux term, we get
[TABLE]
Injecting (30) and (31) into (29), we get the inequality
[TABLE]
Step 3. Passing . From Proposition 1, for any , there is a constant such that for all , we have
[TABLE]
Thus, we can use the dominated convergence theorem to let go to infinity in (32) and we get
[TABLE]
Since from Assumption 1, has polynomial growth, we can bound the second term of the right-hand side: using (2) and (17) with and , we get
[TABLE]
Applying now Lemma 3, we get
[TABLE]
Injecting this last bound in (33), we get the wanted result. ∎
Corollary 2** (Limit of ).**
Under Assumptions 1 and 2, almost surely, and thus almost surely.
Proof.
Let . Writing
[TABLE]
we get from Markov’s inequality,
[TABLE]
We apply now Lemma 4 to get
[TABLE]
Since has been chosen arbitrarily, it follows that almost surely, tends to as . Then, since , we have almost surely. ∎
2.4. Proof of Theorem 1
Under Assumptions 1 and 2, let , and be the maximal solution to Equation (1) given by Lemma 2. By Corollary 2, almost surely. Therefore, is the unique (global) mild solution to Equation (1), and by Proposition 4, it is also the unique (global) strong solution to this equation. It remains to check that this solution depends continuously on .
Lemma 5** (Continuous dependence on initial conditions).**
If is a sequence of satisfying
[TABLE]
then, denoting by the family of associated solutions, for any , we have almost surely
[TABLE]
Proof.
Let us fix a time horizon . Subtracting the mild formulations of and given by Proposition 4 and taking the -norm, we get by the triangle inequality and Proposition 2, for all ,
[TABLE]
Now, for any , we define the stopping times
[TABLE]
and we denote by , according to Proposition 1, the Lipschitz constant of the mapping over the centered ball in of radius . For an arbitrarily fixed , the inequality (34) implies
[TABLE]
In the next step, we iterate this last inequality and apply the Fubini theorem on the double time integral:
[TABLE]
However, by a change of variable, we have
[TABLE]
Hence, Grönwall’s lemma yields the following control
[TABLE]
It follows from this inequality that . Indeed, assuming the opposite, we would have (along a subsequence)
[TABLE]
which would imply
[TABLE]
Hence, necessarily, beyond a certain rank , we have
[TABLE]
Since the solutions of (7) do not explode, the stopping time tends almost surely to as tends to . As a consequence, there exists such that almost surely, so that for all ,
[TABLE]
Hence the result. ∎
3. Invariant measure
This section is dedicated to the proof of Theorem 2. The existence of an invariant measure is proven in Subsection 3.2 using the Krylov-Bogoliubov theorem, whereas the uniqueness is addressed through a coupling argument relying on the -contraction property established in Proposition 5.
The proof of existence of an invariant measure we provide in the next subsection relies plainly on the presence of viscosity. Indeed, the viscous term provides the process with a dissipative – and thus a more stable – behaviour. Still, it has to be borne in mind that when the flux term is nonlinear enough, the presence of a viscous term is not a necessary condition for the stability of the underlying stochastic process. On the physical side, in his theory of turbulent flows [25, 24], Kolmogorov already predicted this idea: the statistical distribution of scales of intermediate size in turbulence are not determined by the viscosity coefficient. On the theoretical side, the same idea was validated theoretically by powerful results on the invariant measure for the inviscid stochastic Burgers’ equation [18] and, quite a few years later, for inviscid stochastic conservation laws with "non-degenerate" flux [14]. However, our framework differs substantially from the inviscid case in the sense that our stability results are driven by regularity issues which cannot be tackled without viscosity.
3.1. Preliminary results
By Definition 2, an invariant measure for Equation (1) is a Borel probability measure on . Our proofs of existence and uniqueness however involve estimates in various spaces, namely , and . In particular, we shall manipulate and identify Borel probability measures on these spaces. We first clarify the relation between the associated Borel -fields thanks to the following result. For any metric space , we respectively denote by and the Borel -field and the set of Borel probability measures on .
Lemma 6** (Borel probability measures on and ).**
For all and , . As a consequence:
- (1)
for any , the mapping defines a Borel probability measure on ; 2. (2)
conversely, for any which gives full weight to , there exists a unique such that for any .
Proof.
Let and . The set defined by
[TABLE]
is a -field on , called the trace -field of in .
(1) We denote by the injection , so that . Since is continuous, and therefore Borel measurable, we have . Thus, for any , the pushforward measure defined by
[TABLE]
is a Borel probability measure on .
(2) Let us first notice that since is separable, the Borel -field is the smallest -field on containing all closed balls. Let be such a ball. Since the -norm is lower semi-continuous on , then is closed in as a level set of a lower semi-continuous function, and thus . It is then clear that , which by the minimality property of entails , and thus .
Now let be a Borel probability measure on which gives full weight to , that is to say such that there exists such that and . Let us define the Borel probability measure on by
[TABLE]
Notice that this definition is not ambiguous, because the identity ensures that any element of writes under the form for some ; besides, if are such that , then because the identity implies that . Finally, the fact that any such that for any needs to coincide with follows again from the identity . ∎
To prove Theorem 2, we will need a standard property of scalar conservation laws, namely the -contraction. In the stochastic setting, we mention that a similar proof of the following proposition is done in [5, Theorem 6.1], but in the case where the flux function is .
Proposition 5** (-contraction).**
Under Assumptions 1 and 2, let and be two strong solutions of (1) starting from different initial conditions and . Then, almost surely and for every , we have
[TABLE]
Proof.
We define a continuous approximation of the sign function by setting for all ,
[TABLE]
which gives rise to the following continuously differentiable approximation of the absolute value function:
[TABLE]
Let . We have
[TABLE]
We fix
[TABLE]
and we denote by a Lipschitz constant of over the interval . Since and belong to almost surely, then is finite almost surely and for all
[TABLE]
with
[TABLE]
Thus, we get from the dominated convergence theorem:
[TABLE]
As for the left-hand side of (35), noticing that increases to as decreases, we have from the monotone convergence theorem
[TABLE]
Hence, (35) yields the wanted result. ∎
3.2. Existence
From the semigroup introduced in Subsection 1.2, we define its time-averaged semigroup by , and for all ,
[TABLE]
[TABLE]
Following the first part of Lemma 6, for any and , we denote by the Borel probability measure on defined by .
Lemma 7**.**
Under Assumptions 1 and 2, for any , there exists an increasing sequence and a probability measure , such that the sequence of measures converges weakly to in .
Proof.
Let . From the inequality (17) with , we can pass to the limit (which we recall implies that almost surely), and we get for all ,
[TABLE]
Applying now the Markov inequality when , we have for all ,
[TABLE]
Setting
[TABLE]
we know from the compact embedding that the set is compact in . Thus, rewriting (37) as
[TABLE]
we deduce that the family of measures is tight in the space . The result is then a consequence of Prokhorov’s theorem [2, Theorem 5.1]. ∎
Lemma 8**.**
Under the assumptions of Lemma 7, for all , if is a random variable in distributed according to , then
[TABLE]
Besides, the probability measure associated with by the second part of Lemma 6 is invariant for the semigroup .
Proof.
We start to show that the measure gives full weight to . Thanks to Lemma 4, since almost surely, we have:
[TABLE]
Let be a sequence of -valued random variables such that and converges in distribution in towards a random variable . From (38) and the definition of , we have
[TABLE]
Now, since is lower semi-continuous on , we get from Portemanteau’s theorem:
[TABLE]
In particular, almost surely, and thus gives full weight to .
We now show that for any , . Let . From Lemma 3, we have for all ,
[TABLE]
Once again, we use Portemanteau’s theorem and the lower semi-continuity, this time of , on :
[TABLE]
and the wanted result follows.
To prove the invariance of the measure with respect to , we wish to apply the Krylov-Bogoliubov theorem [11, Theorem 3.1.1]. However, is a Feller semigroup on the space (Corollary 1) whereas our tightness result (Lemma 7) holds in . To overcome this inconvenience, we use Lemma 6 and we place ourselves at the level of the embedded probability measures in , where we can adapt, thanks to Proposition 5, the proof of [11, Theorem 3.1.1].
Let be associated with by the second part of Lemma 6, and let . In particular, the restriction is bounded and continuous on and we can write
[TABLE]
It follows from the -contraction property that the map is continuous with respect to the -norm. To prove this fact, let and let be a sequence of such that , . Let and , , be the strong solutions of (1) respectively with initial conditions and , . From Proposition 5, we get almost surely and for all ,
[TABLE]
Since is bounded and continuous with respect to the -norm, we have
[TABLE]
so that is continuous with respect to the -norm.
As a consequence, from Lemma 7, we have for all
[TABLE]
For any , gives full weight to and therefore, following the first part of Lemma 6, we can define the associated Borel probability measure on by . From Equation (39) and the above sequence of computations, it follows that for all ,
[TABLE]
Given that has been chosen arbitrarily in , this last equality says that . The second part of Lemma 6 now ensures that . ∎
3.3. Uniqueness
The proof of the uniqueness part of Theorem 2 follows the ideas of the "small-noise" coupling argument from Dirr and Souganidis [15]. On one hand, due to the dissipative nature of the drift, two solutions of (1) perturbed by the same noise and starting from different initial conditions are driven to balls of with small radius whenever this noise is small over sufficiently long time intervals. On the other hand, the -contraction property ensures that when these two solutions get close to one another they stay close forever. Hence, each time the noise gets small enough, the two solutions get closer and closer and eventually, they show the same asymptotical behaviour. This idea allows to show that the law of two solutions have the same limit as the time goes to infinity. Therefore, starting from two invariant measures leads to the equality of these measures. The same kind of argument was used in [14] for the invariant measure of kinetic solutions of inviscid scalar conservation laws and in [12] for the stochastic Navier-Stokes equations.
Let and be two solutions of (1) driven by the same -Wiener process . For all , we define the stopping time:
[TABLE]
Lemma 9**.**
Under Assumptions 1 and 2, there exists such that for any and in , the stopping time is finite almost surely.
Proof.
We can use here, from the statement of Lemma 3, the inequality (17) with . In this case, we get
[TABLE]
from which we deduce, by definition of the stopping time , that
[TABLE]
Taking yields
[TABLE]
from which we derive the wanted result. ∎
The following result asserts that when the coupled processes and start from deterministic initial conditions inside some ball of , then they both attain in finite time any neighbourhood of [math] with positive probability:
Lemma 10**.**
Under Assumptions 1 and 2, for any and any , there exist a time and a value such that for all satisfying ,
[TABLE]
Proof.
Let , be such that , and let us define
[TABLE]
To prove the lemma, we are going to compare the trajectories of and with the trajectories of their noiseless counterparts and , defined by
[TABLE]
Recall that the viscosity yields energy dissipation:
[TABLE]
Applying (3) on the right-hand side, we get
[TABLE]
and we can now apply Grönwall’s lemma:
[TABLE]
With our choice of , the above inequality means that as soon as , we have .
Furthermore, it is a consequence of Lemma 4 that satisfies
[TABLE]
Indeed, when all the noise coefficients are equal to zero, the constant in the statement of Lemma 4 can also be taken equal to zero. Since the same inequality also applies to , we have
[TABLE]
We focus now on the trajectories of the random processes and . We introduce the stopping time
[TABLE]
Following Proposition 4, we may use the expressions of and in the mild sense. From these mild formulations, we write
[TABLE]
where is the stochastic convolution associated with the -Wiener process . According to Proposition 1, we call a local Lipschitz constant of the map over the ball , and we place ourselves in the event
[TABLE]
where has been defined at Proposition 1. Taking , applying the second part of Proposition 2 and Proposition 1 to (40), we get
[TABLE]
Iterating this inequality and using the same arguments as in the proof of Lemma 5, we get for all ,
[TABLE]
Using now Grönwall’s lemma, we deduce
[TABLE]
Since the same arguments apply for the processes and , and given Equation (3), we have shown that for all ,
[TABLE]
We shall prove now that the event is impossible. Indeed, assume for instance that , then we would have
[TABLE]
and thus,
[TABLE]
which is false for too small values of .
We just have proven that for arbitrarily chosen and for all such that , we have
[TABLE]
To conclude the proof, it remains to check that
[TABLE]
We can write where is the closed ball of with radius . Since the process is the mild solution to the stochastic heat equation (i.e. Equation (7) with initial condition and flux ), we can apply the support theorem from [27, Theorem 1.1] which implies , so that (41) is satisfied. ∎
Lemma 11**.**
Under Assumptions 1 and 2, any invariant measure for the process solution to (1) is unique.
Proof.
Step 1. Almost sure confluence. We start by fixing small to which we associate the value defined at Lemma 10, where has been defined at Lemma 9. We define the increasing stopping time sequence
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Lemma 9 and the strong Markov property (Corollary 1) ensure that every is finite almost surely. We claim that
[TABLE]
Indeed, it is true for thanks to the strong Markov property and Lemma 10:
[TABLE]
and the general case follows by induction: assuming that inequality (42) is true for some , we have
[TABLE]
Taking the limit when goes to infinity, we get
[TABLE]
and consequently,
[TABLE]
Since and since the value has been chosen arbitrarily at the beginning of this proof, then Equality (43) means that almost surely,
[TABLE]
Recall however that Proposition 5 states that almost surely, the mapping is non-decreasing. It follows that almost surely,
[TABLE]
Step 2. Uniqueness. Let us now assume that there exist two invariant measures for the solution of (1), and let us take initial conditions and with distributions and respectively. For any test function bounded and Lipschitz continuous, we have for all ,
[TABLE]
Since is Lipschitz continuous, from (44), we have almost surely
[TABLE]
Moreover, for any , we have almost surely . Thus, we may apply the dominated convergence theorem, which yields
[TABLE]
so that , or in other words,
[TABLE]
According to Lemma 6, let and be the probability measures on associated to and respectively. Equation (45) rewrites
[TABLE]
so that and thus, by Lemma 6, . ∎
Proof of Theorem 2.
It follows from Lemmas 8 and 11. ∎
Acknowledgements
The authors would like to thank Sébastien Boyaval for fruitful discussions and for his careful reading of this manuscript.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] D. Alonso-Orán, A. Bethencourt de León, and S. Takao. The Burgers equation with stochastic transport: shock formation, local and global existence of smooth solutions. Nonlinear Differ. Equ. Appl. , page 26:57, 2019.
- 2[2] Patrick Billingsley. Convergence of probability measures . Wiley Series in Probability and Statistics: Probability and Statistics. John Wiley & Sons Inc., New York, second edition, 1999. A Wiley-Interscience Publication.
- 3[3] Alexandre Boritchev. Turbulence de Burgers en 1D : un cas modèle pour la théorie de Kolmogorov. Séminaire Laurent Schwartz — EDP et applications , pages 1–13, 2011-2012.
- 4[4] Alexandre Boritchev. Generalised Burgers equation with random force and small viscosity . Theses, Ecole Polytechnique X, October 2012.
- 5[5] Alexandre Boritchev. Sharp Estimates for Turbulence in White-Forced Generalised Burgers Equation. Geometric and Functional Analysis , 23(6):1730–1771, 2013.
- 6[6] Sébastien Boyaval, Sofiane Martel, and Julien Reygner. Finite-volume approximation of the invariant measure of a viscous stochastic scalar conservation law. Preprint, ar Xiv:1909.08899, 2019.
- 7[7] L. Breiman. Probability . Addison-Wesley series in statistics. Addison-Wesley Publishing Company, 1968.
- 8[8] Gui-Qiang G. Chen and Peter H. C. Pang. Invariant measures for nonlinear conservation laws driven by stochastic forcing. Chinese Annals of Mathematics, Series B , 40(6):967–1004, 2019.
