Yves Achdou
email: achdou@ljll-univ-paris-diderot.fr

Alessio Porretta
email: porretta@mat.uniroma2.it

MEAN FIELD GAMES WITH CONGESTION

des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Recently, an important research activity on mean field games (MFGs for short) has been initiated since the pioneering works [START_REF] Lasry | Jeux à champ moyen. I. Le cas stationnaire[END_REF][START_REF] Lasry | Jeux à champ moyen. II. Horizon fini et contròle optimal[END_REF][START_REF] Lasry | Mean field games[END_REF] of Lasry and Lions: it aims at studying the asymptotic behavior of stochastic differential games (Nash equilibria) as the number n of agents tends to infinity. In these models, it is assumed that the agents are all identical and that an individual agent can hardly influence the outcome of the game. Moreover, each individual strategy is influenced by some averages of functions of the states of the other agents. In the limit when n → +∞, a given agent feels the presence of the others through the statistical distribution of the states. Since perturbations of a single agent's strategy does not influence the statistical states distribution, the latter acts as a parameter in the control problem to be solved by each agent. When the dynamics of the agents are independent stochastic processes, MFGs naturally lead to a coupled system of two partial differential equations (PDEs for short), a forward Kolmogorov or Fokker-Planck equation and a backward Hamilton-Jacobi-Bellman equation, see for example (1.1) below.

The theory of MFGs allows one to model congestion effects, i.e. situations in which the cost of displacement of the agents increases in those regions where the density is large. MFGs models including congestion were introduced and studied in [26]. A typical such model leads to the following system of PDEs:

(1.1)

       -∂ t u -ν∆u + 1 β |Du| β (m+µ) α = F (t, x, m) , (t, x) ∈ (0, T ) × Ω ∂ t m -ν∆m -div(m |Du| β-2 Du (m+µ) α ) = 0 , (t, x) ∈ (0, T ) × Ω m(0, x) = m 0 (x) , u(T, x) = G(x, m(T )) , x ∈ Ω ,
with ν > 0, α > 0, β ∈ (1, 2], µ ∈ R with either µ > 0 or µ = 0. System (1.1) must generally be complemented with suitable boundary conditions on (0, T )×∂Ω, but we will avoid the additional technical difficulties coming from the latter by focusing on the case when Ω is the flat torus, i.e. Ω = T N = R N /Z N and all the data are periodic. Loosely speaking, (1.1) describes the optimization over a stochastic dynamics defined on a standard probability space (X , F, F t , P)

dX t = w t dt + √ 2ν dB t 1
where B t is a N -dimensional Brownian motion, with a cost criterion given by

(1.2) inf β E T 0 c β (m t + µ) γ |w t | β β-1 + F (t, X t , m t ) dt + E(G(X T , m T )) ,
where γ = α β-1 and c β is a suitable normalization constant. In the viewpoint of the generic agent, m t = m(t, X t ) is meant to represent the distribution law of the states, however in the optimization process it is just, a priori, a given frozen density function. The mean field game equilibrium is next given through a fixed point scheme, by requiring, a posteriori, that m t coincides with the probability density in R N of the law of the optimal process X t .

In [26], P-L. Lions put the stress on general stuctural conditions yielding the uniqueness for the following MFG systems with local coupling:

-∂ t u -ν∆u + H(t, x, m, Du) = F (t, x, m) , (t, x) ∈ (0, T ) × Ω (1.3)
∂ t m -ν∆m -div(mH p (t, x, m, Du)) = 0 , (t, x) ∈ (0, T ) × Ω (1.4) m(0, x) = m 0 (x) , u(T, x) = G(x, m(T )) , x ∈ Ω (1.5) namely that F and G be increasing w.r.t. m and that the following matrix be positive semidefinite:

( In the present work, we will show that this hypothesis yields both the existence and the uniqueness of weak solutions. Except for situations in which special tricks may be applied (stationary problems and quadratic Hamiltonian, see [START_REF] Gomes | Existence for stationary mean field games with congestion and quadratic Hamiltonians[END_REF]), the existence of classical solutions of suitable generalizations of (1.1) seems difficult to obtain, because generally neither upper bounds on m nor strict positivity of m are known unless one restricts the growth conditions for the nonlinearities and assumes the time horizon T to be small, see [START_REF] Gomes | Short-time existence of solutions for mean-field games with congestion[END_REF], [START_REF] Graber | Weak solutions for mean field games with congestion[END_REF] (see also [START_REF] Evangelista | On the existence of solutions for stationary mean-field games with congestion[END_REF] for the stationary case). Therefore, in order to get at a sufficiently general result, we aim at proving the existence and uniqueness of suitably defined weak solutions.

For MFG models without congestion, the first results on the existence of weak solutions were supplied in [START_REF] Lasry | Mean field games[END_REF]. Besides, as already observed in [START_REF] Lasry | Jeux à champ moyen. I. Le cas stationnaire[END_REF][START_REF] Lasry | Jeux à champ moyen. II. Horizon fini et contròle optimal[END_REF][START_REF] Lasry | Mean field games[END_REF], in the easiest cases, the system of PDEs can be seen as the optimality conditions of a problem of optimal control driven by a PDE: in such cases, a pair of primal-dual optimization problems can be introduced, leading to a suitable weak formulation for which there exists a unique solution, see [START_REF] Cardaliaguet | Second order mean field games with degenerate diffusion and local coupling[END_REF], where possibly degenerate diffusions are dealt with. A striking fact is that in general, MFGs with congestion cannot be cast into an optimal control problem driven by a PDE, by contrast with simpler cases. For MFG systems (1.3)-(1.5) with H independent of m, a complete analysis is available in [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF], which contains in particular new results on weak solutions of Fokker-Planck equations, and an answer to the delicate question of the uniqueness of weak solutions of MFG systems. Proving uniqueness of weak solutions is difficult for the following reason: to compare two solutions (u, m) and (ũ, m) of (1.1), the main idea is to test the Bellman equations by m -m and the Kolmogorov equations by u-ũ and sum the resulting identities. While this is of course permitted for classical solutions, special care is needed for weak solutions, because the PDEs only hold in the distributional sense. The present work will borrow several ideas and results from [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF].

Let us mention that other kinds of models, which may also include congestion, are obtained by assuming that all the agents use the same distributed feedback strategy and by passing to the limit as N → ∞ before optimizing the common feedback. Given a common feedback strategy, the asymptotics are given by the McKean-Vlasov theory, [START_REF] Mckean | A class of Markov processes associated with nonlinear parabolic equations[END_REF]: the dynamics of a given agent is found by solving a stochastic differential equation with coefficients depending on a mean field, namely the statistical distribution of the states, which may also affect the objective function. Since the feedback strategy is common to all agents, perturbations of the latter affect the mean field. Then, letting each player optimize its objective function amounts to solving a control problem driven by the McKean-Vlasov dynamics. The latter is named control of McKean-Vlasov dynamics by R. Carmona and F. Delarue [START_REF] Carmona | Control of McKean-Vlasov dynamics versus mean field games[END_REF][START_REF] Carmona | Mean field forward-backward stochastic differential equations[END_REF] and mean field type control by A. Bensoussan et al, [START_REF] Bensoussan | Control and Nash games with mean field effect[END_REF][START_REF] Bensoussan | Mean field games and mean field type control theory[END_REF]. By contrast with MFGs, this genuine connection to optimal control problems driven by a PDE makes it possible to use techniques from the calculus of variations and define a suitable notion of weak solutions for which existence and uniqueness can be proved, see [START_REF] Achdou | Mean Field Type Control with Congestion[END_REF]. Finally, let us mention that numerical methods and simulations for MFGs and mean field type control with congestion are dealt with in [START_REF] Achdou | Finite difference methods for mean field games, Hamilton-Jacobi equations: approximations, numerical analysis and applications[END_REF][START_REF] Achdou | On the system of partial differential equations arising in mean field type control[END_REF]. The present paper is devoted to prove the existence and uniqueness of weak solutions to a class of MFG systems which generalize (1.1) including the congestion effects in the structure conditions of the Hamiltonian function. Making reference to the model example (1.1), we will consider both the case when µ > 0 and the case when µ = 0. Let us notice that the congestion effect is essentially contained in the behavior of the Hamiltonian for large values of m, so both cases could be considered as congestion models, since the difference between µ > 0 and µ = 0 will be important in those regions of vanishing density. Hereafter, the case that µ = 0 will be referred to as the case of singular congestion, since in this case the Hamiltonian may not even be defined at m = 0. In the latter case, the Lagrangian vanishes at m = 0, and this vanishing of the cost criterion will lead to a slightly incomplete information on the equation solved by the value function, and therefore to a relaxed formulation of the notion of weak solution.

The work is organized as follows: in Section 2 we set the problem and introduce the assumptions that will be used throughout the paper. In Section 3 we define the notions of weak solutions and state the main existence and uniqueness results: as explained above, we will give a different definition of solution in the case µ = 0, in which special care is required because the Hamiltonian has a singularity at m = 0. In Section 4, we consider non singular Hamiltonians (µ > 0) and give the first steps of the existence proof which consist of studying a sequence of suitably regularized problems. Section 5 contain lemmas which will turn crucial for proving uniqueness but also for identifying the limit of the aforementioned regularized problems. In Section 6, in the case µ > 0, we conclude the proof of existence and obtain uniqueness under suitable additional assumptions. Section 7 is devoted to the singular case µ = 0.

Running assumptions

Let us consider the general MFG system (1.3)- (1.5). Different type of boundary conditions could be considered but for simplicity we neglect this point by assuming that the equation takes place in a standard flat torus Ω = T N = R N /Z N and all functions are Z N -periodic in x. We define Q T = (0, T ) × Ω.

In (1.3)-(1.5), the function H(t, x, m, p) is assumed to be measurable with respect to (t, x), continuous with respect to m and C 1 with respect to p. The notation H p stands for ∂H(t,x,m,p) ∂p . We also assume that H is a convex function of p.

We modify the structure growth conditions introduced in [START_REF] Lasry | Mean field games[END_REF] to take into account the congestion factor. So, we will assume that H satisfies, for some positive β and α H(t, x, m, 0) ≤ 0 , (2.1)

H(t, x, m, p) ≥ c 0 |p| β (m + µ) α -c 1 (1 + m α β-1 ) , (2.2) |H p (t, x, m, p)| ≤ c 2 (1 + |p| β-1 (m + µ) α ) , (2.3) H p (t, x, m, p) • p ≥ (1 + σ) H(t, x, m, p) -c 3 (1 + m α β-1 ) , (2.4)
for a.e. (t, x) ∈ Q T and every (m, p) ∈ R + × R N , where σ, c 0 , . . . , c 3 are positive constants and µ ∈ R, µ ≥ 0. For the sake of clarity, we will separately consider the two cases µ > 0 and µ = 0.

Notice that (2.1) can be alternatively rephrased as: H(t, x, m, 0) is bounded from above; indeed, if H(t, x, m, 0) ≤ K for some K > 0, we can reduce to (2.1) by changing

u into u + K(T -t).
As for the ranges of β, α in the previous assumptions, we will assume that

(2.5) 1 < β ≤ 2 , 0 < α ≤ 4(β -1) β .
The first condition on β is meant to exclude both the case of linear or sub linear growth and the case of super quadratic growth. Actually, if β ≤ 1, then the bound on H p would imply that m is bounded, and this would allow for a huge simplification in the system (at least if µ > 0), in particular the difficult issues coming from the congestion effect would be mostly circumvented and this case would be rather similar to the results which already exist in the literature. The super quadratic case is different in nature since one cannot rely any more on classical tools for parabolic equations and a different approach should be employed, as in the first order case. Since this would bring us too far for the present paper, we decided to assume that β ≤ 2.

As for (2.5), this structural condition seems natural in order to have a completely well-posed formulation in the congestion model. As already mentioned, this limitation is indeed required in order to obtain uniqueness, as pointed out by P.L. Lions [26], when he discussed general conditions for the uniqueness of mean field games systems. In addition, this condition also seems to be useful at other stages of our existence proof. The upper bound required on α can be interpreted as a tolerable growth condition on the cost function which penalizes the motion in the zones where the distribution density is high (since we have γ = α β-1 in the Lagrangian model (1.2)).

Finally, the case α = 0 is excluded in the above conditions since this is already treated in [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF].

Remark 2.1. Straightforward calculus leads to 4/β ≤ β: from this and (2.5), we see that α ≤ β, and that α = β can occur only if

β = α = 2.
Notice that combining (2.1) and the convexity of H in the p variable yields that

(2.6) H p (t, x, m, p) • p -H(t, x, m, p) ≥ -H(t, x, m, 0) ≥ 0.
Moreover, since 1 < β ≤ 2, and from (2.2)-(2.4), we can deduce that there exist nonnegative constants

C 0 , C 1 , C 2 such that (2.7) m α β-1 +1 |H p | 2 -C 0 ≤ C 1 m {H p • p -H + C 2 } .
Using generic non negative constants C 0 , C 1 that can change from line to line, (2.7) is obtained as follows:

m α β-1 |H p | 2 ≤ C 1 m α β-1 + |p| 2(β-1) (m + µ) 2α-α β-1 ≤ C 1 m α β-1 + c 0 |p| β (m + µ) α + C 0 (m + µ) α β-1 ≤ H + C 1 (m α β-1 + 1) ≤ σ -1 (H p • p -H) + C 1 + C 0 m α β-1
where we have used (2.3) in the first line, 1 < β ≤ 2 and Young's inequality in the second line, (2.2) in the third line and (2.4) in the fourth line.

We also stress that, if H satisfies (2.1)-(2.4), then H + b • p satisfies the same conditions for any bounded vector field b. Hence the addition of bounded drift terms in the dynamics is permitted by the above setting of assumptions.

Remark 2.2. Let us compare the above assumptions with previous settings used for congestion models in mean field games. In [START_REF] Gomes | Short-time existence of solutions for mean-field games with congestion[END_REF], congestion models are presented starting from the Lagrangian function

L(t, x, m, q) = m γ L 0 (t, x, q -b(t, x)) + F (t, x, m), with L 0 (t, x, q) = a(t, x)(1 + |q| 2 ) β 2 ,
where β = β/(β -1), a is a smooth positive function and b is a smooth velocity field. This leads to the system

     -∂ t u -ν∆u + m γ H(t, x, Du m γ ) = F (t, x, m) , (t, x) ∈ (0, T ) × Ω ∂ t m -ν∆m -div(mH p (t, x, Du m γ )) = 0 , (t, x) ∈ (0, T ) × Ω m(0, x) = m 0 (x) , u(T, x) = G(x, m(T )) ,
x ∈ Ω where H(t, x, p) = H 0 (t, x, p) + b(t, x) • p, and H 0 (t, x, p) = sup q -q • p -L 0 (t, x, q).

If we set α = γ(β -1), we recover our assumptions if γ ≤ 4/β which is equivalent to saying that ∇ p H 0 (t, x, p) • p -H 0 (t, x, p) ≥ γ 4 p T D pp H 0 (t, x, p)p for all p = 0, as assumed in [START_REF] Gomes | Short-time existence of solutions for mean-field games with congestion[END_REF].

Finally, let us make precise the assumptions on the functions F, G appearing in (1.3)-(1.5). The function F is assumed to be measurable with respect to (t, x) ∈ Q T and continuous with respect to m. In addition, we assume that there exists a nondecreasing function f (s) such that f (s)s is convex and

(2.8) λ f (m) -κ ≤ F (t, x, m) ≤ 1 λ f (m) + κ ∀m ∈ R + ,
for some real numbers λ, κ > 0.

The terminal cost G is assumed to be measurable with respect to x ∈ Ω, continuous with respect to m and such that (2.9) m → G(x, m) is nondecreasing.

In addition, as for F , we assume that there exists a function g(s) such that g(s)s is convex and

(2.10) λ g(m) -κ ≤ G(x, m) ≤ 1 λ g(m) + κ ∀m ∈ R + ,
with for example the same constants λ, κ > 0 as in (2.8).

Let us notice that, as a consequence of assumptions (2.8) ) are the structural conditions under which we will prove the existence of weak solutions to the MFG system with congestion. Eventually, we will need a further condition for H and F in order to have uniqueness. This will be a natural reformulation of condition (1.6). Precisely, we assume the following: for any m ≥ 0, z ≥ -m, p, r ∈ R N such that (z, r) = (0, 0), setting m s = m + sz and p s = p + sr for s ∈ [0, 1], then the function

h : s → -zH(t, x, m s , p s ) + m s H p (t, x, m s , p s ) • r + zF (t, x, m s ) is (strictly) increasing on [0, 1].
(2.13) Remark 2.3. Note that if H(t, x, m, p) = (µ + m) -α |p| β , with (α, β) satisfying (2.5) and if F is nondecreasing with respect to m, then h defined in (2.13) is increasing. Condition (2.5) is thus compatible with the assumption yielding uniqueness.

Finally, we assume m 0 ∈ C(Ω), m 0 ≥ 0 and, in order to be consistent with the interpretation of m, we fix the normalization condition Ω m 0 = 1.

Main results

Case of non singular congestion.

We assume here that conditions (2.1)-(2.4) hold with µ > 0. Let us first make our notion of weak solution precise.

Definition 3.1. A pair (u, m) ∈ L 1 (Q T ) × L 1 (Q T ) is a weak solution to (1.3)-(1.5) if (i) m ∈ C([0, T ]; L 1 (Ω)
), m ≥ 0, u ∈ L q (0, T ; W 1,q (Ω)) for every q < N +2 N +1 , (ii)

F (t, x, m)m ∈ L 1 (Q T ) , G(x, m(T ))m(T ) ∈ L 1 (Ω) , m |Du| β (m + µ) α ∈ L 1 (Q T ) , |Du| β (m + µ) α ∈ L 1 (Q T ) , m ∈ L 1+ α β-1 (Q T ) (iii) u ∈ L ∞ (0, T ; L 1 (Ω))
is bounded below and is a solution of the Bellman equation in the sense of distributions:

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H(t, x, m, Du)ϕ dxdt = T 0 Ω F (t, x, m)ϕ dxdt + Ω G(x, m(T ))ϕ(T ) dx (3.1)
for every ϕ ∈ C ∞ c ((0, T ] × Ω). (iv) m is a solution of the Kolmogorov equation:

(3.2) T 0 Ω m {-ϕ t -ν∆ϕ + H p (t, x, m, Du)Dϕ} dxdt = Ω m 0 ϕ(0) dx for every ϕ ∈ C ∞ c ([0, T ) × Ω).
Remark 3.2. We notice that, on account of conditions (ii), (3.1) implies that u ∈ L 1 (Q T ) solves, in weak sense,

-∂ t u -ν∆u = f (t, x) , (t, x) ∈ (0, T ) × Ω u(T, x) = g(x) , x ∈ Ω for some f ∈ L 1 (Q T ), g ∈ L 1 (Ω). It follows that u ∈ C 0 ([0, T ]; L 1 (Ω))
and, in addition, u is a renormalized solution of the same equation (see e.g. [30, Proposition 2.1]). Indeed, the condition that u belongs to L ∞ (0, T ; L 1 (Ω)) is actually redundant in (iii), since it follows from the above result.

As stated in the previous remark, a solution u of equation (3.1) belongs to C 0 ([0, T ]; L 1 (Ω)). However we also need to work with merely subsolutions of the same equation, for which this kind of continuity may not hold. We recall a lemma, whose proof can be found in [START_REF] Cardaliaguet | Second order mean field games with degenerate diffusion and local coupling[END_REF][START_REF] Achdou | Mean Field Type Control with Congestion[END_REF], which establishes continuity for subsolutions in a weaker sense.

Lemma 3.3. Let u ∈ L ∞ (0, T ; L 1 (Ω)) satisfy, for some function h ∈ L 1 (Q T ) and k ∈ L 1 (Ω), T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt ≤ T 0 Ω hϕ dxdt + Ω k(x)ϕ(T, x) dx,
for every nonnegative function ϕ ∈ C ∞ c ((0, T ] × Ω). Then for any Lipschitz continuous map ξ : Ω → R, the map t → Ω ξ(x)u(t, x)dx has a BV representative on [0, T ]. Moreover, if we note Ω ξ(x)u(t + , x)dx its right limit at t ∈ [0, T ), then the map ξ → Ω ξ(x)u(t + , x)dx can be extended to a bounded linear form on C(Ω).

Remark 3.4. As a consequence of Lemma 3.3, for any subsolution u of (3.1) we can define u(0 + ) as a bounded Radon measure on Ω. For simplicity, we note u(0) = u(0 + ). Remark 3.5. A weaker definition than Definition 3.1 is possible, by replacing (ii) by G(x, m(T

)) ∈ L 1 (Ω), F (t, x, m) ∈ L 1 (Q T ), H(t, x, m, Du) ∈ L 1 (Q T ) and m H p (t, x, m, Du) ∈ L 1 (Q T ).
In this case, we would need to derive firstly the regularity (ii) in order to prove uniqueness and this could be itself a delicate step. It does not seem restrictive to ask directly (ii) in the formulation since such regularity is actually obtained in the existence result.

Finally, let us state below our main result, where we prove existence and uniqueness of weak solutions under the above assumptions. Theorem 3.6. Under Assumptions (2.1)-(2.4) with µ > 0, (2.5) and (2.8)-(2.10), there exists a weak solution of problem (1.3)-(1.5). Furthermore, if the assumption (2.13) is satisfied, then there is a unique weak solution of problem (1.3)-(1.5).

3.2.

Case of singular congestion. Here we consider the limit case µ = 0 in assumptions (2.2)-(2.3), which corresponds to the singular congestion model H = |Du| β m α . Due to the singularity, we will only find a relaxed solution. Before defining weak solutions, we notice that we can always change H(t, x, m, p) into H(t, x, m, p)+ c 1 and u into u -c 1 (T -t); therefore, in view of (2.2), we may assume without restriction that

(3.3) H(t, x, 0, p) ≥ 0.
This allows us to define weak solutions as follows:

Definition 3.7. A pair (u, m) ∈ L 1 (Q T ) × L 1 (Q T ) + is a weak solution to (1.3)-(1.5) if (i) m ∈ C([0, T ]; L 1 (Ω)), u ∈ L q (0, T ; W 1,q (Ω)) for every q < N +2 N +1 , (ii)
F (t, x, m)m ∈ L 1 (Q T ) , G(x, m(T ))m(T ) ∈ L 1 (Ω) , m 1 m>0 |Du| β m α ∈ L 1 (Q T ) , 1 m>0 |Du| β m α ∈ L 1 (Q T ) , m ∈ L 1+ α β-1 (Q T ), Du = 0 a.e. in {m = 0}, (3.4) (iii) u ∈ L ∞ (0, T ; L 1 (Ω))
is bounded below and is a subsolution of the Bellman equation:

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H(t, x, m, Du)1 {m>0} ϕ dxdt ≤ T 0 Ω F (t, x, m)ϕ dxdt + Ω G(x, m(T ))ϕ(T ) dx (3.5) for every nonnegative ϕ ∈ C ∞ c ((0, T ] × Ω).
(iv) m is a solution of the Kolmogorov equation:

(3.6) T 0 Ω m -ϕ t -ν∆ϕ + H p (t, x, m, Du)1 {m>0} Dϕ dxdt = Ω m 0 ϕ(0) dx for every ϕ ∈ C ∞ c ([0, T ) × Ω). (v)
The following identity is satisfied:

Ω m 0 u(0) dx = Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m)m dxdt + T 0 Ω m [H p (t, x, m, Du) • Du -H(t, x, m, Du)] 1 {m>0} dxdt (3.7)
where the first term is understood as the trace of Ω u(t) m 0 dx in BV (0, T ), in view of Lemma 3.3.

Notice that, if α ≤ 1, we actually have m 1-α |Du| β ∈ L 1 (Q T ) and the restriction to the set {m > 0} in (3.6)-(3.7) is not needed. However, for the case α > 1 this restriction applies.

We state below the existence and uniqueness result which we obtain for the case of singular congestion.

Theorem 3.8. Assume that (2.1)-(2.4) hold with µ = 0 and with either β < 2 and 0 < α ≤ 4(β-1) β or β = 2 and 0 < α < 2, and that (2.8)-(2.10) hold true (assuming furthermore the nonrestrictive condition (3.3)). Then there exists a weak solution of problem (1.3)-(1.5). Furthermore, if condition (2.13) holds true for any m > 0, z > -m, p, r ∈ R N such that (z, r) = (0, 0), and if

(3.8) m > 0 ⇒ F (t, x, m) > F (t, x, 0),
then there is a unique weak solution of problem (1.3)-(1.5).

Non singular congestion: approximations of (1.3)-(1.5)

For µ > 0, we consider the following system of PDEs:

-∂ t u -ν∆u + H(t, x, T 1/ m , Du ) = F (t, x, m ) , (t, x) ∈ (0, T ) × Ω (4.1) ∂ t m -ν∆m -div(m H p (t, x, T 1/ m , Du )) = 0 , (t, x) ∈ (0, T ) × Ω (4.2) m (0, x) = m 0 (x) , u (T, x) = G (x, m (T )) , x ∈ Ω (4.3)
where

T 1/ m = min(m, 1/ ), F (t, x, m) = ρ F (t, •, ρ m))(x), G (x, m) = ρ G(•, ρ m))(x), m 0 = ρ m 0 .
Here denotes the convolution with respect to the spatial variable and ρ is a

standard symmetric mollifier, i.e. ρ (x) = 1 N ρ( x ) for a nonnegative function ρ ∈ C ∞ c (R N ) such that R N ρ(x)dx = 1. Lemma 4.1.
There exists a weak solution (u , m ) of (4.1)-( 4.3) such that u ∈ C 1+α,1/2+α/2 ( QT ), m ∈ C α,α/2 ( QT ), and m ≥ 0, Ω m (t, x)dx = Ω m 0 (x)dx, ∀t.

Proof. We set X = m ∈ C( QT ) : m ≥ 0 and Ω m(t, x)dx = Ω m 0 (x)dx, ∀t . For any m ∈ X, the boundary value problem

-∂ t u -ν∆u + H(t, x, T 1/ m, Du) = F (t, x, m) , (t, x) ∈ (0, T ) × Ω (4.4) u(T, x) = G (x, m(T )) ,
x ∈ Ω (4.5) has a unique solution u ∈ L ∞ (0, T ; W 1,∞ (Ω)) (see e.g. [START_REF] Ladyženskaja | Linear and quasilinear equations of parabolic type[END_REF], Chapter 5, Theorem 6.3) and u L ∞ (0,T ;W 1,∞ (Ω)) is bounded independently of m in X. This follows e.g. from the C 1,α estimates in [START_REF] Ladyženskaja | Linear and quasilinear equations of parabolic type[END_REF], Chapter 5, Theorem 3.1, since F is bounded and

|H(t, x, T 1/ m, Du)| ≤ C (1 + |Du| 2 ) for some constant C > 0. In addition, the map m → u is continuous from X to L ∞ (0, T ; W 1,∞ (Ω)).
Thus, the boundary value problem

∂ t m -ν∆ m -div( mH p (t, x, T 1/ m , Du)) = 0 , (t, x) ∈ (0, T ) × Ω m(0, x) = m 0 (x) , x ∈ Ω has a unique solution m ∈ X ∩ C α,α/2 ( QT ), which is bounded in C α,α/2 ( QT )
uniformly with respect to m ∈ X, see e.g. [START_REF] Ladyženskaja | Linear and quasilinear equations of parabolic type[END_REF], Chapter 3, Theorem 10.1. Moreover the map m → m is continuous from X to C α,α/2 ( QT ). Hence, Schauder's theorem implies that the map m → m has a fixed point m .

Lemma 4.2. Let (u , m ) be a solution of system (4.1)-(4.3). Then

u (t, x) ≥ c 4 , (4.6) Ω G (x, m (T ))m (T ) + T 0 Ω F (t, x, m )m dxdt + (T 1/ m ) α β-1 +1 L N +2 N (Q T ) ≤ C, (4.7) T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt ≤ C, (4.8) T 0 Ω |Du | β (T 1/ m + µ) α dxdt + T 0 Ω m |Du | β (T 1/ m + µ) α dxdt ≤ C, (4.9) u L ∞ (0,T ;L 1 (Ω)) ≤ C, (4.10)
for a positive constant C = C(T, H, F, G, m 0 ∞ ) and where c 4 is the constant appearing in (2.11) and (2.12).

Proof. We multiply the equation of u by m , the equation of m by u , integrate by parts some terms and subtract; we get

Ω u (0)m 0 dx - Ω G (x, m (T ))m (T ) - T 0 Ω m H p (t, x, T 1/ m , Du ) • Du dxdt + T 0 Ω m H(t, x, T 1/ m , Du ) dxdt = T 0 Ω F (t, x, m )m dxdt, which implies that Ω G (x, m (T ))m (T ) + T 0 Ω F (t, x, m )m dxdt + Ω u (0) -m 0 dx + T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt ≤ C m 0 L ∞ (Ω) Ω [u (0)] + dx . (4.11)
By using (2.1) and F, G ≥ c 4 , we see by comparison that u (t) ≥ c 4 + c 4 (T -t), and therefore that there exists an absolute constant C such that (4.12)

[u (0)] -≤ C.

Integrating (4.1) and using (2.2) and (4.12), we obtain that

Ω [u (0)] + dx + c 0 T 0 Ω |Du | β (T 1/ m + µ) α dxdt ≤c 1 T 0 Ω (1 + (T 1/ m ) α β-1 )dxdt + T 0 Ω F (t, x, m )dxdt + Ω G (x, m (T ))dx + C =c 1 T 0 Ω (1 + (T 1/ m ) α β-1 )dxdt + T 0 Ω F (t, x, ρ m )dxdt + Ω G(x, ρ m (T ))dx + C (4.13)
From Assumption (2.8) we see that, for any choice of L > 0, (4.14)

F (t, x, m) ≤ f L (t, x) + m L (F (t, x, m) + |c 4 |),
where c 4 is the constant in (2.11) and

f L (t, x) = max 0≤m≤L F (t, x, m) ≤ 1 λ max 0≤m≤L f (m) + κ, from (2.8). Similarly, G(x, m) =≤ g L (x) + m L (G(x, m) + |c 4 |) where g L (x) = max 0≤m≤L G(x, m) ≤ 1 λ max 0≤m≤L g(m)
+ κ, from (2.10). Therefore, (4.13) implies

Ω [u (0)] + dx + c 0 T 0 Ω |Du | β (T 1/ m + µ) α dxdt ≤ c 1 T 0 Ω (1 + (T 1/ m ) α β-1 )dxdt + 1 L T 0 Ω m F (t, x, m )dxdt + Ω m (T )G (x, m (T ))dx + C L , (4.15)
where C L is a constant depending on L. Therefore, by choosing L sufficiently large (depending on m 0 L ∞ (Ω) ) we can combine (4.11) and (4.15) and obtain

Ω G (x, m (T ))m (T ) + T 0 Ω F (t, x, m )m dxdt + T 0 Ω |Du | β (T 1/ m + µ) α dxdt + T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt ≤ C + C T 0 Ω (T 1/ m ) α β-1 dxdt . (4.16)
We now use (4.2) multiplied by (T 1/ m ) α β-1 . We get

∂ t m (T 1/ m ) α β-1 dxdt + να β -1 m <1/ |Dm | 2 (m ) α β-1 -1 dxdt = - α β -1 m <1/ (m ) α β-1 H p (t, x, m , Du ) • Dm dxdt . Calling I (z) = z 0 (T 1/ y) α β-1 dy, we see that t 0 Ω ∂ t m (T 1/ m ) α β-1 dxds = Ω I (m (t))dx - Ω I (m 0 )dx.
Therefore, by Young's inequality sup

t Ω I (m (t))dx - Ω I (m 0 )dx + m <1/ |Dm | 2 (m ) α β-1 -1 dxdt ≤ C m <1/ (m ) α β-1 +1 |H p (t, x, m , Du )| 2 dxdt .
(4.17)

Since I (z) ≥ 1 α β-1 +1 (T 1/ z) α β-1 +1 we use (2.7) to obtain sup t∈(0,T ) Ω (T 1/ m (t)) α β-1 +1 dx + T 0 Ω DT 1/ m 2 (T 1/ m ) α β-1 -1 dxdt ≤ C m <1/ m {H p (t, x, m , Du ) • Du -H(t, x, m , Du )} dxdt + C 1 + m <1/ (m ) α β-1 +1 dxdt + C Ω (m 0 ) α β-1 +1 dx ≤ C T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt + C 1 + T 0 Ω (T 1/ m ) α β-1 +1 dxdt + C Ω (m 0 ) α β-1 +1 dx, (4.18)
where the last inequality is a consequence of (2.6). On the other hand, the left-hand side of (4.18) can be estimated by Gagliardo-Nirenberg interpolation inequality, which implies

(T 1/ m ) α β-1 +1 L N +2 N (Q T ) ≤ C sup t Ω (T 1/ m (t)) α β-1 +1 dx + T 0 Ω |DT 1/ m | 2 (T 1/ m ) α β-1 -1 dxdt .
Then, when comparing with the right-hand side of (4.18), and recalling that m is already bounded in L 1 , we conclude that

(T 1/ m ) α β-1 +1 L N +2 N (Q T ) ≤ C T 0 Ω m H p (t, x, T 1/ m , Du) • Du -H(t, x, T 1/ m , Du ) dxdt + C ,
where C also depends on m 0 . We use now this information in (4.16) and we deduce that

Ω G (x, m (T ))m (T ) + T 0 Ω F (t, x, m )m dxdt + (T 1/ m ) α β-1 +1 L N +2 N (Q T ) ≤ C + C T 0 Ω (T 1/ m ) α β-1 dxdt .
T 0 Ω m |Du | β (µ + T 1/ m ) α dxdt ≤ C,
which completes (4.9). In particular, we also deduce that

T 0 Ω m |Du | β (µ+m ) α dxdt ≤ C. Finally, integrating (4.1) in (t, T ) × Ω, we get that Ω u (t, x)dx = Ω G(x, ρ m (T ))dx + T t Ω F (t, x, ρ m (s))dxds - T t Ω H(s, x, T 1/ m , Du )dxds.
The last term is bounded above by C from (4.7) and (4.9). Let us deal with Ω G(x, ρ m (T ))dx: from (2.12), G(x, ρ m (T ))

≤ g L (x) + 1 L (ρ m (T )G(x, ρ m (T )) + |c 4 |ρ m (T ))
, where L is an arbitrary positive number and

g L (x) = max 0≤m≤L G(x, m) ≤ 1 λ max 0≤m≤L g(m) + κ, from (2.10). Noting that Ω ρ m (T )G(x, ρ m (T ))dx = Ω m G (x, m )dx,

and that

Ω m (x) = Ω m 0 (x), we deduce from (4.7) that Ω G(x, ρ m (T ))dx ≤ C, for some constant C independent of and µ. The same argument can be used for proving that

T t Ω F (t, x, ρ m (s))dxds ≤ C. Thus Ω u (t, x)dx ≤ C.
Combining this with the lower bound on u already obtained, we obtain (4.10).

Remark 4.3. Note that the constant C appearing in Lemma 4.2 does not depend on µ. Lemma 4.4. For each 0 < < 1, there exist two functions w and z defined on

Q T such that m H p (t, x, T 1/ m , Du ) ≤ c 2 m + w + √ m z ,
where the family

(w ) is bounded in L β (Q T ), the family (z ) is bounded in L 2 (Q T ) and, in addition, relatively compact if β < 2.
Proof. From Remark 2.1, we know that α ≤ β. Therefore, if α ≥ 1, then α ≥ β . This yields

(4.21) (1 -α)β ≥ -α, for any α such that 0 < α ≤ 4(β -1) β , because (4.21) is trivial if 0 < α ≤ 1. From (2.3), we see that m H p (t, x, T 1/ m , Du ) ≤ c 2 m (1 + |Du | β-1 (T 1/ m +µ) α ) ≤ c 2 (m + A 1 + A 2 ), where A 1 = 1 {T 1/ m +µ≤1} |Du | β-1 (T 1/ m +µ) α-1 and A 2 = 1 {T 1/ m +µ≥1} m |Du | β-1 (T 1/ m +µ) α . Let us first deal with A 1 : (4.21) implies that A β 1 ≤ 1 {T 1/ m +µ≤1} |Du | β (T 1/ m +µ) α ≤ |Du | β (T 1/ m +µ) α . Therefore, from (4.9), w = c 2 A 1 is bounded in L β (Q T )
uniformly with respect to (and to µ).

Next, we see that

A 2 ≤ √ m z where z = 1 {T 1/ m +µ≥1} √ m |Du | β-1 (T 1/ m +µ) α .
Then, by Hölder inequality, for any E ⊆ Q T we have

E (z ) 2 dxdt = E 1 {T 1/ m +µ≥1} m |Du | 2β-2 (T 1/ m + µ) 2α dxdt ≤ E m |Du | β (T 1/ m + µ) α dxdt 2 β   E 1 {T 1/ m +µ≥1} m (T 1/ m + µ) α 1+ β β -2 dxdt   1-2 β ≤ E m |Du | β (T 1/ m + µ) α dxdt 2 β E m dxdt 1-2 β . (4.22)
We immediately deduce from (4.9) that z is bounded in L 2 (Q T ). In addition, for any k < 1 we have

E m dxdt ≤ k |E| + 1 k α β-1 E m (T 1/ m ) α β-1 ≤ k |E| + C k α β-1
where we used (4.19). This implies that for small enough,

E m dxdt ≤ C |E| α α+β-1
and so we additionally deduce from (4.22

) that z is equi-integrable in L 2 (Q T ) provided β < 2.
We now collect some compactness properties for the family (u , m ), which are mostly borrowed from [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF]. Proposition 4.5.

(1) The family (m ) is relatively compact in L 1 (Q T ) and in C([0, T ]; W -1,r (Ω)) for some r > 1.

(2) There exist m ∈ C([0, T ]; L 1 (Ω)) and u ∈ L ∞ (0, T ; L 1 (Ω)) ∩ L q (0, T ; W 1,q (Ω)) for all q < N +2 N +1 such that after the extraction of a subsequence (not relabeled), m → m in L 1 (Q T ) and almost everywhere, u → u and Du → Du in L 1 (Q T ) and almost everywhere. Moreover, u and m satisfy

T 0 Ω |Du| β (m + µ) α dxdt + T 0 Ω m |Du| β (m + µ) α dxdt ≤ C, (4.23) Ω G(x, m(T ))m(T ) dx + T 0 Ω mF (t, x, m) dxdt + m α β-1 +1 L N +2 N (Q T ) ≤ C, (4.24) for some C = C(T, H, F, G, m 0 ∞ ) (independent of µ). (3) m (t) → m(t) weakly in L 1 (Ω), for every t ∈ [0, T ], and (3.2) holds for any ϕ ∈ C ∞ c ([0, T ) × Ω). If in addition β < 2, then m → m in C([0, T ]; L 1 (Ω)). Proof. Proceeding as in the end of the proof of Lemma 4.2, we see that (F (t, x, m )) is bounded in L 1 (Q T ), (u (t = T )) is bounded in L 1 (Ω) and that H(t, x, T 1/ m, Du ) is bounded in L 1 (Q T ) uniformly in . Therefore, (-∂ t u -ν∆u ) is bounded in L 1 (Q T ) and (u (t = T )) is bounded in L 1 (Ω).
The compactness of u and Du in L 1 (Q T ) (and almost everywhere) follow from classical results on the heat equation with L 1 data, see e.g. [START_REF] Boccardo | Nonlinear parabolic equations with measure data[END_REF], [START_REF] Boccardo | Nonlinear elliptic and parabolic equations involving measure data[END_REF]: there exist u ∈ L ∞ (0, T ; L 1 (Ω)) ∩ L q ((0, T ); W 1,q (Ω)) for every q < N +2 N +1 , such that after the extraction of a subsequence u → u and Du → Du in L 1 (Q T ) and almost everywhere.

As far as m is concerned, thanks to Lemma 4.4 we can apply the estimates and the compactness results in [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF]Theorem 6.1]. In particular, this implies the point (1) in the statement and so, up to extraction of a subsequence, the convergence of m in L 1 (Q T ) and almost everywhere. The bounds (4.23) and (4.24) follow from (4.7), (4.9) and Fatou's lemma, keeping in mind that the constant C in the estimates of Lemma 4.2 does not depend on µ.

To prove that m is a solution of (3.2), we observe that m H p (t, x, T 1/ m , Du ) strongly converges in L 1 (Q T ). Indeed, the decomposition provided by Lemma 4.4 implies that this term is equi-integrable, so using the almost everywhere convergence and Vitali's theorem, we obtain that m H p (t, x, T 1/ m , Du ) → mH p (t, x, m, Du) in L 1 (Q T ). As a consequence, (3.2) holds for any ϕ ∈ C ∞ c ([0, T ) × Ω). Finally, we have that m (t) is bounded in L log L(Ω) uniformly in time (see e.g. [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF]), which means that it is equi-integrable in L 1 (Ω). By Dunford-Pettis theorem, it is weakly relatively compact in L 1 (Ω), and since we already know that it converges to m(t) in W -1,r (Ω), we conclude that it also converges weakly in L 1 (Ω) to the same limit. In fact, if β < 2, from Lemma 4.4 we have the strong convergence of the terms w , z in L 2 (Q T ). This means that we can apply directly Theorem 6.1 in [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF] to deduce the strong convergence of m in C([0, T ]; L 1 (Ω)).

Henceforth, we let (u, m) as well as the convergent subsequence (u , m ) be given by Proposition 4.5. Let us notice that the weak convergence of m (T ) in L 1 (Ω) easily yields that ρ m (T ) also converges to m(T ) weakly in L 1 (Ω). In particular, it is well-known that the sequence ρ m (T ) generates a family of parametrized Young measures {ν x } ∈ P(R), see e.g. [START_REF] Ball | A version of the fundamental theorem for Young measures, PDE's and continuum models of phase transitions[END_REF], [START_REF] Pedregal | Parametrized measures and variational principles[END_REF]. This means that ν x is a probability measure for a.e. x ∈ Ω, the mapping x → ν x is weakly- * measurable and, for a subsequence (not relabeled),

(4.25) f (x, ρ m (T )) R f (x, λ)dν x (λ) weakly in L 1 (Ω) for every Carathéodory function f (x, s) such that f (x, ρ m (T )) is equi-integrable in L 1 (Ω).
Thanks to the Young measures, we can initially identify the limit of u (T ) and give a first description, by now in a relaxed sense, of the equation satisfied by u.

≤ T 0 Ω F (t, x, m)ϕ dxdt + Ω R G(x, λ)dν x (λ) ϕ(T ) dx (4.26)
for every nonnegative function ϕ ∈ C ∞ c ((0, T ] × Ω). In addition, we have that

R G(x, λ)λ dν x (λ) ∈ L 1 (Ω), G(x, m(T ))m(T ) ∈ L 1 (Ω) and (4.27) Ω G(x, m(T ))m(T ) dx ≤ C for some C = C(T, H, F, G, m 0 ∞ ) (independent of µ).
Proof. From (4.7) and the convergence of m , we see that (T 1/ m ) α β-1 converges in L 1 (Q T ). Using this observation and (2.2), we can apply Fatou's lemma and obtain that, for any ϕ ∈

C ∞ c ((0, T ] × Ω), ϕ ≥ 0, T 0 Ω H(t, x, m, Du)ϕ(t, x)dxdt ≤ lim inf →0 T 0 Ω H(t, x, T 1/ m , Du )ϕ(t, x)dxdt.
Moreover, from (4.7) and (2.8), and the definition of F , we see that F (t, x, m ) is equi-integrable in Q T . Therefore, from Vitali's theorem, (4.28) lim

→0 T 0 Ω F (t, x, m )ϕ(t, x)dxdt = T 0 Ω F (t, x, m)ϕ(t, x)dxdt.
Let us deal with the boundary condition at t = T : again from (4.7), the definition of G and (2.10), we deduce that G(x, ρ m (T )) is equi-integrable, so by the properties of Young measures we obtain that lim

→0 Ω u (T, x)ϕ(T, x, )dx = lim →0 Ω G(x, ρ m (T )) ρ ϕ(T, x)dx = Ω R G(x, λ)dν x (λ)ϕ(T, x)dx.
T k (G(x, ρ m (T )))ρ m (T ) R T k (G(x, λ))λ dν x (λ) where T k (s) = min(s, k). Hence Ω R T k (G(x, λ))λ dν x (λ) dx = lim →0 Ω T k (G(x, ρ m (T )))ρ m (T ) dx ≤ C .
By monotone convergence, letting k → ∞ we deduce that R G(x, λ)λ dν x (λ) ∈ L 1 (Ω). Proof. From (4.10), we know that u (0) is bounded in L 1 (Ω). Therefore, there exists some bounded measure χ such that, up to a subsquence, u (0) tends to χ weakly * . Moreover, from Lemma 3.3, Remark 3.4 and Lemma 4.6, u(0) is a well defined Radon measure. Let us prove that χ ≤ u(0). Consider ϕ(t, x) = ψ(x)(1 -t/η) + , where ψ is any nonnegative smooth function defined on Ω. Taking ϕ as a test-function in (4.1), we obtain that

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H (t, x, T 1/ m , Du )ϕ dxdt = T 0 Ω F (t, x, m )ϕ dxdt - Ω ϕ(0, x)u (0, x)dx,
which implies by using Fatou's lemma as in the proof of Proposition 4.5 that

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H(t, x, m, Du)ϕ dxdt ≤ T 0 Ω F (t, x, m)ϕ dxdt -χ, ψ . (4.30)
On the other hand, as η → 0,

T 0 Ω u ϕ t dxdt = -1 η η
0 Ω u ψ(x) dxdt tends to -u(0), ψ . Passing to the limit in (4.30) as η → 0, we obtain that u(0), ψ ≥ χ, ψ , which is the desired result.

Crossed energy inequality

We prove here the main step towards the uniqueness result. Let u ∈ L ∞ (0, T ; L 1 (Ω)) satisfy (4.26) for some family of probability measures {ν x } (weakly- * measurable w.r.t. x) such that

(5.1) R G(x, λ)λ dν x (λ) ∈ L 1 (Ω)
and let m ∈ C 0 ([0, T ]; L 1 (Ω)) be a solution of (1.4), and assume that (u, m) satisfy (ii) in Definition 3.1. Then we have

m0 , u(0) ≤ Ω R G(x, λ) dν x (λ) m(T ) dx + T 0 Ω F (t, x, m) m dxdt + T 0 Ω [ m H p (t, x, m, Dũ) • Du -m H(t, x, m, Du)] dxdt (5.2)
for any couple (ũ, m) satisfying the same conditions as (u, m).

Remark 5.2. The above Lemma applies in particular if ν x = δ m(T,x) , in which case u is simply a subsolution of (1.3) and condition (5.1) is equivalent to G(x, m(T ))m(T ) ∈ L 1 (Ω).

Proof. Let ρ δ (•) be a sequence of standard symmetric mollifiers in R N and set

mδ (t, x) = m(t) ρ δ := R N m(t, y)ρ δ (x -y) dy .
Notice in particular that mδ , D mδ ∈ L ∞ (Q T ) since m ∈ L ∞ (0, T ; L 1 (Ω)). We also take a sequence of 1-d mollifiers ξ ε (t) such that supp(ξ ε ) ⊂ (-ε, 0), and we set mδ,ε := T 0 ξ ε (s -t) mδ (s)ds .

Notice that this function vanishes near t = 0, so we can take it as test function in the inequality satisfied by u. We get

T 0 Ω u [∂ t mδ,ε -ν∆ mδ,ε ] dxdt + T 0 Ω H(t, x, m, Du) mδ,ε dxdt ≤ T 0 Ω F (t, x, m) mδ,ε dxdt + Ω R G(x, λ) dν x (λ) mδ,ε (T ) dx (5.3)
The first integral is equal to

T 0 Ω [-∂ s u δ,ε -ν∆u δ,ε ] m(s, y) dsdy
where u δ,ε (s, y) = T 0 Ω u(t, x)ξ ε (s -t)ρ δ (x -y) dtdx. Notice that this function vanishes near s = T . Thus using the equation of m we have

T 0 Ω [-∂ s u δ,ε -ν∆u δ,ε ] m(s, y) dsdy = - T 0 Ω m(s, y) b(s, y) • D y u δ,ε dsdy + Ω m0 (y)u δ,ε (0)dy ,
where b = H p (t, x, m, Dũ). We shift the convolution kernels from u to m in the right-hand side and we use this equality in (5.3). We get Since now, we distinguish between the cases α ≤ 1 and α > 1. Recall that α ≤ 4(β-1) β is always assumed, hence the latter case only happens for 1 < 4(β-1) β , i.e. β > 4 3 .

- T 0 Ω Du • wδ,ε dtdx + T 0 Ω H(t, x, m, Du) mδ,ε dxdt + Ω ( m0 ρ δ ) T 0 u(t)ξ ε (-t) dt dx ≤ T 0 Ω F (t, x, m) mδ,ε dxdt + Ω R G(x, λ) dν x (λ) mδ,ε (T ) dx
(A) Case when α ≤ 1.

Let us deal with the first two integrals in (5.4). By assumption (2.3) we have

|( m b) ρ δ | ≤ c 2 mδ + c 2 m |Dũ| β ( m + µ) α ρ δ 1 β m ( m + µ) α ρ δ 1 β ≤ c 2 mδ + c 2 m |Dũ| β ( m + µ) α ρ δ 1 β ( m ρ δ ) 1-α β
where we used α ≤ 1 in the latter inequality. Therefore, using also (2.2) we estimate

H(t, x, m, Du) mδ (s) -Du • wδ (s) ≥ c 0 mδ (s) |Du| β (m + µ) α -c 1 mδ (s)(1 + m α β-1 ) -c 2 mδ (s)|Du| -c 2 |Du| m(s) |Dũ(s)| β ( m(s) + µ) α ρ δ 1 β mδ (s) 1-α β
which yields, by Young's inequality,

H(t, x, m, Du) mδ (s) -Du • wδ (s) ≥ c 0 2 mδ (s) |Du| β (m + µ) α -c mδ (s)(1 + m α β-1 ) -σ |Du| β mδ (s) 1-α -C σ m(s) |Dũ(s)| β ( m(s) + µ) α ρ δ .
(5.5)

Since α ≤ 1, choosing σ < c 0 2 we have σ m1-α δ < c 0 2 mδ (m + µ) α + c 0 2 (m + µ) 1-α , while mδ (1 + m α β-1 ) ≤ c (1 + m1+ α β-1 δ + m 1+ α β-1
) .

So we conclude from (5.5) that H(t, x, m, Du) mδ (s) -Du • wδ (s) ≥ -c ( mδ (s)

1+ α β-1 + m 1+ α β-1 ) -c (m + µ) 1-α |Du| β -c m(s) |Dũ(s)| β ( m(s) + µ) α ρ δ .
Notice that condition (ii) in Definition 3.1 implies that (m + µ) 1-α |Du| β belongs to L 1 (Q T ). Due to the summability of all the above terms, we are allowed to use Fatou's lemma and deduce lim inf

ε→0 T 0 Ω [H(t, x, m, Du) mδ,ε -Du • wδ,ε ] dxdt = lim inf ε→0 T 0 Ω T 0 [H(t, x, m, Du) mδ (s) -Du • wδ (s)]ξ ε (s -t) dsdxdt ≥ T 0 Ω [H(t, x, m, Du) mδ -Du • wδ ] dxdt ,
and in the same way

lim inf δ→0 T 0 Ω [H(t, x, m, Du) mδ -Du • wδ ] dxdt ≥ T 0 Ω [H(t, x, m, Du) m -Du • H p (t, x, m, Dũ) m] dxdt .
(5.6)

Now we consider the remaining terms in (5.4), in particular the terms at t = 0 and t = T . For t = T , we have

Ω R G(x, λ) dν x (λ) mδ,ε (T ) dx = Ω R G(x, λ) dν x (λ) mδ (T ) dx + Ω T 0 R G(x, λ) dν x (λ)[ mδ (s) -mδ (T )]ξ ε (s -T ) dsdx .
Since m ∈ C 0 ([0, T ]; L 1 (Ω)) implies that T 0 [ mδ (s) -mδ (T )]ξ ε (s -T ) ds converges to zero uniformly as ε → 0, the last integral will vanish as ε → 0.

For t = 0, we recall that Ω u(t) ( m0 ρ δ )dx has a trace at t = 0 from Lemma 3.3, and this trace is also continuous as δ → 0 since m0 is continuous. Therefore, as ε → 0 we obtain from (5.4)

T 0 Ω [H(t, x, m, Du) mδ -Du • wδ ] dxdt + Ω ( m0 ρ δ )u(0) dx ≤ T 0 Ω F (t, x, m) mδ dxdt + Ω R G(x, λ) dν x (λ) mδ (T ) dx (5.7)
We apply assumptions (2.8) and (2.10) to deal with the last two terms. Indeed, if c 4 is the constant in (2.11), we have

c 4 mδ ≤ F (t, x, m) mδ ≤ (F (m) + |c 4 |) mδ implies that c 4 mδ ≤ F (t, x, m) mδ ≤ (F (t, x, m) + |c 4 |)m1 m≥ mδ + 1 λ f (m) + κ + |c 4 | mδ 1 m< mδ ≤ (F (t, x, m) + |c 4 |)m + 1 λ f ( mδ ) + κ + |c 4 | mδ ≤ F (t, x, m)m + 1 λ f ( mδ ) mδ + C(m + mδ ).
Hence

c 4 mδ ≤ F (t, x, m) mδ ≤ F (t, x, m)m + 1 λ f ( mδ ) mδ + C (m + mδ ) ≤ F (t, x, m)m + 1 λ (f (m)m) ρ δ + C (m + mδ ),
and since F (t, x, m)m ∈ L 1 (Q T ) (and the same holds for f (m)m by (2.8)), we conclude that the sequence F (t, x, m) mδ is bounded above and below by convergent sequences in L 1 (Q T ), which allows us to pass to the limit as δ → 0. We proceed similarly for the term with G. Indeed, by (2.9), R G(x, λ) dν x (λ) mδ (T ) ≤ R G(x, λ)λ dν x (λ) + G(x, mδ (T )) mδ (T ) for a.e. x ∈ Ω, and using (2.10) we deduce

R G(x, λ) dν x (λ) mδ (T ) ≤ R G(x, λ)λ dν x (λ) + 1 λ g( mδ (T )) mδ (T ) + C mδ (T ) ≤ R G(x, λ)λ dν x (λ) + 1 λ (g( m(T )) m(T )) ρ δ + C mδ (T ) .
Thanks to (5.1) and since G(x, m(T )) m(T ) ∈ L 1 (Ω) (and and the same holds for g( m(T )) m(T ) by (2.10)), we can handle this term too. Finally, passing to the limit in (5.7) and using also (5.6), we obtain (5.2).

(B) Case when α > 1.

First recall from Remark 2.1 that α ≤ β. Moreover, since 1 ≥ β -1, we also have α > β -1, i.e. β -α < 1.

We estimate now differently the first two terms in (5.4). First of all, by Young's inequality, we have (omitting to write that mδ , Dũ are evaluated at s and m, Du at t)

Du • wδ ≤ σ max(m + µ, mδ ) m β-α β-1 δ | wδ | β + C σ mβ-α δ max(m + µ, mδ ) β-1 |Du| β
and using β ≥ α and α > 1 this yields (5.8)

Du • wδ ≤ σ m + µ + mδ m β-α β-1 δ | wδ | β + C σ (m + µ) 1-α |Du| β . By (2.3), |( m b) ρ δ | ≤ c 2 mδ + c 2 m |Dũ| β-1 ( m + µ) α ρ δ ≤ c 2 mδ + c 2 |Dũ| β ( m + µ) α ρ δ 1 β mβ ( m + µ) α ρ δ 1 β ≤ c 2 mδ + c 2 |Dũ| β ( m + µ) α ρ δ 1 β m β-α β δ
where we used β -α ≤ 1 in the last inequality. We deduce that (5.9)

| wδ | β m β-α β-1 δ ≤ C m α β-1 δ + C |Dũ| β ( m + µ) α ρ δ .
Combining this information with (5.8), we deduce that

Du • wδ ≤ σ C[m + µ + mδ ] |Dũ| β ( m + µ) α ρ δ + C [1 + m 1+ α β-1 + m1+ α β-1 δ ] + C σ (m + µ) 1-α |Du| β .
(5.10)

We first use this inequality to get

Du • wδ,ε ≤ σ C T 0 [m + µ + mδ ] |Dũ| β ( m + µ) α ρ δ ξ ε (s -t)ds + C[1 + m 1+ α β-1 + T 0 m1+ α β-1 δ ξ ε (s -t)ds] + C σ (m + µ) 1-α |Du| β .
Since Ω m(t, x)dx = Ω m 0 (x)dx and |D ũ| β ( m+µ

) α ∈ L 1 (Q T ), T 0 m(t) |Dũ(s)| β ( m(s) + µ) α ρ δ ξ ε (s -t)ds ε→0 → m(t) |Dũ(t)| β ( m(t) + µ) α ρ δ in L 1 (Q T ).
Using also that mδ (s) is continuous, we deduce that Du • wδ,ε is dominated by a L 1 -convergent sequence, so we are allowed to take ε → 0 in (5.4), obtaining

T 0 Ω [H(t, x, m, Du) mδ -Du • wδ ] dxdt + Ω ( m0 ρ δ )u(0) dx ≤ T 0 Ω F (t, x, m) mδ dxdt + Ω R G(x, λ)dν x (λ) mδ (T ) dx .
(5.11)

Thanks to (5.10), and since the right-hand side is bounded as in the previous case, the above inequality also implies

T 0 Ω H(t, x, m, Du) mδ dxdt ≤ C + σ C T 0 Ω [m + µ + mδ ] |Dũ| β ( m + µ) α ρ δ dxdt + C T 0 Ω [1 + m 1+ α β-1 + m1+ α β-1 δ ] dxdt + C σ T 0 Ω (m + µ) 1-α |Du| β dxdt .
On account of the bounds on m, m in L 1+ α β-1 (Q T ) and using the properties (ii) in Definition 3.1, the last line of the above inequality is uniformly bounded. Therefore, using (2.2) we deduce that

T 0 Ω mδ |Du| β (m + µ) α dxdt ≤ C + σ C T 0 Ω m δ |Dũ| β ( m + µ) α dxdt + σ C T 0 Ω mδ |Dũ| β ( m + µ) α ρ δ dxdt ,
where m δ = m ρ δ . We obtain a similar inequality reversing the roles of (u, m) and (ũ, m), and by addition we get

T 0 Ω m δ |Dũ| β ( m + µ) α dxdt + T 0 Ω mδ |Du| β (m + µ) α dxdt ≤ C + σ C T 0 Ω m δ |Dũ| β ( m + µ) α dxdt + T 0 Ω mδ |Du| β (m + µ) α dxdt + σ C T 0 Ω mδ |Dũ| β ( m + µ) α ρ δ dxdt + T 0 Ω m δ |Du| β (m + µ) α ρ δ dxdt ,
hence choosing σ sufficiently small we conclude that

T 0 Ω m δ |Dũ| β ( m + µ) α dxdt + T 0 Ω mδ |Du| β (m + µ) α dxdt ≤ C + C T 0 Ω mδ |Dũ| β ( m + µ) α ρ δ dxdt + T 0 Ω m δ |Du| β (m + µ) α ρ δ dxdt .
(5.12)

We wish now to estimate last terms in the above inequality. To this purpose, we use (5.11) with m replaced by m δ (and so w by (mH p (t, x, Du)) ρ δ ). Since the terms with F and G can be estimated as before, we get

T 0 Ω m δ [H(t, x, m, Du) ρ δ ] dxdt ≤ T 0 Ω [Du ρ δ ] • [(mH p (t, x, Du)) ρ δ ] dxdt + C .
Therefore, using (2.2) we obtain

c 0 T 0 Ω m δ [ |Du| β (m + µ) α ρ δ ] dxdt ≤ T 0 Ω [Du ρ δ ] • [(mH p (t, x, Du)) ρ δ ] dxdt + C ≤ c 2 T 0 Ω |Du ρ δ | m |Du| β-1 (m + µ) α ρ δ dxdt + C T 0 Ω |Du ρ δ | m δ dxdt ≤ c 2 T 0 Ω |Du ρ δ | |Du| β (m + µ) α ρ δ 1 β m β-α β δ dxdt + C T 0 Ω |Du ρ δ | m δ dxdt ≤ c 0 2 T 0 Ω (m + µ) ρ δ |Du| β (m + µ) α ρ δ dxdt + C T 0 Ω |Du ρ δ | β ((m + µ) ρ δ ) 1-α dxdt + C T 0 Ω 1 + m 1+ α β-1 δ dxdt .
Hence we conclude that

T 0 Ω m δ |Du| β (m + µ) α ρ δ dxdt ≤ C T 0 Ω |Du ρ δ | β ((m + µ) ρ δ ) 1-α dxdt + C .
Now we observe that the function (m, ξ) → (m + µ) 1-α |ξ| β is convex as a function of two variables, therefore

|Du ρ δ | β ((m + µ) ρ δ ) 1-α ≤ (m + µ) 1-α |Du| β ρ δ which is bounded since (m+µ) 1-α |Du| β ∈ L 1 (Q T ).
We deduce the uniform bound (with respect to δ)

T 0 Ω m δ |Du| β (m + µ) α ρ δ dxdt ≤ C ,
and similarly for m. Going back to (5.12), this also implies the bound

T 0 Ω m δ |Dũ| β ( m + µ) α dxdt + T 0 Ω mδ |Du| β (m + µ) α dxdt ≤ C . Now, reasoning as in (5.8)-(5.10), E |Du • wδ | ≤   E (m + µ + mδ ) m β-α β-1 δ | wδ | β   1 β E (m + µ) 1-α |Du| β 1 β ≤ C + E [m + µ + mδ ] |Dũ| β ( m + µ) α ρ δ 1 β E (m + µ) 1-α |Du| β 1 β
.

The bounds previously established yield

E |Du • wδ | ≤ C E (m + µ) 1-α |Du| β 1 β
and since the set E is arbitrary and (m+µ) 1-α |Du| β ∈ L 1 (Q T ), we deduce the equi-integrability of Du • wδ . Finally, this allows us to pass to the limit in (5.11) and to obtain (5.2).

A similar Lemma holds for the case of singular congestion, suitably adapted to the formulation of this case.

Lemma 5.3. Assume that H satisfies assumptions (2.2)-(2.3) with µ = 0, and that (2.8)-(2.10) hold true. Let u ∈ L ∞ (0, T ; L 1 (Ω)) satisfy

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H(t, x, m, Du)1 {m>0} ϕ dxdt ≤ T 0 Ω F (t, x, m)ϕ dxdt + Ω R G(x, λ)dν x (λ) ϕ(T ) dx
for some family of probability measures {ν x } (weakly- * measurable w.r.t. x) such that (5.1) is satisfied, and let m ∈ C 0 ([0, T ]; L 1 (Ω)) be a solution of (3.6). Assume that (u, m) satisfy the conditions (ii) in Definition 3.7. Then

m0 , u(0) ≤ Ω R G(x, λ)dν x (λ) m(T ) dx + T 0 Ω F (t, x, m) m dxdt + T 0 Ω [ m H p (t, x, m, Dũ) • Du -m H(t, x, m, Du)] 1 {m>0 , m>0} dxdt (5.13)
for any couple (ũ, m) satisfying the same conditions as (u, m).

Proof. The proof follows the same lines as that of Lemma 5.1.

(A) Case when α ≤ 1. Proceeding as in the proof of Lemma 5.1, we obtain the counterparts of (5.3)

T 0 Ω u [∂ t mδ,ε -ν∆ mδ,ε ] dxdt + T 0 Ω H(t, x, m, Du) mδ,ε 1 {m>0} dxdt ≤ T 0 Ω F (t, x, m) mδ,ε dxdt + Ω R G(x, λ)dν x (λ) mδ,ε (T ) dx (5.14)
and of (5.4):

- T 0 Ω Du • wδ,ε dtdx + T 0 Ω H(t, x, m, Du) mδ,ε 1 {m>0} dxdt + Ω (m 0 ρ δ ) T 0 u(t)ξ ε (-t) dt dx ≤ T 0 Ω F (t, x, m) mδ,ε dxdt + Ω R G(x, λ)dν x (λ) mδ,ε (T ) dx (5.15)
where we denote wδ = [( m H p (t, x, m, Dũ)

1 { m>0} ) ρ δ ] and wδ,ε = T 0 wδ (s) ξ ε (s -t) ds. Note that wδ = [( m H p (t, x, m, Dũ)) ρ δ ] if α < 1. Set b = H p (t, x, m, Dũ) 1 { m>0} . We deduce from Assumption (2.3) that |( m b) ρ δ | ≤ c 2 mδ + c 2 m |Dũ| β mα 1 { m>0} ρ δ 1 β m mα ρ δ 1 β ≤ c 2 mδ + c 2 m1-α |Dũ| β 1 { m>0} ρ δ 1 β ( m ρ δ ) 1-α β
where we used α ≤ 1 in the latter inequality. Therefore, using also (2.2) we estimate

H(t, x, m, Du)1 {m>0} mδ (s) -Du • wδ (s) ≥ c 0 mδ (s) |Du| β m α 1 {m>0} -c 1 mδ (s)(1 + m α β-1 ) -c 2 mδ (s)|Du| -c 2 |Du| m(s) 1-α |Dũ(s)| β 1 { m(s)>0} ρ δ 1 β mδ (s) 1-α β .
(5.16)

We now use the fact that Du = 0 a.e. in {m = 0}. From Young inequalities, we deduce (5.17)

c 2 mδ (s)|Du| = c 2 mδ (s)|Du|1 {m>0} ≤ c 0 2 mδ (s) |Du| β m α 1 {m>0} + C mδ (s)m α β-1 ,
and

c 2 |Du| m(s) 1-α |Dũ(s)| β 1 { m(s)>0} ρ δ 1 β mδ (s) 1-α β ≤ σ|Du| β 1 {m>0} mδ (s) 1-α + C σ m(s) 1-α |Dũ(s)| β 1 { m(s)>0} ρ δ , (5.18
)

and if 0 < σ ≤ c 0 /2, (5.19
)

σ1 {m>0} mδ (s) 1-α ≤ c 0 2 mδ (s) m α 1 {m>0} + c 0 2 m 1-α 1 {m>0} , because 0 < α ≤ 1. Finally, (5.20) mδ (s)m α β-1 ≤ C 1 + m 1+ α β-1 + mδ (s) 1+ α β-1
.

Combining (5.16) with (5.17)-(5.20), we get

H(t, x, m, Du) mδ (s) -Du • wδ (s) ≥ -C 1 + mδ (s) 1+ α β-1 + m 1+ α β-1 - c 0 2 m 1-α |Du| β 1 {m>0} -c m(s) 1-α |Dũ(s)| β 1 { m(s)>0} ρ δ .
The conclusion of the proof is exactly as for Lemma 5.1.

(B) Case when α > 1.

We start from (5.15), and proceed differently for the first two terms. From (2.3) we have | b| ≤ c 2 (1 + |Dũ| β-1 m-α 1 { m>0} ), so we estimate wδ as in (5.9) of Lemma 5.1 and we get

| wδ | m β-α β-1 δ ≤ C m α β-1 δ + C |Dũ| β mα 1 { m>0} ρ δ .
Since Du = 0 a.e. if m = 0, we can use (5.8) with µ = 0 and, similarly as in Lemma 5.1, we get at

Du • wδ = Du • wδ m β-α β δ m β-α β δ 1 {m>0} ≤ σC( mδ + m1 {m>0} ) |Dũ| β mα 1 { m>0} ρ δ + C 1 + m 1+ α β-1 + m1+ α β-1 δ + C σ m 1-α 1 {m>0} |Du| β .
Thanks to this estimate, we can proceed exactly as in the proof of Lemma 5.1 (when α > 1), and obtain

T 0 Ω [H(t, x, m, Du)1 {m>0} mδ -Du • wδ ] dxdt + Ω ( m0 ρ δ )u(0) dx ≤ T 0 Ω F (t, x, m) mδ dxdt + Ω R G(x, λ)dν x (λ) mδ (T ) dx . (5.21)
and, after some steps,

T 0 Ω m δ |Dũ| β mα 1 { m>0} dxdt + T 0 Ω mδ |Du| β m α 1 {m>0} dxdt ≤ C + C T 0 Ω mδ |Dũ| β mα 1 { m>0} ρ δ dxdt + T 0 Ω m δ |Du| β m α 1 {m>0} ρ δ dxdt .
(

To estimate the last two terms in (5.22), we use (5.21) where we replace m with m δ and proceed as above,

c 0 T 0 Ω m δ |Du| β 1 {m>0} m α ρ δ dxdt ≤ T 0 Ω Du ρ δ • (mH p (t, x, Du)1 {m>0} ) ρ δ dxdt + C ≤ c 2 T 0 Ω |Du ρ δ | m 1-α |Du| β-1 1 {m>0} ρ δ dxdt + C ≤ c 2 T 0 Ω |Du ρ δ | |Du| β m α 1 {m>0} ρ δ 1 β (m ρ δ ) β-α β dxdt + C
Let us observe that m ρ δ = 0 ⇒ Du ρ δ = 0 since Du = 0 at almost every point where m = 0. This means that last integral in the above inequality is restricted to the set {m ρ δ > 0} (even in the limiting case β = α). By Young's inequality, we deduce that (5.23)

T 0 Ω m δ |Du| β m α 1 {m>0} ρ δ dxdt ≤ C T 0 Ω |Du ρ δ | β (m ρ δ ) 1-α 1 {m ρ δ >0} dxdt +C.
Let us define the convex and lower semi-continuous function Ψ on

R N × R by Ψ(p, m) =    m 1-α |p| β if m > 0, 0 if m = 0 and p = 0, +∞ otherwise.
This implies that

(m ρ δ ) 1-α |Du ρ δ | β 1 {m ρ δ >0} = Ψ(Du ρ δ , m ρ δ ) ≤ Ψ(Du, m) ρ δ = m 1-α |Du| β 1 {m>0} ρ δ
and the latter function is bounded in L 1 (Q T ) by assumption. We deduce from (5.23) the uniform bound (with respect to δ)

T 0 Ω m δ |Du| β m α 1 {m>0} ρ δ dxdt ≤ C,
and the same holds replacing m with m and u with ũ. Going back to (5.22), we obtain

T 0 Ω m δ |Dũ| β mα 1 { m>0} dxdt + T 0 Ω mδ |Du| β m α 1 {m>0} dxdt ≤ C
, and we conclude as in the proof of Lemma 5.1.

Existence and uniqueness for non singular congestion

Our first goal is to prove the existence of a weak solution and, to this purpose, we now show the strong convergence of G (x, m (T )) in L 1 (Ω). Lemma 6.1. Consider a subsequence (u , m ) converging to (u, m) as in Proposition 4.5. We have that G (x, m (T )) converges to G(x, m(T )) in L 1 (Ω) and (u, m) satisfies the energy identity:

T 0 Ω m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt + Ω m(T )G(x, m(T ))dx + T 0 Ω mF (t, x, m) dxdt = u(0), m 0 . (6.1)
Proof. We start from the energy identity for the solution of (4.1)-(4.3):

T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt + Ω m (T )G (x, m (T ))dx + T 0 Ω m F (t, x, m ) dxdt = Ω u (0)m 0 dx. (6.2)
Thanks to (2.6) and (2.11), Fatou's lemma implies that lim inf

→0 T 0 Ω m H p (t, x, T 1/ m , Du ) • Du -H(t, x, T 1/ m , Du ) dxdt ≥ T 0 Ω m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt, and lim inf →0 T 0 Ω m F (t, x, m ) dxdt ≥ T 0 Ω mF (t, x, m) dxdt.
From Lemma 4.7 we can assume that u | t=0 converges weakly * to a bounded measure χ on Ω, and χ ≤ u(0). Since m 0 ∈ C(Ω), we deduce that lim

→0 Ω u (0)m 0 dx = χ, m 0 ≤ u(0), m 0 .
Combining the informations above, we get from (6.2)

lim sup →0 Ω m (T )G (x, m (T ))dx ≤ u(0), m 0 - T 0 Ω mF (t, x, m) dxdt - T 0 Ω m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt
We now use (5.2) in Lemma 5.1 with m = m and ũ = m, and we get lim sup

→0 Ω m (T )G (x, m (T ))dx ≤ u(0), m 0 - T 0 Ω mF (t, x, m) dxdt - T 0 Ω m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt ≤ Ω R m(T )G(x, λ)dν x (λ)dx . (6.3)
We now use the monotonicity of G in order to get the strong convergence of G (x, m (T )). Indeed, if we set m := ρ m (T ) and T k (r) = min(r, k),

Ω [T k (G(x, m )) -T k (G(x, m(T )))] [ m -m(T )]dx ≤ Ω G(x, m ) m dx - Ω T k (G(x, m ))m(T ) dx - Ω T k (G(x, m(T ))) [ m -m(T )]dx .
Since m weakly converges to m(T ), the last term vanishes as → 0. The first one is estimated by (6.3), while for the second we can use (4.25) with f (x,

•) = T k (G(x, •)). Therefore, we get lim sup →0 Ω [T k (G(x, m )) -T k (G(x, m(T )))] [ m -m(T )]dx ≤ Ω R G(x, λ)m(T )dν x (λ) dx - Ω R T k (G(x, λ))m(T )dν x (λ) dx .
Letting k → ∞ we conclude that (6.4) lim sup

k→∞ lim sup →0 Ω [T k (G(x, m )) -T k (G(x, m(T )))] [ m -m(T )]dx = 0 .
We claim now that, as a consequence of (6.4), G(x, m ) converges to G(x, m(T )) almost everywhere in Ω, at least for a subsequence. Indeed, we observe that

Ω [G(x, mε )-G(x,m(T ))] [ mε -m(T )] 1+ mε +m(T ) dx ≤ Ω [T k (G(x, m )) -T k (G(x, m(T )))] [ m -m(T )]dx + Ω |G(x, m )| 1 {G(x, m )>k} dx + Ω |G(x, m(T ))| 1 {G(x,m(T ))>k} dx
and last two terms tend to zero as k → ∞ uniformly with respect to ε. So we have lim sup

→0 Ω [G(x, mε )-G(x,m(T ))] [ mε -m(T )] 1+ mε +m(T ) dx ≤ lim sup →0 Ω [T k (G(x, m )) -T k (G(x, m(T )))] [ m -m(T )]dx + sup ε Ω |G(x, m )| 1 {G(x, m )>k} dx + Ω |G(x, m(T ))| 1 {G(x,m(T ))>k} dx
hence, using (6.4) and letting k → ∞, the right-hand side vanishes. We deduce that lim sup

→0 Ω [G(x, mε ) -G(x, m(T ))] [ mε -m(T )] 1 + mε + m(T ) dx = 0
which means, since G is monotone, the L 1 convergence of the integrand function. We deduce that, up to subsequences,

[G(x, mε ) -G(x, m(T ))] [ mε -m(T )] 1 + mε + m(T ) → 0 a.e.
and this readily yields the a.e. convergence of G(x, m ) to G(x, m(T )) in Ω.

Finally, since G(x, m ) is also equi-integrable by (4.7), it is therefore convergent in L 1 (Ω). As a consequence, the L 1 convergence of G ε (x, m (T )) towards G(x, m(T )) is established. In addition, we now deduce that R G(x, λ)dν x (λ) = G(x, m(T )) .

We insert this information in the right-hand side of (6.3), moreover now we can use Fatou's lemma in the left-hand side and we conclude with the identity

Ω G(x, m(T ))m(T )dx = u(0), m 0 - T 0 Ω mF (t, x, m) dxdt - T 0 Ω m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt
which is (6.1). Remark 6.2. The monotonicity of G is only used in the above lemma in order to obtain the strong L 1 convergence of G ε (x, m (T )). If q < 2, this condition would not be required since one already knows from Proposition 4.5 that m (T ) strongly converges in L 1 to m(T ).

We can finally conclude with the existence result. Theorem 6.3. Consider a subsequence (u , m ) converging to (u, m) as in Proposition 4.5. Then (u, m) is a weak solution of (1.3).

Proof. With the results proved in § 4, there only remains to pass to the limit in the Hamilton-Jacobi equation and show that (1.3) holds. By Lemma 6.1, G (x, m (T )) → G(x, m(T )) in L 1 (Ω). We also know that F (t, x, m ) → F (t, x, m) in L 1 (Q T ) (see the proof of Lemma 4.6). Now we observe that assumptions (2.1)-(2.3) imply that H satisfies

c 0 |p| β (m + µ) α -c 1 (1+m α β-1 ) ≤ H(t, x, m, p) ≤ c |p| + |p| β (m + µ) α ≤ c |p| β (m + µ) α + 1 + m α β-1 .
Therefore, up to addition of a L 1 -convergent sequence, the Hamiltonian (t, x, p) → H(t, x, T 1/ m (t, x), p) is nonnegative and has natural growth. Since u is bounded below, without loss of generality we can assume that the solution u is also nonnegative. It is therefore possible to apply (a straightforward adaptation of) the result in [START_REF] Porretta | Existence results for nonlinear parabolic equations via strong convergence of truncations[END_REF]Theorem 3.1], in order to deduce that H(t, x, T 1/ m , Du ) → H(t, x, m, Du) in L 1 (Q T ), and we can pass to the limit in (4.1)-(4.3) and obtain (3.1).

In order to prove the uniqueness of weak solutions, the main step was given by Lemma 5.1. We need however a counterpart which ensures that any weak solution satisfies the energy equality (6.1). To this purpose, we follow an argument developed in [START_REF] Porretta | On the weak theory for mean field games systems[END_REF]. Lemma 6.4. Let (u, m) be any weak solution of (1.3)- (1.5). Then the energy identity (6.1) holds true.

Proof. By Lemma 5.1, we already know that

Ω m 0 u(0) dx ≤ Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m)m dxdt + T 0 Ω [m H p (t, x, m, Du) • Du -m H(t, x, m, Du)] dxdt
where we used the fact that u ∈ C 0 ([0, T ]; L 1 (Ω)). In order to prove the reverse inequality, let u k := min(u, k). Since -u t -ν∆u ∈ L 1 (Q T ), by Kato's inequality,

-(u k ) t -ν∆u k + H(t, x, m, Du)1 {u<k} ≥ F (t, x, m)1 {u<k} .
Since m t -ν∆m -div(m b) = 0 for some b such that m|b| 2 ∈ L 1 (Q T ), by [START_REF] Porretta | Weak solutions to Fokker-Planck equations and mean field games[END_REF]Theorem 3.6] we know that m is also a renormalized solution, hence it satisfies

S n (m) t -ν∆S n (m) -div(S n (m)m H p (t, x, m, Du)) = ω n
where S n (m) is a suitable C 1 truncation (i.e. S n (r) = nS(r/n), for some S compactly supported in [-2, 2] such that S = 1 in [-1, 1]) and where ω n → 0 in L 1 (Q T ). Notice that both u k and S n (m) are bounded functions and belong to L 2 (0, T ; H 1 ) ∩ C 0 ([0, T ]; L 1 (Ω)). By using u k in the equation of S n (m) we get

Ω S n (m 0 )u k (0) dx - Ω u k (T ) S n (m(T )) dx ≥ T 0 Ω F (t, x, m)1 {u<k} S n (m) dxdt + T 0 Ω S n (m)m H p (t, x, m, Du) • Du -S n (m)H(t, x, m, Du) 1 {u<k} dxdt - T 0 Ω ω n u k dxdt .
Letting first n → ∞, the last term vanishes since u k is bounded. The regularity of weak solutions allows us to pass to the limit in the other terms (recall that u(T ) = G(x, m(T ))) and we obtain

Ω m 0 u k (0) dx - Ω u k (T ) m(T ) dx ≥ T 0 Ω F (t, x, m)1 {u<k} m dxdt + T 0 Ω [m H p (t, x, m, Du) • Du -mH(t, x, m, Du)] 1 {u<k} dxdt .
Finally, letting k → ∞ we deduce the desired inequality and we conclude that (6.1) holds true. Lemma 6.5. Under all the assumptions made in Theorem 3.6, there is a unique weak solution of (1.3)-(1.5).

Proof. Let (u, m) and (ũ, m) be two weak solutions of (1.3)- (1.5). By Lemma 6.4, they both satisfy the energy identity (6.1), so using also Lemma 5.1 overall we know that We then deduce that, in a weak sense, (u -ũ) t = 0 and u(T ) = ũ(T ), hence u = ũ a.e. in Q T .

Ω m 0 u(0) dx ≤ Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m) m dxdt + T 0 Ω [ m H p (t, x, m, Dũ) • Du -m H(t, x, m, Du)] dxdt, Ω m 0 ũ(0) dx ≤ Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m)m dxdt + T 0 Ω [m H p (t, x, m, Du) • Dũ -m H(t, x, m, Dũ)] dxdt, Ω m 0 u(0) dx = Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m)m dxdt + T 0 Ω [m H p (t, x, m, Du) • Du -m H(t, x, m, Du)] dxdt, Ω m 0 ũ(0) dx = Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m) m dxdt + T 0 Ω [ m H p (t, x, m, Dũ) • Dũ -m H(t, x, m, Dũ)] dxdt,
and therefore 0 ≥ Ω (G(m(T )) -G( m(T ))) (m(T ) -m(T )) dx + T 0 Ω (F (t, x, m) -F (t, x, m))(m -m) dxdt + T 0 Ω m [H(t,

Existence and uniqueness for singular congestion

We consider the limit case when µ = 0; we aim at proving the existence and uniqueness of a weak solution of (1.3)-(1.5) as defined in Definition 3.7. Recall that it is not restrictive to assume (3.3). To prove existence, we consider the weak solutions (u µ , m µ ) of

-∂ t u µ -ν∆u µ + H(t, x, µ + m µ , Du µ ) = F (t, x, m µ ) , (t, x) ∈ (0, T ) × Ω (7.1) ∂ t m µ -ν∆m µ -div(m µ H p (t, x, µ + m µ , Du µ )) = 0 , (t, x) ∈ (0, T ) × Ω (7.2) m µ (0, x) = m 0 (x) , u µ (T, x) = G(x, m µ (T )) , x ∈ Ω (7.3)
for µ > 0 tending to 0. By (4.23), (4.24), we already know that the following estimates hold

Ω G(x, m µ (T ))m µ (T )dx + T 0 Ω F (t, x, m µ )m µ dxdt + (m µ ) α β-1 +1 L N +2 N (Q T ) ≤ C, (7.4) T 0 Ω m µ {H p (t, x, µ + m µ , Du µ ) • Du µ -H(t, x, µ + m µ , Du µ )} dxdt ≤ C, (7.5) T 0 Ω |Du µ | β (m µ + µ) α dxdt + T 0 Ω m µ |Du µ | β (m µ + µ) α dxdt ≤ C, (7.6)
for some C independent of µ.

With the same proof as for Lemma 4.4, we obtain: so using the bound of z µ in L 2 (Q T ) and the estimate (7.4) for m µ , the right-hand side is bounded as soon as is sufficiently small. Lemma 7.3.

(1) The family (m µ ) is relatively compact in L 1 (Q T ) and in C([0, T ]; W -1,r (Ω)) for some r > 1 (2) There exist m ∈ L 1 (Q T ) and u ∈ L q (0, T ; W 1,q (Ω)) for any q < N +2 N +1 such that, after the extraction of a subsequence (not relabeled), m µ → m in L 1 (Q T ) and almost everywhere, u µ → u and Du µ → Du in L 1 (Ω) and almost everywhere. Moreover,

T 0 Ω m -α |Du| β 1 {m>0} dxdt + T 0 Ω m 1-α |Du| β 1 {m>0} dxdt < +∞, (7.7) T 0 Ω F (x, m)m dxdt + m α β-1 +1 L N +2 N (Q T ) < ∞, (7.8)
and Du = 0 a.e. in {m = 0}. Proof. The compactness of m µ , u µ and Du µ are established exactly as in Proposition 4.5. Then we obtain (7.7) by using (7.6) and Fatou's Lemma in the set {m > 0}. We obtain (7.8) as an easy consequence of (7.4). If Du did not vanish a.e. in {m = 0}, there would exist ρ > 0 and a measurable subset E of {m = 0} in which |Du| > ρ. Then, by Fatou's lemma, E (m µ + µ) -α |Du µ | β dxdt would tend to +∞, in contradiction with (7.6). Therefore Du = 0 a.e. in {m = 0}. To pass to the limit in m µ H p (t, x, µ + m µ , Du µ ), we first observe that this sequence is equiintegrable in L 1 (Q T ) as a consequence of Remark 7.2. Let us now consider the case when α < 1: from

(3) If β ≤ 2 and α < 1, then m µ H p (t, x, µ + m µ , Du µ ) → mH p (t, x, m, Du) in L 1 (Q T ). If β < 2
(2.3) with µ = 0, (t, x, m, p) → mH p (t, x, m, p) is well defined in Q T × [0, +∞) × R N , and m µ H p (t, x, µ + m µ , Du µ ) → mH p (t, x, m, Du) almost everywhere. Then by Vitali's theorem, m µ H p (t, x, µ + m µ , Du µ ) → mH p (t, x, m, Du) in L 1 (Q T ). If β < 2 and 1 ≤ α ≤ 4(β-1)
β or β = 2 and 1 ≤ α < 2, then naming ξ µ = m µ H p (t, x, µ+m µ , Du µ ) for brevity, we proceed through the following steps:

(1) By Remark 7.2, we can extract a subsequence, not relabeled, such that ξ µ ξ in L 1+ (Q T ) weak.

(2) For any bounded and smooth function φ, Lebesgue theorem implies that φ(m µ ) → φ(m) in L p (Q T ) for any p ≥ 1. Therefore, ξ µ φ(m µ ) ξφ(m) in L 1 (Q T ) weak.

(3) If φ is also supported in [δ, +∞) for some positive δ, then using the almost everywhere convergence and Vitali's theorem, ξ µ φ(m µ ) → mH p (t, x, m, Du)φ(m) in L 1 (Q T ). This implies that (ξ -mH p (t, x, m, Du))φ(m) = 0, and that ξ coincides with mH p (t, x, m, Du) at almost every (t, x) such that m(t, x) > 0. (4) We notice that, by definition of ξ µ and assumption (2.3),

T 0 Ω |ξ µ φ(m µ )|dxdt ≤ c 2 T 0 Ω m µ |φ(m µ )|dxdt + c 2 T 0 Ω (m µ ) β-α |φ(m µ )|dxdt 1 β T 0 Ω |Du µ | β (µ + m µ ) α |φ(m µ )|dxdt 1 β and (7.6) implies that T 0 Ω |ξ µ φ(m µ )|dxdt ≤ c 2 T 0 Ω m µ |φ(m µ )|dxdt + C T 0 Ω (m µ ) β-α |φ(m µ )|dxdt 1 β .
Then using the weak convergence of ξ µ φ(m µ ) in L 1 (Q T ) and the fact that 0 < β -α ≤ 1, we can deduce

T 0 Ω |ξφ(m)|dxdt ≤ lim inf µ→0   c 2 T 0 Ω m µ |φ(m µ )|dxdt + C T 0 Ω (m µ ) β-α |φ(m µ )|dxdt 1 β   =c 2 T 0 Ω m|φ(m)|dxdt + C T 0 Ω m β-α |φ(m)|dxdt 1 β
.

Taking now φ(m) = exp(-Km), and letting K tend to +∞ yields that {m=0} |ξ|dxdt = 0, and that ξ = 0 a.e. in {m = 0}. Notice that we used α < β in this step, in particular if β = 2 we needed the restriction α < 2.

(5) Collecting the above results, we obtain that ξ = mH p (t, x, m, Du)1 {m>0} .

Finally, the weak convergence of m µ (t) in L 1 (Q T ) can be justified as in Proposition 4.5 (as well as the strong C 0 ([0, T ]; L 1 (Ω)) convergence if β < 2), and obtaining equation (3.6) is a consequence of the previous points.

We conclude the existence part by showing that u satisfies (3.5) and the energy equality holds true.

Proposition 7.4. Assume that β < 2 and 0 < α ≤ 4(β-1) β or β = 2 and 0 < α < 2. Let (u, m) be given by Lemma 7.3. Then inequality (3.5) holds true for every nonnegative ϕ ∈ C ∞ c ((0, T ] × Ω). Moreover, (u, m) satisfies the energy identity (3.7).

Proof. From (3.1) and (3.3), we deduce that for every nonnegative ϕ

∈ C ∞ c ((0, T ] × Ω), T 0 Ω u µ ϕ t dxdt -ν T 0 Ω u µ ∆ϕ dxdt + T 0 Ω H(t, x, µ + m µ , Du µ )ϕ1 {m>0} dxdt ≤ T 0 Ω F (t, x, m µ )ϕ dxdt + Ω G(x, m µ (T ))ϕ(T ) dx. ( 7
(T )) is equi-integrable, lim µ→0 Ω G(x, m µ (T ))ϕ(T ) dx = Ω R G(x, λ)dν x (λ)ϕ(T ) dx .
We easily pass to the limit in the other terms of (7.9) since u µ and F (t, x, m µ ) converge in L 1 (Ω). So we end up with

T 0 Ω u ϕ t dxdt -ν T 0 Ω u ∆ϕ dxdt + T 0 Ω H(t, x, m, Du)ϕ1 {m>0} dxdt ≤ T 0 Ω F (t, x, m)ϕ dxdt + Ω R G(x, λ)dν x (λ)ϕ(T ) dx . (7.10)
Thanks to (7.10), the conclusion of Lemma 4.7 is also true. Therefore,

lim µ→0 Ω u µ (0)m 0 dx ≤ u(0), m 0 .
Starting now from the energy identity for the solution of (1.3)-(1.5) with µ > 0: Reasoning exactly as in the previous section, from the estimate (7.4) and assumption (2.10) we obtain both G(x, m(T ))m(T ) ∈ L 1 (Ω) and R G(x, λ)λ dν x (λ) ∈ L 1 (Ω). Therefore, we are allowed to use (5.13) Now we proceed as in Lemma 6.1 using the monotonicity of G and we conclude that G(x, m µ (T )) → G(x, m(T )) in L 1 (Ω) and the energy identity holds true. Coming back to (7.10), now we know that R G(x, λ)dν x (λ) = G(x, m(T )) and so (3.5) is proven. We conclude with the uniqueness of the weak solution. To this purpose, since u is now only a subsolution, we will need a stronger version of the energy equality for the uniqueness argument to succeed. for almost every t ∈ (0, T ).

T 0 Ω m µ (H p (t, x, µ + m µ , Du µ ) • Du µ -H(t, x, µ + m µ , Du µ )) dxdt + Ω m µ (T )G(x, m µ (T ))dx + T 0 Ω m µ F (t, x, m µ ) dxdt = Ω u µ (
Proof. Let us set u k := (u -k) + . We first observe that, whenever u is a subsolution, then u k is also a subsolution of a similar problem, namely where E has been defined in the proof of Lemma 6.5. From the last two equations in (7.17 Since u -ũ is only time-dependent, this implies u = ũ a.e. in Q T .

Lemma 4 . 6 .

 46 Let (u, m) be given by Proposition 4.5. Then

(4. 29 )

 29 Combining the latter three points yields(4.26).Finally, we notice that the bound Ω G (x, m (T ))m (T ) ≤ C given by (4.7) implies that both G(x, m(T ))m(T ) ∈ L 1 (Ω) and Ω R G(x, λ)λ dν x (λ) ∈ L 1 (Ω); indeed, on the one hand, thanks to (2.10) we haveΩ g(ρ m (T ))ρ m (T ) dx ≤ Cso by convexity of g(s)s and weak convergence of ρ m (T ) we deduce that g(m(T ))m(T ) ∈ L 1 (Ω), hence G(x, m(T ))m(T ) ∈ L 1 (Ω) again by (2.10), and estimate (4.27) holds. On the other hand, by properties of Young measures we have

Lemma 4 . 7 .

 47 Let (m , u ) be given by Proposition 4.5. Possibly after another extraction of a subsequence, u | t=0 converges weakly * to a bounded measure χ on Ω. Moreover u(0) is a well defined Radon measure on Ω and u(0) ≥ χ.

Lemma 5 . 1 .

 51 Assume that H satisfies assumptions (2.2)-(2.3) with µ > 0, and that (2.8)-(2.10) hold true.

( 5 . 4 )

 54 where we denote wδ = [( b m) ρ δ ] and wδ,ε = T 0 wδ (s) ξ ε (s -t) ds.

Ω

  (G(x, m(T )) -G(x, m(T ))) (m(T ) -m(T )) dx = 0, , x, m, Du, m, Dũ) dxdt = 0.(6.6) From (6.6), we deduce that m = m and Du = Dũ a.e. in Q T , and that G(•, m(T, •)) = G(•, m(T, •)) a.e. in Ω.

and 1 ≤

 1 α ≤ 4(β-1) β or β = 2 and 1 ≤ α < 2, then m µ H p (t, x, µ + m µ , Du µ ) mH p (t, x, m, Du)1 {m>0} weakly in L 1 (Q T ). (4) m ∈ C 0 ([0, T ]; L 1 (Ω)) and m µ (t) converges to m(t) weakly in L 1 (Ω) for any t ∈ [0, T ]. If β < 2, then m µ → m in C([0, T ]; L 1 (Ω)). Finally, (3.6) holds for any ϕ ∈ C ∞ c ([0, T )×Ω).

1

  0)m 0 dx using Fatou's lemma thanks to (2.6) and (2.11), we obtainlim sup µ→0 Ω G(x, m µ (T ))m µ (T )dx ≤ u(0), m 0 -{m>0} m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt.

1

  in Lemma 5.3 with m = m and ũ = u, and we get lim sup µ→0 Ω G(x, m µ (T ))m µ (T )dx ≤ u(0), m 0 -{m>0} m (H p (t, x, m, Du) • Du -H(t, x, m, Du)) dxdt ≤ Ω R G(x, λ)dν x (λ) m(T )dx.

Lemma 7 . 5 .F

 75 Let (u, m) be a weak solution of (1.3)-(1.5). Then um ∈ L 1 (Q T ) and Ω m(t) u(t) dx = Ω G(x, m(T )) m(T ) dx + T t Ω (s, x, m)m dxds + T t Ω [m H p (s, x, m, Du) • Du -m H(s, x, m, Du)] 1 {m>0} dxds(7.11)

F[ 1 F 0 . 0 ΩFFLemma 7 . 6 . 0 ΩF 0 Ω 1 { 1 {

 100760011 (t, x, m)1 {u>k} ϕ dxdt + Ω (G(x, m(T )) -k) + ϕ(T ) dx Henceforth, we can proceed exactly as in Lemma 5.3 with m = m and we obtain (recall that m δ (t, x) = (m(t) ρ δ )(x)) {m>0} H(s, x, m, Du)m δ -Du • w δ ] 1 {u>k} dxds + Ω m δ (t) u k (t) dx ≤ (s, x, m)m δ 1 {u>k} dxds + Ω (G(x, m(T )) -k) + m δ (T ) dx (7.12)which holds for almost every t since m δ ∈ L ∞ (Q T ) (it is indeed even continuous) and u k ∈ L 1 (Q T ) and is bounded below. Starting from (7.12) and following the same steps as in Lemma 5.3, we obtain thatΩ m δ (t) u k (t) dx ≤ T t Ω φ δ 1 {u>k} dxds + Ω χ δ 1 {G(x,m(T ))>k} dxfor some sequences φ δ , χ δ which are convergent, and therefore equi-integrable, in L 1 (Q T ) and in L 1 (Ω) respectively. This implies thatsup δ Ω m δ (t) u k (t) dx k→∞ → Since u -u k is bounded, we know that m δ (u -u k ) strongly converges in L 1 (Q T ) as δ → 0.Therefore, a standard argument allows us to conclude that (7.13)m δ u → mu in L 1 (Q T ).This means that we can repeat the argument of Lemma 5.3 integrating either in (0, t) or in (t, T ) obtaining botht 0 Ω [H(s, x, m, Du)m δ -Du • w δ ]1 {m>0} dxds + m δ (0) u(0) dx ≤ t (s, x, m)m δ dxds + Ω u(t) m δ (t) dx and T t Ω [H(s, x, m, Du)m δ -Du • w δ ]1 {m>0} dxds + Ω m δ (t) u(t) dx ≤ (s, x, m)m δ dxds + Ω G(x, m(T )) m δ (T ) dxfor almost every t ∈ (0, T ). Passing to the limit as δ → 0 is now allowed thanks to (7.13), besides the arguments already given in Lemma 5.3. Finally, as δ → 0 we obtain bothm 0 , u(0) dx ≤ p (s, x, m, Du) • Du -m H(s, x, m, Du)] 1 {m>0} dxds and Ω m(t) u(t) dx ≤ Ω G(x, m(T )) m(T ) dx + p (s, x, m, Du) • Du -m H(s, x, m, Du)] 1 {m>0} dxds .On account of the energy equality (3.7) which holds by definition of weak solution, we conclude that (7.11) holds true, for almost every t ∈ (0, T ). Under all the assumptions made in Theorem 3.8, there is a unique weak solution of (1.3)-(1.5).Proof. Let (u, m) and (ũ, m) be two weak solutions of (1.3)-(1.5). From Lemma 5.3 and Proposition 7.4, we know thatΩ m 0 u(0) dx ≤ Ω G(x, m(T )) m(T ) dx + T (t, x, m) m dxdt + T 0 Ω [ m H p (t, x, m, Dũ) • Du -m H(t, x, m, Du)] 1 {m>0 , m>0} dxdt, Ω m 0 ũ(0) dx ≤ Ω G(x, m(T )) m(T ) dx + T 0 Ω F (t, x, m)m dxdt + T 0 Ω [m H p (t, x, m, Du) • Dũ -m H(t, x, m, Dũ)] 1 {m>0 , m>0} dxdt, Ω m 0 u(0) dx = Ω G(x, m(T )) m(T ) dx + p (t, x, m, Du) • Du -m H(t, x, m, Du)] 1 {m>0} dxdt, Ω m 0 ũ(0) dx = Ω G(x, m(T )) m(T ) dx + p (t, x, m, Dũ) • Dũ -m H(t, x, m, Dũ)] 1 { m>0} dxdt, x, m(T )) -G(x, m(T ))) (m(T ) -m(T )) dx + T 0 Ω (F (t, x, m) -F (t, x, m))(m -m) dxdt + T 0 Ω 1 {m>0 , m>0} m [H(t, x, m, Du) -H(t, x, m, Dũ) -H p (t, x, m, Dũ) • (Du -Dũ)] dxdt + T 0 Ω 1 {m>0 , m>0} m [H(t, x, m, Dũ) -H(t, x, m, Du) -H p (t, x, m, Du) • (Dũ -Du)] dxdt + T 0 Ω 1 {m>0 , m=0} m [-H(t,x, m, Du) + H p (t, x, m, Du) • Du] dxdt + T m>0 , m=0} m [-H(t, x, m, Dũ) + H p (t, x, m, Dũ) • Dũ] dxdt. (7.14) Note that from (2.1) and the convexity of H(t, x, m, p) w.r.t p, m > 0 ⇒ -H(t, x, m, Dũ) + H p (t, x, m, Dũ) • Dũ ≥ -H(t, x, m, 0) ≥ 0, (7.15) m > 0 ⇒ -H(t, x, m, Du) + H p (t, x, m, Du) • Du ≥ -H(t, x, m, 0) ≥ 0. (7.16) Therefore all integrals in (7.14) involve nonnegative functions, which then must be zero almost everywhere. We deduce thatΩ (G(m(T )) -G( m(T ))) (m(T ) -m(T )) dx = 0, T 0 Ω 1 {m>0 , m>0} E(t, x, m, Du, m, Dũ) dxdt = 0, m>0} m(-H(t, x, m, Dũ) + H p (t, x, m, Dũ) • Dũ + F (t, x, m) -F (t, x, 0))dxdt = 0, m=0,m>0} m(-H(t, x, m, Du) + H p (t, x, m, Du) • Du + F (t, x, m) -F (t, x, 0)) = 0, (7.17)

  ),(7.15)-(7.16) and(3.8), we deduce that the sets { m = 0, m > 0} and {m = 0, m > 0} have measure 0. From the second equation in (7.17) and the assumptions made for uniqueness (which apply to the function E), as in Lemma 6.5 we deduce that m = m and Du = Dũ a.e. in {m > 0, m > 0}. Combining all the previous observations, m = m almost everywhere in Q T . Using also that that Du = 0 a.e. in {m = 0} and Dũ = 0 a.e. in { m = 0}, we deduce that Du = Dũ a.e. in Q T . Therefore, u -ũ = Ω (u(t) -ũ(t))dx. By applying Lemma 7.5 to u and ũ, subtracting the two equalities on account of the above informations we finally get Ω m(t)(u -ũ)(t) dx = 0 .

  and (2.10), the functions F, G are bounded below, namely there exists a constant c 4 ∈ R such that (2.11)F (t, x, m) ≥ c 4 , ∀ m ∈ R + , a.e. (t, x) ∈ Q T , ≥ c 4 , ∀ m ∈ R + , a.e. x ∈ Ω .

	and
	(2.12) G(x, m) Assumptions (2.1)-(2.5) and (2.8)-(2.10

  The latter inequality implies that its right hand side is bounded uniformly in , and we obtain (4.7). From (4.16), we also obtain (4.8) and we estimate the first term in (4.9).Finally, from (4.17), we see that sup t Ω I (m (t))dx ≤ C which implies that

	(4.19)	sup t	Ω	m (t)(T 1/ m (t))
	Then using (2.2) and (2.4), (4.8) and (4.19) imply that
	(4.20)

α β-1 dx ≤ C.

  x, m, Du) -H(t, x, m, Dũ) -H p (t, x, m, Dũ)• (Du -Dũ)] dxdt , p 1 , m 2 , p 2 ) = -(H(t, x, m 1 , p 1 ) -H(t, x, m 2 , p 2 ))(m 1 -m 2 ) + (m 1 H p (t, x, m 1 , p 1 ) -m 2 H p (t, x, m 2 , p 2 )) (p 1 -p 2 ) + (F (t, x, m 1 ) -F (t, x, m 2 ))(m 1 -m 2 ) and r = p 2 -p 1 , p s = p 1 + sr, z = m 2 -m 1 , m s = m 1 + sz.From the assumptions, the function h : s → -zH(t, x, m s , p s ) + m s H p (t, x, m s , p s ) • r + zF (t, x, m s ) is increasing on [0, 1]; this yields that E(t, x, m 1 , p 1 , m 2 , p 2 ) ≥ 0 for all m 1 , m 2 ≥ 0 and all p 1 , p 2 ∈ R N and that E(t, x, m 1 , p 1 , m 2 , p 2 ) = 0 if and only if m 1 = m 2 and p 1 = p 2 . Thus, (6.5) and the assumptions imply

	(6.5)
	T
	+	m [H(t, x, m, Dũ) -H(t, x, m, Du) -H p (t, x, m, Du) • (Dũ -Du)] dxdt.
	0	Ω
	Let us define
	E(t, x, m 1

  By properties of Young measures, and since G(x, m µ

	.9)
	By using the a.e. convergence of H(t, x, µ + m µ , Du µ )1 {m>0} to H(t, x, m, Du)1 {m>0} and
	Fatou's lemma,
	T		T
	0 Ω	H(t, x, m, Du)ϕ(t, x)1 {m>0} dxdt ≤ lim inf µ→0	0 Ω	H(t, x, µ + m µ , Du µ )1 {m>0} ϕ(t, x)dxdt.

Acknowledgement.

Most of this work was done during a visit of A. Porretta at the University Paris Diderot (Paris VII). A. Porretta wishes to thank this institution for the invitation and the very kind hospitality given on this occasion. The work was also supported by the Indam Gnampa project 2015 Processi di diffusione degeneri o singolari legati al controllo di dinamiche stocastiche and by ANR projects ANR-12-MONU-0013 and ANR-16-CE40-0015-01.